#include #define rep(i,a,b) for(int i=a;i struct ModInt { static const int Mod = MOD; unsigned x; ModInt() : x(0) { } ModInt(signed sig) { x = sig < 0 ? sig % MOD + MOD : sig % MOD; } ModInt(signed long long sig) { x = sig < 0 ? sig % MOD + MOD : sig % MOD; } int get() const { return (int)x; } ModInt &operator+=(ModInt that) { if ((x += that.x) >= MOD) x -= MOD; return *this; } ModInt &operator-=(ModInt that) { if ((x += MOD - that.x) >= MOD) x -= MOD; return *this; } ModInt &operator*=(ModInt that) { x = (unsigned long long)x * that.x % MOD; return *this; } ModInt &operator/=(ModInt that) { return *this *= that.inverse(); } ModInt operator+(ModInt that) const { return ModInt(*this) += that; } ModInt operator-(ModInt that) const { return ModInt(*this) -= that; } ModInt operator*(ModInt that) const { return ModInt(*this) *= that; } ModInt operator/(ModInt that) const { return ModInt(*this) /= that; } ModInt inverse() const { long long a = x, b = MOD, u = 1, v = 0; while (b) { long long t = a / b; a -= t * b; std::swap(a, b); u -= t * v; std::swap(u, v); } return ModInt(u); } bool operator==(ModInt that) const { return x == that.x; } bool operator!=(ModInt that) const { return x != that.x; } ModInt operator-() const { ModInt t; t.x = x == 0 ? 0 : Mod - x; return t; } }; template ModInt operator^(ModInt a, unsigned long long k) { ModInt r = 1; while (k) { if (k & 1) r *= a; a *= a; k >>= 1; } return r; } typedef ModInt<1000000007> mint; /*---------------------------------------------------------------------------------------------------             ∧_∧       ∧_∧  (´<_` )  Welcome to My Coding Space!      ( ´_ゝ`) /  ⌒i     /   \    | |     /   / ̄ ̄ ̄ ̄/  |   __(__ニつ/  _/ .| .|____      \/____/ (u ⊃ ---------------------------------------------------------------------------------------------------*/ int N; //--------------------------------------------------------------------------------------------------- mint dp[2020][2020]; mint sm[2020][2020]; void dodp() { rep(i, 1, N) { sm[i - 1][0] = 0; rep(j, 1, N + 1) sm[i - 1][j] = sm[i - 1][j - 1] + dp[i - 1][j]; rep(j, 1, N + 1) dp[i][j] = sm[i - 1][min(N, j + max(1, i - 1) - 1)]; } } //--------------------------------------------------------------------------------------------------- void _main() { cin >> N; // dp1: An = any -> ... -> An-1 = any rep(i, 0, 2020) rep(j, 0, 2020) dp[i][j] = 0; rep(j, 1, N + 1) dp[0][j] = 1; dodp(); mint ans = 0; rep(j, 1, N + 1) ans += dp[N - 1][j]; // dp2: An = 1 -> ... -> An-1 = N rep(i, 0, 2020) rep(j, 0, 2020) dp[i][j] = 0; dp[0][1] = 1; dodp(); ans -= dp[N - 1][N]; cout << ans.get() << endl; }