#include #define int long long using namespace std; int n,dp[1005][1005][2][2],mod=1e9+7; signed main(){ ios::sync_with_stdio(0); cin.tie(0),cout.tie(0); cin>>n;dp[1][0][0][0]=1,dp[2][0][0][0]=2; for(int i=3;i<=n;i++){ for(int j=0;j<=i-2;j++){ dp[i][j][0][0]=(dp[i-1][j][0][0]*(i-j-2)%mod+dp[i-1][j][1][0]*(i-j-1)%mod+dp[i-1][j+1][1][0]*j+dp[i-1][j+1][0][0]*(j+1)+dp[i-1][j+1][0][1]+dp[i-1][j+1][1][1])%mod; dp[i][j][1][0]=(dp[i-1][j][0][1]*(i-j-2)%mod+dp[i-1][j][1][1]*(i-j-1)+dp[i-1][j+1][1][1]*max(j-1,0ll)%mod+dp[i-1][j+1][0][1]*j%mod)%mod; dp[i][j][0][1]=(((j)?dp[i-1][j-1][0][0]*2+dp[i-1][j-1][1][0]:0)+dp[i-1][j][1][0])%mod; dp[i][j][1][1]=(((j)?dp[i-1][j-1][0][1]*2+dp[i-1][j-1][1][1]:0)+dp[i-1][j][1][1])%mod; } } cout<