結果

問題 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();
}
0