結果

問題 No.1022 Power Equation
ユーザー QCFiumQCFium
提出日時 2020-04-10 21:54:36
言語 C++14
(gcc 13.3.0 + boost 1.87.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
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#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;
}
0