結果

問題 No.1666 累乗数
コンテスト
ユーザー stng
提出日時 2021-09-04 12:43:26
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,673 bytes
コンパイル時間 663 ms
コンパイル使用メモリ 82,084 KB
実行使用メモリ 77,080 KB
最終ジャッジ日時 2024-12-17 20:53:55
合計ジャッジ時間 12,994 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 1 WA * 18
権限があれば一括ダウンロードができます

ソースコード

diff # 

import math
T = int(input())
cs = [int(input()) for i in range(T)]

def sieve_of_eratosthenes(n):
    prime = [True for i in range(n+1)]
    prime[0] = False
    prime[1] = False

    sqrt_n = math.ceil(math.sqrt(n))
    for i in range(2, sqrt_n):
        if prime[i]:
            for j in range(2*i, n+1, i):
                prime[j] = False

    return prime

def make_divisors(n):
    divisors = []
    for i in range(1, int(n**0.5)+1):
        if n % i == 0:
            divisors.append(i)
            if i != n // i:
                divisors.append(n//i)

    #divisors.sort()
    return divisors

ss = sieve_of_eratosthenes(60)

kb = [0]*61
for i in range(2,61):
    tmp = make_divisors(i)
    t = 0
    for j in range(2,len(tmp)):
        if ss[tmp[j]] == True:
            t += 1
            #if tmp[j] != 2:
            #    ss[tmp[j]] = False
    kb[i] = t
#print(kb)
#exit()
kb[12] = 1

for i in range(T):
    k = cs[i]

    wa = 0#敢えて-1している（ダメな方）
    ac = 10**18#OKな方
    
    while wa+1 < ac:
        tmp = 0
        mid = (wa + ac) //2
        
        #判定書く
        for j in range(2,61):
            #kb = make_divisors(j)
            #kb = kb % 2 + 2*((kb-2)//2)
            val = int(mid**(1/j))
            plus = (val+1)**j
            if plus <= mid:
                tmp += val
                tmp -= (val)*kb[j]
            else:
                tmp += val-1
                tmp -= (val-1)*kb[j]
            #print(val,tmp,j,kb)
        #print(tmp)
        #exit()
        if tmp >= k-1:
            ac = mid
        else:
            wa = mid
        #print(mid,tmp)
            
    #print(i+1,ac)
    print(ac)
0