# 解説を見た。
# 形式的ベキ級数を使うところまではいっていたが、式変形ができなかった。

import sys
input = sys.stdin.readline

N,M,K,L=list(map(int,input().split()))

mod=998244353

FACT=[1]
for i in range(1,10**7+1):
    FACT.append(FACT[-1]*i%mod)

FACT_INV=[pow(FACT[-1],mod-2,mod)]
for i in range(10**7,0,-1):
    FACT_INV.append(FACT_INV[-1]*i%mod)

FACT_INV.reverse()

def Combi(a,b):
    if 0<=b<=a:
        return FACT[a]*FACT_INV[b]%mod*FACT_INV[a-b]%mod
    else:
        return 0

ANS=0
for x in range(L,M+1):
    rest=M*K-x

    score=0
    for i in range(N):
        if i%2==0:
            score=(score+Combi(N-1,i)*(Combi(N-2+rest-(M+1)*i,N-2)-Combi(N-2+rest-(x+1)*i,N-2)))%mod
        else:
            score=(score-Combi(N-1,i)*(Combi(N-2+rest-(M+1)*i,N-2)-Combi(N-2+rest-(x+1)*i,N-2)))%mod

    #print(x,score,rest)

    ANS=(ANS+score)%mod

print(ANS)