結果
問題 | No.1688 Veterinarian |
ユーザー | Hideaki Ito |
提出日時 | 2021-09-24 21:54:53 |
言語 | C++17 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 177 ms / 3,000 ms |
コード長 | 1,426 bytes |
コンパイル時間 | 2,003 ms |
コンパイル使用メモリ | 199,700 KB |
実行使用メモリ | 320,840 KB |
最終ジャッジ日時 | 2023-09-18 21:15:52 |
合計ジャッジ時間 | 4,534 ms |
ジャッジサーバーID (参考情報) |
judge14 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 71 ms
320,720 KB |
testcase_01 | AC | 118 ms
320,704 KB |
testcase_02 | AC | 82 ms
320,668 KB |
testcase_03 | AC | 71 ms
320,668 KB |
testcase_04 | AC | 71 ms
320,616 KB |
testcase_05 | AC | 71 ms
320,684 KB |
testcase_06 | AC | 117 ms
320,732 KB |
testcase_07 | AC | 72 ms
320,668 KB |
testcase_08 | AC | 177 ms
320,740 KB |
testcase_09 | AC | 171 ms
320,616 KB |
testcase_10 | AC | 111 ms
320,672 KB |
testcase_11 | AC | 95 ms
320,616 KB |
testcase_12 | AC | 72 ms
320,764 KB |
testcase_13 | AC | 98 ms
320,840 KB |
testcase_14 | AC | 76 ms
320,736 KB |
testcase_15 | AC | 130 ms
320,680 KB |
testcase_16 | AC | 77 ms
320,660 KB |
ソースコード
#include <bits/stdc++.h> using namespace std; #define ll long long #define ld long double #define REP(i,m,n) for(int i=(int)(m); i<(int)(n); i++) #define rep(i,n) REP(i,0,n) #define RREP(i,m,n) for(int i=(int)(m); i>=(int)(n); i--) #define rrep(i,n) RREP(i,(n)-1,0) #define all(v) v.begin(), v.end() #define endk '\n' const int inf = 1e9+7; const ll longinf = 1LL<<60; const ll mod = 1e9+7; const ll mod2 = 998244353; const ld eps = 1e-10; template<typename T1, typename T2> inline void chmin(T1 &a, T2 b){if(a>b) a=b;} template<typename T1, typename T2> inline void chmax(T1 &a, T2 b){if(a<b) a=b;} ld dp[3][51][51][51][51]; int main() { cin.tie(0); ios::sync_with_stdio(false); int a, b, c, n; cin >> a >> b >> c >> n; auto dfs = [&](auto dfs, int a, int b, int c, int n, int i) -> ld { if(dp[i][a][b][c][n] != -1) return dp[i][a][b][c][n]; if(a == 0 || b == 0 || c == 0 || n == 0) return dp[i][a][b][c][n] = 0; return dp[i][a][b][c][n] = (ld)(dfs(dfs, a, b, c, n-1, i) * (a*b+b*c+c*a) + (dfs(dfs, a-1, b, c, n-1, i)+(i==0)) * a*(a-1)/2 + (dfs(dfs, a, b-1, c, n-1, i)+(i==1)) * b*(b-1)/2 + (dfs(dfs, a, b, c-1, n-1, i)+(i==2)) * c*(c-1)/2) / ((a+b+c)*(a+b+c-1)/2); }; rep(i, 3) { rep(a, 51) rep(b, 51) rep(c, 51) rep(j, n+1) dp[i][a][b][c][j] = -1; cout << fixed << setprecision(10) << dfs(dfs, a, b, c, n, i) << ' '; } cout << endk; return 0; }