#include<bits/stdc++.h>
using namespace std;
#define pb push_back
#define eb emplace_back
//mt19937 rnd(chrono::steady_clock::now().time_since_epoch().count());
#define mod 998244353
#define ll long long
#define int long long
#define inf 0x3f3f3f3f
#define INF 0x3f3f3f3f3f3f3f3f
inline int read()
{
	char ch=getchar(); int nega=1; while(!isdigit(ch)) {if(ch=='-') nega=-1; ch=getchar();}
	int ans=0; while(isdigit(ch)) {ans=ans*10+ch-48;ch=getchar();}
	if(nega==-1) return -ans;
	return ans;
}
void print(vector<int> x){for(int i=0;i<(int)x.size();i++) printf("%d%c",x[i]," \n"[i==(int)x.size()-1]);}
#define N 300005
char _str[5];
int str[3],m;
vector<int> G[N];
int n;
ll f[N][3][3],g[N][3],w12[N],w[N];
int a[N],b[N],dpos[N],reva[N];
vector<int> d[N];
struct Node
{
	int id; ll w;
	bool operator < (const Node &x) const
	{
		if(w!=x.w) return w>x.w;
		return id<x.id;
	}
};
ll solve(int n)
{
	// for(int i=0;i<n;i++) printf("%d%c",g[i][0]," \n"[i==n-1]);
	// for(int i=0;i<n;i++) printf("%d%c",g[i][1]," \n"[i==n-1]);
	// for(int i=0;i<n;i++) printf("%d%c",g[i][2]," \n"[i==n-1]);
	for(int i=0;i<=n;i++) w12[i]=0;
	for(int i=0;i<n;i++)
	{
		int dpos=g[i][1]-g[i][2];
		dpos=min(dpos,n+1);
		dpos=max(dpos,0LL);
		w[i]=g[i][0]-max(g[i][1],g[i][2]);
		w12[dpos+1]++;
		::dpos[i]=dpos;
		d[dpos].push_back(i);
	}
	for(int i=1;i<=n;i++) w12[i]+=w12[i-1];
	for(int i=0;i<n;i++) w12[0]+=max(g[i][1],g[i][2]);
	for(int i=1;i<=n;i++) w12[i]+=w12[i-1];
	for(int i=1;i<=n;i++) a[i]=b[i]=i-1;
	sort(a+1,a+n+1,[&](int x,int y){
		return w[x]>w[y];
	});
	for(int i=1;i<=n;i++) reva[a[i]]=i;
	sort(b+1,b+n+1,[&](int x,int y){
		return w[x]+dpos[x]>w[y]+dpos[y];
	});
	ll ans=w12[0];
	auto gval=[&](int i,int u)
	{
		if(i<=dpos[u]) return w[u];
		else return w[u]-(i-dpos[u]);
	};
	ll suma=0;
	set<Node> ina,inb,outa,outb;
	for(int i=0;i<n;i++)
	{
		if(dpos[i]==0) outb.insert((Node){i,w[i]+dpos[i]});
		else outa.insert((Node){i,w[i]});
	}
	for(int i=1;i<=n;i++)
	{
		while(!inb.empty()&&!outa.empty()&&gval(i,inb.rbegin()->id)<gval(i,outa.begin()->id))
		{
			Node tw=*inb.rbegin();
			inb.erase(inb.find(tw));
			outb.insert(tw);
			tw=*outa.begin();
			outa.erase(outa.begin());
			ina.insert(tw);
			suma++;
		}
		suma-=inb.size();
		for(int j:d[i])
		{
			if(ina.find((Node){j,w[j]})!=ina.end())
			{
				ina.erase(ina.find((Node){j,w[j]}));
				inb.insert((Node){j,w[j]+dpos[j]});
			}
			else
			{
				assert(outa.find((Node){j,w[j]})!=outa.end());
				outa.erase(outa.find((Node){j,w[j]}));
				outb.insert((Node){j,w[j]+dpos[j]});
			}
		}
		if(!outa.empty()&&(outb.empty()||gval(i,outa.begin()->id)>gval(i,outb.begin()->id)))
		{
			Node tw=*outa.begin();
			outa.erase(tw);
			ina.insert(tw);
			suma+=gval(i,tw.id);
		}
		else
		{
			Node tw=*outb.begin();
			outb.erase(tw);
			inb.insert(tw);
			suma+=gval(i,tw.id);
		}
		ans=max(ans,w12[i]+suma);
	}
	for(int i=0;i<=n+1;i++) d[i].clear();
	return ans;
}
void dfs(int u,int fa)
{
	for(int v:G[u]) if(v!=fa) dfs(v,u);
	for(int c=0;c<m;c++) for(int d=0;d<m;d++)
	{
		if(u==1) d=m;
		if(str[1]==c)
		{
			vector<ll> tw[m];
			for(int v:G[u]) if(v!=fa) for(int j=0;j<m;j++) tw[j].push_back(f[v][j][c]);
			for(int j=0;j<m;j++)
			{
				if(str[0]==j&&str[2]==d) for(int i=0;i<(int)tw[j].size();i++) tw[j][i]++;
				if(str[0]==d&&str[2]==j) for(int i=0;i<(int)tw[j].size();i++) tw[j][i]++;
			}
			ll ans=0;
			if(str[0]==str[2])
			{
				if(m==1)
				{
					ans=1LL*(int)tw[0].size()*((int)tw[0].size()-1);
					for(int i=0;i<(int)tw[0].size();i++) ans+=tw[0][i];
				}
				else
				{
					priority_queue<ll> q;
					ll cur=0;
					for(int i=0;i<(int)tw[0].size();i++) q.push(tw[0][i]-tw[1][i]),cur+=tw[1][i];
					ans=cur;
					for(int i=1;i<=(int)tw[0].size();i++)
					{
						cur+=q.top(); q.pop();
						ans=max(ans,cur+1LL*i*(i-1));
					}
				}
			}
			else
			{
				for(int j=0;j<3;j++)
				{
					for(int i=0;i<(int)tw[0].size();i++) g[i][j]=tw[str[j]][i];
				}
				ans=solve(tw[0].size());
			}
			f[u][c][min(d,m-1)]=ans;
		}
		else
		{
			ll w=0;
			for(int v:G[u]) if(v!=fa)
			{
				ll mx=0;
				for(int j=0;j<m;j++) mx=max(mx,f[v][j][c]);
				w+=mx;
			}
			f[u][c][min(d,m-1)]=w;
		}
	}
	// printf("%d :\n",u);
	// for(int i=0;i<m;i++) for(int j=0;j<m;j++) printf("%d%c",f[u][i][j]," \n"[j==m-1]);
}
signed main()
{
	cin>>n; scanf("%s",_str);
	for(int i=1;i<n;i++)
	{
		int u=read(),v=read();
		G[u].pb(v),G[v].pb(u);
	}
	if(strlen(_str)==1) cout<<n<<endl,exit(0);
	if(strlen(_str)==2)
	{
		ll ans=0;
		for(int i=1;i<=n;i++) ans+=G[i].size();
		if(_str[0]==_str[1]) cout<<ans<<endl;
		else cout<<ans/2<<endl;
		return 0;
	}
	for(int i=0;i<3;i++)
	{
		int ok=0;
		for(int j=0;j<i;j++) if(_str[i]==_str[j]) str[i]=str[j],ok=1;
		if(!ok) str[i]=m++;
	}
	dfs(1,0);
	ll ans=0;
	for(int c=0;c<m;c++) ans=max(ans,f[1][c][m-1]);
	cout<<ans<<endl;
	return 0;
}