#include using namespace std; #define repl(i,a,b) for(int i=(int)(a);i<(int)(b);i++) #define rep(i,n) repl(i,0,n) #define mp(a,b) make_pair((a),(b)) #define pb(a) push_back((a)) #define all(x) (x).begin(),(x).end() #define uniq(x) sort(all(x)),(x).erase(unique(all(x)),end(x)) #define fi first #define se second #define dbg(...) _dbg(#__VA_ARGS__, __VA_ARGS__) void _dbg(string){cout< void _dbg(string s,H h,T... t){int l=s.find(',');cout< ostream& operator<<(ostream& o, const pair &p){o<<"("< ostream& operator<<(ostream& o, const vector &v){o<<"[";for(T t:v){o<> mat; // return A*B mat mat_mul(const mat &A, const mat &B){ int n=A.size(), m=B[0].size(), l=B.size(); mat ret(n, vector(m, 0)); rep(i,n) rep(k,l) if(A[i][k]!=0) rep(j,m){ (ret[i][j] += A[i][k] * B[k][j]) %= MOD; } return ret; } // A^p mat mat_pow(const mat &A, long p){ int n = A.size(); mat tmp(A), ret(n, vector(n,0)); rep(i,n) ret[i][i] = 1; while(p>0){ if(p&1) ret = mat_mul(tmp, ret); tmp = mat_mul(tmp, tmp); p /= 2; } return ret; } int main(){ long n; cin>>n; const mat m = { {1,0,1,0,0,0,0,1,0,1}, {0,1,0,0,1,0,0,0,0,0}, {0,1,1,1,1,0,0,0,0,1}, {0,0,1,1,1,1,1,0,0,1}, {1,0,1,1,1,1,1,1,0,1}, {0,0,0,1,1,1,1,0,0,1}, {0,0,0,1,1,1,1,1,0,1}, {0,1,0,0,1,0,1,1,0,1}, {0,1,1,1,1,1,1,1,1,0}, {0,0,0,0,0,0,0,0,0,1} }; mat p = mat_pow(m, n+1); cout << p[8][9] % MOD << endl; return 0; }