#include #include #include #include #include #include // 内部定数 #define D_DICE_MAX 10 // 最大サイコロ数 #define D_VAL_MAX (6 * D_DICE_MAX) // 最大値 // 内部変数 static FILE *szpFpI; // 入力 static int siDCnt; // サイコロ数 static int siCCnt; // イカサマ数 static double sd1P1[D_VAL_MAX + 5]; // P1確率[値] static double sd1P2[D_VAL_MAX + 5]; // P2確率[値] // 内部変数 - テスト用 #ifdef D_TEST static int siRes; static FILE *szpFpA; #endif // 確率 - セット - 1回 int fSetRateOne( double *pdpRate // 確率 , int *pipRngSVal // 範囲 - 開始 , int *pipRngEVal // 範囲 - 終了 , int piAddMin // 加算 - 最小 , int piAddCnt // 加算 - 回数 ) { int i, j; // ワーク - 初期化 double ld1Rate[D_VAL_MAX + 5]; memset(ld1Rate, 0, sizeof(ld1Rate)); // 乗算確率 double ldMulti = 1.0 / (double)piAddCnt; // 範囲でループ for (i = *pipRngSVal; i <= *pipRngEVal; i++) { // 加算回数でループ for (j = 0; j < piAddCnt; j++) { ld1Rate[i + piAddMin + j] += pdpRate[i] * ldMulti; } } // 確率 - コピー memcpy(pdpRate, ld1Rate, sizeof(ld1Rate)); // 範囲 - 更新 *pipRngSVal += piAddMin; *pipRngEVal += piAddMin + piAddCnt - 1; return 0; } // 確率 - セット int fSetRate( double *pdpRate // 確率 , int piDCnt // サイコロ数 , int piCCnt // イカサマ数 ) { int i; // 確率 - 初期値 pdpRate[0] = 1.0; int liSVal = 0; int liEVal = 0; // サイコロ数 for (i = 0; i < piDCnt; i++) { fSetRateOne(pdpRate, &liSVal, &liEVal, 1, 6); } // イカサマ数 for (i = 0; i < piCCnt; i++) { fSetRateOne(pdpRate, &liSVal, &liEVal, 4, 3); } return 0; } // 実行メイン int fMain( int piTNo // テスト番号 1~ ) { int i; char lc1Buf[1024], lc1Out[1024]; // データ - 初期化 memset(sd1P1, 0, sizeof(sd1P1)); // P1確率[値] memset(sd1P2, 0, sizeof(sd1P2)); // P2確率[値] // 入力 - セット #ifdef D_TEST sprintf(lc1Buf, ".\\Test\\T%d.txt", piTNo); szpFpI = fopen(lc1Buf, "r"); sprintf(lc1Buf, ".\\Test\\A%d.txt", piTNo); szpFpA = fopen(lc1Buf, "r"); siRes = 0; #else szpFpI = stdin; #endif // サイコロ数 - 取得 fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d", &siDCnt); // イカサマ数 - 取得 fgets(lc1Buf, sizeof(lc1Buf), szpFpI); sscanf(lc1Buf, "%d", &siCCnt); // 確率 - セット fSetRate(sd1P1, siDCnt - siCCnt, siCCnt); fSetRate(sd1P2, siDCnt, 0); // P2 - 累積和 for (i = 1; i <= siDCnt * 6; i++) { sd1P2[i] += sd1P2[i - 1]; } // 勝率 - 取得 double ldWin = 0.0; for (i = 2; i <= siDCnt * 6; i++) { ldWin += sd1P1[i] * sd1P2[i - 1]; } // 結果 - セット sprintf(lc1Out, "%.5lf\n", ldWin); // 結果 - 表示 #ifdef D_TEST fgets(lc1Buf, sizeof(lc1Buf), szpFpA); if (strcmp(lc1Buf, lc1Out)) { siRes = -1; } #else printf("%s", lc1Out); #endif // 残データ有無 #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", piTNo); } else { printf("NG %d\n", piTNo); } #endif return 0; } int main() { #ifdef D_TEST int i; for (i = D_TEST_SNO; i <= D_TEST_ENO; i++) { fMain(i); } #else fMain(0); #endif return 0; }