結果
問題 | No.144 エラトステネスのざる |
ユーザー |
|
提出日時 | 2017-09-17 17:13:52 |
言語 | C++11 (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 54 ms / 2,000 ms |
コード長 | 1,318 bytes |
コンパイル時間 | 625 ms |
コンパイル使用メモリ | 75,504 KB |
実行使用メモリ | 11,392 KB |
最終ジャッジ日時 | 2024-11-08 01:55:42 |
合計ジャッジ時間 | 1,758 ms |
ジャッジサーバーID (参考情報) |
judge1 / judge4 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 3 |
other | AC * 17 |
ソースコード
#define _CRT_SECURE_NO_WARNINGS#include <iostream>#include <cmath>#include <string>#include <vector>#include <algorithm>#include <queue>#include <map>#include <functional>#include <set>#include <numeric>#include <stack>#include <utility>#include <time.h>//#include "util.h"using namespace std;typedef long long ll;typedef unsigned long long ull;#define PI 3.14159265358979323846#define EPS 1e-6#define MIN(a,b) ((a)<(b)?(a):(b))#define MAX(a,b) ((a)>(b)?(a):(b))#define CHAR_BIT 8#define all(x) (x).begin(), (x).end()template <typename _Ty>ostream& operator << (ostream& ostr, const vector<_Ty>& v) {if (v.empty()) {cout << "{ }";return ostr;}cout << "{" << v.front();for (auto itr = ++v.begin(); itr != v.end(); itr++) {cout << ", " << *itr;}cout << "}";return ostr;}ll GCD(ll m, ll n){ll k;do {k = m % n;m = n; n = k;} while (k != 0);return m;}int yuki0144(){int N;double P, ans = 0;cin >> N >> P;vector<double> p(N + 1, 1.0);p[0] = p[1] = 0;for (int i = 2; i <= N; i++) {ans += p[i];for (int j = 2 * i; j <= N; j += i)p[j] *= 1 - P;}printf("%.10f\n", ans);return 0;}int main(){//clock_t start, end;//start = clock();yuki0144();//end = clock();//printf("%d msec.\n", end - start);return 0;}