#include<bits/stdc++.h>
using namespace std;

#define int long long
typedef pair<int,int>pint;
typedef vector<int>vint;
typedef vector<pint>vpint;
#define pb push_back
#define mp make_pair
#define fi first
#define se second
#define all(v) (v).begin(),(v).end()
#define rep(i,n) for(int i=0;i<(n);i++)
#define reps(i,f,n) for(int i=(f);i<(n);i++)
#define each(it,v) for(__typeof((v).begin()) it=(v).begin();it!=(v).end();it++)
template<class T,class U>inline void chmin(T &t,U f){if(t>f)t=f;}
template<class T,class U>inline void chmax(T &t,U f){if(t<f)t=f;}

int N,K;
int A[1000];

vint prime;
int f[111111];

const int mod=1000000007;
int mpow(int n,int m){
    int ret=1;
    while(m){
        if(m&1)ret=ret*n%mod;
        n=n*n%mod;
        m>>=1;
    }
    return ret;
}


signed main(){
    fill_n(f,111111,true);
    f[0]=f[1]=false;
    for(int i=2;i<111111;i++){
        if(!f[i])continue;
        prime.pb(i);
        for(int j=i+i;j<111111;j+=i)f[j]=false;
    }

    cin>>N>>K;
    map<int,int>hoge;
    rep(i,N){
        cin>>A[i];
        int a=A[i];
        rep(j,prime.size()){
            while(a%prime[j]==0)a/=prime[j];
        }
        if(a!=1)hoge[a]++;
    }

    int ans=1;
    rep(i,prime.size()){
        vint v;
        rep(j,N){
            int a=A[j];
            v.pb(0);
            while(a%prime[i]==0){
                v[v.size()-1]++;
                a/=prime[i];
            }
        }
        sort(all(v));
        int sum=0;
        rep(j,K)sum+=v[N-1-j];
        ans=ans*mpow(prime[i],sum)%mod;
    }
    each(it,hoge)ans=ans*mpow(it->fi,it->se)%mod;
    cout<<ans<<endl;
    return 0;
}