#include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include //#include #pragma warning(disable:4996) typedef long long ll; typedef unsigned long long ull; #define MIN(a, b) ((a)>(b)? (b): (a)) #define MAX(a, b) ((a)<(b)? (b): (a)) #define LINF 9223300000000000000 #define LINF2 1223300000000000000 #define LINF3 1000000000000 #define INF 2140000000 const long long MOD = 1000000007; //const long long MOD = 998244353; using namespace std; ll mpow(ll x, ll n){ //x^n(mod M) ll ans = 1; while(n != 0){ if(n&1) ans = ans*x % MOD; x = x*x % MOD; n = n >> 1; } return ans; } ll minv(ll x){ return mpow( x, MOD-2 ); } void mtxprd( int p, int q, int r, const ll* A, const ll* B, ll* C ) { int i, j, k; for(i=0; i mtx0(dim*dim); vector mtxtmp(dim*dim); mtxcpy( dim, mtx, &mtx0[0] ); mtxuni( dim, mtx2 ); while(n != 0){ if(n&1) { mtxprd( dim, dim, dim, mtx2, &mtx0[0], &mtxtmp[0] ); mtxcpy( dim, &mtxtmp[0], mtx2 ); } mtxprd( dim, dim, dim, &mtx0[0], &mtx0[0], &mtxtmp[0] ); mtxcpy( dim, &mtxtmp[0], &mtx0[0] ); n = n >> 1; } return; } void solve() { ll n; scanf("%lld", &n); vector a(6*6), b(6*6); #define P(i,j) ((i)*6+(j)) int i; for (i = 0; i < 5; i++) { a[P(i + 1, i)] = 1; } for (i = 0; i < 6; i++) { a[P(0, i)] = minv(6); } mtxpow( 6, &a[0], n, &b[0] ); ll ans = 0; printf("%lld\n", b[P(0,0)]); return; } int main(int argc, char* argv[]) { #if 1 solve(); #else int T; scanf("%d", &T); int t; for(t=0; t