10 2 5
1 1 1 2 2 2 2 2 2 2
 10 2 5
1 1 1 1 1 1 2 2 2 2
 3 3 2
1 3 1
 8 8 7
7 2 8 6 8 5 3 8
 7 7 5
2 6 3 6 3 5 1
 5 5 1
5 2 5 1 2
 6 6 6
5 6 5 3 4 2
 5 5 2
2 3 2 2 4
 3 3 3
1 1 1
 111 78 54
24 13 49 22 52 38 38 63 21 10 9 56 60 47 11 78 23 52 55 12 72 33 17 14 33 59 21 32 20 73 1 60 27 64 23 74 36 49 69 14 29 37 19 61 52 33 26 27 10 11 34 25 62 47 36 12 63 45 31 71 78 1 6 26 36 65 56 46 19 28 77 64 56 8 77 38 3 33 63 77 40 7 22 13 60 14 57 43 51 10 30 21 7 67 32 11 12 9 72 24 27 67 25 63 6 67 69 37 26 19 26
 10 2 5
1 1 1 1 1 1 1 1 1 2
 6 6 1
4 5 6 3 5 3
 10 2 5
1 1 1 1 2 2 2 2 2 2
 8 8 3
8 8 2 7 8 2 3 5
 7 7 4
4 6 6 7 2 5 3
 10 2 5
2 2 2 2 2 2 2 2 2 2
 111 67 97
32 67 4 15 17 17 23 6 39 54 19 62 28 59 63 36 44 40 28 67 28 3 55 19 11 26 3 38 22 26 62 54 61 15 22 4 47 53 27 36 31 5 16 59 51 62 51 46 7 55 26 66 50 51 32 17 54 14 18 14 40 30 24 2 61 16 16 27 55 58 35 7 26 3 26 4 13 59 14 32 56 54 5 46 39 30 45 24 37 12 67 38 45 32 41 21 46 26 41 29 15 28 65 5 19 53 13 31 58 18 57
 6 6 2
4 4 2 4 2 3
 10 2 5
1 1 1 1 1 2 2 2 2 2
 10 2 5
1 1 1 1 1 1 1 1 2 2
 6 6 4
5 2 4 3 1 5
 10 2 5
1 1 2 2 2 2 2 2 2 2
 1 1 1
1
 8 8 5
7 5 6 8 2 8 3 4
 10 2 5
1 1 1 1 1 1 1 2 2 2
 111 101 52
54 69 33 81 8 90 48 21 96 35 44 75 68 34 101 82 93 94 86 76 33 91 5 70 22 84 70 71 62 81 50 83 51 43 95 46 19 53 101 34 51 40 6 17 69 60 73 89 65 75 59 68 62 49 89 23 63 83 86 68 62 39 1 18 3 75 44 57 9 24 67 81 83 45 93 66 79 71 43 49 17 2 33 38 98 71 71 21 56 16 31 2 75 61 64 38 69 80 81 12 68 95 26 9 69 82 41 35 75 7 71
 3 3 3
3 3 1
 4 4 2
2 4 4 3
 10 2 5
1 2 2 2 2 2 2 2 2 2
 4 4 4
4 4 4 4
 3 3 2
1 2 3
 2 1 2
1 1
 5 5 3
3 4 1 5 5
 7 7 1
2 6 5 6 2 3 3
 111 9 57
9 9 2 6 8 6 5 8 9 7 3 9 1 2 6 6 5 3 1 1 6 1 8 4 3 9 8 6 8 1 9 8 5 8 6 9 4 5 7 3 1 6 4 3 2 5 7 3 1 7 9 1 7 1 3 5 9 5 6 2 2 7 9 7 1 5 5 5 5 2 6 5 7 3 8 9 7 3 6 5 1 9 2 9 9 8 8 7 9 1 6 2 4 5 7 4 8 9 9 2 7 7 6 9 8 2 1 3 3 1 5
 4 4 1
2 4 1 4
 2 1 1
1 1
 111 16 58
15 6 4 7 7 10 16 2 3 1 4 8 10 1 2 11 6 6 14 16 6 16 1 2 13 2 7 13 1 4 1 4 12 2 4 11 15 12 9 13 10 4 11 7 5 3 12 11 11 3 16 15 9 6 15 11 4 14 8 15 16 3 9 10 9 1 8 4 14 7 13 5 1 8 5 11 8 14 1 8 15 2 5 2 4 5 9 2 10 11 11 5 7 11 11 14 5 16 14 2 4 11 12 10 7 14 2 5 3 2 1
 2 2 2
1 1

4
 5
 2
 6
 12
 3
 1
 4
 1
 648
 2
 4
 5
 19
 25
 1
 214
 5
 6
 3
 11
 3
 1
 41
 4
 881
 1
 4
 2
 1
 3
 1
 7
 4
 446
 3
 1
 552
 1

#include <bits/stdc++.h>
using namespace std;
const long long mod = 1009;
const int MAX = 263005;
template <typename T>
class cmplx {
 private:
  T x, y;

 public:
  cmplx() : x(0.0), y(0.0) {}
  cmplx(T a) : x(a), y(0.0) {}
  cmplx(T a, T b) : x(a), y(b) {}
  T get_real() { return this->x; }
  T get_img() { return this->y; }
  cmplx conj() { return cmplx(this->x, -(this->y)); }
  cmplx operator=(const cmplx& a) {
    this->x = a.x;
    this->y = a.y;
    return *this;
  }
  cmplx operator+(const cmplx& b) {
    return cmplx(this->x + b.x, this->y + b.y);
  }
  cmplx operator-(const cmplx& b) {
    return cmplx(this->x - b.x, this->y - b.y);
  }
  cmplx operator*(const T& num) { return cmplx(this->x * num, this->y * num); }
  cmplx operator/(const T& num) { return cmplx(this->x / num, this->y / num); }
  cmplx operator*(const cmplx& b) {
    return cmplx(this->x * b.x - this->y * b.y, this->y * b.x + this->x * b.y);
  }
  cmplx operator/(const cmplx& b) {
    cmplx temp(b.x, -b.y);
    cmplx n = (*this) * temp;
    return n / (b.x * b.x + b.y * b.y);
  }
};
template <typename T, typename S>
struct FFT {
  S PI;
  int n, L, *rev, last;
  cmplx<S> ONE, ZERO;
  cmplx<S>* omega_powers;
  FFT() {
    PI = acos(-1.0);
    ONE = cmplx<S>(1.0, 0.0);
    ZERO = cmplx<S>(0.0, 0.0);
    last = -1;
    int req = 1 << (32 - __builtin_clz(MAX));
    rev = new int[req];
    omega_powers = new cmplx<S>[req];
  }
  ~FFT() {
    delete rev;
    delete omega_powers;
  }
  void ReverseBits() {
    if (n != last) {
      for (int i = 1, j = 0; i < n; ++i) {
        int bit = n >> 1;
        for (; j >= bit; bit >>= 1) j -= bit;
        j += bit;
        rev[i] = j;
      }
    }
  }
  void DFT(vector<cmplx<S>>& A, bool inverse = false) {
    for (int i = 0; i < n; ++i)
      if (rev[i] > i) swap(A[i], A[rev[i]]);
    for (int s = 1; s <= L; s++) {
      int m = 1 << s, half = m / 2;
      cmplx<S> wm(cos(2.0 * PI / m), sin(2.0 * PI / m));
      if (inverse) wm = ONE / wm;
      omega_powers[0] = ONE;
      for (int k = 1; k < half; ++k) {
        omega_powers[k] = omega_powers[k - 1] * wm;
      }
      for (int k = 0; k < n; k += m) {
        for (int j = 0; j < half; j++) {
          cmplx<S> t = omega_powers[j] * A[k + j + half];
          cmplx<S> u = A[k + j];
          A[k + j] = u + t;
          A[k + j + half] = u - t;
        }
      }
    }
    if (inverse) {
      for (int i = 0; i < n; i++) A[i] = A[i] / n;
    }
  }
  vector<T> Multiply(const vector<T>& a, const vector<T>& b) {
    L = 0;
    int sa = a.size(), sb = b.size(), sc = sa + sb - 1;
    while ((1 << L) < sc) L++;
    n = 1 << L;
    ReverseBits();
    last = n;
    vector<cmplx<S>> aa(n, ZERO), bb(n, ZERO), cc;
    for (int i = 0; i < sa; ++i) aa[i] = cmplx<S>(a[i], 0);
    for (int i = 0; i < sb; ++i) bb[i] = cmplx<S>(b[i], 0);
    DFT(aa, false);
    DFT(bb, false);
    for (int i = 0; i < n; ++i) cc.push_back(aa[i] * bb[i]);
    DFT(cc, true);
    vector<T> ans(sc);
    for (int i = 0; i < sc; ++i) {
      ans[i] = cc[i].get_real() + 0.5;
      ans[i] %= mod;
    }
    return ans;
  }
  vector<T> Multiply(const vector<T>& a) {
    L = 0;
    int sa = a.size(), sc = sa + sa - 1;
    while ((1 << L) < sc) L++;
    n = 1 << L;
    ReverseBits();
    last = n;
    vector<cmplx<S>> aa(n, ZERO), cc;
    for (int i = 0; i < sa; ++i) aa[i] = cmplx<S>(a[i], 0);
    DFT(aa, false);
    for (int i = 0; i < n; ++i) cc.push_back(aa[i] * aa[i]);
    DFT(cc, true);
    vector<T> ans(sc);
    for (int i = 0; i < sc; ++i) ans[i] = cc[i].get_real() + 0.5;
    return ans;
  }
};
FFT<long long, long double> fft;
int n, m, k;
const int N = 2 * 1e5;
vector<int> f;
vector<long long> cal(int st, int en) {
  if (st == en) {
    int tk = min(k, f[st]);
    vector<long long> vec(tk + 1, 1);
    return vec;
  }
  int mid = (st + en) / 2;
  vector<long long> v1 = cal(st, mid), v2 = cal(mid + 1, en);
  v1 = fft.Multiply(v1, v2);
  if (v1.size() > k + 1) v1.resize(k + 1);
  return v1;
}
int main() {
  ios_base::sync_with_stdio(false);
  cin >> n >> m >> k;
  f.resize(m, 0);
  int x;
  for (int i = 0; i < n; ++i) {
    cin >> x;
    --x;
    ++f[x];
  }
  sort(f.begin(), f.end(), greater<int>());
  while (!f.empty() && f.back() == 0) f.pop_back();
  if (f.empty()) {
    cout << 0;
    return 0;
  }
  m = f.size();
  vector<long long> res = cal(0, m - 1);
  if (res.size() < k + 1) {
    cout << 0;
    return 0;
  }
  long long ans = res[k];
  if (ans < 0) ans += mod;
  cout << ans;
  return 0;
}
