結果

問題 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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0