a,b=map(int,input().split())
n,k=map(int,input().split())
m=int(1e9)+7
f=[1]*(2*n+1)
for i in range(2*n):
    f[i+1]=(f[i]*(i+1))%m
def c(i,j):
    return f[i]*(pow(f[j]*f[i-j],m-2,m))
def d(i,j):
    return (c(i-1,j-1)*a+c(i-1,j-2)*b)%m
print(d(n,k))
print((a*d(2*n-1,n)+b*d(2*n-1,n+1))%m)