#include using namespace std; // modint // using mint = MontgomeryModInt64; と定義する // mint::set_mod(M); でmodをMに設定する struct MontgomeryModInt64 { using mint = MontgomeryModInt64; using u64 = uint64_t; using u128 = __uint128_t; // static変数 // R = 2 ^ 64 static inline u64 MOD; static inline u64 INV_MOD; // INV_MOD * MOD ≡ 1 (mod 2 ^ 64) static inline u64 T128; // 2 ^ 128 (mod MOD) u64 val; // コンストラクタ // MODを足す? MontgomeryModInt64(): val(0) {} MontgomeryModInt64(long long v): val(MR((u128(v) + MOD) * T128)) {} // 値を返す u64 get() const { u64 res = MR(val); return res >= MOD ? res - MOD : res; } // static関数 static u64 get_mod() { return MOD; } static void set_mod(u64 mod) { MOD = mod; T128 = -u128(mod) % mod; INV_MOD = get_inv_mod(); } // ニュートン法で逆元を求める static u64 get_inv_mod() { u64 res = MOD; for(int i = 0; i < 5; ++i) res *= 2 - MOD * res; return res; } // モンゴメリリダクション static u64 MR(const u128& v) { return (v + u128(u64(v) * u64(-INV_MOD)) * MOD) >> 64; } // 算術演算子 mint operator - () const { return mint() - mint(*this); } mint operator + (const mint& r) const { return mint(*this) += r; } mint operator - (const mint& r) const { return mint(*this) -= r; } mint operator * (const mint& r) const { return mint(*this) *= r; } mint operator / (const mint& r) const { return mint(*this) /= r; } mint& operator += (const mint& r) { if((val += r.val) >= 2 * MOD) val -= 2 * MOD; return *this; } mint& operator -= (const mint& r) { if((val += 2 * MOD - r.val) >= 2 * MOD) val -= 2 * MOD; return *this; } mint& operator *= (const mint& r) { val = MR(u128(val) * r.val); return *this; } mint& operator /= (const mint& r) { *this *= r.inv(); return *this; } mint inv() const { return pow(MOD - 2); } mint pow(u128 n) const { mint res(1), mul(*this); while(n > 0) { if(n & 1) res *= mul; mul *= mul; n >>= 1; } return res; } // その他演算子 bool operator == (const mint& r) const { return (val >= MOD ? val - MOD : val) == (r.val >= MOD ? r.val - MOD : r.val); } bool operator != (const mint& r) const { return (val >= MOD ? val - MOD : val) != (r.val >= MOD ? r.val - MOD : r.val); } // 入力 friend istream& operator >> (istream& is, mint& x) { long long t; is >> t; x = mint(t); return is; } // 出力 friend ostream& operator << (ostream& os, const mint& x) { return os << x.get(); } friend mint modpow(const mint& r, long long n) { return r.pow(n); } friend mint modinv(const mint& r) { return r.inv(); } }; using mint = MontgomeryModInt64; int main() { mint::set_mod(998244353); int P1, P2, Q1, Q2, T; cin >> P1 >> P2 >> Q1 >> Q2 >> T; mint P = mint(P1) / P2; vector Q(T + 1); Q[0] = 1; for(int i = 0; i < T; i++) { Q[i + 1] = Q[i] * Q1 / Q2; } vector> dp(T + 1, vector(T + 1, 0)); dp[0][0] = 1; vector sm(T + 1, 0); for(int i = 0; i < T; i++) { for(int j = 0; j <= T; j++) { sm[i] += dp[i][j]; } dp[i + 1][0] += sm[i] * P; for(int j = 0; j < T; j++) { dp[i + 1][j + 1] = Q[j + 1] * dp[i][j]; } } mint ans = 0; for(int i = 0; i <= T; i++) { ans += dp[T][i]; } cout << ans << endl; return 0; }