#include <iostream>
#include <algorithm>
#include <cstring>
#include <vector>
using ll=long long;
namespace polynomial
{
	constexpr int P(998244353),N(500010),PR(3),APR((P+1)/PR);
	int ws[2][N<<2],fac[N<<2],ifac[N<<2],fn,fb,lgg[N<<2],_inv[N<<2],rev[N<<2];
	int fast_pow(int a,int b){int ans(1),off(a);while(b){if(b&1) ans=1ll*ans*off%P;off=1ll*off*off%P;b>>=1;}return ans;}
	void init_poly_weights(int len)
	{
		fn=1,fb=0;while(fn<(len<<1)) fn<<=1,++fb;
		_inv[1]=1;for(int i(2);i<=fn;++i) _inv[i]=1ll*(P-P/i)*_inv[P%i]%P;
		ifac[0]=fac[0]=1;for(int i(1);i<=fn;++i) fac[i]=1ll*i*fac[i-1]%P,ifac[i]=1ll*ifac[i-1]*_inv[i]%P;
		for(int i(0);i<fn;++i) rev[i]=(rev[i>>1]>>1)|((i&1)<<(fb-1));
		for(int i(2);i<=fn;++i) lgg[i]=lgg[(i+1)>>1]+1;
		for(int mid((fn>>1)),j0(fast_pow(PR,(P-1)/(mid<<1))),j1(fast_pow(APR,(P-1)/(mid<<1)));mid>=1;mid>>=1,j0=1ll*j0*j0%P,j1=1ll*j1*j1%P)
		for(int i(0),w0(1),w1(1);i<mid;++i,w0=1ll*w0*j0%P,w1=1ll*w1*j1%P) ws[0][i+mid]=w0,ws[1][i+mid]=w1;
	}
	using poly=std::vector<int>;
	int redu(int k){return k>=P?k-P:k;}
	poly operator+(const poly &a,const poly &b)
	{
		poly ret(b);if(ret.size()<a.size()) ret.resize(a.size());
		for(int i(0);i<a.size();++i) ret[i]=redu(ret[i]+a[i]);
		return ret;
	}
	unsigned int tt[N<<2];
	unsigned int W(unsigned int k){return k>=P?k-P:k;}
	void NTT(poly &p,int V)
	{
		int bts(lgg[p.size()]);if(p.size()!=(1<<bts)) p.resize((1<<bts));
		int* w(ws[(V==1)?0:1]),len(1<<bts);for(int i(0);i<len;++i) tt[i]=p[rev[i]>>(fb-bts)];
		unsigned int t1,t2;
		for(int l(2);l<=len;l<<=1)
			for(int j(0),mid=(l>>1);j<len;j+=l)
				for(int i(0);i<mid;++i) t1=tt[j+i],t2=1ull*tt[j+i+mid]*w[mid+i]%P,tt[j+i]=W(t1+t2),tt[j+i+mid]=W(P+t1-t2);
		if(V==1) for(int i(0);i<len;++i) p[i]=tt[i];
		else for(int i(0),j(_inv[len]);i<len;++i) p[i]=1ull*tt[i]*j%P;
	}
	poly operator*(const poly &a,const poly &b)
	{
		static poly p1,p2;poly ret;
		p1=a,p2=b;int len(a.size()+b.size()-1),ff(1<<lgg[len]);
		p1.resize(ff),p2.resize(ff),ret.resize(ff);
		NTT(p1,1);NTT(p2,1);
		for(int i(0);i<ff;++i) ret[i]=1ll*p1[i]*p2[i]%P;
		NTT(ret,-1);ret.resize(len);return ret;
	}
	void print_poly(const poly &a){for(int v:a) printf("%d ",v);printf("\n");}
	poly poly_inv(const poly &a)
	{
		int l(a.size());if(l==1){poly ret(1);ret[0]=fast_pow(a[0],P-2);return ret;}
		poly g0(a);g0.resize((l+1)>>1);g0=poly_inv(g0);
		static poly p1;p1=a;int ff(2<<lgg[l]);poly ret(ff);g0.resize(ff);p1.resize(ff);
		NTT(p1,1);NTT(g0,1);for(int i(0);i<ff;++i) ret[i]=1ll*g0[i]*(2-1ll*g0[i]*p1[i]%P+P)%P;
		NTT(ret,-1);ret.resize(l);return ret;
	}
	poly poly_ln(const poly &a)
	{
		int l(a.size());static poly p1;
		p1.resize(l-1);for(int i(1);i<l;++i) p1[i-1]=1ll*i*a[i]%P;
		p1=p1*poly_inv(a);p1.resize(l-1);poly ret(l);
		for(int i(1);i<l;++i) ret[i]=1ll*_inv[i]*p1[i-1]%P;ret[0]=0;
		return ret;
	}
	poly poly_exp(const poly &a)
	{
		int l(a.size());if(l==1){poly ret(1);ret[0]=1;return ret;}
		poly g0(a);g0.resize((l+1)>>1);g0=poly_exp(g0);static poly g1,g2;g1=g0;g1.resize(l);g1=poly_ln(g1);g2=a;
		int ff(2<<lgg[l]);g0.resize(ff);g1.resize(ff);poly ret(ff);g2.resize(ff);
		NTT(g0,1);NTT(g1,1);NTT(g2,1);
		for(int i(0);i<ff;++i) ret[i]=1ll*g0[i]*(1-g1[i]+g2[i]+P)%P;
		NTT(ret,-1);ret.resize(l);return ret;
	}
	std::pair<int,int> cf[N];// u/(1+dx)
	std::pair<poly,poly> calc(int l,int r)
	{
		if(l==r)
		{
			auto [u,d]=cf[l];
			return {{u},{1,d}};
		}
		int mid((l+r)>>1);
		auto [ls,lt]=calc(l,mid);
		auto [rs,rt]=calc(mid+1,r);
		int T(1<<lgg[r-l+2]);ls.resize(T);lt.resize(T);rs.resize(T);rt.resize(T);
		poly s(T),t(T);
		NTT(ls,1);NTT(lt,1);NTT(rs,1);NTT(rt,1);
		for(int i(0);i<T;++i) s[i]=(1ll*ls[i]*rt[i]+1ll*rs[i]*lt[i])%P,t[i]=1ll*lt[i]*rt[i]%P;
		NTT(s,-1);NTT(t,-1);s.resize(r-l+1);t.resize(r-l+2);
		return {s,t};
	}
	poly solve_fs(const std::vector<std::pair<int,int>> &a,int m)
	{
		int n(a.size());
		std::copy(a.begin(),a.end(),cf+1);
		auto [s,t]=calc(1,n);
		t.resize(m+1);
		t=poly_inv(t);
		s=s*t;
		s.resize(m+1);
		return s;
	}
	const poly unit_poly={1};
}
using polynomial::poly;
using polynomial::operator*;
using polynomial::operator+;
using polynomial::init_poly_weights;
using polynomial::unit_poly;
using polynomial::print_poly;
using polynomial::poly_inv;
using polynomial::poly_exp;
using polynomial::poly_ln;
using polynomial::solve_fs;
constexpr int N(1500010),p(998244353);
int main()
{
	init_poly_weights(N/3);
	int n;std::cin>>n;
	static int ifac[N],fac[N],inv[N];inv[1]=1;
	for(int i(2);i<N;++i) inv[i]=1ll*inv[p%i]*(p-p/i)%p;
	ifac[0]=fac[0]=1;for(int i(1);i<N;++i) ifac[i]=1ll*ifac[i-1]*inv[i]%p,fac[i]=1ll*fac[i-1]*i%p;
	static int s[3][3],d[3],r[3];
	for(int i(0);i<3;++i) for(int j(0);j<3;++j) std::cin>>s[i][j],d[j]+=s[i][j];
	if((d[1]-d[0])%3 || (d[2]-d[0])%3){std::cout<<"0"<<std::endl;return 0;}
	int bs(*std::min_element(d,d+3));
	int mn(1e9),mp(0);
	for(int i(0);i<3;++i)
	{
		r[i]=(d[i]-bs)/3;
		for(int j(0);j<3;++j)
		{
			s[j][i]-=r[i];
			if(s[j][i]<0)
			{
				std::cout<<"0"<<std::endl;
				return 0;
			}
		}
	}
	for(int i(0);i<3;++i){if(s[i][0]<mn) mn=s[i][0],mp=i;}
	if(mp) std::swap(s[0],s[mp]);
	for(int i(0);i<3;++i)
	{
		int diff(s[i][0]+s[i][1]-s[(i+1)%3][2]-s[(i+2)%3][2]);
		if(diff){std::cout<<"0"<<std::endl;return 0;}
	}
	//a:j,b:k,c:i
	int dk(s[1][0]-s[2][2]),di(s[2][0]-s[0][2]+dk);
	poly f(s[0][0]+1),g(s[0][0]+1);
	for(int j(0);j<=s[0][0];++j) if(j+dk>=0 && j+di>=0)
	{
		f[j]=1ll*ifac[j]*ifac[j+dk]%p*ifac[j+di]%p;
	}
	int dbk(s[1][0]-s[0][0]-dk),dci(s[2][0]-s[0][0]-di);
	for(int j(0);j<=s[0][0];++j) if(j+dbk>=0 && j+dci>=0)
	{
		g[j]=1ll*ifac[j]*ifac[j+dbk]%p*ifac[j+dci]%p;
	}
	// std::cerr<<dk<<" "<<di<<" "<<dbk<<" "<<dci<<std::endl;
	f=f*g;f.resize(s[0][0]+1);
	int ans(0);
	for(int i(0);i<=s[0][0];++i)
	{
		int num(3*i+s[1][0]+s[2][0]-s[0][0]*2);
		int cur(f[i]);
		for(int j(0);j<3;++j) cur=1ll*cur*ifac[s[0][0]-i+r[j]]%p,num+=s[0][0]-i+r[j];
		// std::cerr<<num<<std::endl;
		ans=(ans+1ll*cur*fac[n])%p;
	}
	std::cout<<ans<<std::endl;
	return 0;
}