#include <bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=100010;
const int mod=1e9+7;
int ksm(int x,int k){
    int res=1;
    while(k){
        if(k&1)res=res*x%mod;
        x=x*x%mod;
        k/=2;
    }
    return res;
}
int ny(int x){
    return ksm(x,mod-2);
}
void add(int &x,int y){
    if((x+=y)>=mod)x-=mod;
}
void del(int &x,int y){
    if((x-=y)<0)x+=mod;
}
int inv[maxn],fac[maxn];
int C(int n,int m){return n==0?1:fac[n]*inv[n-m]%mod*inv[m]%mod;}
int A(int n,int m){return n==0?1:fac[n]*inv[n-m]%mod;}
void init(){
    inv[0]=fac[0]=1;
    inv[1]=1;
    for(int i=1;i<maxn;i++){
        fac[i]=fac[i-1]*i%mod;
    }
    inv[1]=1;
    for(int i=2;i<maxn;i++){
        inv[i]=(int)(mod-mod/i)*inv[mod%i]%mod;
    }
    inv[0]=1;
    for(int i=1;i<maxn;i++){
        inv[i]=inv[i-1]*inv[i]%mod;
    }
}
int f(int x){
    if(x==1)return 1;
    return ksm(x,x-2);
}
signed main()
{
    init();
    int n,k;cin>>n>>k;
    map<int,int> a;
    for(int i=1;i<=n;i++){
        int x;cin>>x;a[x]++;
    }
    if(a.size()==1){
        if(n==k){
            if(n==1)cout<<1;
            else cout<<ksm(n,n-2);
        }
        else cout<<0;
        return 0;
    }
    int ans=ksm(n,(int)a.size()-2);
    int N=n;
    for(auto [col,n]:a){
        vector<vector<int>>dp(n+1,vector<int>(n+1,0));
        dp[0][0]=1;
        // cout<<"col="<<col<<" n="<<n<<endl;
        for(int i=1;i<=n;i++){
            for(int j=1;j<=n;j++){
                for(int z=max(0ll,j-k);z<j;z++){
                    add(dp[i][j],dp[i-1][z]*C(n-z-1,j-z-1)%mod*f(j-z)%mod*(j-z)%mod);
                }
                // cout<<"dp["<<i<<"]["<<j<<"]="<<dp[i][j]<<endl;
            }
        }
        int res=0;
        for(int i=1;i<=n;i++){
            for(int j=n;j<=n;j++)add(res,dp[i][j]*ksm(N-n,i-1)%mod);
        }
        ans=ans*res%mod;
    }
    cout<<ans<<endl;
}