#include #define all(x) (x).begin(), (x).end() typedef long long ll; #define MOD 1000000007 using namespace std; ll k; long long extGCD(long long a, long long b, long long &x, long long &y) { if(b == 0) { x = 1; y = 0; return a; } long long d = extGCD(b, a % b, y, x); // 再帰的に解く y -= a / b * x; return d; } // 負の数にも対応した mod (a = -11 とかでも OK) inline long long mod(long long a, long long m) { return (a % m + m) % m; } // 逆元計算 (ここでは a と m が互いに素であることが必要) long long modinv(long long a, long long m) { long long x, y; extGCD(a, m, x, y); return mod(x, m); // 気持ち的には x % m だが、x が負かもしれないので } vector> ope(vector> &mat1, vector> &mat2) { vector> res(k, vector(k, 0)); for(int i = 0; i < k; i++) { for(int j = 0; j < k; j++) { for(int l = 0; l < k; l++) { res[i][j] += mat1[i][l] * mat2[l][j] % MOD; res[i][j] %= MOD; } } } return res; } vector> modpow(vector> a, ll n) { vector> res = a; n--; while(n > 0) { if(n & 1) res = ope(res, a); a = ope(a, a); n >>= 1; } return res; } int main() { ll n; cin >> n; n += 5; k = 6; vector init = {1, 0, 0, 0, 0, 0}; vector> mat(k, vector(k, 0)); for(int i = 0; i < k; i++) { mat[0][i] = modinv(6, MOD); } for(int i = 1; i < k; i++) { for(int j = 0; j < k; j++) { mat[i][i - 1] = 1; } } auto ret = modpow(mat, n - 5); ll ans = 0; for(int i = 0; i < k; i++) { ans += ret[0][i] * init[i] % MOD; ans %= MOD; } cout << ans << endl; }