結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2015-09-23 00:30:53 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 71 ms / 2,000 ms |
コード長 | 1,298 bytes |
コンパイル時間 | 611 ms |
コンパイル使用メモリ | 72,100 KB |
実行使用メモリ | 11,776 KB |
最終ジャッジ日時 | 2024-07-19 08:44:58 |
合計ジャッジ時間 | 1,920 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include<iostream>#include<math.h>#include<vector>#include <array>#include <algorithm>#include <queue>#include <numeric>#include <map>using namespace std;#define int long long#define FOR(i, a, b) for(int i=(a);i<(b);i++)#define RFOR(i, a, b) for(int i=(b-1);i>=(a);i--)#define REP(i, n) for(int i=0; i<(n); i++)#define RREP(i, n) for(int i=(n-1); i>=0; i--)#define ALL(a) (a).begin(),(a).end()#define UNIQUE_SORT(l) sort(ALL(l)); l.erase(unique(ALL(l)), l.end());#define CONTAIN(a, b) find(ALL(a), (b)) != (a).end()#define array2(type, x, y) array<array<type, y>, x>#define vector2(type) vector<vector<type> >#define out(...) printf(__VA_ARGS__)typedef pair<int, int> pos;int pos::*x = &pos::first;int pos::*y = &pos::second;int dxy[] = {0, 1, 0, -1, 0};/*================================*/signed main(){int n;float p;cin >> n >> p;int ps[1000001] = {0};double total = 0.0;double q = 1.0 - p;FOR(i,2,n+1) {// 消されずに残っている確率を足す// 1->1.0, 2->(1-p)double k = pow(q, ps[i]);total += k;// 全ての数で篩を行うint j = i;while (j <= n) {ps[j]++;j += i;}}out("%.10lf\n", total);return 0;}