結果

問題 No.1666 累乗数
ユーザー stngstng
提出日時 2021-09-04 12:43:26
言語 PyPy3
(7.3.15)
結果
WA  
実行時間 -
コード長 1,673 bytes
コンパイル時間 237 ms
コンパイル使用メモリ 82,432 KB
実行使用メモリ 77,056 KB
最終ジャッジ日時 2024-05-10 01:38:10
合計ジャッジ時間 12,754 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 66 ms
76,416 KB
testcase_01 AC 299 ms
76,544 KB
testcase_02 WA -
testcase_03 WA -
testcase_04 WA -
testcase_05 WA -
testcase_06 WA -
testcase_07 WA -
testcase_08 WA -
testcase_09 WA -
testcase_10 WA -
testcase_11 WA -
testcase_12 WA -
testcase_13 WA -
testcase_14 WA -
testcase_15 WA -
testcase_16 WA -
testcase_17 WA -
testcase_18 WA -
testcase_19 WA -
権限があれば一括ダウンロードができます

ソースコード

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