#include #include using namespace std; const long long MOD = 998244353; long long modpow(long long a, long long b){ long long ans = 1; while (b > 0){ if (b % 2 == 1){ ans *= a; ans %= MOD; } a *= a; a %= MOD; b /= 2; } return ans; } long long modinv(long long a){ return modpow(a, MOD - 2); } vector inverse(vector f){ int N = f.size(); vector ans(1); ans[0] = 1; for (int i = 1; i < N * 2; i *= 2){ vector ans2 = ans; ans2.resize(i * 4); int N2 = min(N, i * 2); vector f2(i * 4, 0); for (int j = 0; j < N2; j++){ f2[j] = f[j]; } ans2 = atcoder::convolution(ans2, ans2); ans2.resize(i * 2); ans2 = atcoder::convolution(ans2, f2); ans2.resize(i * 2); for (int j = 0; j < i; j++){ ans2[j] = MOD - ans2[j] + ans[j] * 2; ans2[j] %= MOD; } swap(ans, ans2); } ans.resize(N); return ans; } int main(){ long long P1, P2, Q1, Q2; int T; cin >> P1 >> P2 >> Q1 >> Q2 >> T; long long P = P1 * modinv(P2) % MOD; long long Q = Q1 * modinv(Q2) % MOD; vector POW(T + 1); POW[0] = 1; for (int i = 1; i <= T; i++){ POW[i] = POW[i - 1] * Q % MOD; } vector f(T + 1, 0); vector g(T + 1, 0); f[0] = 1; long long p = 1; for (int i = 1; i <= T + 1; i++){ g[i - 1] = p; if (i > T){ break; } f[i] = (MOD - p * P % MOD) % MOD; p *= POW[i]; p %= MOD; } vector finv = inverse(f); vector h = atcoder::convolution(finv, g); cout << h[T] << endl; }