#include <bits/stdc++.h>
using namespace std;
const int inf = 1e9 + 7;
struct node {
  int x;
  int ord;
  int r;
  int vis;
} a[1000005];
int cmp(node m, node n) { return m.x < n.x; }
int cmp1(node m, node n) { return m.ord < n.ord; }
int main() {
  long long n, k;
  long long l;
  cin >> n >> l >> k;
  vector<vector<long long> > dp(n + 3, vector<long long>(k + 2));
  long long i, j;
  for (i = 1; i <= n; i++) {
    scanf("%d", &a[i].x);
    a[i].ord = i;
  }
  long long p = l % n;
  if (l % n) {
    for (i = 1; i <= p; i++) a[i].vis = 1;
  }
  sort(a + 1, a + 1 + n, cmp);
  for (i = 1; i <= n; i++) {
    a[i].r = i;
  }
  for (j = 1; j <= n; j++) dp[j][1] = 1;
  for (i = 2; i <= k; i++) {
    int tot = 1;
    for (j = 1; j <= n; j++) {
      dp[j][i] = dp[j - 1][i];
      while (tot <= n && a[tot].x <= a[j].x) {
        dp[j][i] += dp[tot][i - 1];
        dp[j][i] %= inf;
        tot++;
      }
    }
  }
  long long ans = 0;
  long long len = l / n;
  for (i = 1; i <= n; i++) {
    for (j = 1; j <= k; j++) {
      if (a[i].vis == 1) {
        if (j <= len + 1) {
          ans += ((len + 1) + 1 - j) % inf * dp[i][j];
          ans %= inf;
        }
      } else {
        if (j <= len) {
          ans += (len + 1 - j) % inf * dp[i][j];
          ans %= inf;
        }
      }
    }
  }
  if (ans < 0) ans += inf;
  printf("%I64d\n", ans % inf);
}
