#include using namespace std; using LL = long long; using ULL = unsigned long long; #define rep(i,n) for(int i=0; i<(n); i++) // |--- // -|-- // --|- // ---| // // |LJ- // |L_J // |-LJ // -|LJ // // LJ|- // LJ-| // L_J| // -LJ| // const ULL M = 1000000007; const int matrix_sz = 12; struct Matrix{ ULL X[matrix_sz][matrix_sz] = {}; static Matrix id(){ Matrix res; rep(i, matrix_sz) res[i][i] = 1; return res; } ULL* operator[](int x){ return X[x]; } const ULL* operator[](int x) const { return X[x]; } Matrix pow(ULL N) const; }; Matrix operator*(Matrix l, Matrix r){ Matrix res; rep(i, matrix_sz) rep(j, matrix_sz) rep(k, matrix_sz) res[i][j] = (res[i][j] + l[i][k] * r[k][j]) % M; return res; } Matrix Matrix::pow(ULL N) const{ if (N == 0) return id(); Matrix res = pow(N / 2); res = res * res; if (N % 2 == 1) res = res * *this; return res; } const ULL Gconst[matrix_sz][matrix_sz] = { { 1, 1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0 }, { 1, 1, 1, 1, 0, 0, 1, 1, 0, 0, 0, 0 }, { 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 0, 0 }, { 1, 1, 1, 1, 0, 0, 0, 0, 1, 1, 1, 1 }, { 0, 0, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0 }, { 0, 0, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0 }, { 0, 0, 0, 1, 0, 1, 1, 1, 0, 0, 0, 0 }, { 0, 0, 0, 1, 0, 0, 1, 1, 0, 0, 0, 0 }, { 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0 }, { 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0 }, { 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1 }, { 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1 } }; int main(){ Matrix G; rep(i, matrix_sz) rep(j, matrix_sz) G[i][j] = Gconst[i][j]; ULL N; cin >> N; G = G.pow(N + 1); cout << G[0][3] << endl; return 0; }