#include #define mm 1010 #define mod 1000000007 #define int long long using namespace std; int read() { int x=0,f=1; char ch=getchar(); while(ch<'0' || ch>'9') { if(ch=='-') f=-1; ch=getchar(); } while(ch>='0' && ch<='9') { x=(x<<3)+(x<<1)+(ch^48); ch=getchar(); } return x*f; } int n; int f[mm][mm][2][2]; signed main() { int n=read(); f[1][0][0][0]=1,f[2][0][0][0]=2; for(int i=2;i<=n-1;i++)for(int j=0;j0) { (f[i+1][j-1][0][0]+=f[i][j][0][0]*j+f[i][j][1][0]*(j-1)+f[i][j][0][1]+f[i][j][1][1])%=mod; (f[i+1][j][0][1]+=f[i][j][1][0])%=mod; (f[i+1][j-1][1][0]+=f[i][j][0][1]*(j-1)+f[i][j][1][1]*(j-2))%=mod; } (f[i+1][j][0][0]+=f[i][j][0][0]*(i-j-1)+f[i][j][1][0]*(i-j))%=mod; (f[i+1][j+1][0][1]+=f[i][j][0][0]*2+f[i][j][1][0])%=mod; (f[i+1][j][1][0]+=f[i][j][0][1]*(i-j-1)+f[i][j][1][1]*(i-j))%=mod; (f[i+1][j+1][1][1]+=f[i][j][0][1]*2+f[i][j][1][1])%=mod; (f[i+1][j][1][1]+=f[i][j][1][1])%=mod; } printf("%lld\n",f[n][0][0][0]); return 0; }