結果

問題 No.144 エラトステネスのざる
ユーザー 🍡yurahuna🍡yurahuna
提出日時 2016-02-28 21:41:38
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 54 ms / 2,000 ms
コード長 929 bytes
コンパイル時間 1,449 ms
コンパイル使用メモリ 74,016 KB
実行使用メモリ 8,116 KB
最終ジャッジ日時 2023-10-24 18:47:45
合計ジャッジ時間 1,746 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
4,348 KB
testcase_01 AC 2 ms
4,348 KB
testcase_02 AC 2 ms
4,348 KB
testcase_03 AC 2 ms
4,348 KB
testcase_04 AC 1 ms
4,348 KB
testcase_05 AC 2 ms
4,348 KB
testcase_06 AC 1 ms
4,348 KB
testcase_07 AC 2 ms
4,348 KB
testcase_08 AC 1 ms
4,348 KB
testcase_09 AC 2 ms
4,348 KB
testcase_10 AC 2 ms
4,348 KB
testcase_11 AC 1 ms
4,348 KB
testcase_12 AC 2 ms
4,348 KB
testcase_13 AC 45 ms
8,056 KB
testcase_14 AC 52 ms
8,056 KB
testcase_15 AC 52 ms
8,056 KB
testcase_16 AC 52 ms
8,056 KB
testcase_17 AC 52 ms
8,056 KB
testcase_18 AC 54 ms
8,116 KB
testcase_19 AC 40 ms
7,952 KB
権限があれば一括ダウンロードができます

ソースコード

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