#include<bits/stdc++.h>
using namespace std;
const long long mod=998244353;
const long long g0=3,g1=332748118;
const int N=8010000;
long long Pow(long long x,long long y){
	long long sum=1;
	for(;y;y>>=1,x=x*x%mod)if(y&1)sum=sum*x%mod;
	return sum;
}
namespace Poly
{
	int Pos,L,Rev[N];
	inline void Init(int ll){
		for(Pos=0,L=1;L<ll;L<<=1,++Pos);
		for(register int i=0;i<L;i++)Rev[i]=(Rev[i>>1]>>1)|((i&1)<<(Pos-1));
		return;
	}
	inline void NTT(long long *C,int opt)
	{
		for(register int i=0;i<L;i++)if(i<Rev[i])swap(C[i],C[Rev[i]]);
		for(register int mid=1;(mid<<1)<=L;mid<<=1){
			long long Wn=Pow((opt==0)?g0:g1,(mod-1)/(mid<<1)),W=1;
				for(register int i=0;i<L;i+=(mid<<1)){
					W=1;
					for(register int j=0;j<mid;j++,W=W*Wn%mod){
						long long tmp=C[i+mid+j]*W%mod;
						C[i+mid+j]=(C[i+j]+mod-tmp)%mod;
						C[i+j]=(C[i+j]+tmp)%mod;
					}
				}
			}if(opt==1){
			long long IV=Pow(L,mod-2);
			for(register int i=0;i<L;i++)C[i]=C[i]*IV%mod;
		}return;
	}

	inline void Dot(long long *A,long long *B,long long *C,int ll){
		for(register int i=0;i<ll;i++)C[i]=A[i]*B[i]%mod;return;
	}
	void Inte(long long *A,long long *B,int ll){
		for(int i=ll-1;i>=1;i--)B[i]=A[i-1]*Pow(i,mod-2)%mod;
		B[0]=0;return;
	}
	void Deri(long long *A,long long *B,int ll){
		for(int i=0;i<ll-1;i++)B[i]=A[i+1]*(i+1)%mod;
		B[ll-1]=0;return;
	}
	void Inv(long long *A,long long *B,int ll){
		static long long Tmp[N];for(int i=0;i<(ll<<2);i++)Tmp[i]=0;
		B[0]=Pow(A[0],mod-2);
		for(int now=2;now<(ll<<1);now<<=1){
			Init(now<<1);
			for(int i=0;i<now;i++)Tmp[i]=A[i];for(int i=now;i<(now<<1);i++)Tmp[i]=0;
			NTT(Tmp,0);NTT(B,0);
			for(int i=0;i<(now<<1);i++)B[i]=B[i]*(mod+2-Tmp[i]*B[i]%mod)%mod;
			NTT(B,1);
			for(int i=now;i<(now<<1);i++)B[i]=0;
		}for(int i=ll;i<L;i++)B[i]=0;return;
	}
	void Ln(long long *A,long long *B,int ll){
		static long long Tmp1[N],Tmp2[N];for(int i=0;i<(ll<<2);i++)Tmp1[i]=Tmp2[i]=0;
		Deri(A,Tmp1,ll);Inv(A,Tmp2,ll);
		Init(ll<<1);NTT(Tmp1,0);NTT(Tmp2,0);Dot(Tmp1,Tmp2,Tmp1,L);NTT(Tmp1,1);
		Inte(Tmp1,B,ll);return;
	}
	void Exp(long long *A,long long *B,int ll){
		static long long Tmp1[N],Tmp2[N];for(int i=0;i<(ll<<2);i++)Tmp1[i]=Tmp2[i]=0;
		B[0]=1;
		for(int now=2;now<(ll<<1);now<<=1)
		{
			Ln(B,Tmp1,now);for(int i=now;i<(now<<1);i++)Tmp1[i]=0;
			for(int i=0;i<now;i++)Tmp2[i]=A[i];for(int i=now;i<(now<<1);i++)Tmp2[i]=0;
			Init(now<<1);
			NTT(Tmp1,0);NTT(Tmp2,0);NTT(B,0);
			for(int i=0;i<(now<<1);i++)B[i]=B[i]*(mod+1-Tmp1[i]+Tmp2[i])%mod;
			NTT(B,1);
			for(int i=now;i<(now<<1);i++)B[i]=0;
		}for(int i=ll;i<L;i++)B[i]=0;return;
	}
}
int n;
int A[4];
long long P1[N],P2[N],P3[N];
int main()
{
	for(int i=1;i<=3;i++)cin>>A[i];
	cin>>n;
	P1[0]=Pow(10,mod-2);P1[2]=P1[1]=P1[0];
	Poly::Ln(P1,P2,n*2+1);
	for(int i=0;i<=n*2;i++)P2[i]=P2[i]*n%mod;
	Poly::Exp(P2,P3,n*2+1);
	long long t=Pow(Pow(3,n),mod-2);
	for(int i=0;i<=n*2;i++)P3[i]=P3[i]*t%mod;
	long long res=0;
	for(int i=1;i<=2;i++)
		for(int j=i+1,L;j<=3;j++){
			L=abs(A[i]-A[j]);
			for(int k=0;k<=n*2;k++){
				res+=P3[k]*abs(L+n-k)%mod;
				res%=mod;
			}
		}
	res=res*Pow(2,mod-2)%mod;
	cout<<res;
	return 0;
}