import math
N,K=map(int,raw_input().split())
A=map(int,raw_input().split())

upperBound = 100001
isPrime=[True for i in range((upperBound-3)/2+1)]
for i in range((int(math.sqrt(upperBound-3))-3)/2+1): # k*k = (2i+3)^2 <= upperBound-3 
   if isPrime[i]:
      k = i+i+3 # this is prime number
      # false out non_prime number from k*k
      # k*k = 4i^2+12i+9 = 2*(2i^2+6i+3)+3
      # k*(k+1) = 2*(2i^2+6i+3)+3 + 2i+3 ... even number. skip
      # k*(k+2) = 2*(2i^2+6i+3+k)+3 
      non_prime = k*(i+1)+i
      for j in range(non_prime,(upperBound-3)/2+1,k):
         isPrime[j]=False
primes = [2] + [i+i+3 for i in range((upperBound-3)/2+1) if isPrime[i]==True]
vec=[[] for i in range(101001)]
for a in A:
    for i,p in enumerate(primes):
        cnt=0
        while a%p==0:
            a/=p
            cnt+=1
        if cnt!=0:
            vec[i].append(cnt)
    if a!=1:
        vec[len(primes)].append(1)
        primes.append(a)
ans=1
mod=1000000007
for i,p in enumerate(primes):
    vec[i].sort()
    vec[i].reverse()
    for j,v in enumerate(vec[i]):
        if j>=K:
            break
        ans*=pow(p,v)
        ans%=mod
print ans