結果
問題 | No.1688 Veterinarian |
ユーザー | umezo |
提出日時 | 2021-09-24 22:22:02 |
言語 | C++17 (gcc 12.3.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 207 ms / 3,000 ms |
コード長 | 1,201 bytes |
コンパイル時間 | 2,461 ms |
コンパイル使用メモリ | 201,600 KB |
実行使用メモリ | 58,776 KB |
最終ジャッジ日時 | 2023-09-18 21:40:44 |
合計ジャッジ時間 | 3,688 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
5,556 KB |
testcase_01 | AC | 2 ms
4,380 KB |
testcase_02 | AC | 2 ms
5,776 KB |
testcase_03 | AC | 2 ms
4,376 KB |
testcase_04 | AC | 2 ms
5,616 KB |
testcase_05 | AC | 3 ms
6,104 KB |
testcase_06 | AC | 2 ms
6,164 KB |
testcase_07 | AC | 64 ms
58,648 KB |
testcase_08 | AC | 207 ms
58,776 KB |
testcase_09 | AC | 189 ms
56,052 KB |
testcase_10 | AC | 17 ms
12,816 KB |
testcase_11 | AC | 25 ms
16,732 KB |
testcase_12 | AC | 10 ms
12,952 KB |
testcase_13 | AC | 3 ms
6,580 KB |
testcase_14 | AC | 9 ms
10,740 KB |
testcase_15 | AC | 10 ms
9,496 KB |
testcase_16 | AC | 4 ms
5,048 KB |
ソースコード
#define rep(i,n) for(int i=0;i<(int)(n);i++) #define ALL(v) v.begin(),v.end() typedef long long ll; #include<bits/stdc++.h> using namespace std; double f(ll x){ return (double)x*(x-1)/2; } double dp[52][52][52][52]; int main(){ ios::sync_with_stdio(false); std::cin.tie(nullptr); int A,B,C,N; cin>>A>>B>>C>>N; dp[A][B][C][0]=1.0; for(int i=0;i<N;i++){ for(int a=0;a<=A;a++){ for(int b=0;b<=B;b++){ for(int c=0;c<=C;c++){ if(f(a+b+c)==0){ dp[a][b][c][i+1]+=dp[a][b][c][i]; continue; } if(a>0) dp[a-1][b][c][i+1]+=dp[a][b][c][i]*f(a)/f(a+b+c); if(b>0) dp[a][b-1][c][i+1]+=dp[a][b][c][i]*f(b)/f(a+b+c); if(c>0) dp[a][b][c-1][i+1]+=dp[a][b][c][i]*f(c)/f(a+b+c); dp[a][b][c][i+1]+=dp[a][b][c][i]*(1-(f(a)+f(b)+f(c))/f(a+b+c)); } } } } double asum=0,bsum=0,csum=0; for(int a=0;a<=A;a++){ for(int b=0;b<=B;b++){ for(int c=0;c<=C;c++){ asum+=(A-a)*dp[a][b][c][N]; bsum+=(B-b)*dp[a][b][c][N]; csum+=(C-c)*dp[a][b][c][N]; } } } cout<<fixed<<setprecision(10)<<asum<<" "<<bsum<<" "<<csum<<endl; return 0; }