N = int(input()) A = list(map(int, input().split())) B = [0] * N B[0] = A[0] for i in range(1, N): B[i] = pow(A[i], i + 1, 1000000007) B[i] *= B[i - 1] B[i] %= 1000000007 j = 0 t = 1 a = [] for i in range(N): t *= A[i] while t > 1000000007: t //= A[j] j += 1 if j == 0: a.append(B[i]) else: a.append(B[i] * pow(B[j - 1], 1000000007 - 2, 1000000007) * pow(pow(A[i], j, 1000000007), 1000000007 - 2, 1000000007)) result = 1 for i in range(len(a)): result *= a[i] result %= 1000000007 print(result)