結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー IL_msta
提出日時 2015-07-01 10:19:28
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 83 ms / 2,000 ms
コード長 833 bytes
コンパイル時間 674 ms
コンパイル使用メモリ 84,452 KB
実行使用メモリ 19,200 KB
最終ジャッジ日時 2024-07-07 21:39:29
合計ジャッジ時間 1,910 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define _USE_MATH_DEFINES

#include <iostream>
#include <iomanip>

#include <algorithm>
#include <cmath>

#include <string>
#include <list>
#include <queue>
#include <vector>
#include <complex>
#include <set>

/////////
#define REP(i, x, n) for(int i = x; i < n; i++)
#define rep(i,n) REP(i,0,n)
#define P(p) cout<<(p)<<endl;
/////////
typedef long long LL;
typedef long double LD;
/////////
using namespace::std;
/////////

const int Nmax = (int)1e6;
LD num[Nmax+1];

int main(void){
    std::cin.tie(0); 
    std::ios::sync_with_stdio(false);
    std::cout << std::fixed;//
    //cout << setprecision(7);//
    int N;
	LD p,invp;
	cin>>N>>p;
	invp = 1.0-p;
	LD ans = 0;
	for(int i=2;i<=N;++i){
		num[i] = 1;
	}
	for(int i=2;i<=N;++i){
		ans += num[i];
		for(int k=2;i*k<=N;++k){
			num[i*k] *= invp;
		}
	}
	P(ans);
    return 0;
}
0