結果
| 問題 |
No.1666 累乗数
|
| コンテスト | |
| ユーザー |
|
| 提出日時 | 2021-09-03 22:32:09 |
| 言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
| 結果 |
WA
|
| 実行時間 | - |
| コード長 | 2,326 bytes |
| コンパイル時間 | 1,878 ms |
| コンパイル使用メモリ | 182,512 KB |
| 実行使用メモリ | 74,224 KB |
| 最終ジャッジ日時 | 2024-12-15 15:03:27 |
| 合計ジャッジ時間 | 18,516 ms |
|
ジャッジサーバーID (参考情報) |
judge3 / judge1 |
(要ログイン)
| ファイルパターン | 結果 |
|---|---|
| sample | AC * 1 |
| other | AC * 4 WA * 15 |
ソースコード
#include <bits/stdc++.h>
//#include<atcoder/all>
using namespace std;
//using namespace atcoder;
using ll = long long;
#define all(A) A.begin(),A.end()
using vll = vector<ll>;
#define rep(i, n) for (long long i = 0; i < (long long)(n); i++)
using Graph = vector<vector<ll>>;
vll P = {3,5,7,11,13,17,19,23,29,31,37};
vector<vll> PP(P.size());
set<ll> S;
vll PPP;
bool isOK(ll index, ll key) {
ll k = 0;
ll D = sqrt(index);
D = max(0ll, D - 5);
rep(i, 11) {
if (D * D <= index) {
D++;
}
else {
D--;
break;
}
}
k += D;
auto Ds = upper_bound(all(PPP),index)-PPP.begin();
k += (Ds);
return (k >= key);
}
// 汎用的な二分探索のテンプレ
ll binary_search(ll key) {
ll ng = 0; //「index = 0」が条件を満たすこともあるので、初期値は -1
ll ok = key * key+3; // 「index = a.size()-1」が条件を満たさないこともあるので、初期値は a.size()
/* ok と ng のどちらが大きいかわからないことを考慮 */
while (abs(ok - ng) > 1) {
ll mid = (ok + ng) / 2;
cout << "";
if (isOK(mid, key)) ok = mid;
else ng = mid;
}
return ok;
}
int main() {
//set<ll> S;
S.insert(2e18);
PPP.push_back(2e18);
rep(i, P.size()) {
ll p = P[i];
ll k = 2;
while (1) {
ll d = k;
bool c = true;
rep(u, p - 1) {
if (d <= 2e18 / k) {
d *= k;
}
else {
c = false;
break;
}
}
if (c) {
ll D = sqrt(d);
bool ks = true;
rep(i, 11) {
if ((D - 5 + i) * (D - 5 + i) == d)ks = false;
}
if (ks&&!S.count(d)) {
S.insert(d);
PPP.push_back(d);
}
}
else {
//PP[i].push_back(2e18);
break;
}
k++;
d = 0;
}
}
sort(all(PPP));
ll T;
cin >> T;
rep(t, T) {
ll N;
cin >> N;
cout << binary_search(N) << endl;
}
}