結果

問題 No.574 正多面体サイコロ
コンテスト
ユーザー ryoissy
提出日時 2017-10-06 23:20:07
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 49 ms / 2,000 ms
コード長 775 bytes
コンパイル時間 1,224 ms
コンパイル使用メモリ 158,952 KB
実行使用メモリ 11,676 KB
最終ジャッジ日時 2024-11-17 01:57:30
合計ジャッジ時間 2,175 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:11:14: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   11 |         scanf("%d%d%d",&f,&n,&k);
      |         ~~~~~^~~~~~~~~~~~~~~~~~~

ソースコード

diff # 

#include <bits/stdc++.h>
#define MOD 1000000007LL
using namespace std;
typedef long long ll;
typedef pair<int,int> P;

int f,n,k;
double dp[101][101][101];

int main(void){
	scanf("%d%d%d",&f,&n,&k);
	double ans=0.0;
	for(int num=1;num<=f;num++){
		for(int i=0;i<=n;i++){
			for(int j=0;j<=n;j++){
				for(int l=0;l<=n;l++){
					dp[i][j][l]=0.0;
				}
			}
		}
		dp[0][0][0]=1.0;
		for(int i=0;i<n;i++){
			for(int j=0;j<n;j++){
				for(int l=0;l<n;l++){
					dp[i+1][j+1][l]+=(double)1.0/f*dp[i][j][l];
					dp[i+1][j][l+1]+=(double)(f-num)/f*dp[i][j][l];
					dp[i+1][j][l]+=(double)(num-1)/f*dp[i][j][l];
				}
			}
		}
		for(int i=0;i<=n;i++){
			for(int j=0;j<=n;j++){
				if(j<k && i+j>=k)ans+=(double)num*dp[n][i][j];
			}
		}
	}

	printf("%.10f\n",ans);
	return 0;
}
0