n = int(input()) # ペアをひと固まりと考えると、各グループ、少なくとも一つ以上の塊が必要 # 夫婦で一緒になる組xを固定して、xのスターリング数を求める。 # 残った人は夫婦同氏同じグループにならないように割り振る。 mod = 10 ** 9 + 7 fac=[1]*(n+1) finv=[1]*(n+1) inv=[1]*(n+1) for i in range(2, n + 1): # p = i*p // i + p % i -> i*p//i+p%i=0 mod p # inv[i] = -i*p//i * inv[p%i] inv[i] = mod - mod // i * inv[mod % i] % mod for i in range(1, n + 1): fac[i] = fac[i - 1] * i % mod finv[i] = finv[i - 1] * inv[i] % mod def comb(n, k): if n < 0 or k < 0 or n < k: return 0 return fac[n] * finv[n - k] % mod * finv[k] % mod dp = [[0] * (n + 1) for _ in range(n + 1)] dp[0][0] = 1 for i in range(1, n + 1): for j in range(1, n + 1): dp[i][j] = j * dp[i - 1][j] % mod + dp[i - 1][j - 1] if dp[i][j] >= mod: dp[i][j] -= mod ans = 0 for g in range(1, n + 1): for x in range(g, n + 1): y = n - x z = dp[x][g] tmp = 1 for w in range(y): tmp *= g * (g - 1) % mod tmp %= mod ans += tmp * dp[x][g] % mod * comb(n, x) % mod if ans >= mod: ans -= mod print(ans)