#include <iostream>
#include <cstring>
#include <cstdio>
#include <vector>
#include <cmath>
#include <map>
#define int long long
#define vit vector<int>
#define rep(i, a, b) for (int i = (a); i <= (b); i ++ )
#define rop(i, a, b) for (int i = (a); i < (b); i ++ )

using namespace std;

const int N = 1e7 + 10;
const int mod = 998244353;
int px[N], py[100], pz[100], a[N], b[N], n, x, y, z, ans;
signed main() {
	scanf("%lld%lld%lld%lld", &n, &x, &y, &z);
	if (n > 10000000) return 0;
	px[0] = py[0] = pz[0] = 1;
	rep(i, 1, n) px[i] = px[i - 1] * x % mod;
	rep(i, 1, 50) py[i] = py[i - 1] * y % mod;
	rep(i, 1, 50) pz[i] = pz[i - 1] * z % mod;
	rep(i, 1, n) a[i] = a[i >> 1] + (i & 1);
	rep(i, 1, n) b[i] = b[i / 3] + (i % 3);
	rep(i, 1, n) (ans += px[i] * py[a[i]] % mod * pz[b[i]] % mod) %= mod;
	cout << ans << endl; return 0;
}