結果

問題 No.144 エラトステネスのざる
ユーザー a0171610
提出日時 2020-05-19 21:47:24
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 52 ms / 2,000 ms
コード長 790 bytes
コンパイル時間 1,540 ms
コンパイル使用メモリ 168,120 KB
実行使用メモリ 10,128 KB
最終ジャッジ日時 2024-10-01 23:04:31
合計ジャッジ時間 2,577 ms
ジャッジサーバーID
(参考情報)
judge1 / judge5
このコードへのチャレンジ
(要ログイン)
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function 'int main()':
main.cpp:26:7: warning: 'ans' is used uninitialized [-Wuninitialized]
   26 |   ans += double(p);
      |   ~~~~^~~~~~~~~~~~
main.cpp:25:10: note: 'ans' was declared here
   25 |   double ans;
      |          ^~~

ソースコード

diff #

#include <bits/stdc++.h>
using namespace std;
using lint = long long;
const int MAX_N = 1e6+10;
int prime[MAX_N], dataset[MAX_N];
bool is_prime[MAX_N];
int p = 0;

void sieve(int n){
  for(int i = 0; i <= n; i++) is_prime[i] = true;
  is_prime[0] = false; is_prime[1] = false;
  for(int i = 2; i <= n; i++){
    if(is_prime[i]){
      prime[p++] = i;
      for(int j = 2 * i; j <= n; j += i) is_prime[j] = false;
    }
    for(int j = 2 * i; j <= n; j += i) dataset[j]++;
  }
}

signed main(){
  int n; cin >> n;
  double P; cin >> P;
  sieve(n);
  double ans;
  ans += double(p);
  for(int i = 2; i <= n; i++){
    if(dataset[i]){
      double ans1 = 1;
      for(int j = 0; j < dataset[i]; j++) ans1 *= (1.0 - P);
      ans += ans1;
    }
  }
  cout << setprecision(10) << ans << endl;
}
0