結果
問題 | No.144 エラトステネスのざる |
ユーザー |
|
提出日時 | 2024-06-11 22:03:12 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 43 ms / 2,000 ms |
コード長 | 1,459 bytes |
コンパイル時間 | 2,992 ms |
コンパイル使用メモリ | 139,852 KB |
実行使用メモリ | 11,196 KB |
最終ジャッジ日時 | 2024-06-11 22:03:16 |
合計ジャッジ時間 | 3,938 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#include <algorithm>#include <cassert>#include <cstdio>#include <cstdlib>#include <cstring>#include <cmath>#include <iostream>#include <tuple>#include <vector>#include <map>#include <set>#include <queue>#include <ranges>#include <iomanip>using namespace std;using ll = long long;auto range(int n) { return views::iota(0, n); }template<class T, class U> ostream& operator<<(ostream& os, const pair<T, U>& p){ return os << "{" << p.first << ", " << p.second << "}"; }template<typename T> ostream& operator<<(ostream& os, const vector<T>& obj) { os << "{"; for (const auto& e : obj) os << e << ", "; return os << "}";}template<typename T> ostream& operator<<(ostream& os, const set<T>& obj) { os << "set{"; for (const auto& e : obj) os << e << ", "; return os << "}";}template<typename T, typename U> ostream& operator<<(ostream& os, const map<T, U>& obj) { os << "map{"; for (const auto& e : obj) os << e << ", ";return os << "}"; }#ifdef ONLINE_JUDGE#define dump(expr) ;#else#define dump(expr) { cerr << "\033[33m#L" << __LINE__ << ": " << expr << "\033[39m" << endl; }#endifdouble solve() {int n; double p; cin >> n >> p;vector<double> vs(n + 1, 1.0);double res = 0;for (int i = 2; i <= n; ++i) {res += vs[i];dump(i << " " << vs[i]);for (int j = i * 2; j <= n; j += i) {vs[j] *= 1 - p;}}return res;}int main() {cout << fixed << setprecision(12);cout << solve() << endl;}