結果
問題 | No.144 エラトステネスのざる |
ユーザー |
|
提出日時 | 2015-02-05 23:55:07 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 54 ms / 2,000 ms |
コード長 | 900 bytes |
コンパイル時間 | 1,242 ms |
コンパイル使用メモリ | 158,188 KB |
実行使用メモリ | 11,520 KB |
最終ジャッジ日時 | 2024-06-23 08:54:27 |
合計ジャッジ時間 | 2,442 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <bits/stdc++.h>using namespace std;#define all(c) (c).begin(),(c).end()#define rrep(i,n) for(int i=(int)(n)-1;i>=0;i--)#define REP(i,m,n) for(int i=(int)(m);i<(int)(n);i++)#define rep(i,n) REP(i,0,n)#define iter(c) __typeof((c).begin())#define tr(it,c) for(iter(c) it=(c).begin();it!=(c).end();it++)#define pb(a) push_back(a)#define pr(a) cout<<(a)<<endl#define PR(a,b) cout<<(a)<<" "<<(b)<<endl#define R cin>>#define F first#define S second#define ll long longbool check(int n,int m,int x,int y){return (x<0||x>=n||y<0||y>=m)?false:true;}const ll MAX=1000000007,MAXL=1LL<<60,dx[4]={-1,0,1,0},dy[4]={0,1,0,-1};typedef pair<int,int> P;int main() {int n;double p;cin >> n >> p;double d[n+1];fill(d,d+n+1,1);p=1-p;REP(i,2,n+1) {for(int j=i*2; j<=n; j+=i) d[j]*=p;}double ans=0;REP(i,2,n+1) ans+=d[i];printf("%.10f\n",ans);return 0;}