#pragma GCC target("avx2")
#pragma GCC optimize("Ofast")
#pragma GCC optimize("unroll-loops")
#include <bits/stdc++.h>
using namespace std;
const long long MOD = 998244353;
void dfs(vector<vector<long long>> &H, vector<int> &A, int r, int prv, int x){
  if (r == 0){
    H[prv][x]++;
  } else {
    int K = A.size();
    for (int j = 0; j < K; j++){
      if (A[j] != prv){
        dfs(H, A, r - 1, A[j], x ^ A[j]);
      }
    }
  }
}
vector<long long> fwt(vector<long long> A, int N, bool inv){
	for (int i = 0; i < N; i++){
		for (int j = 0; j < (1 << N); j++){
			if (!(j >> i & 1)){
				long long x = A[j];
				long long y = A[j | (1 << i)];
				A[j] = (x + y) % MOD;
				A[j | (1 << i)] = (x + MOD - y) % MOD;
				if (inv){
					A[j] *= 499122177;
					A[j] %= MOD;
					A[j | (1 << i)] *= 499122177;
					A[j | (1 << i)] %= MOD;
				}
			}
		}
	}
	return A;
}
vector<long long> xor_convolution(vector<long long> A, vector<long long> B, int N){
	A = fwt(A, N, false);
	B = fwt(B, N, false);
	vector<long long> C(1 << N);
	for (int i = 0; i < (1 << N); i++){
		C[i] = A[i] * B[i] % MOD;
	}
	C = fwt(C, N, true);
	return C;
}
int main(){
  int N, K, X, Y;
  cin >> N >> K >> X >> Y;
  Y = min(Y, 1023);
  vector<int> A(K);
  for (int i = 0; i < K; i++){
    cin >> A[i];
  }
  vector<vector<long long>> H(1024, vector<long long>(1024, 0));
  dfs(H, A, N / 2, -1, 0);
  vector<long long> sum(1024, 0);
  for (int i = 0; i < 1024; i++){
    for (int j = 0; j < 1024; j++){
      sum[j] += H[i][j];
      sum[j] %= MOD;
    }
  }
  long long ans = 0;
  for (int i = 0; i < 1024; i++){
    for (int j = 0; j < 1024; j++){
      if (X <= (i ^ j) && (i ^ j) <= Y){
        ans += sum[i] * sum[j];
        ans %= MOD;
      }
    }
  }
  for (int i = 0; i < 1024; i++){
    for (int j = 0; j < 1024; j++){
      for (int k = 0; k < 1024; k++){
        if (X <= (j ^ k) && (j ^ k) <= Y){
          ans += MOD - H[i][j] * H[i][k] % MOD;
          ans %= MOD;
        }
      }
    }
  }
  cout << ans << endl;
}