結果
問題 | No.1255 ハイレーツ・オブ・ボリビアン |
ユーザー |
|
提出日時 | 2020-10-09 22:35:48 |
言語 | C++17(gcc12) (gcc 12.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 30 ms / 2,000 ms |
コード長 | 1,380 bytes |
コンパイル時間 | 3,465 ms |
コンパイル使用メモリ | 134,904 KB |
最終ジャッジ日時 | 2025-01-15 05:18:36 |
ジャッジサーバーID (参考情報) |
judge5 / judge5 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 15 |
ソースコード
#include <iostream>#include <algorithm>#include <map>#include <set>#include <queue>#include <stack>#include <numeric>#include <bitset>#include <cmath>static const int MOD = 1000000007;using ll = long long;using u32 = unsigned;using u64 = unsigned long long;using namespace std;template<class T> constexpr T INF = ::numeric_limits<T>::max()/32*15+208;template<class T>vector<T> divisor(T n){vector<T> ret;for(T i = 1; i * i <= n; i++) {if(n % i == 0) {ret.push_back(i);if(i * i != n) ret.push_back(n / i);}}sort(begin(ret), end(ret));return(ret);}int euler(int n){if(!n) return 0;int ans = n;for (int i = 2; i*i <= n; ++i) {if(n%i == 0){ans -= ans/i;while(n%i == 0) n /= i;}}if(n > 1) ans -= ans/n;return ans;}template <class T>T pow_ (T x, T n, T M){uint64_t u = 1, xx = x;while (n > 0){if (n&1) u = u * xx % M;xx = xx * xx % M;n >>= 1;}return static_cast<T>(u);};int solve(int n){if(n == 1) return 1;for(auto &&d : divisor(euler(2*n-1))){if(pow_(2, d, 2*n-1) == 1) return d;}return -1;}int main() {int t;cin >> t;while(t--){int n;cin >> n;cout << solve(n) << "\n";}return 0;}