#!/usr/bin/env python3.8
# %%
import sys
read = sys.stdin.buffer.read
readline = sys.stdin.buffer.readline
readlines = sys.stdin.buffer.readlines
import numpy as np
from collections import defaultdict

MOD = 10 ** 9 + 7

# %%
N, K, *A = map(int, read().split())


# %%
U = int(10 ** 4.5) + 10
is_prime = np.zeros(U, np.bool)
is_prime[3::2] = 1
is_prime[2] = 1
for p in range(3, U, 2):
    if p * p >= U:
        break
    if is_prime[p]:
        is_prime[p * p:: p + p] = 0
primes = np.where(is_prime)[0]


# %%
def factor(N):
    pf = primes[N % primes == 0].tolist()
    for p in pf:
        e = 0
        while True:
            N //= p
            e += 1
            if N % p != 0:
                break
        yield (p, e)
    if N > 1:
        yield (N, 1)


# %%
multiplicity = defaultdict(list)
for x in A:
    for p, e in factor(x):
        multiplicity[p].append(e)

# %%
answer = 1
for p, E in multiplicity.items():
    if len(E) <= K:
        e = sum(E)
    else:
        E.sort(reverse=True)
        e = sum(E[:K])
    answer *= pow(p, e, MOD)
print(answer % MOD)