結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2019-10-20 18:30:45 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 48 ms / 2,000 ms |
コード長 | 565 bytes |
コンパイル時間 | 1,533 ms |
コンパイル使用メモリ | 166,204 KB |
実行使用メモリ | 7,808 KB |
最終ジャッジ日時 | 2024-07-02 17:28:45 |
合計ジャッジ時間 | 2,826 ms |
ジャッジサーバーID (参考情報) |
judge2 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include<bits/stdc++.h>using namespace std;double m_pow(double x,int n){ //べき乗double ans=1;while(n!=0){if(n&1) ans=ans*x;x=x*x;n=n>>1;}return ans;}int main(){int N; cin >> N;double p; cin >> p;double ans=0;int divisor[N+1]={};for(int i=2;i<=N;i++){for(int j=i*2;j<=N;j+=i){divisor[j]++;}}for(int i=2;i<=N;i++){if(divisor[i]==0) ans+=1;else ans+=m_pow(1-p,divisor[i]);}cout << setprecision(10) << ans << endl;}