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

typedef long long T;
const int mod = 13313;

double PI=3.141592653589793;

complex <double> dft1[2100007]; //poczebne
complex <double> dft2[2100007]; //poczebne

complex <double> * dftn = dft1; //poczebne
complex <double> * dfts = dft2; //poczebne

complex <double> pom[2100007]; //do zwyklego

complex <double> a1[2100007]; //do dokladnego
complex <double> b1[2100007]; //do dokladnego
complex <double> a2[2100007]; //do dokladnego
complex <double> b2[2100007]; //do dokladnego

double cosi[2100007]; //poczebne
complex <double> omega[2100007]; //poczebne

inline int potenga(int v)
{
    for (int i=1; 1; i<<=1)
    {
        if (i>=v)
        {
            return i;
        }
    }
}

inline void dft(int n, int kier)
{
    int n2=n-1;
    int s=0;
    int p;
    int g;

    for (int i=2; i<=n; i<<=1)
    {
        swap(dftn, dfts);

        p=n/i;

        if (kier)
        {
            g=0;
            for (int j=0; j<n; j++)
            {
                dftn[j]=dfts[(2*(j-s)+s)&n2]+omega[g]*dfts[(2*(j-s)+p+s)&n2];
                s++;
                if (s==p)
                {
                    g=j+1;
                    s=0;
                }
            }
        }
        else
        {
            g=n;
            for (int j=0; j<n; j++)
            {
                dftn[j]=dfts[(2*(j-s)+s)&n2]+omega[g]*dfts[(2*(j-s)+p+s)&n2];
                s++;
                if (s==p)
                {
                    g=n-j-1;
                    s=0;
                }
            }
        }
    }
}

void licz_omegi(int n1)
{
    double kat=2.0*PI/n1;
    int n2=n1-1;
    int dod=3*n1/4;

    for (int i=0; i<=n1; i++)
        cosi[i]=cos(kat*i);
    for (int i=0; i<=n1; i++)
        omega[i]=complex <double> (cosi[i], cosi[(i+dod)&n2]);
}

vector<T> brute_mul(const vector<T> & a, const vector<T> & b) {
	if(a.empty() || b.empty()) return vector<T>{};
	vector<T> r(a.size() + b.size() - 1, 0);
	for(int i = 0; i < (int) a.size(); ++i)
		for(int j = 0; j < (int) b.size(); ++j)
			r[i+j] += a[i] * b[j];
	for(auto & x : r) x %= mod;
	return r;
}

vector <long long> fft(vector <long long> jed, vector <long long> dwa)
{
    int n1=potenga(jed.size()+dwa.size());
if(n1 < 500) return brute_mul(jed, dwa);
    licz_omegi(n1);

    for (int i=0; i<(int)jed.size(); i++)
        dftn[i]=jed[i];
    for (int i=jed.size(); i<n1; i++)
        dftn[i]=0;
    dft(n1, 1);
    for (int i=0; i<n1; i++)
        pom[i]=dftn[i];

    for (int i=0; i<(int)dwa.size(); i++)
        dftn[i]=dwa[i];
    for (int i=dwa.size(); i<n1; i++)
        dftn[i]=0;
    dft(n1, 1);
    for (int i=0; i<n1; i++)
        dftn[i]*=pom[i];

    dft(n1, 0);

    vector <long long> ret;

    for (int i=0; i<n1; i++)
        ret.push_back(llround(dftn[i].real()/n1));
    for(auto & x : ret) x %= mod;
    return ret;
}

vector <long long> fft_dokladne(vector <long long> jed, vector <long long> dwa)
{
    int n1=potenga(jed.size()+dwa.size());

    licz_omegi(n1);

    long long M=32000;

    for (int i=0; i<(int)jed.size(); i++)
        dftn[i]=jed[i]/M;
    for (int i=jed.size(); i<n1; i++)
        dftn[i]=0;
    dft(n1, 1);
    for (int i=0; i<n1; i++)
        a1[i]=dftn[i];

    for (int i=0; i<(int)jed.size(); i++)
        dftn[i]=jed[i]%M;
    for (int i=jed.size(); i<n1; i++)
        dftn[i]=0;
    dft(n1, 1);
    for (int i=0; i<n1; i++)
        b1[i]=dftn[i];

    for (int i=0; i<(int)dwa.size(); i++)
        dftn[i]=dwa[i]/M;
    for (int i=dwa.size(); i<n1; i++)
        dftn[i]=0;
    dft(n1, 1);
    for (int i=0; i<n1; i++)
        a2[i]=dftn[i];

    for (int i=0; i<(int)dwa.size(); i++)
        dftn[i]=dwa[i]%M;
    for (int i=dwa.size(); i<n1; i++)
        dftn[i]=0;
    dft(n1, 1);
    for (int i=0; i<n1; i++)
        b2[i]=dftn[i];

    vector <long long> ret;
    for (int i=0; i<n1; i++)
        ret.push_back(0);

    for (int i=0; i<n1; i++)
        dftn[i]=a1[i]*a2[i];
    dft(n1, 0);
    for (int i=0; i<n1; i++)
        ret[i]+=llround(dftn[i].real()/n1)*M*M;

    for (int i=0; i<n1; i++)
        dftn[i]=a1[i]*b2[i]+b1[i]*a2[i];
    dft(n1, 0);
    for (int i=0; i<n1; i++)
        ret[i]+=llround(dftn[i].real()/n1)*M;

    for (int i=0; i<n1; i++)
        dftn[i]=b1[i]*b2[i];
    dft(n1, 0);
    for (int i=0; i<n1; i++)
        ret[i]+=llround(dftn[i].real()/n1);
	for(auto & x : ret) x %= mod;
    return ret;
}






const int nax = 5e5 + 5;
int t[nax];



vector<T> solve(vector<int> freq) {
	if(freq.size() == 1) return vector<T>(freq[0] + 1, 1);
	vector<int> a, b;
	for(int i = 0; i < (int) freq.size(); ++i) {
		if(i < (int) freq.size() / 2) a.push_back(freq[i]);
		else b.push_back(freq[i]);
	}
	return fft(solve(a), solve(b));
}

void te() {
	int n, k;
	scanf("%d%d", &n, &k);
	for(int i = 0; i < n; ++i)
		scanf("%d", &t[i]);
	sort(t, t + n);
	vector<int> freq;
	for(int i = 0; i < n; ++i) {
		int cnt = 1;
		while(i + 1 < n && t[i] == t[i+1]) {
			++i;
			++cnt;
		}
		freq.push_back(cnt);
	}
	printf("%d\n", (int) solve(freq)[k]);
}

int main() {
	te();
}