import sys #input = sys.stdin.readline #input = sys.stdin.buffer.readline #文字列はダメ #sys.setrecursionlimit(1000000) #import bisect #import itertools #import random #from heapq import heapify, heappop, heappush from collections import defaultdict #from collections import deque #import copy import math #from functools import lru_cache #@lru_cache(maxsize=None) #MOD = pow(10,9) + 7 #MOD = 998244353 #dx = [1,0,-1,0] #dy = [0,1,0,-1] #dx8 = [1,1,0,-1,-1,-1,0,1] #dy8 = [0,1,1,1,0,-1,-1,-1] #dx = [1,1,-1,-1] #dy = [1,-1,1,-1] def main(): N = int(input()) X = list(map(int,input().split())) c = 1 p = 1 for i in range(N): dic = defaultdict(int) mn = N for x in X: dic[x] += 1 mn = min(mn, x) c *= dic[mn] p *= N-i g = math.gcd(c,p) c //= g p //= g X.remove(mn) #print(c,p) print(p) if __name__ == '__main__': main()