結果

問題 No.1666 累乗数
ユーザー とりゐとりゐ
提出日時 2021-09-01 23:53:40
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 392 ms / 2,000 ms
コード長 937 bytes
コンパイル時間 446 ms
コンパイル使用メモリ 81,988 KB
実行使用メモリ 238,040 KB
最終ジャッジ日時 2024-12-15 09:02:47
合計ジャッジ時間 9,074 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 329 ms
236,372 KB
testcase_01 AC 371 ms
237,656 KB
testcase_02 AC 390 ms
236,412 KB
testcase_03 AC 385 ms
235,880 KB
testcase_04 AC 376 ms
236,340 KB
testcase_05 AC 365 ms
237,232 KB
testcase_06 AC 372 ms
236,852 KB
testcase_07 AC 372 ms
236,384 KB
testcase_08 AC 373 ms
236,740 KB
testcase_09 AC 372 ms
236,484 KB
testcase_10 AC 378 ms
236,388 KB
testcase_11 AC 392 ms
238,040 KB
testcase_12 AC 379 ms
236,352 KB
testcase_13 AC 379 ms
237,356 KB
testcase_14 AC 380 ms
237,776 KB
testcase_15 AC 375 ms
236,752 KB
testcase_16 AC 377 ms
236,492 KB
testcase_17 AC 372 ms
236,484 KB
testcase_18 AC 374 ms
236,232 KB
testcase_19 AC 371 ms
236,356 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

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=0,10**18
  # f(n)>=k となる最小の n を二分探索で求める
  while abs(ok-ng)>1:
    mid=(ok+ng)//2
    if f(mid)>=n:
      ok=mid
    else:
      ng=mid
  print(ok)
0