結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー toriaezuAC
提出日時 2018-04-04 00:13:07
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 16 ms / 2,000 ms
コード長 572 bytes
コンパイル時間 508 ms
コンパイル使用メモリ 64,440 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-06-26 08:36:40
合計ジャッジ時間 1,538 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <algorithm>
#include <iomanip>
using namespace std;

const int MAX = 1000000;
const int FACTORS_MAX = 240;
int N;
double P;
unsigned char DP[MAX + 1];
int C[FACTORS_MAX];

void Calc() {
	for (int p = 2; p <= MAX; ++p) {
		for (int i = p * 2; i <= MAX; i += p) {
			++DP[i];
		}
	}
}

int main() {
	cin >> N >> P;
	Calc();

	for (int i = 2; i <= N; ++i) {
		++C[DP[i]];
	}

	double sum = 0;
	double prob = 1;
	for (int k = 0; k < FACTORS_MAX; ++k) {
		sum += prob * C[k];
		prob *= 1 - P;
	}
	cout << fixed << setprecision(16) << sum << endl;
}
0