#329249 (C++14) No.144 エラトステネスのざる

• 問題: No.144 エラトステネスのざる
• ユーザー: minami
• 提出日時: 2019-03-27 04:15:44
• 言語: C++14 (gcc 13.2.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 36 ms / 2,000 ms
• コード長: 1,238 bytes
• コンパイル時間: 1,824 ms
• コンパイル使用メモリ: 167,804 KB
• 実行使用メモリ: 7,552 KB
• 最終ジャッジ日時: 2024-04-19 09:07:48
• 合計ジャッジ時間: 2,868 ms
• ジャッジサーバーID （参考情報）: judge4 / judge2

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	1 ms; 5,248 KB
testcase_01	AC	1 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	2 ms; 5,376 KB
testcase_04	AC	2 ms; 5,376 KB
testcase_05	AC	2 ms; 5,376 KB
testcase_06	AC	2 ms; 5,376 KB
testcase_07	AC	2 ms; 5,376 KB
testcase_08	AC	2 ms; 5,376 KB
testcase_09	AC	1 ms; 5,376 KB
testcase_10	AC	2 ms; 5,376 KB
testcase_11	AC	2 ms; 5,376 KB
testcase_12	AC	1 ms; 5,376 KB
testcase_13	AC	33 ms; 7,552 KB
testcase_14	AC	36 ms; 7,552 KB
testcase_15	AC	34 ms; 7,424 KB
testcase_16	AC	33 ms; 7,444 KB
testcase_17	AC	33 ms; 7,296 KB
testcase_18	AC	32 ms; 7,296 KB
testcase_19	AC	33 ms; 7,428 KB

ソースコード

#include "bits/stdc++.h"
using namespace std;
#ifdef _DEBUG
#include "dump.hpp"
#else
#define dump(...)
#endif

//#define int long long
#define rep(i,a,b) for(int i=(a);i<(b);i++)
#define rrep(i,a,b) for(int i=(b)-1;i>=(a);i--)
#define all(c) begin(c),end(c)
const int INF = sizeof(int) == sizeof(long long) ? 0x3f3f3f3f3f3f3f3fLL : 0x3f3f3f3f;
const int MOD = 1'000'000'007;
template<class T> bool chmax(T &a, const T &b) { if (a < b) { a = b; return true; } return false; }
template<class T> bool chmin(T &a, const T &b) { if (b < a) { a = b; return true; } return false; }

// 約数の個数の配列
// 1 から n まで
// O(n log n)
//   n + n / 2 + n / 3 + ... + n / n = n log n
vector<int> numberOfDivisors(int n) {
	vector<int> num(n + 1);
	for (int i = 1; i <= n; i++)
		for (int j = i; j <= n; j += i)
			num[j]++;
	return num;
}

signed main() {
	cin.tie(0);
	ios::sync_with_stdio(false);
	cout << fixed << setprecision(7);
	int N; double p; cin >> N >> p;
	double ans = 0;

	vector<int> num = numberOfDivisors(N);

	vector<double> pw = { 1 };
	int maxi = *max_element(all(num));
	rep(i, 0, maxi)
		pw.push_back(pw.back()*(1 - p));
	dump(pw);
	rep(i, 2, N + 1) {
		ans += pw[num[i] - 2];
	}
	cout << ans << endl;
	return 0;
}