結果

問題 No.144 エラトステネスのざる
ユーザー startcppstartcpp
提出日時 2018-03-16 17:22:18
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 64 ms / 2,000 ms
コード長 569 bytes
コンパイル時間 268 ms
コンパイル使用メモリ 26,740 KB
実行使用メモリ 6,944 KB
最終ジャッジ日時 2023-08-23 11:50:59
合計ジャッジ時間 2,020 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,380 KB
testcase_01 AC 2 ms
4,380 KB
testcase_02 AC 2 ms
4,380 KB
testcase_03 AC 2 ms
4,376 KB
testcase_04 AC 2 ms
4,376 KB
testcase_05 AC 1 ms
4,376 KB
testcase_06 AC 1 ms
4,380 KB
testcase_07 AC 2 ms
4,380 KB
testcase_08 AC 1 ms
4,376 KB
testcase_09 AC 1 ms
4,376 KB
testcase_10 AC 2 ms
4,376 KB
testcase_11 AC 2 ms
4,380 KB
testcase_12 AC 2 ms
4,376 KB
testcase_13 AC 63 ms
6,844 KB
testcase_14 AC 62 ms
6,944 KB
testcase_15 AC 63 ms
6,848 KB
testcase_16 AC 61 ms
6,844 KB
testcase_17 AC 62 ms
6,836 KB
testcase_18 AC 64 ms
6,896 KB
testcase_19 AC 63 ms
6,904 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>

double power(double x, int n) {
	if (n == 0) return 1;
	if (n & 1) return x * power(x, n - 1);
	return power(x * x, n >> 1);
}

int ynum[1000001];

int main() {
	int n; double p;
	int i, j;
	double ans = 0;
	
	scanf("%d%lf", &n, &p);
	for (i = 2; i <= n; i++) ynum[i] = 1;
	for (i = 2; i <= n; i++) {
		if (ynum[i] == 1) {
			for (j = i; j <= n; j += i) {
				int t = j, cnt = 0;
				while (t % i == 0) {
					t /= i;
					cnt++;
				}
				ynum[j] *= (cnt + 1);
			}
		}
		ans += power(1.0 - p, ynum[i] - 2);
	}
	printf("%.14f\n", ans);
	return 0;
}
0