#include<iostream>
#include<iomanip>
using namespace std;
int A,B,C,N;
double dp[2][51][51][51];
main()
{
	cin>>A>>B>>C>>N;
	int now=0;
	dp[now][A][B][C]=1;
	for(int _=0;_<N;_++)
	{
		int nxt=1-now;
		for(int i=1;i<=A;i++)for(int j=1;j<=B;j++)for(int k=1;k<=C;k++)dp[nxt][i][j][k]=0;
		for(int i=1;i<=A;i++)for(int j=1;j<=B;j++)for(int k=1;k<=C;k++)
		{
			int n=i+j+k;
			int all=n*(n-1)/2;
			int usd=all;
			{
				int a=i*(i-1)/2;
				dp[nxt][i-1][j][k]+=dp[now][i][j][k]*a/all;
				usd-=a;
			}
			{
				int b=j*(j-1)/2;
				dp[nxt][i][j-1][k]+=dp[now][i][j][k]*b/all;
				usd-=b;
			}
			{
				int c=k*(k-1)/2;
				dp[nxt][i][j][k-1]+=dp[now][i][j][k]*c/all;
				usd-=c;
			}
			dp[nxt][i][j][k]+=dp[now][i][j][k]*usd/all;
		}
		now=nxt;
	}
	double a=0,b=0,c=0;
	for(int i=1;i<=A;i++)for(int j=1;j<=B;j++)for(int k=1;k<=C;k++)
	{
		a+=dp[now][i][j][k]*(A-i);
		b+=dp[now][i][j][k]*(B-j);
		c+=dp[now][i][j][k]*(C-k);
	}
	cout<<fixed<<setprecision(16)<<a<<" "<<b<<" "<<c<<endl;
}