結果

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

ソースコード

diff # 

#include <bits/stdc++.h>

using namespace std;

#define FOR(i,a,b) for(int i=(a);i<(b);i++)
#define REP(i,n) FOR(i,0,n)
#define ALL(v) (v).begin(),(v).end()

template<typename A, typename B> inline bool chmax(A &a, B b) { if (a<b) { a=b; return 1; } return 0; }
template<typename A, typename B> inline bool chmin(A &a, B b) { if (a>b) { a=b; return 1; } return 0; }

typedef unsigned long long ull;
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<ll, ll> pll;
typedef pair<int, pii> P;

const ll INF = 1ll<<60;
const ll MOD = 1000000007;
const double EPS  = 1e-10;

int n;
double p;
long double pp[10000];
int cnt[1123456];

int main() {	
	cin >> n >> p;
	
	FOR(i, 1, n + 1) {
		for (int j = i; j <= n; j += i) {
			cnt[j]++;
		}
	}
	
	pp[0] = 1;
	FOR(i, 1, 10000) pp[i] = pp[i - 1] * (1 - p);
	
	long double ans = 0.0;
	FOR(i, 2, n + 1) ans += pp[cnt[i] - 2];
	
	printf("%.15Lf\n", ans);
	
	return 0;
}
0