結果

問題 No.144 エラトステネスのざる
ユーザー asugen0402asugen0402
提出日時 2019-03-06 17:28:46
言語 C
(gcc 12.3.0)
結果
AC  
実行時間 49 ms / 2,000 ms
コード長 2,044 bytes
コンパイル時間 122 ms
コンパイル使用メモリ 29,696 KB
実行使用メモリ 10,112 KB
最終ジャッジ日時 2023-09-05 19:11:51
合計ジャッジ時間 1,321 ms
ジャッジサーバーID
(参考情報)
judge14 / judge15
このコードへのチャレンジ
(要ログイン)

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 1 ms
4,376 KB
testcase_01 AC 0 ms
4,376 KB
testcase_02 AC 1 ms
4,376 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 0 ms
4,380 KB
testcase_05 AC 1 ms
4,380 KB
testcase_06 AC 1 ms
4,380 KB
testcase_07 AC 1 ms
4,380 KB
testcase_08 AC 0 ms
4,376 KB
testcase_09 AC 1 ms
4,376 KB
testcase_10 AC 1 ms
4,380 KB
testcase_11 AC 0 ms
4,376 KB
testcase_12 AC 1 ms
4,376 KB
testcase_13 AC 46 ms
9,976 KB
testcase_14 AC 46 ms
10,024 KB
testcase_15 AC 45 ms
10,112 KB
testcase_16 AC 46 ms
10,068 KB
testcase_17 AC 49 ms
10,012 KB
testcase_18 AC 48 ms
10,020 KB
testcase_19 AC 45 ms
10,100 KB
権限があれば一括ダウンロードができます

ソースコード

diff #

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

// 内部定数
#define D_VAL_MAX		2000005									// 最大値

// 内部変数
static FILE *szpFpI;											// 入力
static int siN;													// N
static double sdRate;											// 残る確率
static double sd1Rate[D_VAL_MAX];								// 残る確率[値]

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

// 実行メイン
int
fMain(
	int piTNo					// <I> テスト番号 1~
)
{
	int i, j;
	char lc1Buf[1024], lc1Out[1024];

	// 入力 - セット
#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

	// N・残る確率 - 取得
	fgets(lc1Buf, sizeof(lc1Buf), szpFpI);
	sscanf(lc1Buf, "%d%lf", &siN, &sdRate);
	sdRate = 1.0 - sdRate;

	// 残る確率[値] - 初期化
	for (i = 2; i <= siN; i++) {
		sd1Rate[i] = 1.0;
	}

	// 残る確率[値] - 更新
	for (i = 2; i <= siN; i++) {
		for (j = i * 2; j <= siN; j += i) {
			sd1Rate[j] *= sdRate;
		}
	}

	// 個数 - 取得
	double ldCnt = 0.0;
	for (i = 2; i <= siN; i++) {
		ldCnt += sd1Rate[i];
	}

	// 結果 - セット
	sprintf(lc1Out, "%.7lf\n", ldCnt);

	// 結果 - 表示
#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;
}

0