#include<bits/stdc++.h>
#define int long long
using namespace std;
const int N=5e3+10;
const int mod=1e9+7;
int n,a[N],t[N],dp[N],f[N];
signed main()
{
    ios::sync_with_stdio(0);
    cin.tie(0),cout.tie(0);
    cin>>n;
    for(int i=1;i<=n;i++)
    {
        cin>>a[i];
        a[i]++;
        t[a[i]]++;
    }
    dp[0]=1;
    for(int i=1;i<=n;i++)
    {
        memcpy(f,dp,sizeof(f));
        memset(dp,0,sizeof(dp));
        for(int j=0;j<=n;j++)
        {
            dp[j]=f[j];
            if(j!=0&&t[i]!=0)
            {
                dp[j]+=f[j-1]*t[i];
                dp[j]%=mod;
            }
        }
    }
    for(int i=0;i<=n;i++)
    {
        int us=n-i;
        for(int j=1;j<=us;j++)
        {
            dp[i]*=j;
            dp[i]%=mod;
        }
    }
    int ans=0,f=1;
    for(int i=0;i<=n;i++)
    {
        ans+=dp[i]*f;
        f*=-1;
        ans%=mod;
    }
    ans+=mod;
    ans%=mod;
    cout<<ans<<'\n';
    return 0;
}