結果

問題 No.144 エラトステネスのざる
ユーザー umezoumezo
提出日時 2022-07-24 21:32:07
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 77 ms / 2,000 ms
コード長 536 bytes
コンパイル時間 2,184 ms
コンパイル使用メモリ 203,432 KB
実行使用メモリ 15,176 KB
最終ジャッジ日時 2024-07-06 16:48:09
合計ジャッジ時間 2,949 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 5 ms
11,204 KB
testcase_01 AC 5 ms
11,144 KB
testcase_02 AC 4 ms
11,260 KB
testcase_03 AC 4 ms
11,316 KB
testcase_04 AC 5 ms
11,352 KB
testcase_05 AC 5 ms
11,256 KB
testcase_06 AC 4 ms
11,260 KB
testcase_07 AC 5 ms
11,212 KB
testcase_08 AC 50 ms
11,308 KB
testcase_09 AC 5 ms
11,344 KB
testcase_10 AC 5 ms
11,364 KB
testcase_11 AC 6 ms
11,132 KB
testcase_12 AC 5 ms
11,364 KB
testcase_13 AC 32 ms
15,124 KB
testcase_14 AC 32 ms
15,156 KB
testcase_15 AC 77 ms
15,128 KB
testcase_16 AC 33 ms
15,176 KB
testcase_17 AC 33 ms
15,052 KB
testcase_18 AC 33 ms
15,124 KB
testcase_19 AC 35 ms
15,120 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#define rep(i,n) for(int i=0;i<(int)(n);i++)
#define ALL(v) v.begin(),v.end()
typedef long long ll;
 
#include<bits/stdc++.h>
using namespace std;

int main(){
  ios::sync_with_stdio(false);
  std::cin.tie(nullptr);
  
  int n;
  double p;
  cin>>n>>p;
  
  vector<double> D(1000100);
  D[0]=1.0;
  for(int i=1;i<1000100;i++) D[i]=D[i-1]*(1-p);
  
  vector<int> dp(n+1);
  for(int i=2;2*i<=n;i++){
    for(int j=2*i;j<=n;j+=i) dp[j]++;
  }
  double ans=0;
  for(int i=2;i<=n;i++) ans+=D[dp[i]];
  printf("%.10f\n",ans);
  
  return 0;
}
0