結果

問題 No.144 エラトステネスのざる
ユーザー yuustiyuusti
提出日時 2015-02-13 06:19:09
言語 C++11
(gcc 11.4.0)
結果
WA  
実行時間 -
コード長 1,270 bytes
コンパイル時間 520 ms
コンパイル使用メモリ 75,260 KB
実行使用メモリ 12,016 KB
最終ジャッジ日時 2023-09-06 00:36:54
合計ジャッジ時間 2,332 ms
ジャッジサーバーID
（参考情報） 		judge11 / judge15
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 WA -
testcase_01 WA -
testcase_02 AC 2 ms
4,376 KB
testcase_03 AC 1 ms
4,376 KB
testcase_04 AC 2 ms
4,380 KB
testcase_05 AC 2 ms
4,380 KB
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <vector>
#include <algorithm>
#include <string>
#include <sstream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <numeric>
#include <cctype>
#include <tuple>
#include <array>
#include <climits>
#include <bitset>
#include <cassert>

#ifdef _MSC_VER
#include <agents.h>
#endif

#define FOR(i, a, b) for(int i = (a); i < (int)(b); ++i)
#define rep(i, n) FOR(i, 0, n)
#define ALL(v) v.begin(), v.end()
#define REV(v) v.rbegin(), v.rend()
#define MEMSET(v, s) memset(v, s, sizeof(v))
#define UNIQUE(v) (v).erase(unique(ALL(v)), (v).end())
#define MP make_pair
#define MT make_tuple

using namespace std;

typedef long long ll;
typedef long double ld;
typedef pair<int, int> P;

const int N = 1e6 + 10;
double dp[N];

int main(){
	cin.tie(0);
	ios::sync_with_stdio(false);
	cout.setf(ios::fixed);
	cout.precision(20);

	int n;
	double p;
	cin >> n >> p;	

	rep(i, n + 1) dp[i] = 1;
	for (int i = 2; i <= n; ++i){
		if (dp[i] != 1) continue;
		for (int j = i * 2; j <= n; j += i){
			int cnt = 0, x = j/i;
			while (x%i == 0) x /= i, ++cnt;
			if (!cnt) ++cnt;
			dp[j] *= pow(p, cnt);
		}
	}
	cout << accumulate(dp + 2, dp + n + 1, 0.0) << endl;
}
0