結果

問題 No.574 正多面体サイコロ
コンテスト
ユーザー asugen0402
提出日時 2019-06-14 16:00:34
言語 C
(gcc 13.3.0)
結果
AC  
実行時間 39 ms / 2,000 ms
コード長 2,743 bytes
コンパイル時間 607 ms
コンパイル使用メモリ 32,128 KB
実行使用メモリ 5,248 KB
最終ジャッジ日時 2024-11-07 23:03:37
合計ジャッジ時間 1,655 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 20
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <float.h>
#include <limits.h>
#include <math.h>
#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <time.h>

// 内部変数
static FILE *szpFpI;											// 入力

// 内部変数 - テスト用
#ifdef D_TEST
	static int siRes;
	static FILE *szpFpA;
	static int siTNo;
#endif

// 出力
int
fOut(
	char *pcpLine				// <I> １行
)
{
	char lc1Buf[1024];

#ifdef D_TEST
	fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
	if (strcmp(lc1Buf, pcpLine)) {
		siRes = -1;
	}
#else
	printf("%s", pcpLine);
#endif

	return 0;
}

// 実行メイン
double
fMain(
)
{
	int i, j, k, l;
	char lc1Buf[1024];

	// F, N, K - 取得
	int liF, liN, liK;
	fgets(lc1Buf, sizeof(lc1Buf), szpFpI);
	sscanf(lc1Buf, "%d%d%d", &liF, &liN, &liK);

	// 各目のパターン数でループ
	double ldSumV = 0.0;
	for (i = 1; i <= liF; i++) {
		double ldSumR = 0.0;

		// 各目の範囲でループ
		for (j = 1; j <= liK; j++) {
			for (k = liK; k <= liN; k++) {
				double ldRate = 1.0;

				// パターン数
				for (l = 2; l <= liN; l++) {
					ldRate *= l;
				}

				// 各目より大
				int liLen = j - 1;
				if (liLen > 0) {
					if (i == liF) {
						continue;
					}
					for (l = 0; l < liLen; l++) {
						ldRate *= (double)(liF - i) / (double)liF;
						ldRate /= (l + 1);
					}
				}

				// 各目
				liLen = k - j + 1;
				for (l = 0; l < liLen; l++) {
					ldRate *= 1.0 / (double)liF;
					ldRate /= (l + 1);
				}

				// 各目より小
				liLen = liN - k;
				if (liLen > 0) {
					if (i == 1) {
						continue;
					}
					for (l = 0; l < liLen; l++) {
						ldRate *= (double)(i - 1) / (double)liF;
						ldRate /= (l + 1);
					}
				}

				// 確率合計
				ldSumR += ldRate;
			}
		}

		// 期待値合計
		ldSumV += ldSumR * (double)i;
	}

	return ldSumV;
}

// １回実行
int
fOne(
)
{
	double ldRet;
	char lc1Buf[1024];

	// 入力 - セット
#ifdef D_TEST
	sprintf(lc1Buf, ".\\Test\\T%d.txt", siTNo);
	szpFpI = fopen(lc1Buf, "r");
	sprintf(lc1Buf, ".\\Test\\A%d.txt", siTNo);
	szpFpA = fopen(lc1Buf, "r");
	siRes = 0;
#else
	szpFpI = stdin;
#endif

	// 実行メイン
	ldRet = fMain();

	// 出力
	sprintf(lc1Buf, "%.7lf\n", ldRet);
	fOut(lc1Buf);

	// 残データ有無
#ifdef D_TEST
	lc1Buf[0] = '\0';
	fgets(lc1Buf, sizeof(lc1Buf), szpFpA);
	if (strcmp(lc1Buf, "")) {
		siRes = -1;
	}
#endif

	// テストファイルクローズ
#ifdef D_TEST
	fclose(szpFpI);
	fclose(szpFpA);
#endif

	// テスト結果
#ifdef D_TEST
	if (siRes == 0) {
		printf("OK %d\n", siTNo);
	}
	else {
		printf("NG %d\n", siTNo);
	}
#endif

	return 0;
}

// プログラム開始
int
main()
{

#ifdef D_TEST
	int i;
	for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) {
		siTNo = i;
		fOne();
	}
#else
	fOne();
#endif

	return 0;
}
0