#include<cstdio>
#include<iostream>
using namespace std;
const int o=4e6+10,MOD=998244353;
int T,fac[o],inv[o],a,b,c,d,ans;
inline int C(int x,int y){if(x<y||y<0) return 0;return fac[x]*1ll*inv[y]%MOD*inv[x-y]%MOD;}
int main(){
	inv[1]=1;
	for(int i=2;i<o;++i) inv[i]=MOD-MOD/i*1ll*inv[MOD%i]%MOD;
	for(int i=fac[0]=inv[0]=1;i<o;++i) fac[i]=fac[i-1]*1ll*i%MOD,inv[i]=inv[i-1]*1ll*inv[i]%MOD;
	for(scanf("%d",&T);T--;printf("%d\n",ans),ans=0){
		scanf("%d%d%d%d",&a,&b,&c,&d);
		if(!a&&!b){ans=C(c+d,c);continue;}
		if(!c&&!d){ans=C(a+b,a);continue;}
		if(!a&&!d) continue;
		if(!b&&!c) continue;
		if(a<d) swap(a,d),swap(b,c);
		if(!d){ans=C(a+1-(c-1),c)*1ll*C(a+1-c+b-1,b)%MOD;continue;}
		if(b+c==1){ans=2;continue;}
		if(a==d){ans=2*!b;continue;}
		for(int i=0;i<2;++i){
			b-=i;c-=!i;
			if(b>=0&&c>=0) ans=(ans+C(a-d-1+i-(c-1),c)*2ll*C(a-d+i-c+b-1,b))%MOD;
			b+=i;c+=!i;
		}
		for(int i=0;i<2;++i) for(int j=0;j<2;++j){
			b-=i+j;c-=(!i)+(!j);
			if(b>=0&&c>=0) ans=(ans+C(a-d-3+i+j-(c-1),c)*1ll*C(a-d-1+i+j-c+b-1,b)%MOD*(d-1))%MOD;
			b+=i+j;c+=(!i)+(!j);
		}
	}
	return 0;
}