#include "bits/stdc++.h"
using namespace std;

#define DEBUG(x) cout<<#x<<": "<<x<<endl;
#define DEBUG_VEC(v) cout<<#v<<":";for(int i=0;i<v.size();i++) cout<<" "<<v[i]; cout<<endl

typedef long long ll;
#define vi vector<int>
#define vl vector<ll>
#define vii vector< vector<int> >
#define vll vector< vector<ll> >
#define vs vector<string>
#define pii pair<int,int>
#define pis pair<int,string>
#define psi pair<string,int>
#define pll pair<ll,ll>
const int inf = 1000000001;
const ll INF = 1e18 * 4;
#define MOD 1000000007
#define mod 1000000009
#define pi 3.14159265358979323846
#define Sp(p) cout<<setprecision(15)<<fixed<<p<<endl;
int dx[4] = { 1,0,-1,0 }, dy[4] = { 0,1,0,-1 };
int dx2[8] = { 1,1,0,-1,-1,-1,0,1 }, dy2[8] = { 0,1,1,1,0,-1,-1,-1 };

ll MAXM = 3010;
vl fact(MAXM);

ll mod_pow(ll x, ll p, ll M) {
	ll a = 1;
	while (p) {
		if (p % 2)
			a = a*x%M;
		x = x*x%M;
		p /= 2;
	}
	return a;
}

ll mod_inverse(ll a, ll m) {
	return mod_pow(a, m - 2, m);
}

ll nCr(ll n, ll r) {
	ll ret = fact[n];
	ret = (ret*mod_inverse(fact[r], MOD)) % MOD;
	ret = (ret*mod_inverse(fact[n - r], MOD)) % MOD;
	return ret;
}


int main() {
	int n, r, g, b;
	cin >> n >> r >> g >> b;

	fact[0] = 1;
	for (int i = 1; i <= n; i++) {
		fact[i] = i*fact[i - 1];
	}

	int t1, t2, t3;
	ll ans = 0;
	for (t1 = 0; t1 <= r; t1++) {
		for (t2 = 0; t2 <= min(r, g) && r - t1 - t2 >= 0; t2++) {
			for (t3 = 0; t3 <= min(g, b) && g - t2 - t3 >= 0 && b - r + t1 + t2 - t3 >= 0; t3++) {
				if (n + 1 - 2 * g - 2 * b - r - t1 + t3 > n + 1 || n + 1 - 2 * g - 2 * b - r - t1 + t3 < 0) {
					continue;
				}
				int z = n + 1 - 2 * g - 2 * b - r - t1 + t3;
				int x1 = r - t1 - t2, x2 = t2, x3 = t3;
				int y1 = t1, y2 = b - x1 - t3, y3 = g - t2 - t3;
				int m = z + x1 + x2 + x3 + y1 + y2 + y3;
				ans = (ans + (nCr(m, x1) * mod_pow(2, x1, MOD) % MOD * nCr(m - x1, x2) % MOD * mod_pow(2,x2,MOD) % MOD * nCr(m - x1 - x2, x3) % MOD * mod_pow(2,x3,MOD) % MOD * nCr(m - x1 - x2 - x3, y1) % MOD * nCr(m - x1 - x2 - x3 - y1, y2) % MOD*nCr(m - x1 - x2 - x3 - y1 - y2, y3) % MOD)) % MOD;
			}
		}
	}
	cout << ans << endl;
}