#include<bits/stdc++.h>
using namespace std;
const int mod=998244353,ggg=3,invggg=(mod+1)/ggg;
const int maxn=524288;
int mypow(int x,int y){
	int ret=1;
	while(y){
		if(y&1)ret=1ll*ret*x%mod;
		y>>=1;x=1ll*x*x%mod;
	}return ret;
}
int revs[maxn];
void initrev(int Len){
	for(int i=1;i<Len;++i){
		revs[i]=revs[i>>1]>>1;
		if(i&1)revs[i]|=(Len>>1);
	}
}
void DFT(int *f,int n,int rev){
	int g0=rev==1?ggg:invggg;
	initrev(n);
	for(register int i=0;i<n;++i)
		if(i<revs[i])
			f[i]^=f[revs[i]]^=f[i]^=f[revs[i]];
	for(register int h=2;h<=n;h<<=1){
		int gn=mypow(g0,(mod-1)/h);
		for(register int j=0;j<n;j+=h){
			int gk=1;
			for(register int k=j;k<j+(h>>1);++k){
				int e=f[k],o=1ll*gk*f[k+(h>>1)]%mod;
				f[k]=(e+o)%mod,f[k+(h>>1)]=(e+mod-o)%mod;
				gk=1ll*gk*gn%mod;
			}
		}
	}
	if(rev==-1){
		int invv=mypow(n,mod-2);
		for(register int i=0;i<n;++i)
			f[i]=1ll*f[i]*invv%mod; 
	}
}
int fac[maxn],inv[maxn];
int C(int x,int y){
	if(x<y||y<0)return 0;  
	return 1ll*fac[x]*inv[y]%mod*inv[x-y]%mod;
}
struct node{
	int p,l,r,w;
	void init(){p=0;l=0;r=0;w=1;}
	int calc(int qp,int ql,int qr){
		while(p<qp)w=(2ll*w+C(p,l-1)-C(p,r)+mod)%mod,p++;
		while(l<ql)w=(w-C(p,l++)+mod)%mod;
		while(l>ql)w=(w+C(p,--l))%mod;
		while(r>qr)w=(w-C(p,r--)+mod)%mod;
		while(r<qr)w=(w+C(p,++r))%mod;
		return w;
	}
}w[120003];
const int B=1000;
int dp[2][B+3],d;
void calc(int *f,int n,int p){
	if(n<=B){
		memset(dp,0,sizeof(dp));
		for(int i=0;i<n;i++)
			dp[0][i]=1;
		f[0]=1;
		for(int i=1;i<=d;i++){
			for(int j=0;j<n;j++){
				dp[i&1][j]=(j?dp[i-1&1][j-1]:0)+(j<n-1?dp[i-1&1][j+1]:0);
				if(dp[i&1][j]>=mod)dp[i&1][j]-=mod;
			}f[i]=dp[i&1][p-1];
		}return;
	}else{
		int l=(p-d-n)/(n+1),r=(d+p-1)/(n+1);
		for(int i=l;i<=r;i++)
			w[i-l].init();
		for(int i=0;i<=d;i++)
			for(int j=l;j<=r;j++)
				f[i]=(f[i]+(j%2?mod-1ll:1ll)*w[j-l].calc(i,(i-p+2+j*(n+1))>>1,(i-p+n)+j*(n+1)>>1))%mod;
	}
}
int n,m,k,a,b,c;
int f[maxn],g[maxn],h[maxn];
int main(){
	cin>>d>>n>>m>>k>>a>>b>>c;
	fac[0]=1;
	for(int i=1;i<maxn;i++)
		fac[i]=1ll*fac[i-1]*i%mod;
	for(int i=0;i<maxn;i++)
		inv[i]=mypow(fac[i],mod-2);
	calc(f,n,a);
	calc(g,m,b);
	calc(h,k,c);
	for(int i=0;i<=d;i++)
		f[i]=1ll*f[i]*inv[i]%mod,
		g[i]=1ll*g[i]*inv[i]%mod,
		h[i]=1ll*h[i]*inv[i]%mod;
	DFT(f,maxn,1);
	DFT(g,maxn,1);
	DFT(h,maxn,1);
	for(int i=0;i<maxn;i++)
		f[i]=1ll*f[i]*g[i]%mod*h[i]%mod;
	DFT(f,maxn,-1);
	cout<<1ll*f[d]*fac[d]%mod;
}