結果
問題 | No.1688 Veterinarian |
ユーザー | Manuel1024 |
提出日時 | 2022-12-12 06:36:45 |
言語 | C++14 (gcc 12.3.0 + boost 1.83.0) |
結果 |
MLE
|
実行時間 | - |
コード長 | 1,270 bytes |
コンパイル時間 | 1,045 ms |
コンパイル使用メモリ | 94,344 KB |
実行使用メモリ | 518,656 KB |
最終ジャッジ日時 | 2024-11-06 06:25:30 |
合計ジャッジ時間 | 18,654 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | MLE | - |
testcase_01 | MLE | - |
testcase_02 | MLE | - |
testcase_03 | MLE | - |
testcase_04 | MLE | - |
testcase_05 | MLE | - |
testcase_06 | MLE | - |
testcase_07 | MLE | - |
testcase_08 | MLE | - |
testcase_09 | MLE | - |
testcase_10 | MLE | - |
testcase_11 | MLE | - |
testcase_12 | MLE | - |
testcase_13 | MLE | - |
testcase_14 | MLE | - |
testcase_15 | MLE | - |
testcase_16 | MLE | - |
ソースコード
#include <iostream> #include <vector> #include <cmath> #include <iomanip> using namespace std; vector<vector<vector<vector<vector<double>>>>> mem; vector<double> f(int n, int a, int b, int c){ if(n == 0 || a+b+c <= 1) return {0.0, 0.0, 0.0}; if(mem[n][a][b][c][0] >= 0.0) return mem[n][a][b][c]; mem[n][a][b][c] = {0.0, 0.0, 0.0}; double z = (a+b+c)*(a+b+c-1); double y = z; y -= a*(a-1) + b*(b-1) + c*(c-1); vector<int> v = {a, b, c}; for(int i = 0; i < 3; i++){ if(v[i] < 2) continue; v[i]--; auto res = f(n-1, v[0], v[1], v[2]); res[i] += 1.0; for(int j = 0; j < 3; j++) mem[n][a][b][c][j] += res[j]*(v[i]+1)*v[i]/z; v[i]++; } auto res = f(n-1, a, b, c); for(int i = 0; i < 3; i++) mem[n][a][b][c][i] += res[i]*y/z; return mem[n][a][b][c]; } int main(){ int a, b, c, n; cin >> a >> b >> c >> n; mem = vector<vector<vector<vector<vector<double>>>>>(55, vector<vector<vector<vector<double>>>>(55, vector<vector<vector<double>>>(55, vector<vector<double>>(55, vector<double>(3, -1))))); vector<double> ans = f(n, a, b, c); for(int i = 0; i < 3; i++){ cout << fixed << setprecision(10) << ans[i] << " \n"[i == 2]; } return 0; }