#include<bits/stdc++.h>
using namespace std;
#define int long long
const int maxn=4e6+10;
const int INF=0x3f3f3f3f;
const int mod=998244353;
const int G=3;
int n,sta[3][3],f[maxn],inv[maxn];
int a[6][7],y[maxn],rev[maxn],limit;
int l1[maxn],l2[maxn],l3[maxn];
int qpow(int x,int y){
    int val=1;
    while(y){
        if(y&1)val=val*x%mod;
        x=x*x%mod,y/=2;
    }
    return val;
}
int lcm(int x,int y){
    return x/__gcd(x,y)*y;
}
void init(){//高斯消元
    for(int i=0;i<6;i++){
        int id=0;
        for(int j=i;j<6;j++)if(a[j][i])id=j;        
        swap(a[i],a[id]);
        for(int j=0;j<6;j++){
            if(j==i)continue;
            if(a[j][i]){
                int v=lcm(a[j][i],a[i][i]);
                int val1=v/a[i][i],val2=v/a[j][i];
                for(int k=0;k<7;k++)a[j][k]=a[j][k]*val2,a[i][k]=a[i][k]*val1;
                for(int k=0;k<7;k++)a[j][k]=a[j][k]-a[i][k];
            }
        }
    }
    for(int i=0;i<6;i++){
        if(a[i][6]%a[i][i]){
            cout<<0;
            exit(0);
        }
        y[i]=a[i][6]/a[i][i];
    }
}
void NTT(int dp[],int flag){
    for(int i=0;i<limit;i++)if(i<rev[i])swap(dp[i],dp[rev[i]]);
    for(int mid=1;mid<limit;mid*=2){
        int w=qpow(G,(mod-1)/2/mid);
        if(flag==-1)w=qpow(w,mod-2);
        for(int i=0;i<limit;i=i+2*mid){
            for(int j=i,wn=1;j<i+mid;j++,wn=wn*w%mod){
                int x=dp[j],y=dp[j+mid]*wn%mod;
                dp[j]=(x+y)%mod,dp[j+mid]=(x-y+mod)%mod;
            }
        }
    }
    if(flag==-1){
        int invs=qpow(limit,mod-2);
        for(int i=0;i<limit;i++)dp[i]=dp[i]*invs%mod;
    }
}
void solve(){
    limit=1;
    int v=0,N=n;
    n=sta[0][0];
    while(limit<=2*n)limit*=2,v++;
    for(int i=1;i<limit;i++)rev[i]=((rev[i>>1]>>1)|((i&1)<<(v-1)));
    int v1=max(-y[0],-y[1]),v2=max(-y[2],-y[3]),v3=max(-y[4],-y[5]);
    v1=max(v1,(int)0),v2=max(v2,(int)0),v3=max(v3,(int)0);
    for(int i=v1;i<=n;i++)l1[i]=inv[i]*inv[i+y[0]]%mod*inv[i+y[1]]%mod;
    for(int i=v2;i<=n;i++)l2[i]=inv[i]*inv[i+y[2]]%mod*inv[i+y[3]]%mod;
    for(int i=v3;i<=n;i++)l3[i]=inv[i]*inv[i+y[4]]%mod*inv[i+y[5]]%mod;
    NTT(l1,1);NTT(l2,1);NTT(l3,1);
    for(int i=0;i<limit;i++)l1[i]=l1[i]*l2[i]%mod*l3[i]%mod;
    NTT(l1,-1);
    cout<<l1[n]*f[N]%mod;
}
signed main(){
    std::ios::sync_with_stdio(false);
    cin.tie(0),cout.tie(0);
    cin>>n;
    f[0]=1;
    for(int i=1;i<=n;i++)f[i]=f[i-1]*i%mod;
    inv[n]=qpow(f[n],mod-2);
    for(int i=n-1;i>=0;i--)inv[i]=inv[i+1]*(i+1)%mod;

    for(int i=0;i<3;i++){
        for(int j=0;j<3;j++)cin>>sta[i][j];
    }
    //start
    a[0][0]=a[0][2]=a[0][4]=1;
    a[0][6]=sta[0][1]-sta[0][0];
    
    a[1][1]=a[1][3]=a[1][5]=1;
    a[1][6]=sta[0][2]-sta[0][0];

    a[2][0]=a[2][5]=1,a[2][3]=a[2][4]=-1;
    a[2][6]=sta[1][1]-sta[1][0];

    a[3][1]=a[3][2]=1,a[3][3]=a[3][4]=-1;
    a[3][6]=sta[1][2]-sta[1][0];

    a[4][0]=a[4][3]=1,a[4][2]=a[4][5]=-1;
    a[4][6]=sta[2][1]-sta[2][0];

    a[5][1]=a[5][4]=1,a[5][2]=a[5][5]=-1;
    a[5][6]=sta[2][2]-sta[2][0];
    // for(int i=0;i<6;i++){
    //     for(int j=0;j<6;j++){
    //         cout<<a[i][j]<<" ";
    //     }
    //     cout<<endl;
    // }
    init();
    solve();
    return 0;
}