#699131 (C) No.1688 Veterinarian

• 問題: No.1688 Veterinarian
• ユーザー: ygussany
• 提出日時: 2021-09-12 14:25:04
• 言語: C (gcc 12.3.0)
• 結果: AC
• 実行時間: 120 ms / 3,000 ms
• コード長: 1,142 bytes
• コンパイル時間: 167 ms
• コンパイル使用メモリ: 32,384 KB
• 実行使用メモリ: 160,384 KB
• 最終ジャッジ日時: 2024-07-05 09:52:31
• 合計ジャッジ時間: 1,155 ms
• ジャッジサーバーID （参考情報）: judge1 / judge3

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1 ms; 5,248 KB
testcase_01	AC	1 ms; 5,376 KB
testcase_02	AC	1 ms; 5,376 KB
testcase_03	AC	1 ms; 5,376 KB
testcase_04	AC	1 ms; 5,376 KB
testcase_05	AC	2 ms; 5,376 KB
testcase_06	AC	2 ms; 5,376 KB
testcase_07	AC	70 ms; 160,364 KB
testcase_08	AC	116 ms; 160,384 KB
testcase_09	AC	120 ms; 151,424 KB
testcase_10	AC	17 ms; 22,016 KB
testcase_11	AC	11 ms; 35,968 KB
testcase_12	AC	14 ms; 20,352 KB
testcase_13	AC	4 ms; 5,376 KB
testcase_14	AC	9 ms; 16,896 KB
testcase_15	AC	16 ms; 11,904 KB
testcase_16	AC	3 ms; 5,760 KB

ソースコード

#include <stdio.h>

double memo[51][51][51][51][3];

double recursion(int A, int B, int C, int N, int k)
{
	if (N == 0) return 0.0;
	else if (memo[A][B][C][N][k] >= 0.0) return memo[A][B][C][N][k];
	
	int all = (A + B + C) * (A + B + C - 1) / 2, sum = A * (A - 1) / 2 + B * (B - 1) / 2 + C * (C - 1) / 2;
	memo[A][B][C][N][k] = recursion(A, B, C, N - 1, k) * (all - sum) / all;
	if (A >= 2) memo[A][B][C][N][k] += (recursion(A - 1, B, C, N - 1, k) + ((k == 0)? 1: 0)) * (A * (A - 1) / 2) / all;
	if (B >= 2) memo[A][B][C][N][k] += (recursion(A, B - 1, C, N - 1, k) + ((k == 1)? 1: 0)) * (B * (B - 1) / 2) / all;
	if (C >= 2) memo[A][B][C][N][k] += (recursion(A, B, C - 1, N - 1, k) + ((k == 2)? 1: 0)) * (C * (C - 1) / 2) / all;
	return memo[A][B][C][N][k];
}

int main()
{
	int A, B, C, N;
	scanf("%d %d %d %d", &A, &B, &C, &N);
	
	int a, b, c, n, k;
	for (a = 0; a <= A; a++) for (b = 0; b <= B; b++) for (c = 0; c <= C; c++) for (n = 0; n <= N; n++) for (k = 0; k <= 2; k++) memo[a][b][c][n][k] = -1.0;
	printf("%.7f %.7f %.7f\n", recursion(A, B, C, N, 0), recursion(A, B, C, N, 1), recursion(A, B, C, N, 2));
	fflush(stdout);
	return 0;
}