結果
問題 | 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; }