#include #include #include #include #include #include #include // 内部変数 static FILE *szpFpI; // 入力 // 内部変数 - テスト用 #ifdef D_TEST static int siRes; static FILE *szpFpA; static int siTNo; #endif // 出力 int fOut( char *pcpLine // 1行 ) { 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; } // 1回実行 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; }