結果

問題 No.144 エラトステネスのざる
ユーザー asugen0402
提出日時 2019-03-06 17:25:43
言語 C
(gcc 13.3.0)
結果
RE  
実行時間 -
コード長 2,043 bytes
コンパイル時間 416 ms
コンパイル使用メモリ 31,744 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2024-06-23 14:38:10
合計ジャッジ時間 2,292 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 RE * 7
権限があれば一括ダウンロードができます

ソースコード

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		100005									// 最大値

// 内部変数
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