結果
問題 | No.1666 累乗数 |
ユーザー |
![]() |
提出日時 | 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 |
ソースコード
#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;}