#include using namespace std; typedef signed long long ll; #undef _P #define _P(...) (void)printf(__VA_ARGS__) #define FOR(x,to) for(x=0;x<(to);x++) #define FORR(x,arr) for(auto& x:arr) #define ITR(x,c) for(__typeof(c.begin()) x=c.begin();x!=c.end();x++) #define ALL(a) (a.begin()),(a.end()) #define ZERO(a) memset(a,0,sizeof(a)) #define MINUS(a) memset(a,0xff,sizeof(a)) //------------------------------------------------------- int T; ll mo; const int prime_max = 33000; vector prime; int NP,divp[prime_max]; void cprime() { if(NP) return; for(int i=2;i=i&&j>=1; return r; } int totient(int v) { int ret=v; int i; for(i=0;prime[i]*prime[i]<=v;i++) if(v%prime[i]==0) { ret=ret/prime[i]*(prime[i]-1); while(v%prime[i]==0) v/=prime[i]; } if(v>1) ret=ret/v*(v-1); return ret; } int mod_root(int p,int a) { // x^p=a mod mo vector D; int i; for(i=0;prime[i]*prime[i]<=mo-1;i++) if((mo-1)%prime[i]==0) D.push_back((mo-1)/prime[i]); int g=2; while(1) { int ng=0; FORR(d,D) if(modpow(g,d)==1) { ng=1; break; } if(ng==0) break; g++; } int cur=a; int rg=modpow(g); int mstep=sqrt(mo); unordered_map M; FOR(i,mstep+3) { M[cur]=i; cur=mulmod(cur,rg,mo); } int pg=modpow(g,mstep); int x=-1,step=0; cur=1; while(x==-1) { if(M.count(cur)) x=step+M[cur]; M[cur]=step; cur=mulmod(cur,pg,mo); step+=mstep; } // g^x=aなのでg^(p*q)=g^x=aとしてq=x/p (mod mo-1) mo-1は合成数なのでGCDで割って対応 int tmo=mo-1; int gcd=__gcd(tmo,p); if(x%gcd) return -1; tmo/=gcd; x/=gcd; p/=gcd; return modpow(g,1LL*x*modpow(p,totient(tmo)-1,tmo)%tmo); } void solve() { int i,j,k,l,r,x,y; string s; cprime(); cin>>T; while(T--) { int P,K,A; cin>>P>>K>>A; mo=P; cout<