結果
| 問題 | 
                            No.1022 Power Equation
                             | 
                    
| コンテスト | |
| ユーザー | 
                             nrvft
                         | 
                    
| 提出日時 | 2020-04-11 19:18:34 | 
| 言語 | C++17  (gcc 13.3.0 + boost 1.87.0)  | 
                    
| 結果 | 
                             
                                AC
                                 
                             
                            
                         | 
                    
| 実行時間 | 27 ms / 2,000 ms | 
| コード長 | 1,711 bytes | 
| コンパイル時間 | 1,879 ms | 
| コンパイル使用メモリ | 193,792 KB | 
| 最終ジャッジ日時 | 2025-01-09 17:26:53 | 
| 
                            ジャッジサーバーID (参考情報)  | 
                        judge3 / judge1 | 
(要ログイン)
| ファイルパターン | 結果 | 
|---|---|
| sample | AC * 1 | 
| other | AC * 8 | 
ソースコード
#include <bits/stdc++.h>
using namespace std;
using ll = long long;
using ld = long double;
using vl = vector<ll>;
template<class T> using vc = vector<T>;
template<class T> using vvc = vector<vector<T>>;
const ll MOD = 1e9 + 7;
// const ll MOD = 998244353;
const ll INF = 5e18;
const ld EPS = 1e-11;
const ld PI = acos(-1.0L);
#define eb emplace_back
#define all(x) (x).begin(), (x).end()
#define rep(i, n) for (ll i = 0; i < (n); i++)
#define repr(i, n) for (ll i = (n)-1; i >= 0; i--)
#define repe(i, l, r) for (ll i = (l); i < (r); i++)
#define reper(i, l, r) for (ll i = (r)-1; i >= (l); i--)
template<class T> bool chmax(T &a, const T &b) { if (a<b) { a=b; return 1; } return 0; }
template<class T> bool chmin(T &a, const T &b) { if (b<a) { a=b; return 1; } return 0; }
#ifdef DEBUG
template <class T, class N> void verr(const T& a, const N& n) { rep(i, n) cout << a[i] << " "; cout << "\n"; }
ll dbgt = 1; void err() { cout << "passed " << dbgt++ << "\n"; }
template<class H, class... T> void err(H&& h,T&&... t){ cout<< h << (sizeof...(t)?" ":"\n"); if(sizeof...(t)>0) err(forward<T>(t)...); }
#endif
void init() {cin.tie(0);ios::sync_with_stdio(false);cout << fixed << setprecision(15);}
//--------------------------------------------------------------------------------//
int main() {
    init();
    ll t;
    cin >> t;
    while(t--){
        ll N;
        cin >> N;
        ll ans = N * (2 * N - 1);
        repe(a,2,(ll)sqrt(N)+1){
            ll c = a * a, cp = 2;
            while(c <= N){
                repe(b,1,cp){
                    if (gcd(b, cp) == 1) ans += (N / cp) * 2;
                }
                cp++, c *= a;
            }
        }
        cout << ans << endl;
    }
}
            
            
            
        
            
nrvft