結果

問題 No.1022 Power Equation
コンテスト
ユーザー leaf_1415
提出日時 2020-04-10 22:13:04
言語 C++11（廃止可能性あり）
(gcc 13.3.0)
結果
AC  
実行時間 6 ms / 2,000 ms
コード長 1,376 bytes
コンパイル時間 1,595 ms
コンパイル使用メモリ 81,716 KB
実行使用メモリ 5,376 KB
最終ジャッジ日時 2024-09-15 20:38:38
合計ジャッジ時間 2,142 ms
ジャッジサーバーID
（参考情報） 		judge2 / judge4
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 1
other AC * 8
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <iostream>
#include <cstdio>
#include <cmath>
#include <ctime>
#include <cstdlib>
#include <cassert>
#include <vector>
#include <list>
#include <stack>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <bitset>
#include <string>
#include <algorithm>
#include <utility>
#define llint long long
#define inf 1e18
#define rep(x, s, t) for(llint (x) = (s); (x) < (t); (x)++)
#define Rep(x, s, t) for(llint (x) = (s); (x) <= (t); (x)++)
#define chmin(x, y) (x) = min((x), (y))
#define chmax(x, y) (x) = max((x), (y))

using namespace std;
typedef pair<llint, llint> P;

llint T;
llint n;
vector<llint> vec[35];

llint gcd(llint a, llint b)
{
	if(b == 0) return a;
	return gcd(b, a%b);
}

llint get(llint e, llint n)
{
	if(e >= 2) return upper_bound(vec[e].begin(), vec[e].end(), n) - vec[e].begin();
	else return n;
}

int main(void)
{
	ios::sync_with_stdio(0);
	cin.tie(0);
	
	for(int i = 1; i <= 1000000; i++){
		llint mul = i;
		for(int j = 2; j < 35; j++){
			mul *= i;
			if(mul > 1e9+7) break;
			vec[j].push_back(mul);
		}
	}
	
	cin >> T;
	for(int t = 0; t < T; t++){
		cin >> n;
		
		llint ans = n * n;
		for(int i = 1; i <= min(32LL, n); i++){
			for(int j = 1; j <= min(32LL, n); j++){
				if(gcd(i, j) > 1) continue;
				llint num = n / max(i, j);
				ans += num * (get(max(i, j), n) - 1);
			}
		}
		cout << ans << endl;
	}
	
	return 0;
}
0