#include<bits/stdc++.h>
#define int long long
using namespace std;
#define MOD 998244353
int jc[10200010],T,n,m,l,r,A,x[200010],p,q;
int adeld,deld,H,adm,h0,ny,invj[10200010],mxs[200010],cnt;
int ksm(int p,int q=MOD-2){
    int ans=1;
    while(q){
        if(q&1){ans=ans*p%MOD;}
        p=p*p%MOD;
        q>>=1;
    }
    return ans;
}
map<int,int> mp,mp2;
signed main(){
    ios::sync_with_stdio(false);
    cin.tie();cout.tie();
    jc[0]=invj[0]=jc[1]=1;
    for(int i=2;i<=10200000;i++)
        jc[i]=(jc[i-1]*i)%MOD;
    invj[10200000]=ksm(jc[10200000]);
    for(int i=10199999;i>=1;i--){
        invj[i]=(invj[i+1]*(i+1))%MOD;
        assert(invj[i]*jc[i]%MOD==1);
    }
    cin>>T;
    while(T--){
        mp.clear();mp2.clear();
        cnt=0;H=-1;
        cin>>n>>m>>l>>r;
        A=jc[n+m];
        for(int i=1;i<=n;i++){
            cin>>x[i];deld=r-l+1;
            if(l<=x[i]&&x[i]<=r)
                mp2[x[i]]++;
            else mp[x[i]]++;
        }
        for(auto x:mp2)
            mxs[cnt++]=x.second;
        deld-=cnt;
        // cout<<cnt<<' '<<deld<<'\n';
        sort(mxs,mxs+cnt);
        mxs[cnt]=-1;mxs[cnt+1]=0; 
        for(h0=0;h0<cnt;){
        	adm=0;adeld=0;
            while(mxs[h0+1]==mxs[h0])
                adm+=mxs[h0],adeld++,h0++;
            adm+=mxs[h0],adeld++,h0++;
            if(adm+m>=mxs[h0-1]*(adeld+deld))H=h0,m+=adm,deld+=adeld;
        }
        for(h0=H;h0<cnt;h0++)
            if(mxs[h0])
				mp[l+h0-H]=mxs[h0];
//        if(deld>r-l+1)deld=r-l+1;
//        cout<<A<<' '<<m<<' '<<deld<<'\n';
        p=m/deld;
        q=m%deld;
        if(q&&(p+1)){
            // cout<<"A";
            // cout<<A<<'\n';
            A=(A*ksm(invj[p+1],q))%MOD;
        }
        if((deld-q)&&(p)){
            // cout<<"B";
            // cout<<A<<'\n';
            A=(A*ksm(invj[p],deld-q))%MOD;
        }
        // cout<<A<<'\n';
        for(auto x:mp)
            A=(A*invj[x.second])%MOD;
        cout<<A<<'\n';
    }
    return 0;
}