#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
int N,K,A[1000];
long mod=1e9+7,ans=1;
long power(long a,long b){return b?power(a*a%mod,b/2)*(b%2?a:1)%mod:1;}
main()
{
	cin>>N>>K;
	for(int i=0;i<N;i++)cin>>A[i];
	for(int i=2;i*i<=1000000000;i++)
	{
		vector<int>cnt(N);
		for(int j=0;j<N;j++)
		{
			while(A[j]%i==0)
			{
				A[j]/=i;
				cnt[j]++;
			}
		}
		sort(cnt.begin(),cnt.end());
		int x=0;
		for(int j=0;j<K;j++)x+=cnt[N-j-1];
		ans=ans*power(i,x)%mod;
	}
	for(int i=0;i<N;i++)
	{
		if(A[i]==1)continue;
		int a=A[i];
		int x=0;
		for(int j=0;j<N;j++)
		{
			if(A[j]==a)
			{
				A[j]=1;
				x++;
			}
		}
		ans=ans*power(a,x<K?x:K)%mod;
	}
	cout<<ans<<endl;
}