// Σ[k=1..N]{ C(N,k) * k^(N-k) }
const ll m=1e9+7;
ll N;
mint x,F[2d5],I[2d5];
{
	rd(N);
	*F=1;
	rep(i,1,N+1)F[i]=F[i-1]*i;
	I[N]=1/F[N];
	for(i=N-1;i>=0;--i)I[i]=I[i+1]*(i+1);
	rep(i,1,N+1){
		x+=F[N]*I[i]*I[N-i]*((mint)i).pw(N-i);
	}
	wt(x);
}