結果

問題 No.108 トリプルカードコンプ
コンテスト
ユーザー tkzw_21
提出日時 2014-12-23 16:53:19
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 8 ms / 5,000 ms
コード長 894 bytes
コンパイル時間 672 ms
コンパイル使用メモリ 74,964 KB
実行使用メモリ 11,892 KB
最終ジャッジ日時 2024-06-12 03:56:15
合計ジャッジ時間 1,533 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#include <iostream>
#include <cmath>
#include <vector>
#include <utility>
#include <queue>
#include <map>
#include <string>
#include <cstring>
#include <list>
#include <algorithm>
#define INF 114514

using namespace std;

typedef pair<double,double> P;
int n;

double memo[101][101][101];
double dfs(int n3,int n2,int n1){
	if(memo[n3][n2][n1] >= 0)return memo[n3][n2][n1];

	double ret = n*1.0/(n1+n2+n3);
	if(n1>0)ret += dfs(n3,n2,n1-1) * n1/(n1+n2+n3);
	if(n2>0)ret += dfs(n3,n2-1,n1+1) * n2/(n1+n2+n3);
	if(n3>0)ret += dfs(n3-1,n2+1,n1) * n3/(n1+n2+n3);
	return memo[n3][n2][n1] = ret;
}

int main(void){
	cin >> n;
	vector<int> a(n);
	for(int i=0;i<n;i++)
		cin >> a[i];
	
	int n3,n2,n1;
	n3=n2=n1=0;
	for(int i=0;i<n;i++){
		if(a[i] == 0)n3++;
		if(a[i] == 1)n2++;
		if(a[i] == 2)n1++;
	}
	memset(memo,-1,sizeof(memo));
	memo[0][0][0] = 0;
	printf("%.12lf\n",dfs(n3,n2,n1));

	return 0;
}
0