#include <iostream>

 #include <cstdio>

 #include <cstring>

 #define LL long long

 #define P 1000000007LL

 #define N 100010

 using namespace std;

 int a,b,c,d;

 LL fac[N*];

 LL qpow(LL x,int n)

 {

     LL ans=;

     for(;n;n>>=,x=x*x%P)

         if(n&) ans=ans*x%P;

     return ans;

 }

 LL C(int n,int m)

 {

     if(m==) return 1LL;

     return fac[n]*qpow(fac[n-m],P-)%P*qpow(fac[m],P-)%P;

 } 

 int main()

 {

     fac[]=;

     for(int i=;i<*N;i++) fac[i] = fac[i-]*(LL)i%P;

     while(~scanf("%d%d%d%d",&a,&b,&c,&d))

     {

         if(b== && d==)

         {

             if(a!= && c!=) puts("");

             else if(a!= || c!=) puts("");

             else puts("");

             continue;

         }

         LL ans;

         if(b==d)

         {

             ans = C(c+b-,b-)*C(a+b,b)%P;

             ans += C(c+b,b)*C(a+b-,b-)%P;

             if(ans>=P) ans-=P;

             cout << ans << endl;

         }

         else if(b==d+ || d==b+)

         {

             b = max(b,d);

             ans = C(b+c-,b-)*C(b+a-,b-)%P;

             cout << ans << endl;

         }

         else puts("");

     }

     return ;

 }