#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
int u,v,w,x,n;
ll sum,pre[2001],inv[2001],f[501][501],g[501][501],val[501],ans[501],dp[501];
const int mod=1e9+7;
ll C(int n,int m){
	return n<m||m<0?0:pre[n]*inv[m]%mod*inv[n-m]%mod;
}
ll T(int n){
	return (C(n*2,n)-C(n*2,n-1)+mod)%mod;
}
ll qpow(ll a,ll b){
	ll s=1;
	for(;b;b>>=1,a=a*a%mod) if(b&1) s=s*a%mod;
	return s;
}
int main(){
	scanf("%d%d%d%d",&u,&v,&w,&x),n=u+v+w,pre[0]=inv[0]=f[0][0]=1;
	if(n%2||max({u,v,w})>n/2) return printf("0"),0;
	for(int i=1;i<=2000;i++) pre[i]=pre[i-1]*i%mod,inv[i]=qpow(pre[i],mod-2);
	for(int i=1;i<=n/2;i++) f[i][1]=T(i-1);
	for(int k=2;k<=n/2;k++) for(int i=1;i<=n/2;i++) for(int j=1;i+j<=n/2;j++) (f[i+j][k]+=f[i][k-1]*f[j][1])%=mod;
	for(int i=0;i<=n/2;i++) for(int j=0;j<=n/2;j++) for(int k=i;k>=0;k--) g[i][j]=(g[i][j]*j+f[i][k]*inv[k])%mod;
	for(int i=0;i<=u&&i<=n/2;i++){
		for(int j=0;j<=v&&i+j<=n/2;j++){
			ll l=n/2-i-j,s=0;
			for(int k=0;k<=u-i&&k<=j;k++) (s+=C(i,w-l-j+k)*C(j,k)%mod*C(l,u-i-k))%=mod;
			for(int k=0;k<=n/2;k++) (val[k]+=g[i][k]*g[j][k]%mod*g[l][k]%mod*s)%=mod;
		}
	}
	for(int i=0;i<=n/2;i++){
		memset(dp,0,sizeof(dp)),dp[0]=1;
		for(int j=0;j<=n/2;j++){
			if(i==j) continue;
			ll x=qpow(i-j+mod,mod-2),y=mod-j*x%mod;
			for(int l=n/2;l>=0;l--) dp[l]=(dp[l]*y+(l>0?dp[l-1]:0)*x)%mod;
		}
		for(int j=0;j<=n/2;j++) (ans[j]+=dp[j]*val[i])%=mod;
	}
	for(int i=0;i<=n/2;i++) (sum+=pre[i]*ans[i]%mod*qpow(i,x))%=mod;
	printf("%lld",sum);
	return 0;
}