#include // c #include // io #include #include #include #include // container #include #include #include #include #include #include // other #include #include #include #include #include using namespace std; using i32=int32_t;using i64=int64_t;using ll =i64;using uint=uint32_t;using ull=uint64_t; template using matrix=vector >; #define ALL(c) (begin(c)),(end(c)) #define REP(i,n) FOR(i,0,n) #define REPr(i,n) FORr(i,0,n) #define FOR(i,l,r) for(int i=(int)(l);i<(int)(r);++i) #define FORr(i,l,r) for(int i=(int)(r)-1;i>=(int)(l);--i) #define EACH(it,o) for(auto it = (o).begin(); it != (o).end(); ++it) #define IN(l,v,r) ((l)<=(v) && (v)<(r)) #define UNIQUE(v) v.erase(unique(ALL(v)),v.end()) //debug #define DUMP(x) cerr << #x << " = " << (x) #define LINE() cerr<< " (L" << __LINE__ << ")" class range { private: struct Iter{ int v; int operator*(){return v;} bool operator!=(Iter& itr) {return v < itr.v;} void operator++() {++v;} }; Iter i, n; public: range(int n) : i({0}), n({n}) {} range(int i, int n) : i({i}), n({n}) {} Iter& begin() {return i;} Iter& end() {return n;} }; //output template ostream& operator << (ostream& os, const vector& as){REP(i,as.size()){if(i!=0)os<<" "; os< ostream& operator << (ostream& os, const vector>& as){REP(i,as.size()){if(i!=0)os< ostream& operator << (ostream& os, const set& ss){for(auto a:ss){if(a!=ss.begin())os<<" "; os< ostream& operator << (ostream& os, const pair& p){os< ostream& operator << (ostream& os, const map& m){bool isF=true;for(auto& p:m){if(!isF)os< ostream& operator << (ostream& os, const tuple& t){os << get<0>(t);return os;} template ostream& operator << (ostream& os, const tuple& t){os << get<0>(t)<<" "<(t);return os;} template ostream& operator << (ostream& os, const tuple& t){os << get<0>(t)<<" "<(t)<<" "<(t);return os;} template ostream& operator << (ostream& os, const tuple& t){os << get<0>(t)<<" "<(t)<<" "<(t)<<" "<(t);return os;} template ostream& operator << (ostream& os, const tuple& t){os << get<0>(t)<<" "<(t)<<" "<(t)<<" "<(t)<<" "<(t);return os;} template ostream& operator << (ostream& os, const tuple& t){os << get<0>(t)<<" "<(t)<<" "<(t)<<" "<(t)<<" "<(t)<<" "<(t);return os;} template ostream& operator << (ostream& os, const tuple& t){os << get<0>(t)<<" "<(t)<<" "<(t)<<" "<(t)<<" "<(t)<<" "<(t)<<" "<(t);return os;} //input char tmp[1000]; #define nextInt(n) scanf("%d",&n) #define nextLong(n) scanf("%lld",&n) //I64d #define nextDouble(n) scanf("%lf",&n) #define nextChar(n) scanf("%c",&n) #define nextString(n) scanf("%s",tmp);n=tmp // values template T INF(){assert(false);}; template<> int INF(){return 1<<28;}; template<> ll INF(){return 1LL<<58;}; template<> double INF(){return 1e16;}; template T EPS(){assert(false);}; template<> int EPS(){return 1;}; template<> ll EPS(){return 1LL;}; template<> double EPS(){return 1e-8;}; template<> long double EPS(){return 1e-8;}; template T pmod(T v,U M){return (v%M+M)%M;} ll gcd_positive(ll a,ll b) { return b == 0 ? a : gcd_positive(b,a%b); } ll gcd(ll a,ll b) { return gcd_positive(abs(a), abs(b)); } ll lcm(ll a,ll b){return a/gcd(a,b)*b;} ll MOD=1e9+7; template class modU{ public: ll pmod(ll v){return (v%M+M)%M;} ll ainv(ll a){return pmod(-a);} ll add(ll a,ll b){return pmod(a+b);} ll mul(ll a,ll b){return pmod(a*b);} ll pow(ll x, ll N){ ll res=1; while(N>0){ if(N%2)res=mul(res,x); x=mul(x,x); N/=2; } return res; } //O(logM) ll minv(ll a){return pow(a,M-2);} // O(M) vector minvs(){ vector inv(M + 1);inv[1] = 1; for (int i = 2; i <= M; ++i)inv[i] = mul(inv[M % i],M - M / i); return inv; } //O(N) 条件 N>=0 vector factmemo=vector(2000050,-1); inline ll fact(const int N){ if(factmemo[N]!=-1)return factmemo[N]; return factmemo[N]=N==0?1:mul(N,fact(N-1)); } //条件 N>=0 inline ll nPr(const int N,const int r){ if(r>N)return 0; return mul(fact(N),minv(fact(N-r))); } inline ll nCr(const int N,const int r){ if(N<0){ if(0 <= r) return (r%2?-1:+1)*nCr(-N +r-1,r); else if(r<=N) return ((N-r)%2?-1:+1)*nCr(-r-1,N-r); else return 0; } if(0<= r && r<=N) return mul(mul(fact(N),minv(fact(r))),minv(fact(N-r))); return 0; } }; class Main{ public: modU<> mu; // |A| -> |B| の全射の個数 vector> surjective_count(int A,int B){ vector> dp(A,vector(B)); dp[0][0]=1; REP(a,A)REP(b,B){ if(a-1>=0 && b-1>=0){ // create a new point dp[a][b] += dp[a-1][b-1]; dp[a][b] %= MOD; } if(a-1>=0){// add to an already exist point dp[a][b] += b * dp[a-1][b]; dp[a][b] %= MOD; } } return dp; } void run(){ int N;cin >> N; vector> dp=surjective_count(N+1,N+1); ll res = 0; FOR(i,1,N+1){//create i groups FOR(j,i,N+1){//夫婦 j pairs ll val = mu.mul(mu.nCr(N,j),dp[j][i]); val=mu.mul(val,mu.pow(i*(i-1),N-j)); res=mu.add(res, val); } } cout << res <