結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2019-05-31 05:34:21 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 53 ms / 2,000 ms |
コード長 | 1,097 bytes |
コンパイル時間 | 1,506 ms |
コンパイル使用メモリ | 166,732 KB |
実行使用メモリ | 11,712 KB |
最終ジャッジ日時 | 2024-09-17 16:10:38 |
合計ジャッジ時間 | 2,173 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge6 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include<bits/stdc++.h>using namespace std;using ll = long long;template<class T>bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }template<class T>bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }#define FOR(i,a,b) for(ll i=(a);i<(b);++i)#define ALL(v) (v).begin(), (v).end()#define p(s) cout<<(s)<<endl#define p2(s, t) cout << (s) << " " << (t) << endl#define br() p("")#define pn(s) cout << (#s) << " " << (s) << endl#define p_yes() p("Yes")#define p_no() p("No")const ll mod = 1e9 + 7;const ll inf = 1e18;template < typename T >void vprint(T &V){for(auto v : V){cout << v << " ";}cout << endl;}const int N_MAX = 1000010;ll A[N_MAX];int main(){cin.tie(0);ios::sync_with_stdio(false);// inputll N;cin >> N;double p;cin >> p;FOR(i, 2, N+1){for(int j=2; i*j<=N; j++){A[i*j]++;}}double sum = 0;FOR(i, 2, N+1){ll d = A[i];sum += pow(1-p, d);}cout << setprecision(20);p(sum);return 0;}