結果

問題 No.144 エラトステネスのざる
コンテスト
ユーザー tkmst201
提出日時 2017-05-14 13:12:00
言語 C++11（廃止可能性あり）
(gcc 13.3.0 + boost 1.89.0)
結果
TLE  
実行時間 -
コード長 926 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 1,417 ms
コンパイル使用メモリ 159,164 KB
実行使用メモリ 6,948 KB
最終ジャッジ日時 2024-09-16 05:48:33
合計ジャッジ時間 16,421 ms
ジャッジサーバーID
（参考情報） 		judge6 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 10 TLE * 7
権限があれば一括ダウンロードができます

ソースコード

diff #
raw source code 

#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, pp[100];

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