結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー startcpp
提出日時 2018-03-16 17:22:18
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 34 ms / 2,000 ms
コード長 569 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 356 ms
コンパイル使用メモリ 40,832 KB
実行使用メモリ 7,936 KB
最終ジャッジ日時 2026-05-28 06:30:44
合計ジャッジ時間 1,542 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_0
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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