結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2015-05-11 00:10:36 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 54 ms / 2,000 ms |
コード長 | 799 bytes |
コンパイル時間 | 578 ms |
コンパイル使用メモリ | 83,148 KB |
実行使用メモリ | 7,808 KB |
最終ジャッジ日時 | 2024-07-05 22:09:13 |
合計ジャッジ時間 | 1,683 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <iostream>#include <sstream>#include <string>#include <vector>#include <list>#include <stack>#include <queue>#include <deque>#include <set>#include <map>#include <algorithm>#include <functional>#include <limits>#include <cstdio>#include <cstdlib>#include <cmath>#include <climits>#include <complex>using namespace std;#define INF (INT_MAX/3)#define PI (2*acos(0.0))#define EPS (1e-8)typedef long long ll;typedef unsigned long long ull;int a[1000001];int main(){ios_base::sync_with_stdio(0);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 += 1;else e += pow(1 - P, a[i]);}printf("%.10f\n", e);return 0;}