結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2015-05-29 15:18:33 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 50 ms / 2,000 ms |
コード長 | 1,471 bytes |
コンパイル時間 | 582 ms |
コンパイル使用メモリ | 87,384 KB |
実行使用メモリ | 7,936 KB |
最終ジャッジ日時 | 2024-07-06 10:43:58 |
合計ジャッジ時間 | 1,532 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <iostream>#include <vector>#include <string>#include <stack>#include <queue>#include <deque>#include <set>#include <map>#include <algorithm> // require sort next_permutation count __gcd reverse etc.#include <cstdlib> // require abs exit atof atoi#include <cstdio> // require scanf printf#include <functional>#include <numeric> // require accumulate#include <cmath> // require fabs#include <climits>#include <limits>#include <cfloat>#include <iomanip> // require setw#include <sstream> // require stringstream#include <cstring> // require memset#include <cctype> // require tolower, toupper#include <fstream> // require freopen#include <ctime> // require srand#define rep(i,n) for(int i=0;i<(n);i++)#define ALL(A) A.begin(), A.end()using namespace std;typedef long long ll;typedef pair<int, int> P;const int MAX_N = (int)1e6 + 6;int cnt[MAX_N];// 1と自分自身を除く約数の個数void divisor (int N ){for (int i = 2; i <= N; i++ ){for (int j = 2*i; j <= N; j += i ){cnt[j]++;} // end for} // end for}int main(){memset (cnt, 0, sizeof (cnt ) );ios_base::sync_with_stdio(0);int N; cin >> N;double p; cin >> p;divisor(N );double res = 0.;for (int i = 2; i <= N; i++ ){double curr = pow(1.- p, cnt[i] );// cerr << curr << endl;res += curr;} // end for// for (int i = 2; i <= N; i++ )// cerr << cnt[i] << (i != N ? ' ' : '\n' );printf ("%.5lf\n", res );return 0;}