#include<bits/stdc++.h>
#define L(i, j, k) for(int i = (j); i <= (k); ++i)
#define R(i, j, k) for(int i = (j); i >= (k); --i)
#define ll long long 
#define vi vector < int > 
#define sz(a) ((int) (a).size())
#define ll long long 
#define ull unsigned long long
#define me(a, x) memset(a, x, sizeof(a)) 
using namespace std;
const int N = 1 << 19, mod = 164511353;
int qpow(int x, int y = mod - 2) {
	int res = 1;
	for(; y; x = (ll) x * x % mod, y >>= 1) if(y & 1) res = (ll) res * x % mod;
	return res;
}
int fac[N], ifac[N], inv[N];
void init(int x) {
	fac[0] = ifac[0] = inv[1] = 1;
	L(i, 2, x) inv[i] = (ll) (mod - mod / i) * inv[mod % i] % mod;
	L(i, 1, x) fac[i] = (ll) fac[i - 1] * i % mod, ifac[i] = (ll) ifac[i - 1] * inv[i] % mod;
} 
int C(int x, int y) {
	return x < y || y < 0 ? 0 : (ll) fac[x] * ifac[y] % mod * ifac[x - y] % mod;
}
inline int sgn(int x) {
	return (x & 1) ? mod - 1 : 1;
}

const int r = 41;
int n, k, m, type;
int prm[N], mt[N], rmt[N], cnt[N], top;
int fen[N];
int val[N];

int main() {
	ios :: sync_with_stdio(false); 
	cin.tie(0); cout.tie(0);
	init(N - 10);
	cin >> n >> k >> m;
	cin >> type;
	if(type == 1) {
	
	L(i, 2, k) if(k % i == 0) {
		int c = 0;
		while(k % i == 0) k /= i, ++c;
		++top, prm[top] = i, cnt[top] = c, mt[top] = rmt[top] = 1;
		L(j, 1, c) mt[top] *= i;
		L(j, 1, c - 1) rmt[top] *= i;
	} 
	
	const ll Z = (ll)mod * 41;
	map < ll, int > vc;
	vc[1] = 1; 
	L(i, 1, top) {
		L(cnt, 0, n - 1) {
			int qp = qpow(prm[i], cnt);
			map < ll, int > nvc;
			for(auto&u : vc) {
				ll val = u.first;
				int w = u.second;
				(nvc[mt[i] * val % Z] += w) %= mod;
				(nvc[rmt[i] * val % Z] += mod - (ll) qp * w % mod) %= mod;
			}
			swap(vc, nvc);
		}
	}
	
	int ans = 0;
	for(auto&s : vc) {
		ll x = s.first;
		int y = s.second;
		int pr = (ll)y * qpow(2, x % 41) % mod;
		L(i, 0, m) {
			if(i) pr = (ll) pr * (x % mod - i + 1) % mod * inv[i] % mod;
			(fen[i] += pr) %= mod;
		}
	}
	
	L(i, 0, m) 
		fen[i] = (ll)fen[i] * qpow(qpow(2), i) % mod; 
	L(i, 0, m) {
		L(j, 0, i) {
			(val[j] += (ll)fen[i] * sgn(i - j) % mod * C(i, j) % mod) %= mod;
		}
	}
	
	int ANS = 0;
	L(i, 0, m) {
		(ANS += (ll)val[i] * qpow(i, m) % mod) %= mod;
	}
	cout << ANS << '\n';
	
	}
	return 0;
}