#include int64_t calcGCD(int64_t a, int64_t b) { while (a) { b %= a; std::swap(a, b); } return b; } //////////////////////////// // 拡張ユークリッドの互除法 // //////////////////////////// // x, yは非負整数 // ret[0] * x + ret[1] * y == ret[2] == gcd(x, y) となるようなretを返す std::array extendedEuclidean(const int64_t x, const int64_t y) { if (y == 0) return {1, 0, x}; auto next{extendedEuclidean(y, x % y)}; return {next[1], next[0] - (x / y) * next[1], next[2]}; } ///////////////////////// // Garnerのアルゴリズム // ///////////////////////// int64_t garner(const std::vector& rests, const std::vector& mods, const int64_t ans_mod) { std::vector coefficient(rests); for (int i{}; i < (int)rests.size(); i++) { int64_t mod_multi{1ll}; for (int j{}; j < i; j++) { coefficient[i] = (coefficient[i] + mods[i] - mod_multi * coefficient[j] % mods[i]) % mods[i]; mod_multi = mod_multi * mods[j] % mods[i]; } for (int j{}; j < i; j++) coefficient[i] = coefficient[i] * ((extendedEuclidean(mods[j], mods[i])[0] + mods[i]) % mods[i]) % mods[i]; } int64_t ret{}, mod_multi{1ll}; for (int i{}; i < (int)rests.size(); i++) { ret = (ret + mod_multi * coefficient[i] % ans_mod) % ans_mod; mod_multi = mod_multi * mods[i] % ans_mod; } return ret; } int main() { using namespace std; using ll = int64_t; int n; cin >> n; vector vs(n), ms(n); for (int i = 0; i < n; i++) { cin >> vs[i] >> ms[i]; } vector primes; for (int i = 2; (ll)i * i <= 1e9; i++) { int flag = 1; for (int x : primes) if (i % x == 0) { flag = 0; break; } if (flag) primes.emplace_back(i); } for (int m : ms) { for (int p : primes) while (m % p == 0) m /= p; if (m != 1) primes.emplace_back(m); } for (int p : primes) { vector> ps; for (int i = 0; i < n; i++) { if (ms[i] % p == 0) { int lg = 0; while (ms[i] % p == 0) ms[i] /= p, lg++; ps.emplace_back(lg, i); } } if (ps.size() < 1) continue; sort(begin(ps), end(ps)); int base = vs[ps.back().second]; for (auto d : ps) { if (vs[d.second] % p != base % p) { cout << -1 << endl; return 0; } if (d.second == ps.back().second) { for (int i = 0; i < d.first; i++) ms[d.second] *= p; } } } int allzero = 1; ll ans = 1; for (int i = 0; i < n; i++) { vs[i] = vs[i] % ms[i]; if (vs[i] != 0) allzero = 0; ans = ans * ms[i] % int(1e9 + 7); } if (allzero) { cout << ans << endl; return 0; } cout << garner(vs, ms, 1e9 + 7) << endl; }