結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー onsen_manjuuu
提出日時 2020-09-13 13:18:56
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 1,011 bytes
コンパイル時間 1,551 ms
コンパイル使用メモリ 167,292 KB
実行使用メモリ 11,392 KB
最終ジャッジ日時 2024-06-12 02:32:42
合計ジャッジ時間 2,733 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 14 WA * 3
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include<bits/stdc++.h>

using namespace std;
using ll = long long;
/*********************************/
template<class T> inline ostream& operator<<(ostream& os,vector<T> arr) {os << "["; for(int i = 0; i < (int)arr.size(); i++)os << arr[i] << (i == (int)arr.size() - 1 ? "]" : ",");os << endl;return os;}
template<typename A, typename B> ostream& operator<<(ostream& os, const pair<A,B>& p){os << "{" << p.first << "," << p.second << "}";return os; }
#define prvec(v) cerr<<#v<<": [";for(int __i = 0;__i < (int)(v).size(); __i++)cerr<<((v)[__i])<<(__i+1==(int)(v).size()?"]\n":",");
#define print(x) cerr<<#x<<": "<<(x)<<endl
/*********************************/


int main()
{
	int n;
	double p;
	cin >> n >> p;
	p = 1.0 - p;
	vector<double> ps(n + 1, 1.0);
	ps[0] = ps[1] = 0;
	for(int i = 2; i <= n; i++) {
		for(int j = i + i; j <= n; j += i) {
			ps[j] *= p;
		}
	}
	// print(ps);
	// for(int i = 0; i <= n; i++)cerr << i << ": " << ps[i] << endl;

	cout << accumulate(ps.begin(), ps.end(), 0.0) << endl;
}
0