#include <bits/stdc++.h>
using namespace std;
const int INF = 1 << 29;
const double EPS = 1e-9;
const int MOD = 100000007;
const int dx[] = {1, 0, -1, 0}, dy[] = {0, -1, 0, 1};
const int MAXN = 510;
int dp[2][MAXN][MAXN];
int a[MAXN];
int main() {
  cout.setf(ios::fixed, ios::floatfield);
  cout.precision(8);
  ios_base::sync_with_stdio(false);
  int n, m, b, mod;
  cin >> n >> m >> b >> mod;
  for (int i = 0; i < (n); i++) cin >> a[i];
  for (int i = 0; i < (2); i++)
    for (int j = 0; j < (MAXN); j++)
      for (int k = 0; k < (MAXN); k++) dp[i][j][k] = 0;
  dp[0][0][0] = 1;
  for (int orig_i = 1; orig_i < n + 1; orig_i++) {
    int i = orig_i & 1;
    for (int j = 0; j < (m + 1); j++) {
      for (int k = 0; k < (b + 1); k++) {
        dp[i][j][k] = dp[i ^ 1][j][k];
        if (0 <= j - 1 && 0 <= k - a[orig_i - 1]) {
          dp[i][j][k] += dp[i][j - 1][k - a[orig_i - 1]];
        }
        dp[i][j][k] = (dp[i][j][k] % mod);
      }
    }
  }
  int res = 0;
  for (int i = 0; i < (b + 1); i++) res = (res + dp[n % 2][m][i]) % mod;
  cout << (res % mod) << endl;
  return 0;
}
