MOD = 10 ** 9 + 7

fact_range = 2 * 10**6 + 100
fact = [1] * (fact_range + 1)
for i in range(0, fact_range):
    fact[i+1] = fact[i] * (i + 1) % MOD

ifact = [1] * (fact_range + 1)
ifact[fact_range] = pow(fact[fact_range], MOD - 2, MOD)
for i in range(fact_range, 0, -1):
    ifact[i-1] = ifact[i] * i % MOD

def comb(n, k):
    if k == 0:
        return 1
    if k < 0 or n < k:
        return 0
    else:
        return fact[n] * ifact[n-k] % MOD * ifact[k] % MOD


g, t, x, y = map(int, input().split())
ans = 0
t += 1
cnt_r = t // 4 + (t % 4 >= 1)
cnt_u = t // 4 + (t % 4 >= 2)
cnt_l = t // 4 + (t % 4 >= 3)
cnt_d = t // 4
for k1 in range(g + 1):
    k2 = g - k1
    if (x + k1) % 2 or (y + k2) % 2 or k1 < abs(x) or k2 < abs(y):
        continue
    x_r = (x + k1) // 2
    x_l = (k1 - x) // 2
    y_u = (y + k2) // 2
    y_d = (k2 - y) // 2
    if k1 == 1:
        print(x_r, x_l, y_u, y_d)
    ans += comb(x_r + cnt_r - 1, x_r) * comb(x_l + cnt_l - 1, x_l) % MOD * comb(y_u + cnt_u - 1, y_u) % MOD * comb(y_d + cnt_d - 1, y_d)
    ans %= MOD
print(ans)