結果

問題 No.144 エラトステネスのざる
ユーザー toriaezuACtoriaezuAC
提出日時 2018-04-04 00:13:07
言語 C++11
(gcc 11.4.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
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 15 ms
5,248 KB
testcase_01 AC 15 ms
5,376 KB
testcase_02 AC 16 ms
5,376 KB
testcase_03 AC 15 ms
5,376 KB
testcase_04 AC 16 ms
5,376 KB
testcase_05 AC 15 ms
5,376 KB
testcase_06 AC 15 ms
5,376 KB
testcase_07 AC 14 ms
5,376 KB
testcase_08 AC 16 ms
5,376 KB
testcase_09 AC 15 ms
5,376 KB
testcase_10 AC 15 ms
5,376 KB
testcase_11 AC 15 ms
5,376 KB
testcase_12 AC 15 ms
5,376 KB
testcase_13 AC 16 ms
5,376 KB
testcase_14 AC 15 ms
5,376 KB
testcase_15 AC 16 ms
5,376 KB
testcase_16 AC 16 ms
5,376 KB
testcase_17 AC 16 ms
5,376 KB
testcase_18 AC 15 ms
5,376 KB
testcase_19 AC 16 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

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