結果

問題 No.144 エラトステネスのざる
ユーザー tubo28
提出日時 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
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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);
//     }
// }
0