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