321
社区成员




n = int(input())
mod = int(1e9 + 7)
dp = [0] * max(n,2)
dm = dp + [0]
dp[0] = dm[1] = 1
dp[1] = 2
dm[2] = (mod - mod // 2) % mod
for i in range(2,n):
m = dm[i + 1] = (mod - mod // (i + 1)) * dm[mod % (i + 1)] % mod
dp[i] = (dp[i - 1] * (6 * i - 3) * m % mod - dp[i - 2] * (i - 2) * m % mod) % mod
print(dp[n - 1] % mod)