P = 10 ** 9 + 7
a,b = map(int,input().split())
a %= P
b %= P
n,k = map(int,input().split())

C = n + 5
fact = [1] * C
fact_inv = [1] * C
for i in range(2,C):
    fact[i] = fact[i-1] * i % P
fact_inv[-1] = pow(fact[-1],P - 2,P)
for i in range(C - 2,1,-1):
    fact_inv[i] = fact_inv[i + 1] * (i + 1) % P
def calc(n,k):
    return fact[n] * fact_inv[k] % P * fact_inv[n - k] % P

if k == 1:
    ans = a
elif k == n + 1:
    ans = b
else:
    ans = calc(n - 1,k - 1) * a + calc(n - 1,n - k + 1) * b
print(ans % P)
ans = 0
for k in range(1,n + 1 + 1):
    if k == 1:
        tmp = a
    elif k == n + 1:
        tmp = b
    else:
        tmp = calc(n - 1,k - 1) * a + calc(n - 1,n - k + 1) * b
    ans += (tmp % P) ** 2
    ans %= P
print(ans)