from copy import deepcopy

N,K,X,Y=map(int,input().split())
A=list(map(int,input().split()))
Mod=998244353
T=1023

DP=[[0]*(T+1) for _ in range(K)]

for i in range(K):
    DP[i][A[i]]=1

for _ in range(N-1):
    ODP,DP=DP,[[0]*(T+1) for _ in range(K)]

    S=[0]*(T+1)
    for k in range(K):
        E=ODP[k]
        for x in range(T+1):
            S[x]+=E[x]
            S[x]%=Mod

    for i in range(K):
        E=DP[i]
        F=ODP[i]
        for x in range(T+1):
            E[x]=S[x^A[i]]-F[x^A[i]]

B=0
for L in DP:
    B+=sum(L[X:Y+1])
    B%=Mod
print(B)