結果

問題 No.1688 Veterinarian
ユーザー hourenhouren
提出日時 2021-09-24 22:32:17
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 110 ms / 3,000 ms
コード長 1,499 bytes
コンパイル時間 1,662 ms
コンパイル使用メモリ 177,112 KB
実行使用メモリ 61,696 KB
最終ジャッジ日時 2024-07-05 10:58:33
合計ジャッジ時間 2,568 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
using ll = long long;
const ll MOD = 1000000007;
using P = pair<int,int>;
#define rep(i, n) for(int i = 0; i < n; i++)
#define all(x) (x).begin(),(x).end()

int main(){
    int n,a,b,c;
    cin >> a >> b >> c >> n;
    vector<vector<vector<vector<double>>>> dp(n+1,vector<vector<vector<double>>>(a+1,vector<vector<double>>(b+1,vector<double>(c+1,0))));
    dp[0][a][b][c] = 1;
    rep(i,n){
        for(double j=1;j<=a;j++){
            for(double k=1;k<=b;k++){
                for(double l=1;l<=c;l++){
                    double na,nb,nc;
                    na = dp[i][j][k][l] * (j/(j+k+l)) * ((j-1)/(j+k+l-1));
                    nb = dp[i][j][k][l] * (k/(j+k+l)) * ((k-1)/(j+k+l-1));
                    nc = dp[i][j][k][l] * (l/(j+k+l)) * ((l-1)/(j+k+l-1));
                    //printf("%.6lf %.6lf %.6lf\n",na,nb,nc);
                    dp[i+1][j-1][k][l] += na;
                    dp[i+1][j][k-1][l] += nb;
                    dp[i+1][j][k][l-1] += nc;
                    dp[i+1][j][k][l] += dp[i][j][k][l] - na - nb - nc;
                }
            }
        }
    }
    double ansa = 0, ansb = 0, ansc = 0;
    rep(i,a+1){
        rep(j,b+1){
            rep(k,c+1){
                ansa += (a-(double)i) * dp[n][i][j][k];
                ansb += (b-(double)j) * dp[n][i][j][k];
                ansc += (c-(double)k) * dp[n][i][j][k];
            }
        }
    }
    printf("%.12lf %.12lf %.12lf\n",ansa,ansb,ansc);
    return 0;
}
0