結果

問題 No.144 エラトステネスのざる
ユーザー tossy
提出日時 2017-06-06 15:23:36
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 45 ms / 2,000 ms
コード長 891 bytes
コンパイル時間 1,512 ms
コンパイル使用メモリ 168,284 KB
実行使用メモリ 11,408 KB
最終ジャッジ日時 2024-09-22 11:44:44
合計ジャッジ時間 2,572 ms
ジャッジサーバーID
（参考情報） 		judge3 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
#define repl(i,a,b) for(int i=(int)(a);i<(int)(b);i++)
#define rep(i,n) repl(i,0,n)
#define mp(a,b) make_pair((a),(b))
#define pb(a) push_back((a))
#define all(x) (x).begin(),(x).end()
#define uniq(x) sort(all(x)),(x).erase(unique(all(x)),end(x))
#define fi first
#define se second
#define dbg(x) cout<<#x" = "<<((x))<<endl
template<class T,class U> ostream& operator<<(ostream& o, const pair<T,U> &p){o<<"("<<p.fi<<","<<p.se<<")";return o;}
template<class T> ostream& operator<<(ostream& o, const vector<T> &v){o<<"[";for(T t:v){o<<t<<",";}o<<"]";return o;}

#define INF 2147483600

int main(){
  int n; double p;
  cin>>n>>p;
  vector<double> vec(n+1, 1);
  // 残っている確率

  double ans = 0;
  for(int i=2; i<=n; i++){
    ans += vec[i];
    for(int j=2*i; j<=n; j+=i) vec[j] *= (1.0-p);
  }

  printf("%.8f\n", ans);

  return 0;
}
0