# 包除原理の典型問題
# 解説を見てしまったが、これを解けないのはダメ。

N=int(input())
A=[int(input()) for i in range(N)]
C=[0]*N
for a in A:
    if a<N:
        C[a]+=1

mod=10**9+7

DP=[0]*(N+1)
DP[0]=1

for v in C:
    if v==0:
        continue
    for i in range(N-1,-1,-1):
        DP[i+1]+=DP[i]*v
        DP[i+1]%=mod      
ANS=0
fac=1
for i in range(N+1):
    if i%2==0:
        ANS+=fac*DP[N-i]
    else:
        ANS-=fac*DP[N-i]

    fac=fac*(i+1)%mod

    ANS%=mod

if N%2==1:
    print((-ANS)%mod)
else:
    print(ANS)