#include using namespace std; #define ALL(x) x.begin(),x.end() #define rep(i,n) for(int i=0;i<(n);i++) #define debug(v) cout<<#v<<":";for(auto x:v){cout<bool chmax(T &a,const T &b){if(abool chmin(T &a,const T &b){if(b ostream &operator<<(ostream &os,const pair&p){ os< ostream &operator<<(ostream &os,const vector&v){ for(int i=0;i<(int)v.size();i++) os<0){ res.push_back('0'+num%10); num/=10; } if(neg) res.push_back('-'); reverse(begin(res),end(res)); os<>(istream &is,__int128 &n){ string s; is>>s; int idx=0; bool neg=false; if(s[0]=='-')neg=true,idx++; n=0; for(;idx<(int)s.size();idx++){ n=10*n+s[idx]-'0'; } if(neg) n=-n; return is; } __int128 abs(__int128 x){ return x<0?-x:x; } __int128 gcd(__int128 a,__int128 b){ if(a==0) return b; if(b==0) return a; return b>0?gcd(b,a%b):a; } namespace fastprime{ int128 pow_mod(int128 a,int128 n,int128 m){ int128 ret=1; while(n){ if(n&1)ret=ret*a%m; a=a*a%m; n>>=1; } return ret; } bool miller_rabin(int128 n){ if(n<=1) return false; if(n==2) return true; if(n%2==0) return false; int128 d=n-1,s=0; while(!(d&1))d>>=1,s++;// n-1=2^s*d for(int128 a:{2,3,5,7,11,13,17,19,23,29,31,37}){ if(n<=a) break; int128 y=pow_mod(a,d,n);// y=a^d (mod n) int128 t; for(t=d;t pollard_rho(int128 n){ if(n==1) return {}; int128 x=pollard_rho_single(n); if(x==n) return {x}; vector l=pollard_rho(x); vector r=pollard_rho(n/x); l.insert(l.end(),r.begin(),r.end()); return l; } } void solve(){ int128 a,p;cin>>a>>p; if(!fastprime::miller_rabin(p)){ cout<<-1<>q; while(q--) solve(); return 0; }