結果
問題 | No.574 正多面体サイコロ |
ユーザー |
![]() |
提出日時 | 2017-10-07 17:35:06 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 196 ms / 2,000 ms |
コード長 | 874 bytes |
コンパイル時間 | 1,874 ms |
コンパイル使用メモリ | 159,748 KB |
実行使用メモリ | 164,608 KB |
最終ジャッジ日時 | 2024-11-17 04:33:37 |
合計ジャッジ時間 | 3,940 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#include <bits/stdc++.h> #define REP(i, a, n) for(ll i = ((ll) a); i < ((ll) n); i++) #define EPS 1e-10 using namespace std; typedef long long ll; ll F, N, K; double dp[20][101][101][101]; int main(void) { cin >> F >> N >> K; REP(f, 0, F) REP(i, 0, 101) REP(j, 0, 101) REP(k, 0, 101) dp[f][i][j][k] = 0; REP(f, 0, F) dp[f][0][0][0] = 1; REP(f, 0, F) { REP(i, 0, N) REP(j, 0, N) REP(k, 0, N) if(dp[f][i][j][k] > EPS) { REP(g, 0, F) { if(g > f) dp[f][i + 1][j + 1][k] += dp[f][i][j][k]; if(g == f) dp[f][i + 1][j][k + 1] += dp[f][i][j][k]; if(g < f) dp[f][i + 1][j][k] += dp[f][i][j][k]; } } } double p = 1; REP(i, 0, N) p /= F; double ans = 0; REP(f, 0, F) { REP(j, 0, N + 1) REP(k, 0, N + 1) { if(j < K && K <= j + k) ans += (f + 1) * dp[f][N][j][k] * p; } } printf("%.15lf\n", ans); }