結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2015-07-01 10:19:28 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 83 ms / 2,000 ms |
コード長 | 833 bytes |
コンパイル時間 | 674 ms |
コンパイル使用メモリ | 84,452 KB |
実行使用メモリ | 19,200 KB |
最終ジャッジ日時 | 2024-07-07 21:39:29 |
合計ジャッジ時間 | 1,910 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#define _USE_MATH_DEFINES#include <iostream>#include <iomanip>#include <algorithm>#include <cmath>#include <string>#include <list>#include <queue>#include <vector>#include <complex>#include <set>/////////#define REP(i, x, n) for(int i = x; i < n; i++)#define rep(i,n) REP(i,0,n)#define P(p) cout<<(p)<<endl;/////////typedef long long LL;typedef long double LD;/////////using namespace::std;/////////const int Nmax = (int)1e6;LD num[Nmax+1];int main(void){std::cin.tie(0);std::ios::sync_with_stdio(false);std::cout << std::fixed;////cout << setprecision(7);//int N;LD p,invp;cin>>N>>p;invp = 1.0-p;LD ans = 0;for(int i=2;i<=N;++i){num[i] = 1;}for(int i=2;i<=N;++i){ans += num[i];for(int k=2;i*k<=N;++k){num[i*k] *= invp;}}P(ans);return 0;}