結果
問題 | No.108 トリプルカードコンプ |
ユーザー |
|
提出日時 | 2015-08-21 11:08:26 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 11 ms / 5,000 ms |
コード長 | 1,115 bytes |
コンパイル時間 | 1,390 ms |
コンパイル使用メモリ | 159,632 KB |
実行使用メモリ | 11,776 KB |
最終ジャッジ日時 | 2024-07-18 11:12:36 |
合計ジャッジ時間 | 2,158 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#include <bits/stdc++.h> #define rep(i, a) for (int i = 0; i < (a); i++) #define rep2(i, a, b) for (int i = (a); i < (b); i++) #define repr(i, a) for (int i = (a) - 1; i >= 0; i--) #define repr2(i, a, b) for (int i = (b) - 1; i >= (a); i--) using namespace std; typedef long long ll; const ll inf = 1e9; const ll mod = 1e9 + 7; int N; double memo[101][101][101]; // num 0, num 1, num 2 double E(int n0, int n1, int n2) { if (n0 == 0 && n1 == 0 && n2 == 0) return 0; if (memo[n0][n1][n2] >= 0) return memo[n0][n1][n2]; int n3 = N - n0 - n1 - n2; double res = (double)n3 / N; if (n0 > 0) res += (E(n0 - 1, n1 + 1, n2) + 1) * ((double)n0 / N); if (n1 > 0) res += (E(n0, n1 - 1, n2 + 1) + 1) * ((double)n1 / N); if (n2 > 0) res += (E(n0, n1, n2 - 1) + 1) * ((double)n2 / N); res /= 1 - (double)n3 / N; return memo[n0][n1][n2] = res; } int main() { cin >> N; vector<int> A(N), num(4); rep (i, N) { cin >> A[i]; A[i] = min(3, A[i]); num[A[i]]++; } fill_n((double *)memo, sizeof(memo) / sizeof(memo[0][0][0]), -1); double ans = E(num[0], num[1], num[2]); printf("%.20f\n", ans); return 0; }