結果

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

ソースコード

diff #
プレゼンテーションモードにする

#define _CRT_SECURE_NO_WARNINGS
#pragma GCC target("avx")
#pragma GCC optimize("O3")
#pragma GCC optimize("unroll-loops")
#include <algorithm>
#include <bitset>
#include <cassert>
#include <cfloat>
#include <climits>
#include <cmath>
#include <complex>
#include <ctime>
#include <deque>
#include <fstream>
#include <functional>
#include <iomanip>
#include <iostream>
#include <list>
#include <map>
#include <queue>
#include <random>
#include <set>
#include <stack>
#include <string>
#include <unordered_map>
#include <unordered_set>
#include <vector>
#define rep(i, n) for(int i=0;i<(n);i++)
#define REP(i, n) for(int i=1;i<=(n);i++)
#define all(V) V.begin(),V.end()
typedef long long lint;
typedef std::pair<lint, lint> P;
constexpr int INF = INT_MAX / 10;
constexpr lint LINF = LLONG_MAX / 10;
constexpr double eps = 1e-9;
template<class T>
class prique :public std::priority_queue<T, std::vector<T>, std::greater<T>> {};
template <class T, class U>
inline bool chmax(T& lhs, const U& rhs) {
    if (lhs < rhs) {
        lhs = rhs;
        return 1;
    }
    return 0;
}
template <class T, class U>
inline bool chmin(T& lhs, const U& rhs) {
    if (lhs > rhs) {
        lhs = rhs;
        return 1;
    }
    return 0;
}
inline lint gcd(lint a, lint b) {
    while (b) {
        lint c = a;
        a = b; b = c % b;
    }
    return a;
}
inline lint lcm(lint a, lint b) {
    return a / gcd(a, b) * b;
}
bool isprime(lint n) {
    if (n == 1)return false;
    for (int i = 2; i * i <= n; i++) {
        if (n % i == 0)return false;
    }
    return true;
}
lint mypow(lint a, lint b) {
    if (!b)return 1;
    if (b & 1)return mypow(a, b - 1) * a;
    lint memo = mypow(a, b >> 1);
    return memo * memo;
}
lint modpow(lint a, lint b, lint m) {
    if (!b)return 1;
    if (b & 1)return modpow(a, b - 1, m) * a % m;
    lint memo = modpow(a, b >> 1, m);
    return memo * memo % m;
}
void printArray(std::vector<int>& vec) {
    rep(i, vec.size() - 1)std::cout << vec[i] << " ";
    std::cout << vec.back() << std::endl;
}
template<typename T>
void printArray(T l, T r) {
    T rprev = r;
    rprev--;
    for (T i = l; i != rprev; i++) {
        std::cout << *i << " ";
    }
    std::cout << *rprev << std::endl;
}
std::string to_string(std::vector<int>& vec) {
    std::string res = "[";
    rep(i, vec.size() - 1)res += std::to_string(vec[i]) + ", ";
    res += std::to_string(vec.back()) + "]";
    return res;
}
int t;
lint n;
bool used[100000];
int main() {
    std::cin >> t;
    rep(i, t) {
        std::cin >> n;
        lint ans = n * n + (n - 1) * n;
        std::fill(used, used + 100000, false);
        for (lint j = 2; j * j <= n; j++) {
            if (used[j])continue;
            lint m = j;
            while (m * m <= n) {
                used[m] = true;
                m *= j;
            }
            for (lint k = 1; mypow(j, k) <= n; k++) {
                for (lint l = 1; mypow(j, l) <= n; l++) {
                    if (k == l)continue;
                    ans += n / (std::max(k, l) / gcd(k, l));
                }
            }
        }
        std::cout << ans << std::endl;
    }
    return 0;
}
הההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההההה
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
0