結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2015-02-06 00:07:13 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 59 ms / 2,000 ms |
コード長 | 2,733 bytes |
コンパイル時間 | 1,299 ms |
コンパイル使用メモリ | 157,540 KB |
実行使用メモリ | 11,520 KB |
最終ジャッジ日時 | 2024-06-23 09:01:34 |
合計ジャッジ時間 | 2,430 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#define _CRT_SECURE_NO_WARNINGS //#define _GLIBCXX_DEBUG #include <bits/stdc++.h> using namespace std; typedef long long ll; typedef vector<int> vi; typedef vector<vi> vvi; typedef pair<int,int> pii; #define all(c) (c).begin(), (c).end() #define loop(i,a,b) for(ll i=a; i<ll(b); i++) #define rep(i,b) loop(i,0,b) #define pb push_back #define eb emplace_back #define mp make_pair #define mt make_tuple template<class T> istream& operator>>(istream& is, vector<T>& v); template<class T> ostream& operator<<(ostream& os, vector<T> const& v); template<int n, class...T> typename enable_if<(n>=sizeof...(T))>::type _ot(ostream&, const tuple<T...>&){} template<int n, class...T> typename enable_if<(n< sizeof...(T))>::type _ot(ostream& os, const tuple<T...>& t){ os << (n==0?"":" ") << get<n>(t); _ot<n+1>(os, t); } template<class...T> ostream& operator<<(ostream& os, const tuple<T...>& t){ _ot<0>(os, t); return os; } template<int n, class...T> typename enable_if<(n>=sizeof...(T))>::type _it(istream&, tuple<T...>&){} template<int n, class...T> typename enable_if<(n< sizeof...(T))>::type _it(istream& is, tuple<T...>& t){ is >> get<n>(t); _it<n+1>(is, t); } template<class...T> istream& operator>>(istream& is, tuple<T...>& t){ _it<0>(is, t); return is; } template<class T> istream& operator>>(istream& is, vector<T>& v){ rep(i,v.size()) is >> v[i]; return is; } template<class T> ostream& operator<<(ostream& os, vector<T> const& v){ rep(i,v.size()) os << v[i] << (i+1==(int)v.size()?"":" "); return os; } #ifdef DEBUG #define dump(...) (cerr<<#__VA_ARGS__<<" = "<<mt(__VA_ARGS__)<<" ["<<__LINE__<<"]"<<endl) #else #define dump(...) #endif int N; double p; int cnt[1000010]; double dp[1000010]; int main(){ while(cin>>N>>p){ // for(int i=2;i<=N;i++){ // dp[i] = 1; // for(int j=2;j*j<=i;j++){ // if(i%j==0){ // dp[i] *= (1-p)*dp[j]; // if(j*j!=i) dp[i] *= (1-p)*dp[i/j]; // } // } // } for(int i=2;i<=N;i++) dp[i] = 1; for(int i=2;i<=N;i++){ for(int j=i+i; j<=N; j+=i){ dp[j] *= 1-p; } } double ans = 0; for(int i=2;i<=N;i++){ ans += dp[i]; } printf("%.10lf\n", ans); } } // memset(cnt,0,sizeof(cnt)); // for(int i=2;i<=1000000;i++){ // for(int j=i; j<=1000000; j+=i){ // cnt[j]++; // } // } // dump(vi(cnt, cnt+20)); // while(cin>>N>>p){ // double ans = 0; // for(int i=2;i<=N;i++){ // dump(pow(p,cnt[i])); // ans += 1-pow(1-p,cnt[i]-1); // } // printf("%.10lf\n", ans); // } // }