#include using namespace std; using ll = long long; constexpr char newl = '\n'; // https://noshi91.hatenablog.com/entry/2019/03/31/174006 template struct ModInt { using u64 = std::uint_fast64_t; static constexpr u64 MOD = Modulus; u64 val; constexpr ModInt(const u64 x = 0) noexcept : val(x % MOD) {} constexpr ModInt operator+() const noexcept { return ModInt(*this); } constexpr ModInt operator-() const noexcept { ModInt res(*this); if (res.val != 0) res.val = MOD - res.val; return res; } constexpr bool operator==(const ModInt& rhs) const noexcept { return val == rhs.val; } constexpr bool operator!=(const ModInt& rhs) const noexcept { return val != rhs.val; } // prefix increment/decrement constexpr ModInt& operator++() noexcept { return *this += ModInt(1); } constexpr ModInt& operator--() noexcept { return *this -= ModInt(1); } // postfix increment/decrement constexpr ModInt operator++(int) noexcept { ModInt tmp(*this); ++*this; return tmp; } constexpr ModInt operator--(int) noexcept { ModInt tmp(*this); --*this; return tmp; } constexpr ModInt operator+(const ModInt& rhs) const noexcept { return ModInt(*this) += rhs; } constexpr ModInt operator-(const ModInt& rhs) const noexcept { return ModInt(*this) -= rhs; } constexpr ModInt operator*(const ModInt& rhs) const noexcept { return ModInt(*this) *= rhs; } constexpr ModInt operator/(const ModInt& rhs) const noexcept { return ModInt(*this) /= rhs; } constexpr ModInt& operator+=(const ModInt& rhs) noexcept { val += rhs.val; if (val >= MOD) val -= MOD; return *this; } constexpr ModInt& operator-=(const ModInt& rhs) noexcept { if (val < rhs.val) val += MOD; val -= rhs.val; return *this; } constexpr ModInt& operator*=(const ModInt& rhs) noexcept { val = val * rhs.val % MOD; return *this; } // prime Modulus only constexpr ModInt& operator/=(const ModInt& rhs) noexcept { return *this *= rhs.inv(); } // prime Modulus only constexpr ModInt inv() const noexcept { return pow(*this, MOD - 2); } }; template constexpr ModInt pow(ModInt x, std::uint_fast64_t n) { ModInt res(1); while (n) { if (n & 1) res *= x; x *= x; n >>= 1; } return res; } template istream& operator>>(istream& is, ModInt& x) { std::uint_fast64_t val; is >> val; x = ModInt(val); return is; } template ostream& operator<<(ostream& os, const ModInt& x) { return os << x.val; } using mint = ModInt<1000000007>; using vec = vector; using mat = vector; mat mul(const mat& A, const mat& B) { mat C(A.size(), vec(B[0].size(), 0)); for (int i = 0; i < A.size(); i++) { for (int k = 0; k < B.size(); k++) { for (int j = 0; j < B[0].size(); j++) { C[i][j] += A[i][k] * B[k][j]; } } } return C; } mat pow(mat A, ll n) { mat B(A.size(), vec(A.size(), 0)); for (int i = 0; i < A.size(); i++) { B[i][i] = 1; } while (n > 0) { if (n & 1) B = mul(B, A); A = mul(A, A); n >>= 1; } return B; } int main() { cin.tie(nullptr); ios::sync_with_stdio(false); mat A(6, vec(6, 0)); for (int i = 1; i < 6; i++) { A[i - 1][i] = 1; } A.back().assign(6, 166666668); ll n; cin >> n; A = pow(A, n); cout << A[5][5] << newl; return 0; }