#include using namespace std; #define pb push_back #define mp make_pair #define all(x) (x).begin(),(x).end() #define rep(i,n) for(int i=0;i<(n);i++) constexpr int MOD = 1000000007; typedef long long ll; typedef unsigned long long ull; typedef pair pii; constexpr int dx[] = {1, 0, -1, 0, 1, 1, -1, -1}; constexpr int dy[] = {0, -1, 0, 1, 1, -1, -1, 1}; template ostream &operator<<(ostream &os, const vector &vec){os << "["; for (const auto &v : vec) {os << v << ","; } os << "]"; return os; } const int mod = 1e9 + 7; ll mod_pow(ll x, ll n, ll mod) { ll res = 1; while(n > 0) { if(n & 1) res = res * x % mod; x = x * x % mod; n >>= 1; } return res; } ll a(int n) { ll fact2N = 1; for(int i = 1; i <= 2 * n; i++) { fact2N *= i; fact2N %= mod; } ll power2 = 1; for(int i = 0; i < n; i++) { power2 *= 2; power2 %= mod; } return (fact2N * mod_pow(power2, mod - 2, mod)) % mod; } void solve() { /* // The experiment code -> 1, 6, 90, 2520, 113400... // found: A000680 (OEIS) // The answer is fact(N) * A000680 int N; cin >> N; if (N > 10) { exit(1); return; } vector> dp(1 << N, vector(1 << N)); dp[0][0] = 1; for(int i = 0; i < 1 << N; i++) { for(int j = 0; j < 1 << N; j++) { for(int k = 0; k < N; k++) { if (!(i >> k & 1)) { dp[i | (1 << k)][j] += dp[i][j]; dp[i | (1 << k)][j] %= mod; } } for(int k = 0; k < N; k++) { if ((i >> k & 1) && !(j >> k & 1)) { dp[i][j | (1 << k)] += dp[i][j]; dp[i | (1 << k)][j] %= mod; } } } } cout << (dp.back().back()) % mod << endl; */ int N; cin >> N; ll factN = 1; for(int i = 1; i <= N; i++) { factN *= i; factN %= mod; } cout << (a(N) * factN) % mod << endl; } int main() { std::cin.tie(0); std::ios::sync_with_stdio(false); cout.setf(ios::fixed); cout.precision(16); solve(); return 0; }