import sys
input = sys.stdin.readline

A,B,C,N=map(int,input().split())

def c2(x):
    return x*(x-1)//2

def plus(A,B):
    ANS=[0,0,0]
    for i in range(3):
        ANS[i]=A[i]+B[i]

    return ANS

def mul(k,A):
    ANS=[0,0,0]
    for i in range(3):
        ANS[i]=k*A[i]
    return ANS
        
from functools import lru_cache
@lru_cache(maxsize=None)
def calc(A,B,C,N):
    if N==0:
        return [0,0,0]
    ALL=c2(A+B+C)

    if ALL==0:
        return [0,0,0]

    ANS=[0,0,0]

    ANS=plus(ANS,mul(c2(A)/ALL,plus([1,0,0],calc(A-1,B,C,N-1))))
    ANS=plus(ANS,mul(c2(B)/ALL,plus([0,1,0],calc(A,B-1,C,N-1))))
    ANS=plus(ANS,mul(c2(C)/ALL,plus([0,0,1],calc(A,B,C-1,N-1))))
    ANS=plus(ANS,mul((ALL-c2(A)-c2(B)-c2(C))/ALL,calc(A,B,C,N-1)))

    return ANS
    
print(calc(A,B,C,N))