結果

問題 No.1666 累乗数
ユーザー 沙耶花沙耶花
提出日時 2021-09-03 22:02:42
言語 C++17
(gcc 12.3.0 + boost 1.83.0)
結果
AC  
実行時間 1,951 ms / 2,000 ms
コード長 1,458 bytes
コンパイル時間 4,157 ms
コンパイル使用メモリ 261,456 KB
実行使用メモリ 4,380 KB
最終ジャッジ日時 2023-08-21 21:52:55
合計ジャッジ時間 42,099 ms
ジャッジサーバーID
（参考情報） 		judge15 / judge13
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 49 ms
4,376 KB
testcase_01 AC 1,923 ms
4,380 KB
testcase_02 AC 1,924 ms
4,380 KB
testcase_03 AC 1,926 ms
4,376 KB
testcase_04 AC 1,925 ms
4,376 KB
testcase_05 AC 1,944 ms
4,380 KB
testcase_06 AC 1,946 ms
4,380 KB
testcase_07 AC 1,943 ms
4,376 KB
testcase_08 AC 1,946 ms
4,376 KB
testcase_09 AC 1,926 ms
4,376 KB
testcase_10 AC 1,935 ms
4,376 KB
testcase_11 AC 1,935 ms
4,376 KB
testcase_12 AC 1,934 ms
4,380 KB
testcase_13 AC 1,950 ms
4,376 KB
testcase_14 AC 1,951 ms
4,376 KB
testcase_15 AC 1,951 ms
4,376 KB
testcase_16 AC 1,950 ms
4,376 KB
testcase_17 AC 1,943 ms
4,376 KB
testcase_18 AC 1,944 ms
4,380 KB
testcase_19 AC 1,948 ms
4,376 KB
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint998244353;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000001

long long get(long long n,long long b){
	
	long long ok = 1,ng = Inf;
	if(b>=20){
		ng = 10;
	}
	else if(b>=10){
		ng = 70;
	}
	while(ng-ok>1LL){
		long long mid = (ok+ng)/2;
		long long t = n;
		rep(i,b){
			t /= mid;
		}
		if(t==0)ng = mid;
		else ok = mid;
	}
	return ok;
}

int main(){
	
	vector<long long> p;
	
	for(int i=2;i<60;i++){
		int c = 0;
		for(int j=1;j<=i;j++){
			if((i%j==0))c++;
		}
		if(c==2)p.push_back(i);
	}
	int m = p.size();
	vector<int> mul(1<<m,0);
	vector<int> si(1<<m,0);
	for(int i=1;i<(1<<m);i++){
		mul[i] = 1;
		si[i] = -1;
		rep(j,m){
			if((i>>j)&1){
				si[i] *= -1;
				mul[i] = min((long long)60,mul[i]*p[j]);
			}
		}
	}
	
	vector<int> mm,ss;
	rep(i,mul.size()){
		if(si[i]!=0 && mul[i]!=60){
			mm.push_back(mul[i]);
			ss.push_back(si[i]);
		}
	}
	//cout<<mm.size()<<endl;
	int Sum = 0;
	rep(i,mul.size()){
		if(mul[i]==60)Sum += si[i];
	}
	
	int _t;
	cin>>_t;
	
	rep(_,_t){
		long long K;
		cin>>K;
		
		
		
		long long ok = 0,ng = 1000000000000000000;
		
		while(ng-ok>1LL){
			long long mid = (ok+ng)/2;
			long long sum = 0LL;
			sum += Sum;
			rep(i,mm.size()){
				
				sum += get(mid,mm[i])*ss[i];
				
			}
			
			if(sum >= K)ng = mid;
			else ok = mid;
		}
		cout<<ok+1<<endl;
	}
	
	return 0;
}
0