import sys
input = sys.stdin.readline
sys.setrecursionlimit(10**7)

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

dp = [[[(0,0,0)]*(C+1) for _ in range(B+1)] for _ in range(A+1)]
for _ in range(N):
    new_dp = [[[(0,0,0)]*(C+1) for _ in range(B+1)] for _ in range(A+1)]
    for i in range(1,A+1):
        for j in range(1,B+1):
            for k in range(1,C+1):
                s = i+j+k
                ss = s*(s-1)
                ww = i*(i-1)/ss
                bb = j*(j-1)/ss
                cc = k*(k-1)/ss
                q = 1-ww-bb-cc
                prob = [ww,bb,cc,q]
                white,black,brown = ww,bb,cc
                X = [(1,0,0),(0,1,0),(0,0,1),(0,0,0)]
                for l in range(4):
                    p = prob[l]
                    di,dj,dk = X[l]
                    white += p*dp[i-di][j-dj][k-dk][0]
                    black += p*dp[i-di][j-dj][k-dk][1]
                    brown += p*dp[i-di][j-dj][k-dk][2]
                new_dp[i][j][k] = (white,black,brown)
    dp = new_dp

w,b,c = dp[-1][-1][-1]
print(w,b,c)