#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #define X first #define Y second #define pb push_back #define rep(X,Y) for (int (X) = 0;(X) < (Y);++(X)) #define rrep(X,Y) for (int (X) = (Y-1);(X) >=0;--(X)) #define repe(X,Y) for ((X) = 0;(X) < (Y);++(X)) #define peat(X,Y) for (;(X) < (Y);++(X)) #define all(X) (X).begin(),(X).end() #define rall(X) (X).rbegin(),(X).rend() using namespace std; typedef long long ll; typedef pair pii; typedef pair pll; template using vv=vector>; template ostream& operator<<(ostream &os, const vector &t) { os<<"{"; rep(i,t.size()) {os< ostream& operator<<(ostream &os, const pair &t) { return os<<"("< primes; int prime[PRIME_MAX]; void findPrime(){ for(int i=2;i a,vector m,vector p,int f,ll mod=MOD){ int n=a.size(); m.pb(mod); p.pb(mod); vector d(n+1,1),sum(n+1); rep(i,n){ // sum[i]+x*d[i]==a[i] mod m[i]; ll x=(a[i]-sum[i])%m[i]*modpow(d[i],m[i]/p[i]*(p[i]-1)-1,m[i])%m[i]; //cout< &a,const vector &m){ ll p=1,re=0,pp,x,y; rep(i,m.size()) p*=m[i]; //cout< &mp,ll a,ll m,ll p){ //cout< &a,const vector &m,int f=0){ map mp; //PRIME -> REM * MOD pll tmp; ll num,p=1; rep(i,a.size()){ num=m[i]; rep(j,primes.size())if(m[i]%primes[j]==0){ p=1; while(num%primes[j]==0){ p*=primes[j]; num/=primes[j]; } if(updChina(mp,a[i],p,primes[j])) return -1; } if(num>1) if(updChina(mp,a[i],num,num)) return -1; } //for(auto t:mp)cout< aa,mm,pp; for(auto temp:mp){ num=temp.Y.X; p=temp.Y.Y; aa.pb((num%p+p)%p); mm.pb(p); pp.pb(temp.X); } //cout<>n; vector a(n),m(n); long long te,mp; rep(i,n){ cin>>te>>mp; a[i]=te; m[i]=mp; } ll re=china(a,m),lcm=1; //cout<