#include <bits/stdc++.h>

using namespace std;

const int MAXN = 1e3 + 10, P = 1e9 + 7;

int n;
long long dp[MAXN][MAXN][2][2];

int main(){
  cin >> n;
  dp[1][0][0][0] = 1, dp[2][0][0][0] = 2;
  for(int i = 2; i <= n; i++){
    for(int j = 0; j <= i; j++){
      (dp[i + 1][j + 1][0][1] += dp[i][j][0][0] * 2) %= P;
      (dp[i + 1][j - 1][0][0] += j * dp[i][j][0][0]) %= P;
      (dp[i + 1][j][0][0] += (i + 1 - j - 2) * dp[i][j][0][0]) %= P;

      (dp[i + 1][j + 1][1][1] += dp[i][j][0][1] * 2) %= P;
      (dp[i + 1][j - 1][1][0] += (j - 1) * dp[i][j][0][1]) %= P;
      (dp[i + 1][j - 1][0][0] += dp[i][j][0][1]) %= P;
      (dp[i + 1][j][1][0] += (i + 1 - (j - 1) - 1 - 2) * dp[i][j][0][1]) %= P;

      (dp[i + 1][j + 1][1][1] += dp[i][j][1][1]) %= P;
      (dp[i + 1][j + 1][1][1] += dp[i][j][1][1]) %= P;
      (dp[i + 1][j - 1][1][1] += (j - 1) * dp[i][j][1][1]) %= P;
      (dp[i + 1][j - 1][0][1] += dp[i][j][1][1]) %= P;
      (dp[i + 1][j][1][0] += (i + 1 - (j - 1) - 1 - 2) * dp[i][j][1][1]) %= P;

      (dp[i + 1][j + 1][0][1] += dp[i][j][1][0] * 2) %= P;
      (dp[i + 1][j - 1][0][0] += j * dp[i][j][1][0]) %= P;
      (dp[i + 1][j][1][0] += (i + 1 - j - 2) * dp[i][j][1][0]) %= P;
    }
  }
  cout << dp[n][0][0][0];
  return 0;
}