結果

問題 No.144 エラトステネスのざる
ユーザー HachimoriHachimori
提出日時 2015-02-05 23:52:22
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 264 ms / 2,000 ms
コード長 1,195 bytes
コンパイル時間 1,316 ms
コンパイル使用メモリ 59,556 KB
実行使用メモリ 5,380 KB
最終ジャッジ日時 2023-09-05 13:15:33
合計ジャッジ時間 3,157 ms
ジャッジサーバーID
（参考情報） 		judge13 / judge15
このコードへのチャレンジ（β）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 7 ms
4,976 KB
testcase_01 AC 6 ms
5,112 KB
testcase_02 AC 6 ms
5,020 KB
testcase_03 AC 7 ms
5,188 KB
testcase_04 AC 6 ms
5,108 KB
testcase_05 AC 7 ms
5,020 KB
testcase_06 AC 7 ms
5,160 KB
testcase_07 AC 6 ms
5,280 KB
testcase_08 AC 7 ms
5,344 KB
testcase_09 AC 6 ms
5,332 KB
testcase_10 AC 7 ms
5,188 KB
testcase_11 AC 6 ms
5,380 KB
testcase_12 AC 7 ms
5,160 KB
testcase_13 AC 238 ms
5,168 KB
testcase_14 AC 261 ms
5,204 KB
testcase_15 AC 261 ms
5,356 KB
testcase_16 AC 263 ms
5,180 KB
testcase_17 AC 263 ms
5,280 KB
testcase_18 AC 264 ms
5,376 KB
testcase_19 AC 225 ms
5,092 KB
権限があれば一括ダウンロードができます

ソースコード

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