#include #include using mint=atcoder::static_modint<998244353>; using namespace std; template struct comb{ vector dat,idat; comb(int mx=3000000):dat(mx+1,1),idat(mx+1,1){ for(int i=1;i<=mx;++i){ dat[i]=dat[i-1]*i; } idat[mx]/=dat[mx]; for(int i=mx;i>0;--i){ idat[i-1]=idat[i]*i; } } T operator()(int n,int k){ if(n==k){ return 1; } if(n<0||k<0||n>S>>A>>B>>C; int N=S.size(); int P=count(S.begin(),S.end(),'a'),Q=count(S.begin(),S.end(),'b'),R=count(S.begin(),S.end(),'c'); int T=A-P,U=B-Q,V=C-R; comb Cb; auto f=[&](int n,int m){ return Cb(n+m-1,m-1); }; mint ans=0; for(int k=0;k<=T+U;k++){ mint cnt1=0,cnt2=0; for(int i=0;i<=k;i++){ cnt1+=f(i,R)*f(V,P+Q+T+U-i+1); } for(int j=max(0,k-U);j<=min(T,k);j++){ cnt2+=Cb(k,j)*f(T-j,Q)*f(U-(k-j),P); } ans+=cnt1*cnt2; } cout<