結果
問題 | No.1255 ハイレーツ・オブ・ボリビアン |
ユーザー |
![]() |
提出日時 | 2020-10-09 22:12:57 |
言語 | C++11(廃止可能性あり) (gcc 13.3.0) |
結果 |
AC
|
実行時間 | 98 ms / 2,000 ms |
コード長 | 1,556 bytes |
コンパイル時間 | 1,214 ms |
コンパイル使用メモリ | 83,380 KB |
実行使用メモリ | 6,944 KB |
最終ジャッジ日時 | 2024-07-23 08:43:02 |
合計ジャッジ時間 | 2,303 ms |
ジャッジサーバーID (参考情報) |
judge4 / judge1 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 15 |
ソースコード
#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;llint mod;llint modpow(llint a, llint n){if(n == 0) return 1;if(n % 2){return ((a%mod) * (modpow(a, n-1)%mod)) % mod;}else{return modpow((a*a)%mod, n/2) % mod;}}int main(void){ios::sync_with_stdio(0);cin.tie(0);cin >> T;for(int t = 0; t < T; t++){cin >> n;if(n == 1){cout << 1 << endl;continue;}if(n == 2){cout << 2 << endl;continue;}if(n == 3){cout << 4 << endl;continue;}n = 2*n-1;mod = n;llint res = n;for(int i = 2; i * i <= n; i++){if(n % i == 0){res = res / i * (i-1);for(; n%i == 0; n/=i);}}if(n > 1) res = res / n * (n-1);llint ans = res; vector<llint> vec;for(int i = 1; i * i <= res; i++){if(res % i == 0){vec.push_back(i);vec.push_back(res/i);}}for(int i = 0; i < vec.size(); i++){if(modpow(2, vec[i]) == 1) ans = min(ans, vec[i]);}cout << ans << endl;}return 0;}