結果

問題 No.1022 Power Equation
ユーザー ayaoniayaoni
提出日時 2020-11-21 13:44:00
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 229 ms / 2,000 ms
コード長 872 bytes
コンパイル時間 159 ms
コンパイル使用メモリ 82,128 KB
実行使用メモリ 80,404 KB
最終ジャッジ日時 2024-07-23 15:36:32
合計ジャッジ時間 1,903 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 41 ms
53,520 KB
testcase_01 AC 42 ms
53,884 KB
testcase_02 AC 54 ms
64,120 KB
testcase_03 AC 66 ms
70,672 KB
testcase_04 AC 144 ms
80,028 KB
testcase_05 AC 216 ms
79,892 KB
testcase_06 AC 211 ms
80,068 KB
testcase_07 AC 229 ms
80,404 KB
testcase_08 AC 105 ms
79,544 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

import sys
from math import gcd
def I(): return int(sys.stdin.readline().rstrip())


def lcm(a,b):
    return a*b//gcd(a,b)


def f(N):
    res = 2*N**2-N
    M = int(N**.5)
    flag = [1]*(M+1)  # flag[i] == 1 ⇔ iが冪乗数でない
    for i in range(2,M+1):
        if flag[i] == 0:
            continue
        x = i**2
        while x <= M:
            flag[x] = 0
            x *= i
    for a in range(2,M+1):
        if flag[a] == 0:
            continue
        e = 2
        while a**e <= N:
            res += 2*(N//e)
            e += 1
    for aa in range(2,M+1):
        if flag[aa] == 0:
            continue
        m = 2
        while aa**m <= N:
            e = m+1
            while aa**e <= N:
                res += 2*(N*m//lcm(m,e))
                e += 1
            m += 1
    return res


T = I()
for _ in range(T):
    N = I()
    print(f(N))
0