#include #include #include #include using namespace std; long P,R; long power(long a,long b){return b?power(a*a%P,b/2)*(b%2?a:1)%P:1;} main() { int q;cin>>P>>R>>q; int m=sqrt(q*P); vector >table(m); long e=1; for(int i=0;i>a>>b>>c; long d=(b*b%P-4*a*c%P+P)%P*power(4*a*a%P,P-2)%P; long t=(P-b*power(2*a%P,P-2)%P)%P; if(d==0) { cout<p=*lower_bound(table.begin(),table.end(),make_pair((int)now,-1)); if(p.first==now) { k=i*m+p.second; break; } now=now*inv%P; } if(k%2==0) { long x=(power(R,k/2)+t)%P,y=(power(R,(k/2+P/2)%(P-1))+t)%P; if(x