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

typedef long long ll;
const ll mo=1e9+7;

ll fac[305], ifac[305], c[305], val[305];
ll F[305], G[305], f[305];

ll qpow(ll a,ll p) {
	ll res=1;
	while(p) {
		if(p&1) res=res*a%mo;
		a=a*a%mo, p>>=1;
	}
	return res;
}

ll C(ll n,ll m) {
	if(n<0||m<0||n<m) return 0;
	return fac[n]*ifac[n-m]%mo*ifac[m]%mo;
}

int main() {
	ios::sync_with_stdio(false);
	cin.tie(0), cout.tie(0);
	ll n,K; cin>>n>>K;
	fac[0]=1;
	for(ll i=1;i<=n;i++) fac[i]=fac[i-1]*i%mo, val[i]=qpow(i,i-1);
	ifac[n]=qpow(fac[n],mo-2);
	for(ll i=n-1;i>=0;i--) ifac[i]=ifac[i+1]*(i+1)%mo;
	for(ll i=1,x;i<=n;i++) cin>>x, c[x]++;
	F[0]=1;
	for(ll i=1;i<=n;i++) {
		f[0]=1;
		for(ll j=1;j<=i;j++) {
			f[j]=0;
			for(ll k=1;k<=min(j,K);k++) f[j]=(f[j]-C(j-1,k-1)*f[j-k]%mo*val[k]%mo*i+mo)%mo;
		}
		F[i]=(mo-f[i]*qpow(i,mo-2)%mo)%mo;
	}
	ll ans=1;
	for(ll i=1;i<=n;i++) {
		if(!c[i]) continue;
		G[0]=1;
		for(ll j=1;j<=c[i];j++) {
			G[j]=0;
			for(ll k=1;k<=j;k++) G[j]=(G[j]+C(j-1,k-1)*G[j-k]%mo*F[k]%mo*n%mo)%mo;
		}
		ans=ans*G[c[i]]%mo;
	}
	cout<<ans*qpow(n,mo-3)%mo<<'\n';
	return 0;
}