#include <bits/stdc++.h>
using namespace std;
const int N = 3e5 + 5;
typedef long long ll;
template<typename T> inline bool ckmax(T &x,const T &y) { return (x < y) && (x = y,true);}
int n;
vector<int> G[N];
char s[5];

namespace SolveA { 
	inline void Main() { printf("%d\n",n);exit(0);}
}
namespace SolveAB {
	inline void Main() { printf("%d\n",n - 1);exit(0);}
}
namespace SolveAA {
	inline void Main() { printf("%d\n",n + n - 2);exit(0);}
}
namespace SolveABC {
	ll f[N][2],g[N]; // g 表示 u 点上方还有一个 1/3，当前是 2 的代价
	void dfs0(int x,int fa) {
		for(auto y : G[x]) if(y != fa) dfs0(y,x);
		vector<ll> V;ll sum1 = 0;
		for(auto y : G[x])
			if(y != fa) f[x][0] += max({g[y],f[y][0]}),V.push_back(f[y][0] - f[y][1]),sum1 += f[y][1];
		sort(V.begin(),V.end(),greater<ll>());
		f[x][1] = sum1;
		for(int i = 0;i < V.size();i++) {
			sum1 += V[i];
			ll val = sum1 + 1ll * (i + 2 >> 1) * (i + 1 >> 1);
			ckmax(f[x][1],val);ckmax(g[x],val + (i + 2 >> 1));
		}
		// printf("%d %d %d\n",f[x][0],f[x][1],g[x]);
	}
	void Main() { dfs0(1,0); printf("%lld\n",max(f[1][0],f[1][1]));exit(0);}
}
namespace SolveABA {
	ll f[N][2],g[N];
	void dfs0(int x,int fa) {
		for(auto y : G[x]) if(y != fa) dfs0(y,x);
		vector<ll> V;ll sum1 = 0;
		for(auto y : G[x])
			if(y != fa) f[x][0] += max({g[y],f[y][0]}),V.push_back(f[y][0] - f[y][1]),sum1 += f[y][1];
		sort(V.begin(),V.end(),greater<ll>());
		f[x][1] = sum1;
		for(int i = 0;i < V.size();i++) {
			sum1 += V[i];
			ll val = sum1 + 1ll * (i + 1) * i / 2;
			ckmax(f[x][1],val);ckmax(g[x],val + (i + 1));
		}
	}
	void Main() { dfs0(1,0); printf("%lld\n",2 * max(f[1][0],f[1][1]));exit(0);}
}
namespace SolveAAB {
	ll f[N][2],g[N][2]; // g 表示 x = A，上面还有一个 A/B
	void dfs0(int x,int fa) {
		for(auto y : G[x]) if(y != fa) dfs0(y,x);
		vector<ll> V;ll sum1 = 0;
		for(auto y : G[x]) if(y != fa) f[x][0] += g[y][0],V.push_back(f[y][1] - g[y][0]),sum1 += g[y][0];
		sort(V.begin(),V.end(),greater<ll>());
		for(int i = 0;i < V.size();i++) {
			sum1 += V[i];
			ll val = sum1 + 1ll * (i + 1) * ((int)V.size() - i - 1);
			ckmax(f[x][0],val);ckmax(g[x][0],val + (i + 1));ckmax(g[x][1],val + ((int)V.size() - i - 1));
		}
		for(auto y : G[x]) if(y != fa) f[x][1] += max({f[y][1],g[y][1]});
	}
	void Main() { dfs0(1,0);printf("%lld\n",max(f[1][0],f[1][1]));exit(0);}
}
int main() {
	cin >> n;cin >> (s + 1);
	for(int i = 1,u,v;i < n;i++) cin >> u >> v,G[u].push_back(v),G[v].push_back(u);
	
	int len = strlen(s + 1);
	if(len == 1) SolveA::Main();
	if(len == 2 && s[1] == s[2]) SolveAA::Main();
	if(len == 2 && s[1] != s[2]) SolveAB::Main();
	if(s[1] != s[3] && s[1] != s[2] && s[2] != s[3]) SolveABC::Main();
	if(s[1] == s[2] || s[2] == s[3]) SolveAAB::Main();
	if(s[1] == s[3]) SolveABA::Main();
	return 0;
}