結果
問題 |
No.1666 累乗数
|
ユーザー |
|
提出日時 | 2021-09-03 23:19:41 |
言語 | C++14 (gcc 13.3.0 + boost 1.87.0) |
結果 |
WA
|
実行時間 | - |
コード長 | 2,017 bytes |
コンパイル時間 | 1,450 ms |
コンパイル使用メモリ | 178,308 KB |
実行使用メモリ | 80,940 KB |
最終ジャッジ日時 | 2024-12-15 17:58:18 |
合計ジャッジ時間 | 11,461 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge3 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
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 = ll(sqrt(index)); 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); 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 = round(sqrt(d)); if (D*D!=d) { S.insert(d); } } else { //PP[i].push_back(2e18); break; } k++; d = 0; } } for (auto l : S) { PPP.push_back(l); } ll T; cin >> T; rep(t, T) { ll N; cin >> N; cout << binary_search(N) << endl; } }