結果
問題 | No.1022 Power Equation |
ユーザー | QCFium |
提出日時 | 2020-04-10 21:54:36 |
言語 | C++14 (gcc 13.2.0 + boost 1.83.0) |
結果 |
AC
|
実行時間 | 5 ms / 2,000 ms |
コード長 | 897 bytes |
コンパイル時間 | 1,795 ms |
コンパイル使用メモリ | 166,168 KB |
実行使用メモリ | 4,376 KB |
最終ジャッジ日時 | 2023-10-14 00:25:38 |
合計ジャッジ時間 | 2,524 ms |
ジャッジサーバーID (参考情報) |
judge12 / judge13 |
テストケース
テストケース表示入力 | 結果 | 実行時間 実行使用メモリ |
---|---|---|
testcase_00 | AC | 2 ms
4,372 KB |
testcase_01 | AC | 2 ms
4,376 KB |
testcase_02 | AC | 2 ms
4,372 KB |
testcase_03 | AC | 2 ms
4,372 KB |
testcase_04 | AC | 4 ms
4,376 KB |
testcase_05 | AC | 5 ms
4,372 KB |
testcase_06 | AC | 5 ms
4,372 KB |
testcase_07 | AC | 5 ms
4,372 KB |
testcase_08 | AC | 4 ms
4,372 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; }