#include <bits/stdc++.h>
using namespace std;

#define M 1000003
#define m(t)r=r*fac[(t)%M]%M*fac[(t)/M%M]%M*fac[(t)/M/M]%M

char T[M];
long S[26],x[M],y[M];
int inv[M],fac[M];
long r=1;

main(){
	fac[0]=1;
	for(int i=1;i<M;++i){
		fac[i]=1ll*fac[i-1]*i%M;
		long e,s=M,t=i,u=1,v=0;
		while(s){
			e=t/s;
			t-=e*s;
			u-=e*v;
			swap(t,s);
			swap(u,v);
		}
		inv[i]=(u+M)%M;
	}
	for(int c=0;c<26;++c)cin>>S[c];
	cin>>T;
	for(int c=0;c<26;++c){
		priority_queue<pair<double,int>> q;
		int n=0,s=0,t=-1;
		for(int i=0;!i||T[i-1];++i){
			if(T[i]==c+'a'){
				if(t<0)t=i;
			}else{
				if(t>=0){
					int xl=i-t;
					x[n]=xl;
					++n;
					s+=xl;
					if(s>S[c])r=0;
					t=-1;
				}
			}
		}
		if(s) {
			long d=S[c]/s;
			for(int j=0;j<n;++j){
				y[j]=x[j]*d;
				q.push(make_pair(1.*~y[j]/x[j],j));
			}
			for(S[c]%=s;S[c]--;){
				int j=q.top().second;q.pop();
				++y[j];
				q.push(make_pair(1.*~y[j]/x[j],j));
			}
			for(int j=0;j<n;++j){
				m(y[j]);
				r=inv[r];
				m(y[j]-x[j]);
				m(x[j]);
				r=inv[r];
			}
		}
	}
	cout<<r;
}