結果

問題 No.1514 Squared Matching
ユーザー nok0nok0
提出日時 2021-04-29 20:18:23
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 2,193 ms / 4,000 ms
コード長 1,587 bytes
コンパイル時間 2,318 ms
コンパイル使用メモリ 208,340 KB
実行使用メモリ 216,764 KB
最終ジャッジ日時 2024-07-17 21:38:05
合計ジャッジ時間 38,564 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
6,812 KB
testcase_01 AC 452 ms
216,764 KB
testcase_02 AC 2 ms
6,940 KB
testcase_03 AC 2 ms
6,940 KB
testcase_04 AC 2 ms
6,940 KB
testcase_05 AC 3 ms
6,940 KB
testcase_06 AC 103 ms
6,940 KB
testcase_07 AC 1,483 ms
45,704 KB
testcase_08 AC 1,803 ms
209,180 KB
testcase_09 AC 734 ms
184,672 KB
testcase_10 AC 1,374 ms
199,332 KB
testcase_11 AC 1,705 ms
206,648 KB
testcase_12 AC 1,947 ms
210,932 KB
testcase_13 AC 2,038 ms
213,952 KB
testcase_14 AC 2,082 ms
215,380 KB
testcase_15 AC 2,151 ms
214,980 KB
testcase_16 AC 2,193 ms
215,720 KB
testcase_17 AC 2,131 ms
215,648 KB
testcase_18 AC 2,128 ms
215,232 KB
testcase_19 AC 2,107 ms
215,396 KB
testcase_20 AC 2,099 ms
216,320 KB
testcase_21 AC 2,137 ms
216,064 KB
testcase_22 AC 1,510 ms
46,824 KB
testcase_23 AC 1,737 ms
90,536 KB
testcase_24 AC 1,896 ms
129,312 KB
testcase_25 AC 2,003 ms
177,416 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;

#pragma region Math Osa_k
struct osa_k {
private:
	std::vector<int> spf, pr;

public:
	osa_k() = default;

	osa_k(int MAX) : spf(MAX + 1) {
		for(int i = 2; i <= MAX; i++) {
			if(spf[i] == 0) {
				spf[i] = i;
				pr.push_back(i);
			}
			for(int j = 0; j < pr.size() and pr[j] <= spf[i] and i * pr[j] <= MAX; j++)
				spf[i * pr[j]] = pr[j];
		}
	}

	std::vector<std::pair<int, int>> prime_factorization(int n) {
		std::vector<std::pair<int, int>> divisor;
		if(n == 1) return divisor;
		int before = spf[n], cnt = 0;
		while(n > 1) {
			if(spf[n] == before) {
				cnt++;
				n /= spf[n];
			} else {
				divisor.emplace_back(before, cnt);
				before = spf[n];
				cnt = 1;
				n /= spf[n];
			}
		}
		divisor.emplace_back(before, cnt);
		return divisor;
	}

	int smallestprimefactor(const int n) const {
		return spf[n];
	}

	bool is_prime(const int n) const {
		return n == spf[n];
	}
};
#pragma endregion

int embed[6] = {0, 99371160, 207162592, 318140852, 431159896, 545761220};

int n, res;
int main() {
	scanf("%d", &n);
	osa_k pf(n);

	/*
	愚直
	for(int a = 1; a <= n; a++) {
		auto v = pf.prime_factorization(a);
		int b = 1;
		for(auto [val, cnt] : v)
			if(cnt & 1) b *= val;
		res += sqrt(n / b);
	}
	printf("%d\n", res);
	*/

	int base = n / 10000000 * 10000000;
	res = embed[n / 10000000];
	for(int a = base + 1; a <= n; a++) {
		auto v = pf.prime_factorization(a);
		int b = 1;
		for(auto [val, cnt] : v)
			if(cnt & 1) b *= val;
		int d = sqrt(n / b);
		res += d + min(d, (int)sqrt(base / b));
	}

	printf("%d\n", res);
}
0