#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)) //------------------------------------------------------- ll mo,P,R; int Q; ll A,B,C; ll modpow(ll a, ll n = mo-2) { ll r=1; while(n) r=r*((n%2)?a:1)%mo,a=a*a%mo,n>>=1; return r; } inline int mulmod(int a,int b,int mo) { int d,r; if(a==0 || b==0) return 0; if(a==1 || b==1) return max(a,b); __asm__("mull %4;" "divl %2" : "=d" (r), "=a" (d) : "r" (mo), "a" (a), "d" (b)); return r; } ll mod_log(ll a,ll b,ll mo) { // a^x=b static ll pre_a=-1; static unordered_map M; static ll a_sq=1; static int sq=0; int i,j; if(pre_a!=a) M.clear(); if(M.empty()) { pre_a=a; ll ap=1,ar=modpow(a),as=1; for(i=0;i*i>P>>R; mo=P; cin>>Q; while(Q--) { cin>>A>>B>>C; ll r=modpow(A); A = A*r%mo; B = B*r%mo; C = C*r%mo; if(B%2==1) B+=P; B/=2; C=(P-(C-B*B%P+P)%P)%P; if(C==0) { cout<<(P-B)%mo<Y) swap(X,Y); cout<