結果

問題 No.1666 累乗数
ユーザー tamatotamato
提出日時 2021-09-03 23:24:13
言語 PyPy3
(7.3.15)
結果
AC  
実行時間 998 ms / 2,000 ms
コード長 1,267 bytes
コンパイル時間 188 ms
コンパイル使用メモリ 82,696 KB
実行使用メモリ 259,904 KB
最終ジャッジ日時 2024-05-09 07:23:52
合計ジャッジ時間 21,212 ms
ジャッジサーバーID
（参考情報） 		judge4 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 911 ms
259,412 KB
testcase_01 AC 966 ms
259,668 KB
testcase_02 AC 960 ms
259,448 KB
testcase_03 AC 964 ms
259,408 KB
testcase_04 AC 998 ms
259,648 KB
testcase_05 AC 975 ms
259,904 KB
testcase_06 AC 972 ms
259,836 KB
testcase_07 AC 966 ms
259,492 KB
testcase_08 AC 952 ms
259,404 KB
testcase_09 AC 959 ms
259,336 KB
testcase_10 AC 976 ms
259,460 KB
testcase_11 AC 964 ms
259,340 KB
testcase_12 AC 940 ms
259,464 KB
testcase_13 AC 944 ms
259,676 KB
testcase_14 AC 933 ms
259,344 KB
testcase_15 AC 934 ms
259,740 KB
testcase_16 AC 924 ms
259,876 KB
testcase_17 AC 944 ms
259,900 KB
testcase_18 AC 974 ms
259,620 KB
testcase_19 AC 965 ms
259,904 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

mod = 1000000007
eps = 10**-9


def main():
    import sys
    from bisect import bisect_left
    from math import sqrt, floor
    input = sys.stdin.readline

    def square_root(n):
        r = floor(sqrt(n))
        for x in range(3, -4, -1):
            if (r+x) ** 2 <= n:
                return r+x

    M = 15 * 10 ** 5
    MAX = 2 * 10 ** 18
    S_set = {1}
    for a in range(2, M+1):
        cur = a
        for b in range(2, 70):
            cur *= a
            if cur > MAX:
                break
            S_set.add(cur)
    S = sorted(list(S_set))
    kk = bisect_left(S, M+1)

    for _ in range(int(input())):
        K = int(input())
        ok = MAX
        ng = 0
        mid = (ok + ng) // 2
        while ok - ng > 1:
            r = square_root(mid)
            if r <= M:
                j = bisect_left(S, mid+1)
                if j >= K:
                    ok = mid
                else:
                    ng = mid
            else:
                k = bisect_left(S, r+1)
                j = bisect_left(S, mid + 1) + (r - M) - (k - kk)
                if j >= K:
                    ok = mid
                else:
                    ng = mid
            mid = (ok + ng) // 2
        print(ok)


if __name__ == '__main__':
    main()
0