#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef pair<int,int> pr;
const int mod=998244353,N=2097152,g=3,invg=(mod+1)/g;
ll n,pw3[40]={1};
int X,Y,Z,pwY[105]={1},pwZ[105]={1},cnt3[50000005],ans=0;
int tr[2100005],wk[2100005],a[2100005],b[2100005],c[2100005];
int Power(int x,ll y){
	int r=1;
	while(y){
		if(y&1)r=1ll*r*x%mod;
		x=1ll*x*x%mod,y>>=1; 
	}
	return r;
}
void GetTr(int l){
	for(int i=0;i<l;i++)tr[i]=((tr[i>>1]>>1)|((i&1)?(l>>1):0));
}
void NTT(int a[],int n,int fl){
	for(int i=0;i<n;i++)if(i<tr[i])swap(a[i],a[tr[i]]);
	for(int i=1;i<n;i<<=1){
		int w=Power(fl==1?g:invg,(mod-1)/(i<<1));
		wk[0]=1;
		for(int j=1;j<i;j++)wk[j]=1ll*wk[j-1]*w%mod;
		for(int j=0;j<n;j+=(i<<1)){
			int *p=a+j,*q=a+i+j,*t=wk;
			for(int k=0;k<i;k++,++p,++q,++t){
				int x=*p,y=1ll*(*q)*(*t)%mod,z=x;
				x+=y,((x>=mod)&&(x-=mod)),*p=x;
				x=z-y,((x<0)&&(x+=mod)),*q=x;
			}
		}
	}
	if(fl==-1){
		int inv=Power(n,mod-2);
		for(int i=0;i<n;i++)a[i]=1ll*a[i]*inv%mod;
	}
}
void upd(int &x,int y){
	x=(x+y)%mod;
}
int main(){
	cin>>n>>X>>Y>>Z;
	for(int i=1;i<=100;i++)pwY[i]=1ll*pwY[i-1]*Y%mod,pwZ[i]=1ll*pwZ[i-1]*Z%mod;
	for(int i=1;i<=30;i++)pw3[i]=3*pw3[i-1];
	for(int i=0;i<pw3[16];i++)cnt3[i]=cnt3[i/3]+(i%3);
	int pwX12=Power(X,pw3[12]);
	for(int i=0,t=1;i<pw3[16];i++,t=1ll*t*pwX12%mod){
		if(i*pw3[12]>n)break;
		if((i+1)*pw3[12]-1>n){
			for(int j=0;j<pw3[12];j++){
				ll x=i*pw3[12]+j;
				if(x>n)break;
				upd(ans,1ll*pwZ[cnt3[i]+cnt3[j]]*pwY[__builtin_popcountll(x)]%mod*Power(X,x)%mod);
			}
			continue;
		}
		ll x=i*pw3[12];
		//不进位 
		upd(a[x&((N-1)>>1)],1ll*pwZ[cnt3[i]]*pwY[__builtin_popcountll(x>>20)]%mod*t%mod);
		//进位 
		upd(b[x&((N-1)>>1)],1ll*pwZ[cnt3[i]]*pwY[__builtin_popcountll((x>>20)+1)]%mod*t%mod);
	}
	for(int i=0,t=1;i<pw3[12];i++,t=1ll*t*X%mod)upd(c[i],1ll*pwZ[cnt3[i]]*t%mod);
	GetTr(N),NTT(a,N,1),NTT(b,N,1),NTT(c,N,1);
	for(int i=0;i<N;i++)a[i]=1ll*a[i]*c[i]%mod,b[i]=1ll*b[i]*c[i]%mod;
	NTT(a,N,-1),NTT(b,N,-1);
	for(int i=0;i<N/2;i++){
		upd(ans,1ll*a[i]*pwY[__builtin_popcountll(i)]%mod);
		upd(ans,1ll*b[i+N/2]*pwY[__builtin_popcountll(i)]%mod);
	}
	cout<<(ans-1+mod)%mod;
}