結果
問題 |
No.144 エラトステネスのざる
|
ユーザー |
![]() |
提出日時 | 2018-08-14 21:41:37 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 41 ms / 2,000 ms |
コード長 | 850 bytes |
コンパイル時間 | 1,656 ms |
コンパイル使用メモリ | 158,508 KB |
実行使用メモリ | 11,556 KB |
最終ジャッジ日時 | 2024-09-24 09:06:39 |
合計ジャッジ時間 | 2,470 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include "bits/stdc++.h" using namespace std; #define ll long long int #define rep(i,n) for( int i = 0; i < n; i++ ) #define rrep(i,n) for( int i = n; i >= 0; i-- ) #define REP(i,s,t) for( int i = s; i <= t; i++ ) #define RREP(i,s,t) for( int i = s; i >= t; i-- ) #define dump(x) cerr << #x << " = " << (x) << endl; #define INF 2000000000 #define mod 1000000007 #define INF2 1000000000000000000 double L[1000001] = {}; int main(void) { cin.tie(0); ios::sync_with_stdio(false); int N; cin >> N; double p; cin >> p; rep(i, N + 1) L[i] = 1; L[0] = 0; L[1] = 0; for(int i = 2; i <= N; i++) { int j = 2; while (i * j <= N) { L[i * j] *= 1 - p; j++; } } double ans = 0; rep(i, N + 1) ans += L[i]; //rep(i, N + 1) cout << i << " " << L[i] << endl; cout << fixed << setprecision(14); cout << ans << endl; return 0; }