#462027 (C++14) No.1022 Power Equation

• 問題: No.1022 Power Equation
• ユーザー: QCFium
• 提出日時: 2020-04-10 21:54:36
• 言語: C++14 (gcc 12.3.0 + boost 1.83.0)
• 結果: AC
• 実行時間: 6 ms / 2,000 ms
• コード長: 897 bytes
• コンパイル時間: 1,588 ms
• コンパイル使用メモリ: 168,388 KB
• 実行使用メモリ: 5,376 KB
• 最終ジャッジ日時: 2024-09-15 20:15:21
• 合計ジャッジ時間: 2,138 ms
• ジャッジサーバーID （参考情報）: judge3 / judge2

入力	結果	実行時間; 実行使用メモリ
testcase_00	AC	2 ms; 5,248 KB
testcase_01	AC	2 ms; 5,376 KB
testcase_02	AC	2 ms; 5,376 KB
testcase_03	AC	3 ms; 5,376 KB
testcase_04	AC	4 ms; 5,376 KB
testcase_05	AC	6 ms; 5,376 KB
testcase_06	AC	6 ms; 5,376 KB
testcase_07	AC	5 ms; 5,376 KB
testcase_08	AC	4 ms; 5,376 KB

ソースコード

#include <bits/stdc++.h>

int ri() {
	int n;
	scanf("%d", &n);
	return n;
}

std::vector<int> num;

void init(int n) {
	num.assign(30, 0);
	num[1] = n;
	for (int i = 2; i < 30; i++) {
		for (int j = 2; ; j++) {
			int cur = 1;
			bool ok = true;
			for (int k = 0; k < i; k++) {
				if ((int64_t) cur * j > n) {
					ok = false;
					break;
				}
				cur *= j;
			}
			if (ok) num[i]++;
			else break;
		}
	}
}

int gcd(int a, int b) {
	while (a && b) {
		if (a > b) a %= b;
		else b %= a;
	}
	return a + b;
}

void solve() {
	int n = ri();
	init(n);
	int64_t res = 0;
	res += (int64_t) n * n; // x, x
	res += (int64_t) n * n;
	res -= n;
	for (int i = 2; i < 30 && i <= n; i++) for (int j = 1; j < i; j++) {
		if (gcd(i, j) != 1) continue;
		res += (int64_t) num[i] * (n / i) * 2;
	}
	printf("%" PRId64 "\n", res);
}

int main() {
	int t = ri();
	for (int i = 0; i < t; i++) solve();
	return 0;
}