結果

問題 No.574 正多面体サイコロ
コンテスト
ユーザー chocorusk
提出日時 2019-07-19 23:53:49
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 2 ms / 2,000 ms
コード長 1,212 bytes
コンパイル時間 846 ms
コンパイル使用メモリ 97,096 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-12-26 05:08:02
合計ジャッジ時間 1,836 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <unordered_set>
#include <random>
#include <cassert>
#include <fstream>
#define popcount __builtin_popcount
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;

int main()
{
    int f, n, k;
    cin>>f>>n>>k;
    double p[110];
    p[f+1]=0;
    for(int i=1; i<=f; i++){
        p[i]=0;
        double c=1;
        for(int j=1; j<k; j++){
            c*=(double)(n-j+1);
            c/=(double)j;
        }
        for(int j=k; j<=n; j++){
            c*=(double)(n-j+1);
            c/=(double)j;
            double a=c;
            for(int l=0; l<j; l++){
                a*=(double)(f-i+1);
                a/=(double)f;
            }
            for(int l=0; l<n-j; l++){
                a*=(double)(i-1);
                a/=(double)f;
            }
            p[i]+=a;
        }
    }
    double ans=0;
    for(int i=1; i<=f; i++){
        ans+=(double)i*(p[i]-p[i+1]);
    }
    printf("%.7lf\n", ans);
    return 0;
}
0