結果

問題 No.1666 累乗数
ユーザー KoDKoD
提出日時 2021-09-03 21:53:31
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,470 ms / 2,000 ms
コード長 2,966 bytes
コンパイル時間 3,694 ms
コンパイル使用メモリ 243,716 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-05-09 03:31:49
合計ジャッジ時間 27,063 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 17 ms
5,248 KB
testcase_01 AC 801 ms
5,376 KB
testcase_02 AC 821 ms
5,376 KB
testcase_03 AC 836 ms
5,376 KB
testcase_04 AC 844 ms
5,376 KB
testcase_05 AC 1,470 ms
5,376 KB
testcase_06 AC 1,467 ms
5,376 KB
testcase_07 AC 1,462 ms
5,376 KB
testcase_08 AC 1,465 ms
5,376 KB
testcase_09 AC 1,007 ms
5,376 KB
testcase_10 AC 1,027 ms
5,376 KB
testcase_11 AC 1,037 ms
5,376 KB
testcase_12 AC 1,037 ms
5,376 KB
testcase_13 AC 1,404 ms
5,376 KB
testcase_14 AC 1,401 ms
5,376 KB
testcase_15 AC 1,400 ms
5,376 KB
testcase_16 AC 1,398 ms
5,376 KB
testcase_17 AC 1,172 ms
5,376 KB
testcase_18 AC 1,247 ms
5,376 KB
testcase_19 AC 1,215 ms
5,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>

using i32 = std::int32_t;
using u32 = std::uint32_t;
using i64 = std::int64_t;
using u64 = std::uint64_t;
using i128 = __int128_t;
using u128 = __uint128_t;
using isize = std::ptrdiff_t;
using usize = std::size_t;

class rep {
    struct Iter {
        usize itr;
        constexpr Iter(const usize pos) noexcept : itr(pos) {}
        constexpr void operator++() noexcept { ++itr; }
        constexpr bool operator!=(const Iter& other) const noexcept { return itr != other.itr; }
        constexpr usize operator*() const noexcept { return itr; }
    };
    const Iter first, last;

  public:
    explicit constexpr rep(const usize first, const usize last) noexcept : first(first), last(std::max(first, last)) {}
    constexpr Iter begin() const noexcept { return first; }
    constexpr Iter end() const noexcept { return last; }
};

class revrep {
    struct Iter {
        usize itr;
        constexpr Iter(const usize pos) noexcept : itr(pos) {}
        constexpr void operator++() noexcept { --itr; }
        constexpr bool operator!=(const Iter& other) const noexcept { return itr != other.itr; }
        constexpr usize operator*() const noexcept { return itr; }
    };
    const Iter first, last;

  public:
    explicit constexpr revrep(const usize first, const usize last) noexcept
        : first(last - 1), last(std::min(first, last) - 1) {}
    constexpr Iter begin() const noexcept { return first; }
    constexpr Iter end() const noexcept { return last; }
};

template <class T, T Div = 2> constexpr T INFTY = std::numeric_limits<T>::max() / Div;

template <class T> using Vec = std::vector<T>;

void main_() {
    usize K;
    std::cin >> K;
    const auto check = [&](const u64 thres) {
        std::array<u64, 60> cnt{};
        for (const auto e : rep(2, 60)) {
            u64 ok = 1, ng = thres + 1;
            while (ng - ok > 1) {
                const auto md = (ok + ng) / 2;
                u64 cur = 1;
                bool f = true;
                for (const auto k : rep(0, e)) {
                    if (cur > thres / md) {
                        f = false;
                        break;
                    } else {
                        cur *= md;
                    }
                }
                (f ? ok : ng) = md;
            }
            cnt[e] = ok - 1;
        }
        for (const auto i : revrep(2, 60)) {
            for (usize j = i + i; j < 60; j += i) {
                cnt[i] -= cnt[j];
            }
        }
        u64 sum = 1;
        for (const auto i : rep(2, 60)) {
            sum += cnt[i];
        }
        return sum >= K;
    };
    u64 ok = INFTY<u64>, ng = 0;
    while (ok - ng > 1) {
        const auto md = (ok + ng) / 2;
        (check(md) ? ok : ng) = md;
    }
    std::cout << ok << '\n';
}

int main() {
    std::ios_base::sync_with_stdio(false);
    std::cin.tie(nullptr);
    usize cases;
    std::cin >> cases;
    while (cases--) main_();
    return 0;
}
0