結果
問題 | 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 1000000000000000000double 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;}