結果
問題 | No.1688 Veterinarian |
ユーザー |
![]() |
提出日時 | 2021-09-24 22:13:31 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 122 ms / 3,000 ms |
コード長 | 1,719 bytes |
コンパイル時間 | 2,441 ms |
コンパイル使用メモリ | 193,852 KB |
最終ジャッジ日時 | 2025-01-24 17:22:56 |
ジャッジサーバーID (参考情報) |
judge5 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 14 |
ソースコード
#include <bits/stdc++.h> using namespace std; double dp[51][51][51][51]; int main() { int a,b,c,d; cin >> a >> b >> c >> d; dp[0][a][b][c] = 1; for(int i = 0; i < d; i++) { for(int j = 0; j <= a; j++) { for(int k = 0; k <= b; k++) { for(int l = 0; l <= c; l++) { if(j+k+l <= 1) { continue; } double pp = 1; if(j) { double p = 1; p *= (double)(j*(j-1))/((j+k+l)*(j+k+l-1)); dp[i+1][j-1][k][l] += dp[i][j][k][l]*p; pp -= p; } if(k) { double p = 1; p *= (double)(k*(k-1))/((j+k+l)*(j+k+l-1)); dp[i+1][j][k-1][l] += dp[i][j][k][l]*p; pp -= p; } if(l) { double p = 1; p *= (double)(l*(l-1))/((j+k+l)*(j+k+l-1)); dp[i+1][j][k][l-1] += dp[i][j][k][l]*p; pp -= p; } dp[i+1][j][k][l] += dp[i][j][k][l]*pp; } } } } double ans1 = 0,ans2 = 0,ans3 = 0; for(int i = 0; i <= a; i++) { for(int j = 0; j <= b; j++) { for(int k = 0; k <= c; k++) { ans1 += (a-i)*dp[d][i][j][k]; ans2 += (b-j)*dp[d][i][j][k]; ans3 += (c-k)*dp[d][i][j][k]; } } } cout << fixed << setprecision(18) << ans1 << " " << ans2 << " " << ans3 << endl; }