結果

問題 No.1666 累乗数
コンテスト
ユーザー KoD
提出日時 2021-09-03 21:53:31
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 1,310 ms / 2,000 ms
コード長 2,966 bytes
コンパイル時間 3,078 ms
コンパイル使用メモリ 244,360 KB
実行使用メモリ 6,820 KB
最終ジャッジ日時 2024-12-15 12:48:38
合計ジャッジ時間 23,767 ms
ジャッジサーバーID
（参考情報） 		judge1 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 19
権限があれば一括ダウンロードができます

ソースコード

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