結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー Hachimori
提出日時 2015-02-05 23:52:22
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 226 ms / 2,000 ms
コード長 1,195 bytes
コンパイル時間 487 ms
コンパイル使用メモリ 58,848 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-23 08:54:06
合計ジャッジ時間 2,691 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<iostream>
#include<cstdio>
#include<cmath>
using namespace std;
const int BUF = 1000005;


int nPrime;
int prime[BUF];

void makePrime() {
    bool isPrime[BUF];
    for (int i = 0; i < BUF; ++i)
        isPrime[i] = true;

    for (int i = 2; i * i < BUF; ++i)
        if (isPrime[i])
            for (int j = i * 2; j < BUF; j += i)
                isPrime[j] = false;

    nPrime = 0;
    for (int i = 2; i < BUF; ++i)
        if (isPrime[i])
            prime[nPrime++] = i;
}


int N;
double p;

void read() {
    cin >> N >> p;
}


int calcNumDivisor(int n) {
    int nDivisor = 1;
    for (int i = 0; prime[i] * prime[i] <= n; ++i) {
        if (n % prime[i] == 0) {
            int cnt = 0;
            while (n % prime[i] == 0) {
                n /= prime[i];
                ++cnt;
            }
            nDivisor *= cnt + 1;
        }
    }
    if (n > 1)
        nDivisor *= 2;
    return nDivisor;
}


void work() {
    double sum = 0;

    for (int i = 2; i <= N; ++i) {
        int nDivisor = calcNumDivisor(i);
        sum += pow(1 - p, nDivisor - 2);
    }
    
    printf("%.10lf\n", sum);
}


int main() {
    makePrime();
    read();
    work();
    return 0;
}
0