#661827 (C++17) No.108 トリプルカードコンプ

• 問題: No.108 トリプルカードコンプ
• ユーザー: kusainu7
• 提出日時: 2021-05-29 19:08:59
• 言語: C++17 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 18 ms / 5,000 ms
• コード長: 705 bytes
• コンパイル時間: 616 ms
• コンパイル使用メモリ: 68,864 KB
• 実行使用メモリ: 12,800 KB
• 最終ジャッジ日時: 2024-11-08 01:58:07
• 合計ジャッジ時間: 1,610 ms
• ジャッジサーバーID （参考情報）: judge4 / judge3

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,248 KB
testcase_02	AC	2 ms; 5,248 KB
testcase_03	AC	2 ms; 5,248 KB
testcase_04	AC	2 ms; 5,248 KB
testcase_05	AC	2 ms; 5,248 KB
testcase_06	AC	2 ms; 5,248 KB
testcase_07	AC	18 ms; 12,800 KB
testcase_08	AC	18 ms; 12,800 KB
testcase_09	AC	18 ms; 12,800 KB
testcase_10	AC	18 ms; 12,800 KB
testcase_11	AC	17 ms; 12,800 KB
testcase_12	AC	2 ms; 5,248 KB
testcase_13	AC	4 ms; 5,760 KB
testcase_14	AC	4 ms; 5,632 KB
testcase_15	AC	4 ms; 5,632 KB
testcase_16	AC	4 ms; 5,888 KB
testcase_17	AC	3 ms; 5,504 KB
testcase_18	AC	16 ms; 12,160 KB
testcase_19	AC	16 ms; 11,904 KB
testcase_20	AC	18 ms; 12,672 KB
testcase_21	AC	18 ms; 12,544 KB
testcase_22	AC	15 ms; 12,032 KB

ソースコード

#include<iostream>
using namespace std;

double c[110][110][110] = {};

void dp(int n, int *a) {
	c[n][n][n] = 0;
	for (int i=n; i>=0; i--) {
		for (int j=n; j>=0; j--) {
			for (int k=n; k>=0; k--) {
				if (k>=n) {
					continue;
				}
				c[i][j][k] = 1.0*n/(n-k);
				if (i!=n) c[i][j][k] += c[i+1][j][k]*(n-i)/(n-k);
				if (i!=j) c[i][j][k] += c[i][j+1][k]*(i-j)/(n-k);
				if (j!=k) c[i][j][k] += c[i][j][k+1]*(j-k)/(n-k);
			}
		}
	}
}

int main() {
	int n, x;
	cin >> n;
	int a[4] = {};

	for (int i=0; i<n; i++) {
		cin >> x;
		if (x>3) {
			a[3]++;
		} else if (x!=0) {
			a[x]++;
		}
	}
	
	for (int i=2; i>0; i--) {
		a[i] += a[i+1];
	}
	
	dp(n, a);

	printf("%.10f\n", c[a[1]][a[2]][a[3]]);
}