from copy import deepcopy

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

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

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

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

    S=[sum(G)%Mod for G in ODP]

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

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