# include "bits/stdc++.h" using namespace std; using LL = long long; using ULL = unsigned long long; const double PI = acos(-1); templateconstexpr T INF() { return ::std::numeric_limits::max(); } templateconstexpr T HINF() { return INF() / 2; } template T_char TL(T_char cX) { return tolower(cX); }; template T_char TU(T_char cX) { return toupper(cX); }; const int vy[] = { -1, -1, -1, 0, 1, 1, 1, 0 }, vx[] = { -1, 0, 1, 1, 1, 0, -1, -1 }; const int dx[4] = { 0,1,0,-1 }, dy[4] = { 1,0,-1,0 }; int popcnt(unsigned long long n) { int cnt = 0; for (int i = 0; i < 64; i++)if ((n >> i) & 1)cnt++; return cnt; } int d_sum(LL n) { int ret = 0; while (n > 0) { ret += n % 10; n /= 10; }return ret; } int d_cnt(LL n) { int ret = 0; while (n > 0) { ret++; n /= 10; }return ret; } LL gcd(LL a, LL b) { if (b == 0)return a; return gcd(b, a%b); }; LL lcm(LL a, LL b) { LL g = gcd(a, b); return a / g*b; }; # define ALL(qpqpq) (qpqpq).begin(),(qpqpq).end() # define UNIQUE(wpwpw) sort(ALL((wpwpw)));(wpwpw).erase(unique(ALL((wpwpw))),(wpwpw).end()) # define LOWER(epepe) transform(ALL((epepe)),(epepe).begin(),TL) # define UPPER(rprpr) transform(ALL((rprpr)),(rprpr).begin(),TU) # define FOR(i,tptpt,ypypy) for(LL i=(tptpt);i<(ypypy);i++) # define REP(i,upupu) FOR(i,0,upupu) # define INIT std::ios::sync_with_stdio(false);std::cin.tie(0) //行列は二次元vectorを用いることにする typedef vector vec; typedef vector mat; //A*Bの計算 mat mul(mat &A, mat&B, int M) { mat C((int)A.size(), vec(B[0].size())); for (int i = 0; i < (int)A.size(); i++) { for (int k = 0; k < (int)B.size(); k++) { for (int j = 0; j < (int)B[0].size(); j++) { C[i][j] = (C[i][j] + A[i][k] * B[k][j]) % M; } } } return C; } //A^nの計算 mat pow(mat A, LL n, int M) { mat B((int)A.size(), vec(A.size())); for (int i = 0; i < (int)A.size(); i++) { B[i][i] = 1; } while (n > 0) { if (n & 1)B = mul(B, A, M); A = mul(A, A, M); n >>= 1; } return B; } LL n; const LL mod = 1e9 + 7; int main() { cin >> n; mat A(2, vec(1)); A[0][0] = 1, A[1][0] = 1; mat B(2, vec(2)); B[0][0] = 1, B[0][1] = 1; B[1][0] = 1, B[1][1] = 0; B = pow(B, n - 1, mod); A = mul(B, A, mod); cout << A[0][0] * A[1][0] % mod << endl; }