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

using VI = vector<int>;
using VVI = vector<VI>;
using PII = pair<int, int>;
using LL = long long;
using VL = vector<LL>;
using VVL = vector<VL>;
using PLL = pair<LL, LL>;
using VS = vector<string>;

#define ALL(a)  begin((a)),end((a))
#define RALL(a) (a).rbegin(), (a).rend()
#define PB push_back
#define EB emplace_back
#define MP make_pair
#define SZ(a) int((a).size())
#define SORT(c) sort(ALL((c)))
#define RSORT(c) sort(RALL((c)))
#define UNIQ(c) (c).erase(unique(ALL((c))), end((c)))

#define FOR(i,a,b) for(int i=(a);i<(b);++i)
#define REP(i,n)  FOR(i,0,n)

#define FF first
#define SS second
template<class S, class T>
istream& operator>>(istream& is, pair<S,T>& p){
  return is >> p.FF >> p.SS;
}
template<class S, class T>
ostream& operator<<(ostream& os, const pair<S,T>& p){
  return os << p.FF << " " << p.SS;
}
template<class T>
void maxi(T& x, T y){
  if(x < y) x = y;
}
template<class T>
void mini(T& x, T y){
  if(x > y) x = y;
}


const double EPS = 1e-10;
const double PI  = acos(-1.0);
const LL MOD = 1e6+3;

using LLL = __int128;

LLL powm(LLL a, LLL b){
  return (b==0? 1: powm(a*a%MOD,b/2) * (b%2?a:1) % MOD);
}
LLL inv(LLL x){
  return powm(x, MOD-2);
}

LLL solve(LLL rest, LLL sum, VL& vs){
  if(rest == 0 || vs.empty()) return 1;

  LLL n = SZ(vs);
  if(rest >= n*(MOD-1)+1)
	return 0;
  
  using FRAC = pair<LLL,LLL>; // (first / second)
  auto lt = [](const FRAC& l, const FRAC& r){
	return l.first * r.second < r.first * l.second;
  };
  priority_queue<FRAC,vector<FRAC>,decltype(lt)> pq(lt);
  for(LL v: vs)
	pq.push(MP(v+1, 1));
  rest -= sum;

  LLL res = 1;
  while(rest){
	--rest;
	auto p = pq.top();
	pq.pop();

	if(p.FF >= MOD || res == 0) return 0;
	(res = res * p.first % MOD * inv(p.second)) %= MOD;

	pq.push(MP(p.first+1, p.second+1));
  }
  
  return res;
}

int main(){
  cin.tie(0);
  ios_base::sync_with_stdio(false);

  VL S(26);
  REP(i,26) cin >> S[i];

  string T;
  cin >> T;
  VL sum(26);
  VVL vs(26);

  REP(i,SZ(T)){
	int c = T[i] - 'a';
	sum[c]++;
	if(!i || T[i] != T[i-1])
	  vs[c].PB(1);
	else
	  vs[c].back()++;
  }

  REP(i,26) if(S[i] < sum[i]){
	cout << 0 << endl;
	return 0;
  }

  LLL ans = 1;
  REP(c,26){
	LLL res = solve(S[c], sum[c], vs[c]);
	(ans *= res) %= MOD;
  }
  cout << (LL)(ans) << endl;

  return 0;
}