結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2019-10-12 07:41:38 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
TLE
(最新)
AC
(最初)
|
実行時間 | - |
コード長 | 1,344 bytes |
コンパイル時間 | 2,110 ms |
コンパイル使用メモリ | 171,652 KB |
実行使用メモリ | 6,824 KB |
最終ジャッジ日時 | 2024-11-27 07:02:53 |
合計ジャッジ時間 | 18,425 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 10 TLE * 7 |
ソースコード
#include <bits/stdc++.h>#define rep(i,n) for(int i = 0; i < (n); ++i)#define rrep(i,n) for(int i = 1; i <= (n); ++i)#define drep(i,n) for(int i = (n)-1; i >= 0; --i)#define srep(i,s,t) for (int i = s; i < t; ++i)#define rng(a) a.begin(),a.end()using namespace std;typedef long long int ll;typedef pair<int,int> P;typedef vector<int> vi;typedef vector<vi> vvi;typedef vector<ll> vl;typedef vector<P> vp;#define dame { puts("-1"); return 0;}#define yn {puts("Yes");}else{puts("No");}#define MAX_N 200005int yakusuu(int x){int x_ = x;vector<P> v;srep(i,2,x){if(x_ == 1)break;if(i*i>x)break;if(x_ % i == 0){int num = i;int cnt = 0;while(x_ % i == 0){cnt++;x_ /= i;}v.push_back(P(num, cnt));}}if(x_ > 1){v.push_back(P(x_, 1));}int res = 1;rep(i,v.size()){res *= v[i].second + 1;}return res;}int main() {int n;cin >> n;double p; cin >> p;p = 1-p;double ans = 0;srep(i,2,n+1){int y = yakusuu(i);y -= 2;double tmp = 1;rep(j,y){tmp *= p;}ans += tmp;}cout << fixed << setprecision(10) << ans << endl;return 0;}