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
                wb = 2*i*j/ss
                bc = 2*j*k/ss
                cw = 2*i*k/ss
                white,black,brown = 0,0,0
                white += (wb+bc+cw)*dp[i][j][k][0]
                black += (wb+bc+cw)*dp[i][j][k][1]
                brown += (wb+bc+cw)*dp[i][j][k][2]
                white += ww*(1+dp[i-1][j][k][0])
                black += ww*dp[i-1][j][k][1]
                brown += ww*dp[i-1][j][k][2]
                white += bb*dp[i][j-1][k][0]
                black += bb*(1+dp[i][j-1][k][1])
                brown += bb*dp[i][j-1][k][2]
                white += cc*dp[i][j][k-1][0]
                black += cc*dp[i][j][k-1][1]
                brown += cc*(1+dp[i][j][k-1][2])
                new_dp[i][j][k] = (white,black,brown)
    dp = new_dp

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