#include #include #include #include using namespace std; typedef long long int64; int64 v[210000],tx[210000],ty[210000],rk[210000],a,b,c,d,s,t,P,pp,x,y,sp,gg,tot,rs,ans,cnt; int64 gcd(int64 a,int64 b) { x=1,y=0; if(!b) return a; int64 d=gcd(b,a%b); int64 xx=x; x=y; y=xx-a/b*y; return d; } int64 qpow(int64 a,int64 b,int64 c) { int64 res=1; for(;b;b>>=1,a=a*a%c) if(b&1) res=res*a%c; return res; } int64 dt2(int64 a) { int l=1,r=cnt,mid; while(l+1>1; if(tx[rk[mid]]>1; if(v[mid]=a) return v[l]; return v[r]; } int64 solve(int64 a,int64 b,int64 P,int64 t) { a=(a+P)%P; b=(b+P)%P; tot=0; int64 sp=0,ans=b; for(;(sp+1)*(sp+1)<=t;++sp); for(int i=0;sp*sp+i<=t;++i) ans=min(ans,(a*(sp*sp+i)+b)%P); if(!sp) return ans; for(int i=1;i<=sp;++i) v[++tot]=(a*sp*i+b)%P; sort(v+1,v+tot+1); for(int i=0;i<=sp;++i) { int64 x=a*i%P; ans=min(ans,(v[1]+P-x)%P); ans=min(ans,(t2(x)-x+P)%P); } return ans; } bool cmp(int a,int b) { return tx[a]1) v[++tot]=pp; pp=P-1; gg=getpg(P); cnt=0; for(int i=1;i<=sp;++i) tx[++cnt]=qpow(gg,sp*i,P),ty[cnt]=sp*i,rk[cnt]=cnt; sort(rk+1,rk+cnt+1,cmp); a=dlog(a); b=dlog(b); c=dlog(c); a=(pp-a)%pp; c=(c-b*s%pp+pp)%pp; t-=s; d=gcd(pp,gcd(a,b)); if(c%d) { printf("impossible\n"); continue; } a/=d; b/=d; c/=d; pp/=d; d=gcd(pp,a); gcd(b,d); x=(x%d+d)%d; rs=c*x%d; if(t