結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2015-02-06 00:44:06 |
言語 | D (dmd 2.109.1) |
結果 |
AC
|
実行時間 | 49 ms / 2,000 ms |
コード長 | 587 bytes |
コンパイル時間 | 2,693 ms |
コンパイル使用メモリ | 151,296 KB |
実行使用メモリ | 10,240 KB |
最終ジャッジ日時 | 2024-06-12 02:10:20 |
合計ジャッジ時間 | 4,062 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
import std.stdio, std.conv, std.math, std.string, std.range, std.array,std.algorithm;double[] is_prime;void Eratosthenes(long lim, double q){is_prime[] = 1.0;is_prime[0] = 0.0;is_prime[1] = 0.0;for(long p=2; p<=lim; p++) {for(long k=2*p; k<=lim; k+=p) is_prime[k] *= q;}}void main(){auto buf = readln().strip().split();immutable N = buf[0].to!int;immutable p = buf[1].to!double;is_prime.length = N + 5;Eratosthenes(N, 1-p);double ans = 0;foreach(immutable i; 2 .. N+1) {ans += is_prime[i];}writefln("%.10f", ans);}