#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
#define REP(i, n)           for(int(i)=0;(i)<(n);++(i))

const int MOD = 1000000007;
ll dp[1010][1010][2][2];

int main(){
    int N;
    cin >> N;

    dp[0][0][0][0] = 1;
    REP(n,N){
        const int mm = n+1; // 置ける場所の数
        REP(k,n+1){
            int m = mm;
            if(n<2){
                // 現在の置いている数が2個未満ならどこに置いてもいい
                dp[n+1][max(k-1,0)][0][0] += dp[n][k][0][0] * k;        // 減少
                dp[n+1][k         ][0][0] += dp[n][k][0][0] * (m-k);    // 維持
                continue;
            }
            m = mm;
            dp[n+1][k+1       ][0][1] += dp[n][k][0][0] * 2; m-=2;  // 増加 [n-1]の両隣2個
            dp[n+1][max(k-1,0)][0][0] += dp[n][k][0][0] * k; m-=k;  // 減少 (k個)
            dp[n+1][k         ][0][0] += dp[n][k][0][0] * m;        // 維持

            if(k == 0) continue;
            // [n-2,n]がある場合
            m = mm;
            dp[n+1][k+1][1][1] += dp[n][k][0][1] * 2;     m-=2;     // 増加 [n-1]の両隣2個
            dp[n+1][k-1][0][0] += dp[n][k][0][1] * 1;     m-=1;     // 減少 [n-2,n]の間
            dp[n+1][k-1][1][0] += dp[n][k][0][1] * (k-1); m-=k-1;   // 減少 [n-2,n]以外の間(k-1個)
            dp[n+1][k  ][1][0] += dp[n][k][0][1] * m;               // 維持

            // [n-3,n-1]がある場合
            m = mm;
            dp[n+1][k+1][0][1] += dp[n][k][1][0] * 1;     m-=1;     // 増加 [n-1]の片隣1個
            dp[n+1][k  ][0][1] += dp[n][k][1][0] * 1;     m-=1;     // 維持 [n-1]の片隣かつ[n-1,n-3]の間1個
            dp[n+1][k-1][0][0] += dp[n][k][1][0] * (k-1); m-=k-1;   // 減少 [n-3,n-1]以外の間(k-1個)
            dp[n+1][k  ][0][0] += dp[n][k][1][0] * m;               // 維持

            // [n-2,n],[n-3,n-1]がある場合
            m = mm;
            dp[n+1][k+1][1][1] += dp[n][k][1][1] * 1;     m-=1;     // 増加 [n-1]の片隣1個
            dp[n+1][k  ][1][1] += dp[n][k][1][1] * 1;     m-=1;     // 維持 [n-1]の片隣かつ[n-1,n-3]の間1個
            dp[n+1][k-1][0][0] += dp[n][k][1][1] * 1;     m-=1;     // 減少 [n-2,n]の間
            dp[n+1][k-1][1][0] += dp[n][k][1][1] * (k-2); m-=k-2;   // 減少 [n-2,n],[n-3,n-1]以外の間(k-2個)
            dp[n+1][k  ][1][0] += dp[n][k][1][1] * m;               // 維持
        }

        REP(k,n+1){
            REP(a,2)REP(b,2) dp[n+1][k][a][b] %= MOD;
        }
    }

    cout << dp[N][0][0][0] << endl;
    return 0;
}