結果

問題 No.144 エラトステネスのざる
ユーザー Kyutatsu
提出日時 2025-10-15 23:17:48
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 84 ms / 2,000 ms
コード長 445 bytes
コンパイル時間 989 ms
コンパイル使用メモリ 100,520 KB
実行使用メモリ 15,236 KB
最終ジャッジ日時 2025-10-15 23:17:51
合計ジャッジ時間 2,195 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<vector>
#include<iomanip>

using namespace std;

int main() {
    int N; cin >> N;
    double p; cin >> p;
    vector<double> P(N+1, 1);
    for (int i=1;i<=N;i++) P[i] = P[i-1] * (1.0 - p);

    double E = 0;
    vector<int> cnt(N+1);
    for (int i=2;i<=N;i++) {
        E += P[cnt[i]];
        for (int j=i+i;j<=N;j+=i) {
            cnt[j]++;
        }
    }
    cout << fixed << setprecision(12) << E << endl;
}
0