結果

問題 No.144 エラトステネスのざる
ユーザー ferin
提出日時 2017-03-03 16:09:54
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 55 ms / 2,000 ms
コード長 742 bytes
コンパイル時間 1,717 ms
コンパイル使用メモリ 158,468 KB
実行使用メモリ 7,856 KB
最終ジャッジ日時 2024-06-22 05:05:23
合計ジャッジ時間 2,860 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge2
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 17
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
typedef vector<int> VI;
typedef vector<VI> VVI;
typedef vector<ll> VL;
typedef vector<VL> VVL;
typedef pair<int, int> PII;

#define FOR(i, a, n) for (ll i = (ll)a; i < (ll)n; ++i)
#define REP(i, n) FOR(i, 0, n)
#define ALL(x) x.begin(), x.end()
#define MOD 1000000007
#define INF (1LL<<25)     //33554432
#define PI 3.14159265359
#define EPS 1e-12
//#define int ll

int a[1000010];
signed main(void)
{
  int n;
  double p;
  cin >> n >> p;

  FOR(i, 2, n+1)
    for(int j=i*2; j<=n; j+=i)
      a[j]++;
  double ans = 0;
  FOR(i, 2, n+1) {
    if(a[i] == 0) ans += 1;
    else ans += pow((1-p), a[i]);
  }

  cout << fixed << setprecision(15) << ans << endl;
  return 0;
}
0