結果
問題 | No.574 正多面体サイコロ |
ユーザー |
![]() |
提出日時 | 2017-10-07 10:09:54 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 39 ms / 2,000 ms |
コード長 | 1,857 bytes |
コンパイル時間 | 1,807 ms |
コンパイル使用メモリ | 166,452 KB |
実行使用メモリ | 5,248 KB |
最終ジャッジ日時 | 2024-11-17 04:22:30 |
合計ジャッジ時間 | 2,822 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#include<bits/stdc++.h> #define rep(i,a,b) for(int i=a;i<b;i++) #define rrep(i,a,b) for(int i=a;i>=b;i--) #define fore(i,a) for(auto &i:a) #pragma GCC optimize ("-O3") using namespace std; void _main(); int main() { cin.tie(0); ios::sync_with_stdio(false); _main(); } //--------------------------------------------------------------------------------------------------- /*--------------------------------------------------------------------------------------------------- ∧_∧ ∧_∧ (´<_` ) Welcome to My Coding Space! ( ´_ゝ`) / ⌒i / \ | | / / ̄ ̄ ̄ ̄/ | __(__ニつ/ _/ .| .|____ \/____/ (u ⊃ ---------------------------------------------------------------------------------------------------*/ int F, N, K; //--------------------------------------------------------------------------------------------------- void _main() { cin >> F >> N >> K; double ans = 0; rep(x, 1, F + 1) { double p1 = (double)(F - x) / (double)F; double p2 = 1.0 / (double)F; double p3 = 1 - (p1 + p2); double p = 0; rep(i, 0, K) rep(j, 0, N - i + 1) { if (i + j < K) continue; int k = N - (i + j); // より大きい数がi回、同じ数がj回、より小さい数がk回 double pp = 1; rep(a, 0, i) pp *= p1; rep(a, 0, j) pp *= p2; rep(a, 0, k) pp *= p3; rep(a, 0, N) pp *= (a + 1); rep(a, 0, i) pp /= (a + 1); rep(a, 0, j) pp /= (a + 1); rep(a, 0, k) pp /= (a + 1); p += pp; } ans += x * p; } printf("%.10f\n", ans); }