結果
問題 | No.108 トリプルカードコンプ |
ユーザー | Kyutatsu |
提出日時 | 2025-01-01 04:38:50 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 15 ms / 5,000 ms |
コード長 | 1,270 bytes |
コンパイル時間 | 7,194 ms |
コンパイル使用メモリ | 283,368 KB |
実行使用メモリ | 11,904 KB |
最終ジャッジ日時 | 2025-01-01 04:38:59 |
合計ジャッジ時間 | 5,903 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
6,816 KB |
testcase_01 | AC | 2 ms
6,820 KB |
testcase_02 | AC | 2 ms
6,820 KB |
testcase_03 | AC | 2 ms
6,820 KB |
testcase_04 | AC | 1 ms
6,816 KB |
testcase_05 | AC | 2 ms
6,816 KB |
testcase_06 | AC | 2 ms
6,820 KB |
testcase_07 | AC | 15 ms
11,904 KB |
testcase_08 | AC | 9 ms
11,904 KB |
testcase_09 | AC | 9 ms
11,904 KB |
testcase_10 | AC | 9 ms
11,904 KB |
testcase_11 | AC | 10 ms
11,904 KB |
testcase_12 | AC | 2 ms
6,820 KB |
testcase_13 | AC | 3 ms
6,816 KB |
testcase_14 | AC | 3 ms
6,820 KB |
testcase_15 | AC | 3 ms
6,816 KB |
testcase_16 | AC | 3 ms
6,820 KB |
testcase_17 | AC | 3 ms
6,820 KB |
testcase_18 | AC | 12 ms
11,008 KB |
testcase_19 | AC | 10 ms
10,496 KB |
testcase_20 | AC | 10 ms
11,776 KB |
testcase_21 | AC | 12 ms
11,776 KB |
testcase_22 | AC | 8 ms
10,880 KB |
ソースコード
#include <bits/stdc++.h> #include <iomanip> #include <functional> #include <sstream> #include <bit> using namespace std; int main() { int N; cin >> N; vector<int> A(4); int a; vector DP(N+1, vector(N+1, vector<double>(N+1, -1))); DP[0][0][0] = 0; for (int i=0;i<N;i++) { cin >> a; A[3 - min(3, a)]++; } for (int i=3;1<=i;i--) A[i-1] += A[i]; function<double(int,int,int)> dfs = [&](int a1, int a2, int a3) { if (DP[a1][a2][a3] != -1) return DP[a1][a2][a3]; double& dp = DP[a1][a2][a3]; dp = 1.0; // 1パック買った! // 残り0枚必要なカードを引いた double p0 = (double)(N - a1) / N; // 残り1枚必要なカードを引いた double p1 = (double)(a1 - a2) / N; if (0<a1-a2) dp += dfs(a1-1, a2, a3) * p1; // 残り2枚必要なカードを引いた double p2 = (double)(a2 - a3) / N; if (0<a2-a3) dp += dfs(a1, a2-1, a3) * p2; // 残り3枚必要なカードを引いた double p3 = (double)a3 / N; if (0<a3) dp += dfs(a1, a2, a3-1) * p3; dp /= (1.0 - p0); return dp; }; cout << fixed << setprecision(12) << dfs(A[1], A[2], A[3]) << endl; }