結果
| 問題 |
No.1666 累乗数
|
| コンテスト | |
| ユーザー |
とりゐ
|
| 提出日時 | 2021-09-01 21:50:39 |
| 言語 | PyPy3 (7.3.15) |
| 結果 |
AC
|
| 実行時間 | 391 ms / 2,000 ms |
| コード長 | 875 bytes |
| コンパイル時間 | 397 ms |
| コンパイル使用メモリ | 82,040 KB |
| 実行使用メモリ | 238,468 KB |
| 最終ジャッジ日時 | 2024-12-15 09:03:19 |
| 合計ジャッジ時間 | 8,911 ms |
|
ジャッジサーバーID (参考情報) |
judge1 / judge5 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 19 |
ソースコード
import bisect
#floor((10**18)**(1/n)) を求めておく
lim=[0]*60
for i in range(2,60):
lim[i]=int((10**18+10000)**(1/i))
#平方数を記録
square=[0]*1000001
for i in range(1001):
square[i*i]=1
#10^18 以下の平方数でない累乗数を列挙する
s=[]
for i in range(3,60,2):
j=1
while j<=lim[i]:
if square[j]:
pass
else:
s.append(j**i)
j+=1
s=sorted(list(set(s)))
def f(n):
ng,ok=10**9+1,0
#floor(sqrt(n)) を二分探索で求める
while abs(ok-ng)>1:
mid=(ok+ng)//2
if mid**2<=n:
ok=mid
else:
ng=mid
#s の要素のうち n 以下のものの個数を二分探索で求める
m=bisect.bisect(s,n)
return ok+m
t=int(input())
for _ in range(t):
n=int(input())
ng,ok=10**18+1,0
while abs(ok-ng)>1:
mid=(ok+ng)//2
if f(mid)>=n:
ng=mid
else:
ok=mid
print(ok+1)
とりゐ