#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
#define int ll
#define inf 0x3f3f3f3f
#define mod 1000000007
#define pb push_back
ll qow(ll a,ll p) {ll ans=1;for(;p;a=a*a%mod,p>>=1) if(p&1) ans=ans*a%mod;return ans;}
ll inv(ll a) {return qow(a,mod-2);}
const int N=2000010;
int T,n,a,b,c;
int ans;
int fac[N],invfac[N];
int C(int x,int y)
{
	if(x<0||y<0||y>x) return 0;
	return fac[x]*invfac[y]%mod*invfac[x-y]%mod;
}
int cal(int x)
{
	int a1=a,b1=b,c1=c;
	int ret=C(n,x)*C(n-x,a1-2*x)%mod*qow(2,a1-2*x)%mod*inv(C(2*n,a1))%mod;//概率

	int bo=x;//两个均被占
	int si=a1-2*x;//单个座位
	int tot=a1;//贡献
	if(si>=b1){
		tot+=b1;
		bo+=b1,si-=b1;
		//剩下的是1
		if(c1<=n-si-bo) tot+=c1;
		else{
			c1-=(n-si-bo),tot+=n-si-bo;
			if(c1<=si) ;
			else{

				c1-=si;
				c1=min(c1,n-si-bo);
				tot-=c1;//原来c满意，现在不满意了
			}
		}
	}
	else{
		//先跟a匹配
		tot+=si,b1-=si,bo+=si,si=0;
		//然后自己匹配
		tot+=b1/2*2,bo+=b1/2,si=b1%2;
		if(c1<=n-si-bo) tot+=c1;
		else{
			c1-=n-si-bo,tot+=n-si-bo;
			if(c1<=si) tot+=c1;//c不满意，但是b满意了
			else{
				c1-=si;
				tot+=si;//b满意了
				c1=min(c1,n-si-bo);
				tot-=c1;
			}
		}
	}
	//cout<<"|||"<<x<<' '<<ret<<' '<<tot<<'\n';
	return ret*tot%mod;

}
signed main()
{
	fac[0]=invfac[0]=1;
	for(int i=1;i<N;i++) fac[i]=fac[i-1]*i%mod,invfac[i]=inv(fac[i]);
	cin>>n>>a>>c>>b;
	if(a+b>=2*n){
		cout<<2*n<<'\n';
		return 0;
	}
	for(int i=0;i<=a/2;i++){//枚举有几个成对的
		if(i+a-2*i>n) continue;
		ans=(ans+cal(i))%mod;
	}

	int tot=0;
	for(int i=0;i<=a/2;i++){
		if(i+a-2*i>n) continue;
		tot=(tot+C(n,i)*C(n-i,a-2*i)%mod*qow(2,a-2*i)%mod*inv(C(2*n,a))%mod)%mod;
	}
	assert(tot==1);
	cout<<ans<<'\n';
}
/*
11 10 15 6

10 5 13 11

*/