#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int N = 5000005, mod = 998244353;
int n, a, b, c;
int fac[N], inv[N];
int ans;
int qpow(int x, int y) {
int res = 1;
while (y) {
if (y & 1) res = 1ll * res * x % mod;
x = 1ll * x * x % mod;
y >>= 1;
}
return res;
}
void init(int maxn) {
fac[0] = 1;
for (int i = 1; i <= maxn; ++i) fac[i] = 1ll * fac[i - 1] * i % mod;
inv[maxn] = qpow(fac[maxn], mod - 2);
for (int i = maxn - 1; ~i; --i) inv[i] = 1ll * inv[i + 1] * (i + 1) % mod;
}
int C(int n, int m) {
if (n < 0 || m < 0 || n < m) return 0;
return 1ll * fac[n] * inv[m] % mod * inv[n - m] % mod;
}
void add(int &a, int b) {
a += b;
if (a >= mod) a -= mod;
if (a < 0) a += mod;
}
void mul(int &a, int b) {
a = 1ll * a * b % mod;
}
int main() {
scanf("%d%d%d%d", &n, &a, &b, &c);
init(5000000);
int x = 1, y = 1, z = 1, f;
if (n & 1) f = -1; else f = 1;
for (int i = 0; i <= n; ++i) {
int tmp = C(n, i);
mul(tmp, x), mul(tmp, y), mul(tmp, z);
mul(x, 2), add(x, -C(i, a));
mul(y, 2), add(y, -C(i, b));
mul(z, 2), add(z, -C(i, c));
add(ans, tmp * f);
f *= -1;
}
printf("%d", ans);
return 0;
}