結果
| 問題 |
No.144 エラトステネスのざる
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2017-11-22 23:08:03 |
| 言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
| 結果 |
AC
|
| 実行時間 | 50 ms / 2,000 ms |
| コード長 | 544 bytes |
| コンパイル時間 | 1,110 ms |
| コンパイル使用メモリ | 157,672 KB |
| 実行使用メモリ | 7,936 KB |
| 最終ジャッジ日時 | 2024-11-26 11:27:13 |
| 合計ジャッジ時間 | 2,054 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 3 |
| other | AC * 17 |
ソースコード
#include "bits/stdc++.h"
using namespace std;
#define print(x) cout<<x<<endl;
#define rep(i,a,b) for(int i=a;i<b;i++)
#define REP(i,a) for(int i=0;i<a;i++)
typedef long long ll;
typedef pair<int, int> PI;
typedef pair<int, PI> V;
typedef vector<int> VE;
const ll mod = 1000000007;
int a[1000001];
int main() {
int n;
double p;
cin>>n>>p;
for(int i=2;i<=n;i++){
for(int j=i*2;j<=n;j+=i){
a[j]++;
}
}
double e=0;
for(int i=2;i<=n;i++){
if(a[i]==0)e++;
else e += pow(1 - p, a[i]);
}
printf("%.10f\n",e);
}