from _collections import defaultdict
mod=998244353
def sol1(n,a,b,c):
    def op(x,y):
        return (a*(x&y)+b*(x|y)+c*(x^y))%4
    ans=0
    for bit in range(4**n):
        A=[]
        for i in range(n):
            A.append(bit%4)
            bit//=4
        res1=A[0]
        res2=A[n-1]
        for i in range(n-1):
            res1=op(res1,A[i+1])
            res2=op(res2,A[n-2-i])
        if res1==res2:
            ans+=1
    return ans%mod
def sol2(n,a,b,c):
    def op(x,y):
        return (a*(x&y)+b*(x|y)+c*(x^y))%4
    def mp(f,g):
        ff=[]
        gg=[]
        for i in range(4):
            ff.append(f%4)
            gg.append(g%4)
            f//=4
            g//=4
        res=[0]*4
        for i in range(4):
            res[i]=ff[gg[i]]
        return res[0]+4*res[1]+16*res[2]+64*res[3]
    def tof(f):
        ff=[]
        for i in range(4):
            ff.append(f%4)
            f//=4
        return ff
    func=[0]*4
    for i in range(4):
        fs=[0]*4
        for x in range(4):fs[x]=op(x,i)
        func[i]=fs[0]+4*fs[1]+16*fs[2]+64*fs[3]
    dp=[defaultdict(int) for i in range(n+1)]
    for i in range(4):
        dp[1][func[i],func[i]]+=1
    for i in range(1,n):
        for f,g in dp[i]:
            dp[i][f,g]%=mod
            for j in range(4):
                dp[i+1][mp(func[j],f),mp(g,func[j])]+=dp[i][f,g]
    ans=0
    for i in range(4):
        for j in range(4):
            for f,g in dp[n-2]:
                res1=tof(f)[i]
                res1=op(res1,j)
                res2=tof(g)[j]
                res2=op(res2,i)
                if res1==res2:
                    ans+=dp[n-2][f,g]
    return ans%mod


n,a,b,c=map(int,input().split())
print(sol2(n,a,b,c))
exit()
for n in range(6,7):
    for a in range(4):
        for b in range(4):
            for c in range(4):
                ans=sol2(n,a,b,c)
                if ans!= 4**n and ans!=(4**n)//2:
                    print(a,b,c)
                    print(ans,4**n)