#pragma GCC optimize("Ofast") #include using namespace std; #define rep(i, n) for(int i = 0; i < (n); ++i) #define all(x) (x).begin(),(x).end() #define ln '\n' const long long MOD = 1000000007LL; //const long long MOD = 998244353LL; typedef long long ll; typedef unsigned long long ull; typedef pair pii; typedef pair pll; template inline bool chmax(T &a, T b) { if (a < b) { a = b; return true;} return false; } template inline bool chmin(T &a, T b) { if (a > b) { a = b; return true;} return false; } /////////////////////////////////////////////////////////////////////////////////////////////////// template class modint { using u64 = std::uint_fast64_t; public: u64 a; constexpr modint(const u64 x = 0) noexcept : a(x % Modulus) {} constexpr u64 &value() noexcept { return a; } constexpr const u64 &value() const noexcept { return a; } 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 { a += rhs.a; if (a >= Modulus) { a -= Modulus; } return *this; } constexpr modint &operator-=(const modint rhs) noexcept { if (a < rhs.a) { a += Modulus; } a -= rhs.a; return *this; } constexpr modint &operator*=(const modint rhs) noexcept { a = a * rhs.a % Modulus; return *this; } constexpr modint &operator/=(modint rhs) noexcept { u64 exp = Modulus - 2; while (exp) { if (exp % 2) { *this *= rhs; } rhs *= rhs; exp /= 2; } return *this; } }; using mint = modint; template struct Matrix { array, 500> A; int N; Matrix() {} Matrix(int N) : N(N) {} inline const array &operator[](int k) const {return A[k];} inline array &operator[](int k) {return A[k];} static Matrix I(int N) { Matrix mat(N); for(int i = 0; i < N; ++i) mat[i][i] = 1; return mat; } Matrix &operator+=(const Matrix &B) { for(int i = 0; i < N; ++i) { for(int j = 0; j < N; ++j) { (*this)[i][j] += B[i][j]; } } return (*this); } Matrix &operator-=(const Matrix &B) { for(int i = 0; i < N; ++i) { for(int j = 0; j < N; ++j) { (*this)[i][j] -= B[i][j]; } } return (*this); } Matrix &operator*=(const Matrix &B) { Matrix C(N); for(int i = 0; i < N; ++i) { for(int j = 0; j < N; ++j) { for(int k = 0; k < N; ++k) { C[i][j] = (C[i][j] + (*this)[i][k] * B[k][j]); } } } for (int i = 0; i < N; ++i) { for (int j = 0; j < N; ++j) { swap(A[i][j],C.A[i][j]); } } return *this; } Matrix &operator^=(long long k) { Matrix B = Matrix::I(N); while(k > 0) { if(k & 1) B *= *this; *this *= *this; k >>= 1LL; } for (int i = 0; i < N; ++i) { for (int j = 0; j < N; ++j) { swap(A[i][j],B.A[i][j]); } } return *this; } Matrix operator+(const Matrix &B) const {return (Matrix(*this) += B);} Matrix operator-(const Matrix &B) const {return (Matrix(*this) -= B);} Matrix operator*(const Matrix &B) const {return (Matrix(*this) *= B);} Matrix operator^(const long long k) const {return (Matrix(*this) ^= k);} }; int main() { ios::sync_with_stdio(false); cin.tie(nullptr); ll N; cin >> N; Matrix mat(2); mat[0][0] = mat[0][1] = mat[1][0] = 1; mat[1][1] = 0; mat ^= N; cout << (mat[0][0]*mat[1][0]).a << ln; }