結果
| 問題 |
No.1666 累乗数
|
| コンテスト | |
| ユーザー |
ygd.
|
| 提出日時 | 2021-09-03 23:18:13 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 1,479 bytes |
| コンパイル時間 | 433 ms |
| コンパイル使用メモリ | 82,248 KB |
| 実行使用メモリ | 85,528 KB |
| 最終ジャッジ日時 | 2024-12-15 17:53:17 |
| 合計ジャッジ時間 | 39,208 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | WA * 12 TLE * 7 |
ソースコード
import sys
input = sys.stdin.readline
import math
def main():
B = set([])
MAX = pow(10,6) + 10
L = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59]
for i in L:
for x in range(2,pow(10,3) + 10):
v = pow(x, i)
if v <= MAX:
B.add(v)
#print(B)
B = list(B); B.sort()
#A[x]: x以下の累乗数の値(ただし1を除く)
A = [0]*MAX
#A[1] = 1
for i in range(2,MAX):
if i in B:
A[i] = A[i-1] + 1
else:
A[i] = A[i-1]
#print(A[:100])
#print(check(10))
T = int(input())
for _ in range(T):
K = int(input())
ok = pow(10,18)
ng = 0
while abs(ok - ng) > 1:
mid = (ok + ng)//2
if check(mid,A) >= K :
ok = mid
else:
ng = mid
print(ok)
def check(x,A):
ret = 1 #1の分
L = [2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, 41, 43, 47, 53, 59]
for i in L: #i乗根以下の素数
if pow(2,i) > x: break
ok = 2
ng = math.ceil(x ** (1/i)) + 10
while abs(ok - ng) > 1:
mid = (ok + ng) // 2
if pow(mid,i) <= x:
ok = mid
else:
ng = mid
#print(i,ok)
if i == 2:
ret += ok - 1
else:
ret += ok - 1 - A[ok]
return ret
if __name__ == "__main__":
main()
ygd.