結果
問題 | No.144 エラトステネスのざる |
ユーザー |
![]() |
提出日時 | 2021-01-24 15:22:17 |
言語 | C++17 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 73 ms / 2,000 ms |
コード長 | 1,086 bytes |
コンパイル時間 | 3,064 ms |
コンパイル使用メモリ | 193,088 KB |
最終ジャッジ日時 | 2025-01-18 07:51:08 |
ジャッジサーバーID (参考情報) |
judge1 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <bits/stdc++.h> #define M_PI 3.14159265358979323846 // pi using namespace std; typedef long long ll; typedef unsigned long long ull; typedef vector<ll> VI; typedef pair<ll, ll> P; typedef tuple<ll, ll, ll> t3; typedef tuple<ll, ll, ll, ll> t4; #define rep(a,n) for(ll a = 0;a < n;a++) #define repi(a,b,n) for(ll a = b;a < n;a++) #include <bits/stdc++.h> using namespace std; template<typename T> void chmax(T& reference, T value) { reference = max(reference, value); } template<typename T> void chmin(T& reference, T value) { reference = min(reference, value); } constexpr double mpow(double x, ll n) { double ans = 1; while (n != 0) { if (n & 1) ans = ans * x; x = x * x; n = n >> 1; } return ans; } int main() { ll n; double p; cin >> n >> p; vector<int> table(n + 1, 0); for (int i = 2; i <= n; i++) { for (int j = 1; i * j <= n; j++) { table[i * j]++; } } double ex = 0; for (int i = 2; i <= n; i++) { int c = table[i]; double q = mpow(1-p, c - 1); ex += q; } cout << setprecision(10) << fixed << ex << endl; return 0; }