#include <bits/stdc++.h>

#define _overload(_1,_2,_3,name,...) name
#define _rep(i,n) _range(i,0,n)
#define _range(i,a,b) for(int i=int(a);i<int(b);++i)
#define rep(...) _overload(__VA_ARGS__,_range,_rep,)(__VA_ARGS__)

#define _rrep(i,n) _rrange(i,n,0)
#define _rrange(i,a,b) for(int i=int(a)-1;i>=int(b);--i)
#define rrep(...) _overload(__VA_ARGS__,_rrange,_rrep,)(__VA_ARGS__)

#define _all(arg) begin(arg),end(arg)
#define uniq(arg) sort(_all(arg)),(arg).erase(unique(_all(arg)),end(arg))
#define getidx(ary,key) lower_bound(_all(ary),key)-begin(ary)
#define clr(a,b) memset((a),(b),sizeof(a))
#define bit(n) (1LL<<(n))
#define popcount(n) (__builtin_popcountll(n))

using namespace std;

template<class T>bool chmax(T &a, const T &b) { return (a<b)?(a=b,1):0;}
template<class T>bool chmin(T &a, const T &b) { return (b<a)?(a=b,1):0;}

using ll=long long;
using R=long double;
const R EPS=1e-9; // [-1000,1000]->EPS=1e-8 [-10000,10000]->EPS=1e-7
inline int sgn(const R& r){return(r > EPS)-(r < -EPS);}
inline R sq(R x){return sqrt(max<R>(x,0.0));}

const int dx[8]={1,0,-1,0,1,-1,-1,1};
const int dy[8]={0,1,0,-1,1,1,-1,-1};

// Problem Specific Parameter:

const int nmax=100010;
string s;

vector<int> graph[nmax];
vector<int> inv[nmax];

vector<ll> dpc[nmax];
vector<ll> dpw[nmax];


void dfs(int v,int i){
	int p=graph[v][i];
	if(dpc[v][i]>=0) return;
	dpc[v][i]=(s[v]=='c');
	dpw[v][i]=(s[v]=='w');	
	rep(j,graph[v].size()){
		const int v2=graph[v][j];
		if(v2==p) continue;
		const int k=inv[v][j];
		dfs(v2,k);
		dpc[v][i]+=dpc[v2][k];
		dpw[v][i]+=dpw[v2][k];
	}
}

int main(void){
	int n;
	cin >> n >> s;
	
	rep(i,n-1){
		int a,b;
		cin >> a >> b;
		graph[a-1].push_back(b-1);
		inv[a-1].push_back(inv[b-1].size());
		graph[b-1].push_back(a-1);
		inv[b-1].push_back(inv[a-1].size()-1);
	}

	rep(v,n) dpc[v].assign(graph[v].size(),-1),dpw[v].assign(graph[v].size(),-1);
	rep(v,n) rep(i,graph[v].size()) dfs(v,i);

	ll ans=0LL;

	rep(v,n){
		if(s[v]=='c') continue;
		ll tc=0LL,tw=0LL;

		rep(i,graph[v].size()){	
			const int v2=graph[v][i];
			const int j=inv[v][i];
			tc+=dpc[v2][j];
			tw+=dpw[v2][j];
			ans-=1LL*dpc[v2][j]*dpw[v2][j];
		}
		ans+=1LL*tc*tw;
	}

	cout << ans << endl;
	return 0;
}