結果
| 問題 |
No.574 正多面体サイコロ
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2025-10-04 13:18:00 |
| 言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
AC
|
| 実行時間 | 4 ms / 2,000 ms |
| コード長 | 1,025 bytes |
| コンパイル時間 | 1,810 ms |
| コンパイル使用メモリ | 95,440 KB |
| 実行使用メモリ | 7,716 KB |
| 最終ジャッジ日時 | 2025-10-04 13:18:04 |
| 合計ジャッジ時間 | 2,315 ms |
|
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 20 |
ソースコード
#include <iostream>
#include <algorithm>
#include <vector>
#include <string>
#include <set>
#include <unordered_set>
#include <cstring>
#include <queue>
#include <bitset>
#include <map>
#include <stack>
using namespace std;
long double dp[101][101];
int f,n,k;
long double C[101][101];
int main()
{
cin.tie(0);
ios::sync_with_stdio(false);
cin >> f >> n >> k;
long double p = 1;
p/=f;
C[0][0] = 1;
for(int i=1;i<=n;i++)
{
C[i][0] = 1;
C[i][i] = 1;
for(int j=1;j<i;j++)
{
C[i][j] = C[i-1][j-1] + C[i-1][j];
}
}
k = (n+1) - k;
dp[0][0] = 1;
for(int i=0;i<=n;i++)
{
for(int j=0;j<f;j++)
{
if(dp[i][j]==0) continue;
//cout << i << ' ' << j << ' ' << dp[i][j] << '\n';
dp[i][j+1] += dp[i][j];
double P = 1;
for(int a=i+1;a<=n;a++)
{
P*=p;
long double value = P;
value*=C[(n-i)][a-i];
value*=dp[i][j];
if(i<k && a>=k) value*=(j+1);
dp[a][j+1] += value;
}
}
}
cout << fixed;
cout.precision(16);
cout << dp[n][f] << '\n';
return 0;
}