結果

問題 No.108 トリプルカードコンプ
コンテスト
ユーザー はまやんはまやん
提出日時 2016-11-15 03:50:51
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 12 ms / 5,000 ms
コード長 895 bytes
コンパイル時間 1,982 ms
コンパイル使用メモリ 166,496 KB
実行使用メモリ 11,776 KB
最終ジャッジ日時 2024-11-26 01:37:31
合計ジャッジ時間 2,553 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include<bits/stdc++.h>
using namespace std;
#define rep(i,a,b) for(int i=a;i<b;i++)



int N;
int A[101];
//-----------------------------------------------------------------
double memo[101][101][101];
double rec(int a, int b, int c) {
	if (0 <= memo[a][b][c]) return memo[a][b][c];
	int sm = a + b + c;
	if (sm == 0) return 0;

	double ret = ((double)N / (double)sm);
	
	if (0 < a) ret += rec(a - 1, b + 1, c) * ((double)a / (double)sm);
	if (0 < b) ret += rec(a, b - 1, c + 1) * ((double)b / (double)sm);
	if (0 < c) ret += rec(a, b, c - 1) * ((double)c / (double)sm);

	return memo[a][b][c] = ret;
}

int main() {
	cin >> N;
	rep(i, 0, N) cin >> A[i];

	int a = 0, b = 0, c = 0;
	rep(i, 0, N) {
		if (A[i] == 0)
			a++;
		else if (A[i] == 1)
			b++;
		else if (A[i] == 2)
			c++;
	}

	rep(i, 0, N + 1) rep(j, 0, N + 1) rep(k, 0, N + 1) memo[i][j][k] = -1;
	printf("%.10f\n", rec(a, b, c));
}
0