結果

問題 No.144 エラトステネスのざる
ユーザー 🍡yurahuna
提出日時 2016-02-28 21:41:38
言語 C++11
(gcc 13.3.0)
結果
AC  
実行時間 53 ms / 2,000 ms
コード長 929 bytes
コンパイル時間 719 ms
コンパイル使用メモリ 74,072 KB
実行使用メモリ 7,936 KB
最終ジャッジ日時 2024-09-24 12:09:08
合計ジャッジ時間 1,814 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff #
プレゼンテーションモードにする

#include <cstdio>
#include <iostream>
#include <algorithm>
#include <complex>
#include <queue>
#include <map>
#include <string>
using namespace std;
#define FOR(i,a,b) for (int i=(a);i<(b);i++)
#define FORR(i,a,b) for (int i=(b)-1;i>=(a);i--)
#define REP(i,n) for (int i=0;i<(n);i++)
#define RREP(i,n) for (int i=(n)-1;i>=0;i--)
#define pb push_back
#define ALL(a) (a).begin(),(a).end()
#define PI 3.1415926535
typedef long long ll;
typedef pair<int, int> P;
//typedef complex<double> C;
const int MAX_N = 1000000;
int N;
double p;
int num_div[MAX_N + 1]; // 1と自身を除く約数の個数
void input() {
    cin >> N >> p;
}
void solve() {
    FOR(i, 2, N + 1) {
        for (int j = i * 2; j <= N; j += i) {
            num_div[j]++;
        }
    }
    double sum = 0;
    FOR(i, 2, N + 1) {
        if (num_div[i] == 0) {
            sum += 1;
        } else {
            sum += pow(1 - p , num_div[i]);
        }
    }
    printf("%.10f\n", sum);
}
int main() {
    input();
    solve();
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0