結果
| 問題 | No.1666 累乗数 |
| コンテスト | |
| ユーザー |
とりゐ
|
| 提出日時 | 2021-09-01 21:50:39 |
| 言語 | PyPy3 (7.3.17) |
| 結果 |
AC
|
| 実行時間 | 254 ms / 2,000 ms |
| コード長 | 875 bytes |
| 記録 | |
| コンパイル時間 | 287 ms |
| コンパイル使用メモリ | 85,556 KB |
| 実行使用メモリ | 280,832 KB |
| 最終ジャッジ日時 | 2026-05-25 22:02:17 |
| 合計ジャッジ時間 | 5,727 ms |
|
ジャッジサーバーID (参考情報) |
judge1_0 / judge3_0 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| 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)
とりゐ