#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 > V, ll mo=1000000007) { const int prime_max = 40000; static int NP,prime[10000],divp[prime_max+3]; int i,j,k,l,r,x,y; string s; int N=V.size(); set P; vector num; if(NP==0) { for(int i=2;ii&&j1) P.insert(v); } num=vector(N,0); FORR(r,P) { y=0; FOR(x,N) { num[x]=1; j=V[x].second; while(j%r==0) j/=r, num[x]*=r; if(num[y] > V2; FOR(x,N) if(V[x].second>1) V2.emplace_back(V[x].second,V[x].first); sort(V2.begin(),V2.end()); V=V2; N=V.size(); // garner FOR(y,N) FOR(x,y) { ll k=V[y].second-V[x].second; if(k<0) k+=V[y].first; V[y].second = k*inv(V[x].first,V[y].first)%V[y].first; } ll ret=0; for(x=N-1;x>=0;x--) ret = (ret*V[x].first+V[x].second)%mo; return ret; } void solve() { int N,i,x,y,nz=0; vector > V; cin>>N; FOR(i,N) { cin>>x>>y; V.emplace_back(x,y); nz+=x>0; } if(nz==0) { FOR(y,N) FOR(x,y) V[y].second/=__gcd(V[x].second,V[y].second); ll ret=1; FOR(x,N) ret=ret*V[x].second%1000000007; cout<