#include using i64=long long; using u64=unsigned long long; using namespace std; template struct modint{ u64 val; modint(i64 val_=0):val((val_%i64(mod)+i64(mod))%i64(mod)){} modint operator-(){ return (val==0)?0:mod-val; } modint operator+(modint rhs){ return modint(*this)+=rhs; } modint operator-(modint rhs){ return modint(*this)-=rhs; } modint operator*(modint rhs){ return modint(*this)*=rhs; } modint operator/(modint rhs){ return modint(*this)/=rhs; } modint pow(i64 rhs){ modint res=1,now=(*this); while(rhs){ res*=((rhs&1)?now:1),now*=now,rhs>>=1; } return res; } modint &operator+=(modint rhs){ val+=rhs.val,val-=((val>=mod)?mod:0); return (*this); } modint &operator-=(modint rhs){ val+=((val>(std::istream& is,modint& x){ u64 t; is>>t,x=t; return is; } }; 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<0||k<0||n mint; string S; int A,B,C; cin>>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; if(T<0 || U<0 || V<0){ cout<<0< Cb; mint ans=0; for(int i=0;i<=min(P,T);++i){ for(int j=0;j<=min(Q,U);++j){ for(int k=0;k<=min(R,V);++k){ mint tmp=1; tmp*=Cb(P,i)*Cb(Q,j)*Cb(R,k); tmp*=Cb(N+(T-i)+j+k,T-i); tmp*=Cb(N+T+(U-j)+k,U-j); tmp*=Cb(N+T+U+(V-k),V-k); ans+=tmp*((i+j+k)%2?-1:1); } } } cout<