/* -*- coding: utf-8 -*- * * 1073.cc: No.1073 無限すごろく - yukicoder */ #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; /* constant */ const int N = 6; const int MOD = 1000000007; /* typedef */ typedef long long ll; typedef int vec[N]; typedef vec mat[N]; /* global variables */ mat ma, mb; /* subroutines */ int powmod(int a, int n) { // a^n % MOD int pm = 1; while (n > 0) { if (n & 1) pm = (ll)pm * a % MOD; a = (ll)a * a % MOD; n >>= 1; } return pm; } inline void initvec(vec a) { memset(a, 0, sizeof(vec)); } inline void initmat(mat a) { memset(a, 0, sizeof(mat)); } inline void unitmat(mat a) { initmat(a); for (int i = 0; i < N; i++) a[i][i] = 1; } inline void copymat(const mat a, mat b) { memcpy(b, a, sizeof(mat)); } inline void mulmat(const mat a, const mat b, mat c) { for (int i = 0; i < N; i++) for (int j = 0; j < N; j++) { c[i][j] = 0; for (int k = 0; k < N; k++) c[i][j] = (c[i][j] + (ll)a[i][k] * b[k][j] % MOD) % MOD; } } inline void powmat(const mat a, ll b, mat c) { mat s, t; copymat(a, s); unitmat(c); while (b > 0) { if (b & 1) { mulmat(c, s, t); copymat(t, c); } mulmat(s, s, t); copymat(t, s); b >>= 1; } } /* main */ int main() { ll n; scanf("%lld", &n); int p = powmod(N, MOD - 2); ma[0][0] = ma[0][1] = ma[0][2] = ma[0][3] = ma[0][4] = ma[0][5] = p; ma[1][0] = ma[2][1] = ma[3][2] = ma[4][3] = ma[5][4] = 1; powmat(ma, n, mb); printf("%d\n", mb[0][0]); return 0; }