// // S i r o t a n w a t c h e s o v e r y o u . // // ...Jggg+J+JJJg@NQmgJa....., // ....gH@@@@HHB""""7"YYYYHMMMMMMM@@@@@@@Hma,. // ...JH@@MMHY"! ? __""YH@@@@N&... // ..JH@@HY"~ _TW@@@Mme. // .Jg@@#"= _TTM@@N.. // .Jg@@MY! ?T@@@h,. // .-@@@B! (TM@@@L // .(H@MY ?W@@@+ // .W@@@ .T@@@[ // .d@@M! .T@@N, // .d@M@' -W@@o. // (@@M\ ?M@N, // -d@M% .., .., j@@b // d@@M= TM9 ?MD W@@[ // .H@M: .W@H, // H@Ht ,@@# // (@@M~ .@@h. // .@@M% ..gmHHJ. .JdHga. .H@@e // j@@#_ .@@@@@@@b J@@@@@@@h. .H@@\ // d@@@ .4@@@@@@MF (@@@@@@@@ H@@b // d@@[ ?"BMY"= .d@@@@H, ?TWHHY"! d@@e // J@@b .JJJ-.., ,@@@@@@M% ......... -@@@M. // ?@@M\ ?YHHM@@@@b .. .W@@HP X@@HMWY"= d@@@# // ,@@@L. ?H@Ng&+gd@@#H@@NHaJ+gH@[ J@@@] // X@@@[ ...... ?"YYYYY"" ?"YHHHB"^ ..... (@@@# // WH@N+. .W@@@@MHB= .TWH@M@Hmc .H@@M~ // .H@@@@N, _!~ .d@@@N, // .J@@#T@@@N, .d@@@@@@@b. // (@@@@! .T@@@n, .(H@@@H>.W@@@x // (@@@F 4@@@@MaJ. .d@@@@Y77 4@@@r //.H@@P ?TM@@@@N... .-JH@HMY= d@@N, //(@@@F ?"WM@@@MQa-,. .(J(JN@@M#" Z@@@L // d@@H, (M@@@@@@@Ng&maJ.... .. ...J.J+W@@@@@@HY! .dH@b // ?M@@@N&. ..(JW@@@MM"?7""TYHMH@@HH@@@@@HHHgkHagHa(mggdmmagH@H@@Q@@HMMMHY"7!TMM@@@HaJ,. ..H@@@M^ // ?"W@@@@MN@@@@@H@@MMY" _???!"= ?WMMBYYTMH=7""Y@""?"~^ _"YM@@@@@@@@MH@@@@@#"^ // ?77WMMMYB""! _7"WWMMMY"7= // need #include #include // data structure #include #include #include #include #include #include #include #include #include //#include #include #include #include #include // etc #include #include #include #include #include #include #include #include // input #define INIT std::ios::sync_with_stdio(false);std::cin.tie(0); #define VAR(type, ...)type __VA_ARGS__;MACRO_VAR_Scan(__VA_ARGS__); template void MACRO_VAR_Scan(T& t) { std::cin >> t; } templatevoid MACRO_VAR_Scan(First& first, Rest& ...rest) { std::cin >> first; MACRO_VAR_Scan(rest...); } #define VEC_ROW(type, n, ...)std::vector __VA_ARGS__;MACRO_VEC_ROW_Init(n, __VA_ARGS__); for(int w_=0; w_ void MACRO_VEC_ROW_Init(int n, T& t) { t.resize(n); } templatevoid MACRO_VEC_ROW_Init(int n, First& first, Rest& ...rest) { first.resize(n); MACRO_VEC_ROW_Init(n, rest...); } template void MACRO_VEC_ROW_Scan(int p, T& t) { std::cin >> t[p]; } templatevoid MACRO_VEC_ROW_Scan(int p, First& first, Rest& ...rest) { std::cin >> first[p]; MACRO_VEC_ROW_Scan(p, rest...); } #define VEC(type, c, n) std::vector c(n);for(auto& i:c)std::cin>>i; #define MAT(type, c, m, n) std::vector> c(m, std::vector(n));for(auto& R:c)for(auto& w:R)std::cin>>w; // output templatevoid MACRO_OUT(const T t) { std::cout << t; } templatevoid MACRO_OUT(const First first, const Rest...rest) { std::cout << first << " "; MACRO_OUT(rest...); } #define OUT(...) MACRO_OUT(__VA_ARGS__); #define FOUT(n, dist) std::cout<=(a);--w) #define RREP(w, n) RFOR(w, 0, n) template bool IN(S a, T x, U b) { return a <= x && x < b; } template inline T CHMAX(T& a, const T b) { return a = (a < b) ? b : a; } template inline T CHMIN(T& a, const T b) { return a = (a > b) ? b : a; } // test template using V = std::vector; template using VV = V>; template std::ostream& operator<<(std::ostream& os, std::pair p) { os << "(" << p.first << ", " << p.second << ")"; return os; } // type/const using i64 = std::int_fast64_t; using u64 = std::uint_fast64_t; using ll = long long; using ull = unsigned long long; using ld = long double; using PAIR = std::pair; using PAIRLL = std::pair; constexpr int INFINT = (1 << 30) - 1; // 1.07x10^ 9 constexpr int INFINT_LIM = (1LL << 31) - 1; // 2.15x10^ 9 constexpr ll INFLL = 1LL << 60; // 1.15x10^18 constexpr ll INFLL_LIM = (1LL << 62) - 1 + (1LL << 62); // 9.22x10^18 constexpr double EPS = 1e-6; constexpr i64 MOD = 1000000007; constexpr double PI = 3.141592653589793238462643383279; template void FILL(T(&a)[N], const T & val) { for (auto& x : a) x = val; } template void FILL(ARY(&a)[N][M], const T & val) { for (auto& b : a) FILL(b, val); } template void FILL(std::vector & a, const T & val) { for (auto& x : a) x = val; } template void FILL(std::vector> & a, const T & val) { for (auto& b : a) FILL(b, val); } // ------------>8--------- template class MatrixMOD { private: std::valarray> mat; public: MatrixMOD(size_t m = 0, size_t n = 0, T init = 0) { if (n == 0) n = m; mat.resize(m); for (size_t i = 0; i < m; ++i) mat[i].resize(n, init); } MatrixMOD(std::valarray> a) { mat = a; } MatrixMOD init(size_t m = 0, size_t n = 0, T init = 0) { if (n == 0) n = m; mat.resize(m); for (size_t i = 0; i < m; ++i) mat[i].resize(n, init); return *this; } std::valarray & operator[](size_t i) { return mat[i]; } const std::valarray& operator[](size_t i) const { return mat[i]; } MatrixMOD& operator=(const MatrixMOD & r) { for (size_t i = 0; i < mat.size(); ++i) mat[i] = r[i]; return *this; } MatrixMOD operator+() const { return mat; } MatrixMOD operator-() const { MatrixMOD res(mat.size()); for (size_t i = 0; i < mat.size(); ++i) res[i] = (MOD - mat[i]) %= MOD; return res; } MatrixMOD & operator+=(const MatrixMOD & r) { for (size_t i = 0; i < mat.size(); ++i) (mat[i] += r[i]) %= MOD; return *this; } MatrixMOD & operator+=(const T & x) { for (size_t i = 0; i < mat.size(); ++i) (mat[i] += x) %= MOD; return *this; } MatrixMOD & operator-=(const MatrixMOD & r) { return *this += MOD - r; } MatrixMOD& operator-=(const T & x) { return *this += MOD - x; } MatrixMOD& operator*=(const MatrixMOD & r) { // O(N^3) MatrixMOD res(*this); for (size_t i = 0; i < mat.size(); ++i) { for (size_t j = 0; j < r[0].size(); ++j) { res[i][j] = 0; for (size_t k = 0; k < mat[0].size(); ++k) { (res[i][j] += mat[i][k] * r[k][j] % MOD) %= MOD; } } } return *this = res; } MatrixMOD& operator*=(const T & x) { for (size_t i = 0; i < mat.size(); ++i) (mat[i] *= x) %= MOD; return *this; } MatrixMOD & operator%=(const T & mod) { for (size_t i = 0; i < mat.size(); ++i) mat[i] %= MOD; return *this; } MatrixMOD& operator^=(ll p) { // O(N^3 logP) MatrixMOD res(*this); for (size_t i = 0; i < mat.size(); ++i) { for (size_t j = 0; j < mat[0].size(); ++j) { res[i][j] = i == j; } } while (p) { if (p & 1) (res *= (*this)) %= MOD; ((*this) *= (*this)) %= MOD; p >>= 1; } for (size_t i = 0; i < mat.size(); ++i) mat[i] = res[i]; return *this; } MatrixMOD operator+(const MatrixMOD & r) const { MatrixMOD res(mat); return res += r; } MatrixMOD operator-(const MatrixMOD & r) const { MatrixMOD res(mat); return res -= r; } MatrixMOD operator*(const MatrixMOD & r) const { MatrixMOD res(mat); return res *= r; } MatrixMOD operator*(const T & r) const { MatrixMOD res(mat); return res *= r; } MatrixMOD operator^(const int& p) const { MatrixMOD res(mat); return res ^= p; } MatrixMOD t() const { MatrixMOD res(mat[0].size(), mat.size(), 0); for (size_t i = 0; i < mat[0].size(); ++i) { for (size_t j = 0; j < mat.size(); ++j) { res[i][j] = mat[j][i]; } } return res; } static MatrixMOD getUnit(size_t n) { MatrixMOD res(n, n, 0); for (size_t i = 0; i < n; ++i) res[i][i] = 1; return res; } void show() const { for (const auto& r : mat) { for (const auto& x : r) { std::cerr << x << "\t"; } std::cerr << std::endl; } } }; ll powMod(ll n, ll p) { ll res = 1; while (p) { if (p & 1) (res *= n) %= MOD; (n *= n) %= MOD; p >>= 1; } return res; } std::map inv; ll invMod(ll n) { if (inv.count(n)) return inv[n]; return inv[n] = powMod(n, MOD - 2); } signed main() { INIT; VAR(i64, n); MatrixMOD A(6, 6); i64 inv6 = invMod(6); REP(j, 6) A[0][j] = inv6; REP(i, 5) A[i + 1][i] = 1; A ^= n; OUT(A[0][0])BR; return 0; }