#include #include #include using namespace atcoder; using mint = modint998244353; using namespace std; #define rep(i,n) for (int i = 0; i < (n); ++i) #define Inf 100000000000000000 template struct matrix{ F0 func0; F1 func1; vector> value; int height,width; T init_value0,init_value1; matrix(vector> X,F0 f0,F1 f1,T iv0,T iv1):func0(f0),func1(f1){ height = X.size(); width = X[0].size(); value = X; init_value0 = iv0,init_value1 = iv1; } matrix(int h,int w,F0 f0,F1 f1,T iv0,T iv1):func0(f0),func1(f1){ vector> d(h,vector(w,iv0)); height = h,width = w; value = d; init_value0 = iv0,init_value1 = iv1; } void set_1(){ for(int i=0;i R(height,another.width,func0,func1,init_value0,init_value1); for(int i=0;i R(height,width,func0,func1,init_value0,init_value1); R.set_1(); auto temp = *this; while(cnt!=0LL){ if(cnt%2==1){ R *= temp; } temp *= temp; cnt/=2; } return R; } }; int main(){ mint Pa,Pb; int S,T; { int x,y; cin>>x>>y; Pa = x; Pa /= y; cin>>S; } { int x,y; cin>>x>>y; Pb= x; Pb /= y; cin>>T; } //cout<> A(105,vector(105,0)); int c = 52; S = c+S; T = c-T; rep(i,105){ int cur = i; if(cur == S || cur == T){ A[cur][cur] = 1; continue; } if(cur > S || cur < T)continue; rep(j,105){ int c2 = cur + j; if(c2 > S)break; mint p = Pa.pow(j); if(c2 != S){ p *= mint(1)-Pa; } else{ A[c2][cur] = p; continue; } rep(k,105){ int c3 = c2 - k; if(c3 < T)continue; mint pp = p * Pb.pow(k); if(c3 != T){ pp *= mint(1)-Pb; } /* if(i==c&&j==0&&k==0){ cout< AA(A,f0,f1,0,1); long long K; cin>>K; AA = AA.beki(K); vector> B(105,vector(1,0)); B[c][0] = 1; matrix BB(B,f0,f1,0,1); AA = AA * BB; mint ans0 = AA.value[S][0]; mint ans1 = AA.value[T][0]; //cout<<(ans0*81).val()<