結果
問題 | No.1688 Veterinarian |
ユーザー |
![]() |
提出日時 | 2022-01-14 16:48:49 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,665 bytes |
コンパイル時間 | 1,139 ms |
コンパイル使用メモリ | 87,024 KB |
実行使用メモリ | 159,232 KB |
最終ジャッジ日時 | 2024-11-19 10:28:29 |
合計ジャッジ時間 | 2,204 ms |
ジャッジサーバーID (参考情報) |
judge3 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 WA * 2 |
other | AC * 1 WA * 13 |
ソースコード
#include <iostream> #include <cstdio> #include <string> #include <algorithm> #include <utility> #include <cmath> #include <vector> #include <stack> #include <queue> #include <deque> #include <set> #include <map> #include <tuple> #include <numeric> #include <functional> using namespace std; typedef long long ll; typedef vector<ll> vl; typedef vector<vector<ll>> vvl; typedef pair<ll, ll> P; #define rep(i, n) for(ll i = 0; i < n; i++) #define exrep(i, a, b) for(ll i = a; i <= b; i++) #define out(x) cout << x << endl #define exout(x) printf("%.10f\n", x) #define chmax(x, y) x = max(x, y) #define chmin(x, y) x = min(x, y) #define all(a) a.begin(), a.end() #define rall(a) a.rbegin(), a.rend() #define pb push_back #define re0 return 0 const ll mod = 998244353; const ll INF = 1e16; double dp1[51][51][51][51]; double dp2[51][51][51][51]; double dp3[51][51][51][51]; int c2(int x) { return x * (x - 1) / 2; } int main() { int a, b, c, n; cin >> a >> b >> c >> n; rep(i, n) { for(int j1 = a; j1 >= 0; j1--) { for(int j2 = b; j2 >= 0; j2--) { for(int j3 = c; j3 >= 0; j3--) { if(j1 + j2 + j3 <= 1) { continue; } double x1 = 0.0; if(j1 >= 2) { x1 = (double)c2(j1) / (double)c2(j1 + j2 + j3); } double x2 = 0.0; if(j2 >= 2) { x2 = (double)c2(j2) / (double)c2(j1 + j2 + j3); } double x3 = 0.0; if(j3 >= 2) { x3 = (double)c2(j3) / (double)c2(j1 + j2 + j3); } dp1[i+1][j1 - 1][j2][j3] += x1 * (dp1[i][j1][j2][j3] + 1.0); dp2[i+1][j1][j2 - 1][j3] += x2 * (dp2[i][j1][j2][j3] + 1.0); dp3[i+1][j1][j2][j3 - 1] += x3 * (dp3[i][j1][j2][j3] + 1.0); dp1[i+1][j1][j2][j3] += (1.0 - (x1 + x2 + x3)) * dp1[i][j1][j2][j3]; dp2[i+1][j1][j2][j3] += (1.0 - (x1 + x2 + x3)) * dp2[i][j1][j2][j3]; dp3[i+1][j1][j2][j3] += (1.0 - (x1 + x2 + x3)) * dp3[i][j1][j2][j3]; } } } } double ans1 = 0.0; double ans2 = 0.0; double ans3 = 0.0; exrep(j1, 0, a) { exrep(j2, 0, b) { exrep(j3, 0, c) { ans1 += dp1[n][j1][j2][j3]; ans2 += dp2[n][j1][j2][j3]; ans3 += dp3[n][j1][j2][j3]; } } } printf("%.10f %.10f %.10f\n", ans1, ans2, ans3); re0; }