#include using namespace std; constexpr int mod = 1000000007; int dp[5050][5050]; void mpl(int &x,int y) { x += y; if(x >= mod) x -= mod; } int main() { int N; cin >> N; vectorcnt(5000); for(int i = 0; i < N; i++) { int a; cin >> a; cnt[a]++; } dp[0][0] = 1; for(int i = 0; i < N; i++) { for(int j = 0; j < N; j++) { mpl(dp[i+1][j],dp[i][j]); mpl(dp[i+1][j+1],1ll*dp[i][j]*cnt[i]%mod); } } vectorf(N+1); f[0] = 1; for(int i = 1; i <= N; i++) { f[i] = 1ll*f[i-1]*i%mod; } int ans = 0; for(int i = 0; i <= N; i++) { if(i%2 == 0) { mpl(ans,1ll*dp[N][i]*f[N-i]%mod); } else { mpl(ans,mod-1ll*dp[N][i]*f[N-i]%mod); } } cout << ans << endl; }