import sys
import numpy as np
input=lambda:sys.stdin.readline().rstrip()
mod=10**9+7
N=int(input())
if N<=3:
	print([1,2,2][N-1])
	sys.exit()
dp=np.zeros((N+1,N,2,2),dtype=int)
dp[4,2,1,1]=8
dp[4,1,0,1]=4
dp[4,1,1,0]=4
dp[4,0,0,0]=8
for i in range(4,N):
	for j in range(i):
		for k in range(2):
			for l in range(2):
				dp[i,j,k,l]%=mod
				dp[i+1,j+1,l,1]+=(2-k)*dp[i,j,k,l]
				dp[i+1,j-1,l,0]+=(j-k-l)*dp[i,j,k,l]
				dp[i+1,j,l,0]+=(i-j-1+k)*dp[i,j,k,l]
				dp[i+1,j-1,0,0]+=l*dp[i,j,k,l]
				dp[i+1,j,l,1]+=k*dp[i,j,k,l]
print(dp[N,0,0,0]%mod)