#include <bits/stdc++.h>
using namespace std;
template <class T>
ostream &operator<<(ostream &os, const vector<T> &vec) {
  os << '{';
  for (int i = (0), __i = (((int)(vec).size())); i < __i; i++) {
    os << vec[i];
    if (i + 1 != ((int)(vec).size())) os << ',';
  }
  os << '}';
  return os;
}
template <class T1, class T2>
ostream &operator<<(ostream &os, const pair<T1, T2> &par) {
  os << '(' << par.first << ',' << par.second << ')';
  return os;
}
int gcd(int x, int y) { return y ? gcd(y, x % y) : abs(x); }
template <class T>
T sqr(T x) {
  return x * x;
}
int mod1(int x, int p) { return x >= p ? x - p : x; }
int a[550];
int d[2][550][550];
int main() {
  ios::sync_with_stdio(false);
  int n, m, b, p;
  cin >> n >> m >> b >> p;
  for (int i = (0), __i = (n); i < __i; i++) cin >> a[i];
  for (int i = (0), __i = (n + 1); i < __i; i++)
    for (int j = (0), __j = (m + 1); j < __j; j++)
      for (int k = (0), __k = (b + 1); k < __k; k++)
        d[i & 1][j][k] = mod1(
            i ? d[~i & 1][j][k] +
                    (j && k >= a[i - 1] ? d[i & 1][j - 1][k - a[i - 1]] : 0)
              : j == 0,
            p);
  cout << d[n & 1][m][b] << "\n";
  return 0;
}
