A,B,C,N=map(int,input().split()) DP=[[[[0]*(C+2) for ___ in range(B+2)] for _ in range(A+2)] for _ in range(N+1)] DP[0][A][B][C]=1 for n in range(1,N+1): for a in range(A+1): for b in range(B+1): for c in range(C+1): K=a+b+c if K<=1: continue s=(K*(K-1)-(a*(a-1)+b*(b-1)+c*(c-1)))/(K*(K-1)) t=((a+1)*a)/((K+1)*K) u=((b+1)*b)/((K+1)*K) v=((c+1)*c)/((K+1)*K) DP[n][a][b][c]=s*DP[n-1][a][b][c]+t*DP[n-1][a+1][b][c]+u*DP[n-1][a][b+1][c]+v*DP[n-1][a][b][c+1] S=T=U=0 for a in range(A+1): for b in range(B+1): for c in range(C+1): S+=(A-a)*DP[N][a][b][c] T+=(B-b)*DP[N][a][b][c] U+=(C-c)*DP[N][a][b][c] print(S,T,U)