結果

問題 No.1022 Power Equation
コンテスト
ユーザー kcvlex
提出日時 2020-04-10 22:21:19
言語 C++14
(gcc 13.3.0 + boost 1.87.0)
結果
AC  
実行時間 10 ms / 2,000 ms
コード長 3,578 bytes
コンパイル時間 1,110 ms
コンパイル使用メモリ 124,128 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-15 20:52:27
合計ジャッジ時間 1,661 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge3
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <limits>
#include <initializer_list>
#include <utility>
#include <bitset>
#include <tuple>
#include <type_traits>
#include <functional>
#include <string>
#include <array>
#include <deque>
#include <list>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <set>
#include <unordered_map>
#include <unordered_set>
#include <iterator>
#include <algorithm>
#include <complex>
#include <random>
#include <numeric>
#include <iostream>
#include <iomanip>
#include <sstream>
#include <regex>
#include <cassert>
#include <cstddef>
#define endl codeforces
#define ALL(v) std::begin(v), std::end(v)
#define ALLR(v) std::rbegin(v), std::rend(v)
using ll = std::int64_t;
using ull = std::uint64_t;
using pii = std::pair<int, int>;
using tii = std::tuple<int, int, int>;
using pll = std::pair<ll, ll>;
using tll = std::tuple<ll, ll, ll>;
template <typename T> using vec = std::vector<T>;
template <typename T> using vvec = vec<vec<T>>;
template <typename T> const T& var_min(const T &t) { return t; }
template <typename T> const T& var_max(const T &t) { return t; }
template <typename T, typename... Tail> const T& var_min(const T &t, const Tail&... tail) { return std::min(t, var_min(tail...)); }
template <typename T, typename... Tail> const T& var_max(const T &t, const Tail&... tail) { return std::max(t, var_max(tail...)); }
template <typename T, typename... Tail> void chmin(T &t, const Tail&... tail) { t = var_min(t, tail...); }
template <typename T, typename... Tail> void chmax(T &t, const Tail&... tail) { t = var_max(t, tail...); }
template <typename T> const T& clamp(const T &t, const T &low, const T &high) { return std::max(low, std::min(high, t)); }
template <typename T> void chclamp(T &t, const T &low, const T &high) { return t = clamp(t, low, high); }
template <typename T> T make_v(T init) { return init; }
template <typename T, typename... Tail> auto make_v(T init, std::size_t s, Tail... tail) { auto v = std::move(make_v(init, tail...)); return vec<decltype(v)>(s, v); }
template <typename T, std::size_t Head, std::size_t ...Tail> struct multi_dem_array { using type = std::array<typename multi_dem_array<T, Tail...>::type, Head>; };
template <typename T, std::size_t Head> struct multi_dem_array<T, Head> { using type = std::array<T, Head>; };
template <typename T, std::size_t ...Args> using mdarray = typename multi_dem_array<T, Args...>::type;
namespace init__ { struct InitIO { InitIO() { std::cin.tie(nullptr); std::ios_base::sync_with_stdio(false); std::cout << std::fixed << std::setprecision(30); } } init_io; }

ll gcd(ll a, ll b) { return b ? gcd(b, a % b) : a; }

ll lcm(ll a, ll b) { return a / gcd(a, b) * b; }

ll solve() {
    ll n;
    std::cin >> n;
    ll ans = n * n;  // a = c = 1

    auto check = [&](ll a, ll b) {
        ll val = 1, cur = a;
        bool ok = true;
        while (b && ok) {
            if (b & 1) val *= cur;
            cur *= cur;
            b /= 2;
            if (b == 0) break;
            if (n < val) ok = false;
            if (n < cur) ok = false;
        }
        return (b == 0 && val <= n);
    };

    for (ll x = 1; x <= 32; x++) for (ll y = 1; y <= 32; y++) {
        ll g = gcd(x, y);
        if (2 <= g) continue;
        ll k = std::max(x, y);
        ll ok = 1, ng = n + 1;
        while (1 < std::abs(ok - ng)) {
            ll mid = (ok + ng) / 2;
            (check(mid, k) ? ok : ng) = mid;
        }
        ans += (ok - 1) * (n / k);
    }
    return ans;
}

int main() {
    ll t;
    std::cin >> t;
    while (t--) std::cout << solve() << '\n';
    return 0;
}
0