import sys
from collections import defaultdict

sys.setrecursionlimit(10 ** 6)
int1 = lambda x: int(x) - 1
p2D = lambda x: print(*x, sep="\n")
def II(): return int(sys.stdin.readline())
def MI(): return map(int, sys.stdin.readline().split())
def LI(): return list(map(int, sys.stdin.readline().split()))
def LLI(rows_number): return [LI() for _ in range(rows_number)]
def SI(): return sys.stdin.readline()[:-1]

def pf(x):
    bb=[]
    ee=[]
    cnt2=(x&-x).bit_length()-1
    if cnt2:
        bb.append(2)
        ee.append(cnt2)
        x>>=cnt2
    b=3
    while b**2<=x:
        e=0
        while x%b==0:
            x//=b
            e+=1
        if e:
            bb.append(b)
            ee.append(e)
        b+=2
    if x>1:
        bb.append(x)
        ee.append(1)
    return bb,ee

def main():
    md=10**9+7
    n,k=MI()
    aa=LI()
    be=defaultdict(list)
    for a in aa:
        bb,ee=pf(a)
        for b,e in zip(bb,ee):
            be[b].append(e)
    ans=1
    for b,ee in be.items():
        ee.sort(reverse=True)
        s=sum(ee[:k])
        ans*=pow(b,s,md)
        ans%=md
    print(ans)

main()