結果

問題 No.144 エラトステネスのざる
ユーザー NrKmDevNrKmDev
提出日時 2020-04-02 20:00:34
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 55 ms / 2,000 ms
コード長 1,072 bytes
コンパイル時間 1,526 ms
コンパイル使用メモリ 166,232 KB
実行使用メモリ 7,936 KB
最終ジャッジ日時 2024-06-28 15:16:48
合計ジャッジ時間 2,582 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef long double ld;
typedef pair<int,int> P;
typedef pair<ll,ll> l_l;
typedef vector<int> vi;
typedef vector<vi> vvi;
typedef vector<ll> vl;
typedef vector<P> vp;
#define fi first
#define se second
#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)
const int INF=1001001000;
const int mINF=-1001001000;
const ll LINF=1010010010010010000;
template<class T> inline bool chmin(T& a, T b) {
if (a > b) {
a = b;
return true;
}
return false;
}
template<class T> inline bool chmax(T& a, T b) {
if (a < b) {
a = b;
return true;
}
return false;
}
//グリッド：（典型）dp,dfs,bfs,最短経路,その他
int divisor[1100000];
int main(){
    int n;cin >> n;
    double p;cin>>p;
    double ans=0;
    for(int k=2;k<=n;k++){
        for(int i=k;i<=n;i+=k){
            divisor[i]++;
        }
        ans+=pow(1-p,divisor[k]-1);
    }
    cout <<fixed<<setprecision(10);
    cout << ans << endl;
    return 0;
}
0