結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2018-06-30 09:02:45 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 55 ms / 2,000 ms |
コード長 | 2,115 bytes |
コンパイル時間 | 853 ms |
コンパイル使用メモリ | 100,064 KB |
実行使用メモリ | 7,680 KB |
最終ジャッジ日時 | 2024-07-01 00:50:03 |
合計ジャッジ時間 | 2,008 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <iostream>#include <cstdio>#include <cstdlib>#include <cstring>#include <climits>#include <cmath>#include <complex>#include <map>#include <set>#include <vector>#include <stack>#include <queue>#include <bitset>#include <algorithm>#include <numeric>#include <functional>#include <cassert>#include <iomanip>using namespace std;#define Rep(b, e, i) for(int i = b; i <= e; i++)#define Repr(e, b, i) for(int i = e; i >= b; i--)#define rep(n, i) Rep(0, n-1, i)#define repr(n, i) Repr(n-1, 0, i)#define all(v) (v).begin(), (v).end()#define pb(v) push_back(v)#define uniq(v) (v).erase(unique(all(v)),(v).end())#define bitcnt(x) __builtin_popcount(x)#define fst first#define snd second#define Pqaz(T) priority_queue<T,vector<T>,greater<T>>#define Pqza(T) priority_queue<T>#define put(x) cout << x;#define putsp(x) cout << x << ' ';#define putln(x) cout << x << endl;#define ENJYU std::ios::sync_with_stdio(false);std::cin.tie(0);typedef long long ll;typedef pair<ll, ll> llP;typedef pair<int, int> intP;typedef complex<double> comp;typedef vector <int> vec;typedef vector <ll> vecll;typedef vector <double> vecd;typedef vector <vec> mat;typedef vector <vecll> matll;typedef vector <vecd> matd;//vector の中身を出力template <class T>ostream &operator<<(ostream &o,const vector<T>&v){o<<"{";for(int i=0;i<(int)v.size();i++)o<<(i>0?", ":"")<<v[i];o<<"}";return o;}const int MAX = 200020;const double PI = 3.14159265358979323846;const double EPS = 1e-12;const int INF = 1<<29;const ll INFL = 1e18;const int MOD = 1000000007;int dx4[4]={1,0,-1,0};int dy4[4]={0,1,0,-1};int dx8[8]={1,0,-1,1,-1,1,0,-1};int dy8[8]={1,1,1,0,0,-1,-1,-1};//********************************template END****************************************//void solve(void){int N; double p;cin >> N >> p;double ans = 0;vec fact(N+1, 0);Rep(2, N, n){for (int i = n; i <= N; i += n){fact[i]++;}ans += pow(1-p, fact[n]-1);}cout << fixed << setprecision(6) << ans << endl;}int main(void){solve();//cout << "yui(*-v・)yui" << endl;return 0;}