#include <bits/stdc++.h>
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;
    vector<int>cnt(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);
        }
    }
    vector<int>f(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;
}