#pragma GCC target("avx") #pragma GCC optimize("O3") #pragma GCC optimize("unroll-loops") #include using namespace std; using ll=long long; using datas=pair; using ddatas=pair; using tdata=pair; using vec=vector; using mat=vector; using pvec=vector; using pmat=vector; #define For(i,a,b) for(i=a;i<(ll)b;++i) #define bFor(i,b,a) for(i=b,--i;i>=(ll)a;--i) #define rep(i,N) For(i,0,N) #define rep1(i,N) For(i,1,N) #define brep(i,N) bFor(i,N,0) #define brep1(i,N) bFor(i,N,1) #define all(v) (v).begin(),(v).end() #define allr(v) (v).rbegin(),(v).rend() #define vsort(v) sort(all(v)) #define vrsort(v) sort(allr(v)) #define uniq(v) vsort(v);(v).erase(unique(all(v)),(v).end()) #define endl "\n" #define popcount __builtin_popcountll #define eb emplace_back #define print(x) cout< ostream& operator<<(ostream& os,const pair& p){return os<<"("< ostream& operator<<(ostream& os,const vector& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const set& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const multiset& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< ostream& operator<<(ostream& os,const map& v){ os<<"{";bool f=false; for(auto& x:v){if(f)os<<",";os< inline bool chmax(T& a,const T b){bool x=a inline bool chmin(T& a,const T b){bool x=a>b;if(x)a=b;return x;} #ifdef DEBUG void debugg(){cout<void debugg(const T& x,const Args&... args){cout<<" "<H)rep(j,this->W){ (*this)[i][j]+=b[i][j]; if((*this)[i][j]>=mod)(*this)[i][j]-=mod; } return (*this); } modmatrix operator -=(modmatrix b){ ll i,j; rep(i,this->H)rep(j,this->W){ (*this)[i][j]-=b[i][j]; if((*this)[i][j]<0)(*this)[i][j]+=mod; } return (*this); } modmatrix operator *=(modmatrix b){ ll i,j,k; assert(this->W==b.H); modmatrix c(this->H,b.W); rep(i,this->H){ rep(k,this->W)rep(j,b.W){ c[i][j]+=(*this)[i][k]*b[k][j]%mod; } rep(j,b.W)c[i][j]%=mod; } (*this)=c; return (*this); } modmatrix operator ^=(ll K){ assert(this->H==this->W); modmatrix c(this->H); ll i; rep(i,this->H)c[i][i]=1; if(K&1)c*=(*this); while(K){ K>>=1; (*this)*=(*this); if(K&1)c*=(*this); } this->a.swap(c.a); return (*this); } modmatrix operator +(modmatrix c){ return modmatrix(*this)+=c; } modmatrix operator -(modmatrix c){ return modmatrix(*this)-=c; } modmatrix operator *(modmatrix c){ return modmatrix(*this)*=c; } modmatrix operator ^(ll K){ return modmatrix(*this)^=K; } modmatrix del(ll eh,ll ew){ ll i,j; mat res; rep(i,H){ if(i==eh)continue; res.resize(res.size()+1); rep(j,W){ if(j==ew)continue; res.back().eb(a[i][j]); } } return res; } ll determinant(){ assert(H==W); ll i,j,k,ans=1; auto b(a); rep(i,H){ if(!b[i][i]){ For(j,i+1,H)if(b[j][i]){ swap(b[i],b[j]); ans*=-1; break; } if(j==H)return 0; } (ans*=b[i][i])%=mod; For(j,i+1,W)(b[i][j]*=modinv(b[i][i]))%=mod; For(j,i+1,H)if(b[j][i]){ ll x=mod-b[j][i]; b[j][i]=0; For(k,i+1,W){ b[j][k]+=x*b[i][k]%mod; if(b[j][k]>=mod)b[j][k]-=mod; } } } if(ans<0)ans+=mod; return ans; } modmatrix inv(){ assert(H==W); ll i,j,k; modmatrix b(a),c(H,W); rep(i,H)c[i][i]=1; rep(i,H){ if(!b[i][i]){ For(j,i+1,H)if(b[j][i]){ swap(b[i],b[j]); swap(c[i],c[j]); break; } if(j==H)assert(false); } ll x=modinv(b[i][i]); rep(j,W){ (b[i][j]*=x)%=mod; (c[i][j]*=x)%=mod; } rep(j,H){ if(i==j)continue; if(b[j][i]){ x=mod-b[j][i]; rep(k,W){ b[j][k]+=x*b[i][k]%mod; if(b[j][k]>=mod)b[j][k]-=mod; c[j][k]+=x*c[i][k]%mod; if(c[j][k]>=mod)c[j][k]-=mod; } } } } return c; } void out(){ for(auto x:a)output(x); } }; ll p,q; ll S,T,K; int main(){ startupcpp(); ll i,j; cin>>p>>j>>S; (p*=modinv(j))%=mod; cin>>q>>j>>T>>K; (q*=modinv(j))%=mod; mat A(S+T+1,vec(S+T+1,0)); A[0][0]=A[S+T][S+T]=1; mat B=A; rep1(i,S+T){ A[i][i]=(mod+1-p)%mod; A[S+T][S+T-i]=A[S+T][S+T+1-i]*p%mod; B[i][i]=(mod+1-q)%mod; B[0][i]=B[0][i-1]*q%mod; } rep1(i,S+T){ for(j=i;j0;--j)A[i][j-1]=A[i][j]*p%mod; } debug(A); debug(B); modmatrix a(A),b(B),first(S+T+1,1); first[T][0]=1; a*=b; a^=K; a*=first; print(a[S+T][0]); print(a[0][0]); }