結果

問題 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;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0