#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; ll modpow(ll a, ll n = mo-2, ll m=mo) { ll r=1; while(n) r=r*((n%2)?a:1)%m,a=a*a%m,n>>=1; return r; } int totient(int v) { int ret=v; for(int i=2;i*i<=v;i++) if(v%i==0) { ret=ret/i*(i-1); while(v%i==0) v/=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; for(int i=2;i*i<=mo-1;i++) if((mo-1)%i==0) D.push_back(i),D.push_back((mo-1)/i); int g=2; while(1) { int ng=0; FORR(d,D) if(modpow(g,d)==1) ng=1; if(ng==0) break; g++; } ll cur=a; int rg=modpow(g); int mstep=sqrt(mo); map M; int i; FOR(i,mstep+3) { M[cur]=i; cur=cur*rg%mo; } ll 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=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; cin>>T; while(T--) { int P,K,A; cin>>P>>K>>A; mo=P; cout<