結果
問題 |
No.108 トリプルカードコンプ
|
ユーザー |
![]() |
提出日時 | 2018-12-27 20:59:54 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 22 ms / 5,000 ms |
コード長 | 882 bytes |
コンパイル時間 | 2,254 ms |
コンパイル使用メモリ | 200,896 KB |
最終ジャッジ日時 | 2025-01-06 20:05:05 |
ジャッジサーバーID (参考情報) |
judge1 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 20 |
ソースコード
#include <bits/stdc++.h> using namespace std; using ld = long double; int main() { cout << fixed << setprecision(8); int N; cin >> N; int X = 0, Y = 0, Z = 0; for (int i = 0; i < N; i++) { int A; cin >> A; if (A == 0) ++X; if (A == 1) ++Y; if (A == 2) ++Z; } vector<vector<vector<ld>>> dp(N + 1, vector<vector<ld>>(N + 1, vector<ld>(N + 1))); for (int i = 0; i <= N; i++) { for (int j = 0; j <= N; j++) { for (int k = 0; k <= N; k++) if (i + j + k <= N) { if (0 <= i - 1) { dp[i][j][k] += (dp[i - 1][j + 1][k] + (ld)N / (i + j + k)) * ((ld)i / (i + j + k)); } if (0 <= j - 1) { dp[i][j][k] += (dp[i][j - 1][k + 1] + (ld)N / (i + j + k)) * ((ld)j / (i + j + k)); } if (0 <= k - 1) { dp[i][j][k] += (dp[i][j][k - 1] + (ld)N / (i + j + k)) * ((ld)k / (i + j + k)); } } } } cout << dp[X][Y][Z] << endl; return 0; }