結果

問題 No.144 エラトステネスのざる
ユーザー toriaezuACtoriaezuAC
提出日時 2018-04-04 00:13:07
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 17 ms / 2,000 ms
コード長 572 bytes
コンパイル時間 711 ms
コンパイル使用メモリ 65,384 KB
実行使用メモリ 4,648 KB
最終ジャッジ日時 2023-09-08 15:35:00
合計ジャッジ時間 1,918 ms
ジャッジサーバーID
（参考情報） 		judge12 / judge14
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 16 ms
4,476 KB
testcase_01 AC 15 ms
4,476 KB
testcase_02 AC 15 ms
4,484 KB
testcase_03 AC 15 ms
4,596 KB
testcase_04 AC 15 ms
4,512 KB
testcase_05 AC 16 ms
4,616 KB
testcase_06 AC 16 ms
4,572 KB
testcase_07 AC 16 ms
4,512 KB
testcase_08 AC 15 ms
4,604 KB
testcase_09 AC 15 ms
4,636 KB
testcase_10 AC 16 ms
4,588 KB
testcase_11 AC 15 ms
4,628 KB
testcase_12 AC 16 ms
4,516 KB
testcase_13 AC 16 ms
4,648 KB
testcase_14 AC 17 ms
4,604 KB
testcase_15 AC 16 ms
4,600 KB
testcase_16 AC 17 ms
4,604 KB
testcase_17 AC 17 ms
4,476 KB
testcase_18 AC 17 ms
4,604 KB
testcase_19 AC 16 ms
4,592 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