結果

問題 No.108 トリプルカードコンプ
コンテスト
ユーザー kazuma
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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;
}
0