#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 i = max(0,k-U); i <= T && i <= k; i++) { cnt2+=Cb(k,i)*f(T-i,Q)*f(U-(k-i),P); } ans += cnt1*cnt2; } cout<