結果
| 問題 |
No.574 正多面体サイコロ
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2017-10-06 23:36:03 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 32 ms / 2,000 ms |
| コード長 | 896 bytes |
| コンパイル時間 | 702 ms |
| コンパイル使用メモリ | 80,256 KB |
| 実行使用メモリ | 6,820 KB |
| 最終ジャッジ日時 | 2024-11-17 02:12:04 |
| 合計ジャッジ時間 | 1,757 ms |
|
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 20 |
ソースコード
#include <iostream>
#include <algorithm>
#include <vector>
#include <map>
#include <cmath>
using namespace std;
int main() {
int F, N, K;
cin >> F >> N >> K;
double ans = 0;
for (int t = 1; t <= F; t++) {
static double dp0[101][101];
static double dp1[101][101];
fill_n(*dp0, 101 * 101, 0.0);
dp0[0][0] = 1;
for (int k = 0; k < N; k++) {
fill_n(*dp1, 101 * 101, 0.0);
for (int i = 0; i <= k; i++) {
for (int j = 0; i + j <= k; j++) {
dp1[i + 1][j] += 1.0 * (F - t) / F * dp0[i][j];
dp1[i][j + 1] += 1.0 / F * dp0[i][j];
dp1[i][j] += 1.0 * (t - 1) / F * dp0[i][j];
}
}
swap(dp0, dp1);
}
for (int i = 0; i <= N; i++) {
for (int j = 0; i + j <= N; j++) {
if (i + 1 <= K && K <= i + j) {
ans += t * dp0[i][j];
}
}
}
}
printf("%.20f\n", ans);
}