mod = 998244353


def main():
    import sys
    input = sys.stdin.readline

    N, a, b, c = map(int, input().split())
    Z = [[0] * 4 for _ in range(4)]
    for x in range(4):
        for y in range(4):
            Z[x][y] = (a * (x & y) + b * (x | y) + c * (x ^ y)) % 4

    ans = 0
    dp = [[0] * 256 for _ in range(4)]
    for i in range(4):
        state = 0
        for j in range(4):
            r = Z[i][j]
            state += r << (2 * j)
        dp[i][state] += 1
    for t in range(N - 1):
        dp_new = [[0] * 256 for _ in range(4)]
        for l in range(4):
            for state in range(256):
                A = []
                s = state
                for _ in range(4):
                    A.append(s % 4)
                    s >>= 2
                for x in range(4):
                    l_new = Z[l][x]
                    state_new = 0
                    for j in range(4):
                        r = A[Z[j][x]]
                        state_new += r << (2 * j)
                    dp_new[l_new][state_new] = (dp_new[l_new][state_new] + dp[l][state]) % mod
                    if t == N - 2:
                        if l_new == A[x]:
                            ans = (ans + dp[l][state]) % mod
        dp = dp_new

    print(ans)


if __name__ == '__main__':
    main()