#include #include #include #include #include #include #include #include #include #include static const int MOD = 1000000007; using ll = long long; using u32 = uint32_t; using namespace std; template constexpr T INF = ::numeric_limits::max()/32*15+208; template struct modint{ ll val; modint(): val(0){} template explicit modint(T t){val = t%M; if(val < 0) val += M;} modint pow(ll k){ modint res(1), x(val); if(!val) return modint(0); while(k){ if(k&1) res *= x; x *= x; k >>= 1; } return res; } template modint& operator=(T a){ val = a%M; if(val < 0) val += M; return *this; } modint inv() {return pow(M-2);} modint& operator+=(modint a){ val += a.val; if(val >= M) val -= M; return *this;} modint& operator-=(modint a){ val += M-a.val; if(val >= M) val -= M; return *this;} modint& operator*=(modint a){ val = 1LL*val*a.val%M; return *this;} modint& operator/=(modint a){ return (*this) *= a.inv();} modint operator+(modint a) const {return modint(val) +=a;} modint operator-(modint a) const {return modint(val) -=a;} modint operator*(modint a) const {return modint(val) *=a;} modint operator/(modint a) const {return modint(val) /=a;} modint operator-(){ return modint(-val);} bool operator==(const modint a) const {return val == a.val;} bool operator!=(const modint a) const {return val != a.val;} bool operator<(const modint a) const {return val < a.val;} }; using mint = modint; template struct SquareMatrix { using ar = array; using mat = array; mat A; SquareMatrix() = default; static SquareMatrix I(T e){ SquareMatrix X; for (int i = 0; i < SIZE; ++i) { X[i][i] = e; } return X; } friend ar operator*=(ar &x, const SquareMatrix &Y) { ar ans; for (int i = 0; i < SIZE; ++i) { for (int j = 0; j < SIZE; ++j) { ans[j] += x[i]*Y[i][j]; } } x.swap(ans); return x; } friend ar operator*(ar x, const SquareMatrix &Y) { return x *= Y; } inline const ar &operator[](int k) const{ return (A.at(k)); } inline ar &operator[](int k) { return (A.at(k)); } SquareMatrix &operator+= (const SquareMatrix &B){ for (int i = 0; i < SIZE; ++i) { for (int j = 0; j < SIZE; ++j) { (*this)[i][j] += B[i][j]; } } return (*this); } SquareMatrix &operator-= (const SquareMatrix &B){ for (int i = 0; i < SIZE; ++i) { for (int j = 0; j < SIZE; ++j) { (*this)[i][j] -= B[i][j]; } } return (*this); } SquareMatrix &operator*=(const SquareMatrix &B) { SquareMatrix C; for (int i = 0; i < SIZE; ++i) { for (int j = 0; j < SIZE; ++j) { for (int k = 0; k < SIZE; ++k) { C[i][j] += ((*this)[i][k] * B[k][j]); } } } A.swap(C.A); return (*this); } SquareMatrix pow(ll n) const { SquareMatrix a = (*this), res = I(T(1)); while(n > 0){ if(n & 1) res *= a; a *= a; n >>= 1; } return res; } SquareMatrix operator+(const SquareMatrix &B) const {return SquareMatrix(*this) += B;} SquareMatrix operator-(const SquareMatrix &B) const {return SquareMatrix(*this) -= B;} SquareMatrix operator*(const SquareMatrix &B) const {return SquareMatrix(*this) *= B;} }; using ar = array; using mat = SquareMatrix; int main() { int n; cin >> n; vector> v; mat A; A[0][1] = 1; A[1][0] = 1; A[1][1] = 1; mint ans(1); for (int i = 0; i < n; ++i) { ll x; string y; cin >> x >> y; ll d = 0; for (auto &&j : y) { d = (d*10+j-'0')%(MOD-1); } mat X = A.pow(x+1); ans *= X[1][1].pow(d); } cout << ans.val << "\n"; return 0; }