結果

問題 No.144 エラトステネスのざる
ユーザー GyuutoGyuuto
提出日時 2015-11-24 23:38:43
言語 C++11
(gcc 11.4.0)
結果
AC  
実行時間 52 ms / 2,000 ms
コード長 2,176 bytes
コンパイル時間 520 ms
コンパイル使用メモリ 75,612 KB
実行使用メモリ 12,376 KB
最終ジャッジ日時 2023-10-11 18:28:43
合計ジャッジ時間 1,876 ms
ジャッジサーバーID
（参考情報） 		judge14 / judge12
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 8 ms
6,236 KB
testcase_01 AC 7 ms
6,172 KB
testcase_02 AC 8 ms
6,104 KB
testcase_03 AC 8 ms
6,296 KB
testcase_04 AC 7 ms
6,288 KB
testcase_05 AC 8 ms
6,104 KB
testcase_06 AC 8 ms
6,236 KB
testcase_07 AC 7 ms
6,124 KB
testcase_08 AC 8 ms
6,244 KB
testcase_09 AC 8 ms
6,296 KB
testcase_10 AC 7 ms
6,204 KB
testcase_11 AC 7 ms
6,168 KB
testcase_12 AC 8 ms
6,164 KB
testcase_13 AC 51 ms
12,364 KB
testcase_14 AC 51 ms
12,364 KB
testcase_15 AC 50 ms
12,256 KB
testcase_16 AC 50 ms
12,376 KB
testcase_17 AC 50 ms
12,372 KB
testcase_18 AC 52 ms
12,332 KB
testcase_19 AC 50 ms
12,256 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

// Standard I/O
#include <iostream>
#include <sstream>
#include <cstdio>
// Standard Library
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <cmath>
// Template Class
#include <complex>
#include <string>
#include <vector>
#include <list>
#include <set>
#include <map>
#include <queue>
#include <stack>
// Container Control
#include <algorithm>

using namespace std;

#define rep( i, n ) for( int i = 0; i < n; ++i )
#define irep( i, n ) for( int i = n-1; i >= 0; --i )
#define reep( i, s, n ) for ( int i = s; i < n; ++i )
#define ireep( i, n, s ) for ( int i = n-1; i >= s; --i )
#define foreach(itr, x) for( typeof(x.begin()) itr = x.begin(); itr != x.end(); ++itr)

#define mp( a, b ) make_pair( a, b )
#define pb( a ) push_back( a )
#define all( v ) v.begin(), v.end()
#define fs first
#define sc second
#define vc vector

// for visualizer.html
double SCALE = 1.0;
double OFFSET_X = 0.0;
double OFFSET_Y = 0.0;
#define LINE(x,y,a,b) cerr << "line(" << SCALE*(x) + OFFSET_X << ","	\
	<< SCALE*(y) + OFFSET_Y << ","										\
	<< SCALE*(a) + OFFSET_X << ","										\
	<< SCALE*(b) + OFFSET_Y << ")" << endl;
#define CIRCLE(x,y,r) cerr << "circle(" << SCALE*(x) + OFFSET_X << ","	\
	<< SCALE*(y) + OFFSET_Y << ","										\
	<< SCALE*(r) << ")" << endl;

typedef long long ll;
typedef complex<double> Point;

typedef pair<int, int> pii;
typedef pair<int, pii> ipii;
typedef vector<int> vi;
typedef vector<double> vd;
typedef vector< vector<int> > vii;
typedef vector< vector<double> > vdd;

typedef vector<int>::iterator vi_itr;

const int IINF = 1 << 28;
const double INF = 1e30;
const double EPS = 1e-10;
const double PI = acos(-1.0);

// Direction : L U R D
const int dx[] = { -1, 0, 1, 0};
const int dy[] = { 0, -1, 0, 1 };

bool isprime[1000010];
double dp[1000010];
int main()
{
	int N;
	double p;
	reep(i, 2, 1000010) isprime[i] = true;
	reep(i, 2, 1000010) if( isprime[i] ) for( int j = 2*i; j < 1000010; j+=i ) isprime[j] = false;

	cin >> N >> p;

	reep(i, 2, N+1) dp[i] = 1.0;
	reep(i, 2, N+1){
		for( int j = 2*i; j <= N; j+=i ){
			dp[j] *= (1 - p);
		}
	}

	double ans = 0.0;
	reep(i, 2, N+1){
		ans += dp[i];
	}
	printf("%.16f\n", ans);
}
0