結果
問題 | No.574 正多面体サイコロ |
ユーザー |
|
提出日時 | 2017-10-06 23:31:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 53 ms / 2,000 ms |
コード長 | 1,282 bytes |
コンパイル時間 | 1,152 ms |
コンパイル使用メモリ | 106,004 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-17 02:00:12 |
合計ジャッジ時間 | 1,943 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#define _USE_MATH_DEFINES#include <cstdio>#include <iostream>#include <sstream>#include <fstream>#include <iomanip>#include <algorithm>#include <cmath>#include <complex>#include <string>#include <vector>#include <list>#include <queue>#include <stack>#include <set>#include <map>#include <bitset>#include <numeric>#include <limits>#include <climits>#include <cfloat>#include <functional>#include <iterator>using namespace std;int main(){int f, n, k;cin >> f >> n >> k;double ans = 0.0;for(int x=1; x<=f; ++x){vector<vector<double> > dp(n+1, vector<double>(n+1, 0.0));dp[0][0] = 1.0;for(int i=0; i<n; ++i){vector<vector<double> > nextDp(n+1, vector<double>(n+1, 0.0));for(int a=0; a<n; ++a){for(int b=0; b<n; ++b){nextDp[a][b] += dp[a][b] * (1.0 / f);nextDp[a+1][b] += dp[a][b] * ((x - 1.0) / f);nextDp[a][b+1] += dp[a][b] * ((f - x) / (double)f);}}dp.swap(nextDp);}for(int a=0; a<=n-k; ++a){for(int b=0; b<k; ++b){ans += dp[a][b] * x;}}}printf("%.10f\n", ans);return 0;}