#ifndef ONLINE_JUDGE
#define _GLIBCXX_DEBUG
#endif
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace std;
using ll=long long;
using ld=long double;
using st=string;
using P=pair<ll,ll>;
typedef atcoder::modint mint;
ll inf=9e18;
template<class... T>
constexpr auto max(T... a){return max(initializer_list<common_type_t<T...>>{a...});}
template<class... T>
constexpr auto min(T... a){return min(initializer_list<common_type_t<T...>>{a...});}
template<typename T, int s, int i = 0>
auto vec(const ll (&sizes)[s], const T& init = T()){
  if constexpr(i < s) return vector(sizes[i], vec<T, s, i+1>(sizes, init));
  else return init;
}



int main(){
  cin.tie(nullptr);
  ios_base::sync_with_stdio(false);
  
  auto v=vec<mint>({1001},0);
  mint x=1;
  for(ll i=1;i<=1000;i++){
    x*=i;
    v[i]=x;
  }

  ll n;
  cin>>n;
  auto dp=vec<ll>({n+1,n+1,3},0);
  dp[0][0][1]=1;
  for(ll i=1;i<=n;i++){
    for(ll j=0;j<=n;j++){
      if(j>0) dp[i][j][0]=dp[i-1][j-1][1]+dp[i-1][j-1][2];
      dp[i][j][1]=dp[i-1][j][0];
      dp[i][j][2]=dp[i-1][j][1];
    }
  }

  mint ans=0;
  for(ll j=0;j<=n;j++){
    ans+=v[j]*v[n-j]*(dp[n][j][0]+dp[n][j][1]);
  }
  if(n!=1) cout<<ans.val()<<"\n";
  else cout<<"1\n";

  // for(ll j=0;j<n+1;j++){
  //   for(ll i=0;i<n+1;i++){
  //     cout<<dp[i][j][0]<<dp[i][j][1]<<dp[i][j][2]<<" ";
  //   }
  //   cout<<"\n";
  // }
}