結果

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

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 924 ms
259,752 KB
testcase_01 AC 949 ms
259,748 KB
testcase_02 AC 962 ms
259,500 KB
testcase_03 AC 944 ms
259,624 KB
testcase_04 AC 901 ms
259,068 KB
testcase_05 AC 898 ms
259,596 KB
testcase_06 AC 915 ms
259,496 KB
testcase_07 AC 898 ms
259,372 KB
testcase_08 AC 891 ms
259,244 KB
testcase_09 AC 929 ms
259,364 KB
testcase_10 AC 1,005 ms
259,500 KB
testcase_11 AC 960 ms
259,752 KB
testcase_12 AC 923 ms
259,500 KB
testcase_13 AC 903 ms
259,732 KB
testcase_14 AC 966 ms
259,624 KB
testcase_15 AC 922 ms
259,236 KB
testcase_16 AC 932 ms
259,624 KB
testcase_17 AC 953 ms
259,476 KB
testcase_18 AC 965 ms
259,244 KB
testcase_19 AC 1,041 ms
259,616 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