Codeforces 571E

首先将所有的$a_i$，$b_i$均分解质因数，由于$n$较小，采用根号试除法即可。

假设所有的等比数列共有$m$个不同的质因子，从小到大依次为$p_1$，$p_2$，$...$，$p_m$。于是我们可以知道第$i$个等比数列的第$k$项为${\prod }_{j=1}^m{p}_j^{c_{i,j}+d_{i,j}\cdot k}$，其中$c_{i,j}$是$a_i$中$p_j$的指数，$d_{i,j}$是$b_i$中$p_j$的幂次。

若$\exists j$，使得$\exists u$有$d_{u,j}=0$。我们发现当$\forall 1\le i\le n$都有$d_{i,j}=0$时，显然当$c_{i,j}$均相等时可以忽略掉这个质因子（最后答案乘上即可），否则直接无解。若$\exists v$有$d_{v,j}>0$，那么合法的$x$至多只有一个，直接尝试解出来判定是否合法即可。

经过上面的处理后，我们可以认为所有的$d_{i,j}$均不为$0$。此时如果仍有$m>1$，我们尝试消到只剩一个质因子，具体是对于某个$1<j\le m$，我们枚举$1<i\le n$，若有$\frac{d_{1,1}}{d_{i,1}}\ne \frac{d_{1,j}}{d_{i,j}}$，那么同样合法的$x$至多只有一个，尝试解出来判定是否合法即可，否则先看$c_{i,j}$是否成比例，如果不成比例也无解，成的话就可以忽略掉$p_j$了。

这样我们终于转化为了只有一个质因子$p_1$的情况，那么我们设最终$x$中$p_1$的幂次为$y$，显然有$\forall 1\le i\le n$，$y\equiv c_{i,1}(\mathrm{m}\mathrm{o}\mathrm{d}{d}_{i,1})$。这样问题转化为了给出$n$个线性同余方程，要求出最小的非负整数解。这个问题有经典的扩展CRT做法，即每次尝试合并两个同余方程（可能无解）。求出了最小的$y$后就容易求出答案了。

注意有一个坑点，刚刚的同余方程其实忽略了一些信息，有可能这样求出来的最小的$y$会比某个$c_{i,1}$小，因此需要加上$lc{m}_{i=1}^n{d}_{i,1}$的倍数。

时间复杂度是$\mathcal{O}(n(\sqrt{V}+m))$。

#include <bits/stdc++.h>
#define MOD 1000000007

using namespace std;

typedef long long ll;

inline ll lcm(ll x,ll y) {
  return x/__gcd(x,y)*y;
}

ll pow_mod(ll x,ll k) {
  ll ans=1;
  while (k) {
  	if (k&1) ans=ans*x%MOD;
  	x=x*x%MOD;
  	k>>=1;
  }
  return ans;
}

void exgcd(ll a,ll b,ll &x,ll &y) {
  if (!b) {
  	x=1;y=0;
  }
  else {
    exgcd(b,a%b,y,x);
    y-=(a/b)*x;
  }
}

ll getinv(ll x,ll mod) {
  ll t1,t2;
  exgcd(x,mod,t1,t2);
  return (t1%mod+mod)%mod;
}

map <int,int> mp;
int prime[2005],cnt;
bool vis[2005];

int a[105][2005],b[105][2005];

void check_sol(int n,int x,ll d) {
  int ans=1;
  for(int i=1;i<=cnt;i++) {
  	ll v=a[x][i]+b[x][i]*d;
    for(int j=1;j<=n;j++)
      if (b[j][i]) {
      	if (v<a[j][i]||(v-a[j][i])%b[j][i]!=0) return;
	  }
	  else if (a[j][i]!=v) return;
	ans=ans*pow_mod(prime[i],v)%MOD;
  }
  printf("%d\n",ans);
  exit(0);
}

bool check(int n) {
  for(int i=1;i<=cnt;i++) {
  	int id=0;
  	for(int j=1;j<=n;j++)
  	  if (!b[j][i]) id=j;
  	if (!id) continue;
	vis[i]=1;
	for(int j=1;j<=n;j++)
	  if (b[j][i]) {
	  	int v=a[id][i]-a[j][i];
	  	if (v%b[j][i]!=0) return 0;
	  	v/=b[j][i];
	  	if (v<0) return 0;
		check_sol(n,j,v);
		return 0; 
	  }
	  else if (a[j][i]!=a[id][i]) return 0;
  }
  int id=1;
  while (id<=cnt&&vis[id]) id++;
  if (id>cnt) check_sol(n,1,0);
  for(int i=id+1;i<=cnt;i++)
    if (!vis[i]) {
    	for(int j=2;j<=n;j++) {
    		int u=a[j][id]-a[1][id],v=a[j][i]-a[1][i];
    		if (b[1][id]*b[j][i]!=b[1][i]*b[j][id]) {
    			int t1=b[1][id]*b[j][i]-b[1][i]*b[j][id],t2=u*b[j][i]-v*b[j][id];
    			if (t2%t1!=0) return 0;
    			t2/=t1;
    			if (t2<0) return 0;
    			check_sol(n,1,t2);
			}
			else if (b[1][id]*v!=b[1][i]*u) return 0;
		}
	}
  return 1;
}

bool solve(int n) {
  int id=1;
  while (id<=cnt&&vis[id]) id++;
  ll s=a[1][id]%b[1][id],w=b[1][id];
  for(int i=2;i<=n;i++) {
  	ll u=a[i][id]%b[i][id],v=b[i][id];
  	ll d=__gcd(w,v);
  	if ((u-s)%d!=0) return 0;
  	ll ww=w/d,vv=v/d,cc=(u-s)/d;
  	ll x=(getinv(ww,vv)*cc%vv+vv)%vv;
  	s=w*x+s;w=lcm(w,v);
  }
  for(int i=1;i<=n;i++) {
  	while (s<a[i][id]) s+=w;
  }
  check_sol(n,1,(s-a[1][id])/b[1][id]);
}

void getfact(int x,int *p) {
  int t=x;
  for(int i=2;i*i<=t;i++)
    if (x%i==0) {
    	if (!mp.count(i)) {
		  mp[i]=++cnt;
		  prime[cnt]=i;
	    }
    	int u=mp[i];
    	while (x%i==0) {
    		p[u]++;
    		x/=i;
		} 
	}
  if (x>1) {
  	if (!mp.count(x)) {
	  mp[x]=++cnt;
	  prime[cnt]=x;
    }
  	int u=mp[x];
  	p[u]=1;
  }
}

int main() {
  int n;
  scanf("%d",&n);
  for(int i=1;i<=n;i++) {
  	int x,y;
  	scanf("%d%d",&x,&y);
  	getfact(x,a[i]);
  	getfact(y,b[i]);
  }
  if (!check(n)) {
  	puts("-1");
  	return 0;
  }
  if (!solve(n)) {
  	puts("-1");
  	return 0;
  }
  return 0;
}