結果
問題 | No.1140 EXPotentiaLLL! |
ユーザー |
|
提出日時 | 2025-01-29 12:31:08 |
言語 | C++23 (gcc 13.3.0 + boost 1.87.0) |
結果 |
AC
|
実行時間 | 947 ms / 2,000 ms |
コード長 | 1,541 bytes |
コンパイル時間 | 5,525 ms |
コンパイル使用メモリ | 333,568 KB |
実行使用メモリ | 6,820 KB |
最終ジャッジ日時 | 2025-01-29 12:31:24 |
合計ジャッジ時間 | 14,661 ms |
ジャッジサーバーID (参考情報) |
judge5 / judge2 |
(要ログイン)
ファイルパターン | 結果 |
---|---|
sample | AC * 1 |
other | AC * 12 |
ソースコード
#include <bits/stdc++.h>#include <atcoder/all>using mint = atcoder::static_modint<998244353>;// using mint = atcoder::static_modint<1000000007>;using namespace std;using namespace atcoder;using ld = long double;using ll = long long;#define mp(a,b) make_pair(a,b)#define rep(i,s,n) for(int i=s; i<n; i++)const vector<int> dx{1,0,-1,0},dy{0,1,0,-1};bool is_prime(ll x){// miller rabinauto pow_mod=[&](__int128_t a,ll k,__int128_t mod){__int128_t output=1;a%=mod;while(k){if(k&1)output=(output*a)%mod;a=(a*a)%mod;k>>=1;}return output;};if(x<=1)return false;vector<ll> test;if(x<4759123141)test={2,7,61};else test={2,325,9375,28178,450775,9780504,1795265022};ll K=x-1;while(!(K&1))K>>=1;for(ll a:test){if(a>=x)return true;bool ok=false;ll k=K;__int128_t c=pow_mod(a,k,x);if(c==1)ok=true;else{while(k<=x-1){if(c==x-1){ok=true;break;}c=c*c%x;k<<=1;}}if(!ok)return false;}return true;}int main(){ios::sync_with_stdio(false);cin.tie(nullptr);int T;cin >> T;rep(i,0,T){ll a,p;cin >> a >> p;a%=p;if(is_prime(p)){if(a)cout << 1 << "\n";else cout << 0 << "\n";}else cout << -1 << "\n";}}