#include #define rep(i, n) for (int (i) = 0; (i) < (int)(n); (i)++) const int dx[] = {1, 0, -1, 0}; const int dy[] = {0, 1, 0, -1}; using namespace std; typedef long long ll; const int MOD = 1e9+7; const int MM = MOD-1; //typedef double number; //const number eps = 1e-8; typedef long long number; typedef vector vec; typedef vector matrix; // O( n ) matrix identity(int n) { matrix A(n, vec(n)); for (int i = 0; i < n; ++i) A[i][i] = 1; return A; } // O( n ) number inner_product(const vec &a, const vec &b) { number ans = 0; for (int i = 0; i < a.size(); ++i) ans += a[i] * b[i]; return ans; } // O( n^2 ) vec mul(const matrix &A, const vec &x) { vec y(A.size()); for (int i = 0; i < A.size(); ++i) for (int j = 0; j < A[0].size(); ++j) y[i] = A[i][j] * x[j]; return y; } // O( n^3 ) matrix mul(const matrix &A, const matrix &B) { matrix C(A.size(), vec(B[0].size())); for (int i = 0; i < C.size(); ++i) for (int j = 0; j < C[i].size(); ++j) for (int k = 0; k < A[i].size(); ++k) { C[i][j] += A[i][k] * B[k][j]; C[i][j] %= MOD; } return C; } // O( n^3 log e ) matrix pow(const matrix &A, ll e) { return e == 0 ? identity(A.size()) : e % 2 == 0 ? pow(mul(A, A), e/2) : mul(A, pow(A, e-1)); } // O( n^3 ) //number det(matrix A) { // const int n = A.size(); // number D = 1; // for (int i = 0; i < n; ++i) { // int pivot = i; // for (int j = i+1; j < n; ++j) // if (abs(A[j][i]) > abs(A[pivot][i])) pivot = j; // swap(A[pivot], A[i]); // D *= A[i][i] * (i != pivot ? -1 : 1); // if (abs(A[i][i]) < eps) break; // for(int j = i+1; j < n; ++j) // for(int k = n-1; k >= i; --k) // A[j][k] -= A[i][k] * A[j][i] / A[i][i]; // } // return D; //} // O(n) number tr(const matrix &A) { number ans = 0; for (int i = 0; i < A.size(); ++i) ans += A[i][i]; return ans; } // O( n^3 ). //int rank(matrix A) { // const int n = A.size(), m = A[0].size(); // int r = 0; // for (int i = 0; r < n && i < m; ++i) { // int pivot = r; // for (int j = r+1; j < n; ++j) // if (abs(A[j][i]) > abs(A[pivot][i])) pivot = j; // swap(A[pivot], A[r]); // if (abs(A[r][i]) < eps) continue; // for (int k = m-1; k >= i; --k) // A[r][k] /= A[r][i]; // for(int j = r+1; j < n; ++j) // for(int k = i; k < m; ++k) // A[j][k] -= A[r][k] * A[j][i]; // ++r; // } // return r; //} ll getMod(const string& s) { int n = s.size(); ll ans = 0; for (int i = 0; i < n; i++) { ans *= 10; ans += (s[i]-'0'); ans %= MM; } return ans; } // x^p ll powmod(ll x, ll p, ll m = MOD) { if (x == 0) return 0; if (p == 0) return 1; if (p == 1) return x; if (p % 2 == 0) { ll tmp = powmod(x, p/2, m); return (tmp*tmp)%m; } else { ll tmp = powmod(x, p-1, m); return (tmp*x)%m; } } int main() { int N; cin >> N; ll ans = 1; matrix A(2, vector(2)); A[0][1] = A[1][0] = A[1][1] = 1; while (N--) { ll C; string D; cin >> C >> D; ll d = getMod(D); matrix B = pow(A, C-1); ll tmp = B[1][0] + B[1][1]*2; ans *= powmod(tmp, d); ans %= MOD; } cout << ans << endl; return 0; }