結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー IL_msta
提出日時 2015-07-01 10:19:28
言語 C++11
(gcc 15.2.0 + boost 1.89.0)
コンパイル:
g++-15 -O2 -lm -std=gnu++11 -Wuninitialized -DONLINE_JUDGE -o a.out _filename_
実行:
./a.out
結果
AC  
実行時間 61 ms / 2,000 ms
コード長 833 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 469 ms
コンパイル使用メモリ 106,484 KB
実行使用メモリ 19,584 KB
最終ジャッジ日時 2026-03-29 12:47:08
合計ジャッジ時間 1,717 ms
ジャッジサーバーID
（参考情報） 		judge3_0 / judge2_1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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