#include #include #include #include long long pow(long long a, long long n, long long p) { long long r=1; for (;n>0;n>>=1,a=a*a%p)if(n%2==1)r=r*a%p; return r; } int cnt(long long a,long long base,long long p) { int ret=0; while (a!=1) { a=pow(a,base,p); ++ret; } return ret; } long long inv(long long a, long long p) { a%=p; long long u=1, v=0; long long b=p; while (b>0) { long long q=a/b; a%=b; u-=v*q%p; u=(u%p+p)%p; { u^=v;v^=u;u^=v; a^=b;b^=a;a^=b; } } return u<0?u+p:u; } long long gcd(long long a,long long b) { return a==0?b:gcd(b%a,a); } long long peth_root(long long a,long long p,int e,long long mod) { long long q=mod-1; int s=0; while (q%p==0) { q/=p; ++s; } long long pe=pow(p,e,mod); long long ans=pow(a,((pe-1)*inv(q,pe)%pe*q+1)/pe,mod); long long c=2; while (pow(c,(mod-1)/p,mod)==1) ++c; c=pow(c,q,mod); std::map map; long long add=1; int v=(int)std::sqrt((double)(s-e)*p)+1; long long mul=pow(c,v*pow(p,s-1,mod-1)%(mod-1),mod); for (int i=0;i<=v;++i) { map[add]=i; add=add*mul%mod; } mul=inv(pow(c,pow(p,s-1,mod-1),mod),mod); for (int i=e;i0) { long long b=peth_root(a,div,sz,p); a=b; } } if (g>1) a=peth_root(a,g,1,p); return a; } int main() { int T; scanf("%d",&T); for(int t=0;t