#include #define rep(i,n) for(int i=0;i<(int)(n);i++) #define rep1(i,n) for(int i=1;i<=(int)(n);i++) #define all(c) c.begin(),c.end() #define pb push_back #define fs first #define sc second #define show(x) cout << #x << " = " << x << endl #define chmin(x,y) x=min(x,y) #define chmax(x,y) x=max(x,y) using namespace std; typedef long long ll; const ll c=32000; bool isp[c+1]; vector ps; void makeprime(){ for(int i=2;i<=c;i++) isp[i]=1; for(int i=2;i<=c;i++) if(isp[i]){ for(int j=i*2;j<=c;j+=i) isp[j]=0; ps.pb(i); } } ll inv(ll a,ll m){ ll b=m,u=1,v=0; while(b){ ll t=a/b; swap(a-=t*b,b); swap(u-=t*v,v); } u%=m; if(u<0) u+=m; return u; } typedef pair P; //val,mod; map > mp; vector

vp; bool all0=1; ll garner(vector

rm,ll mod){ //remainder,mod rm.pb(P(0,mod)); int N=rm.size(); vector as(N,1); //coefficients vector bs(N,0); //constants rep(i,N-1){ ll v=(rm[i].fs-bs[i])*inv(as[i],rm[i].sc)%rm[i].sc; if(v<0) v+=rm[i].sc; for(int j=i+1;j>N; rep(i,N){ ll x,m; cin>>x>>m; if(x) all0=0; for(int p:ps){ if(m%p==0){ int k=0; int a=1; while(m%p==0) m/=p,a*=p,k++; mp[p].pb(P(x%a,a)); } } if(m>1){ mp[m].pb(P(x%m,m)); } } for(auto it:mp){ vector

& vc=it.sc; ll q=0,x=0; for(P p:vc){ if(q