#include<bits/stdc++.h>
typedef long long ll;
using namespace std;
const int N = 1e6 + 5, inf = 0x3f3f3f3f;
int n, m, f[N], x[N], y[N], fd[N * 20], Tp[N * 20], fa[N], col[N], pos[N];
ll ans, w[N];
struct node {
	pair<int, int> fi, se;
	node(int a = inf, int b = 0, int c = inf, int d = 0) : fi({a, b}), se({c, d}) {}
	void chkmin(node tmp) {
		if(fi.first > tmp.fi.first) swap(fi, tmp.fi);
		if(se.first > tmp.fi.first && col[fi.second] != col[tmp.fi.second]) swap(se, tmp.fi);
		if(se.first > tmp.se.first && col[fi.second] != col[tmp.se.second]) swap(se, tmp.se);
	}
};
void chkmin(int p, node tmp, int dat) {
	// cerr<<p<<' '<<col[p]<<" "<<tmp.fi.second<<' '<<tmp.se.second<<' '<<dat<<' '<<tmp.fi.first<<'\n';
	p = col[p];
	if(tmp.fi.first + dat < w[p] && p != col[tmp.fi.second]) w[p] = tmp.fi.first + dat, pos[p] = tmp.fi.second;
	if(tmp.se.first + dat < w[p] && p != col[tmp.se.second]) w[p] = tmp.se.first + dat, pos[p] = tmp.se.second;
}
int FindSet(int x) {return x == fa[x] ? x : fa[x] = FindSet(fa[x]);}
vector<int> E[N], vec[N], Vec[N];
int dfn[N], siz[N], cnt, tcnt;
void dfs(int x) {
	dfn[x] = ++cnt, fd[++tcnt] = x, siz[x] = 1;
	for(auto y : E[x]) dfs(y), siz[x] += siz[y];
	fd[++tcnt] = -x;
}
int rt, rot[N], tot, lc[N * 20], rc[N * 20];
node tr[N * 20];
void pushup(int p) {
	tr[p] = node();
	if(lc[p]) tr[p].chkmin(tr[lc[p]]);
	if(rc[p]) tr[p].chkmin(tr[rc[p]]);
}
void change1(int &p, int l, int r, int x, node val) {
	int mid = l + r >> 1;
	if(!p) p = ++tot, tr[p] = node(), lc[p] = rc[p] = 0;
	if(l == r) return tr[p].chkmin(val), void();
	if(x <= mid) change1(lc[p], l, mid, x, val);
	else change1(rc[p], mid + 1, r, x, val); 
	pushup(p);
}
node query1(int p, int l, int r, int ql, int qr) {
	int mid = l + r >> 1;
	if(!p) return node();
	if(ql <= l && r <= qr) return tr[p];
	node ans;
	if(ql <= mid) ans.chkmin(query1(lc[p], l, mid, ql, qr));
	if(mid < qr) ans.chkmin(query1(rc[p], mid + 1, r, ql, qr));
	return ans;
}
node res, *st[N * 20], Val[N * 20]; int tp;
void change2(int &p, int l, int r, int ql, int qr, node val) {
	int mid = l + r >> 1;
	if(!p) p = ++tot, tr[p] = node(), lc[p] = rc[p] = 0;
	if(ql <= l && r <= qr) return st[++tp] = tr + p, Val[tp] = *st[tp], tr[p].chkmin(val), void();
	if(ql <= mid) change2(lc[p], l, mid, ql, qr, val);
	if(mid < qr) change2(rc[p], mid + 1, r, ql, qr, val);
}
void query2(int p, int l, int r, int x) {
	int mid = l + r >> 1;
	if(!p) return;
	res.chkmin(tr[p]);
	if(l == r) return;
	if(x <= mid) query2(lc[p], l, mid, x);
	else query2(rc[p], mid + 1, r, x);
}
int merge(int pl, int pr) {
	if(!pl || !pr) return pl | pr;
	tr[pl].chkmin(tr[pr]);
	lc[pl] = merge(lc[pl], lc[pr]), rc[pl] = merge(rc[pl], rc[pr]);
	return pl;
}
void solve1() {//Case1: sxi+syi+sxj+syj
	node tmp;
	for(int i = 1; i <= m; i++) tmp.chkmin(node(siz[x[i]] + siz[y[i]], i));
	for(int i = 1; i <= m; i++) chkmin(i, tmp, siz[x[i]] + siz[y[i]]);
}
void solve2() {//Case2: sxi+syi+sxj-syj -> yj in yi
	rt = tot = tp = 0;
	for(int i = 1; i <= m; i++) 
		change1(rt, 1, n, dfn[y[i]], node(siz[x[i]] - siz[y[i]], i));
	for(int i = 1; i <= m; i++) 
		chkmin(i, query1(rt, 1, n, dfn[y[i]], dfn[y[i]] + siz[y[i]] - 1), siz[x[i]] + siz[y[i]]);
}
void solve3() {//Case3: sxi-syi+sxj+syj -> yi in yj
	rt = tot = tp = 0;
	for(int i = 1; i <= m; i++) 
		change2(rt, 1, n, dfn[y[i]], dfn[y[i]] + siz[y[i]] - 1, node(siz[x[i]] + siz[y[i]], i));
	for(int i = 1; i <= m; i++)
		res = node(), query2(1, 1, n, dfn[y[i]]), chkmin(i, res, siz[x[i]] - siz[y[i]]);
}
void solve4() {//Case4: sxi+syi-sxi+syi -> xj in xi
	rt = tot = tp = 0;
	for(int i = 1; i <= m; i++) 
		change1(rt, 1, n, dfn[x[i]], node(-siz[x[i]] + siz[y[i]], i));
	for(int i = 1; i <= m; i++) 
		chkmin(i, query1(rt, 1, n, dfn[x[i]], dfn[x[i]] + siz[x[i]] - 1), siz[x[i]] + siz[y[i]]);
}
void solve5() {//Case5: -sxi+syi+sxj+syj -> xi in xj
	rt = tot = tp = 0;
	for(int i = 1; i <= m; i++) 
		change2(rt, 1, n, dfn[x[i]], dfn[x[i]] + siz[x[i]] - 1, node(siz[x[i]] + siz[y[i]], i));
	for(int i = 1; i <= m; i++)
		res = node(), query2(1, 1, n, dfn[x[i]]), chkmin(i, res, -siz[x[i]] + siz[y[i]]);
}
void solve6() {//Case6: sxi+syi-sxj-syj -> xj in xi, yj in yi
	tot = tp = 0;
	for(int i = 1; i <= n; i++) rot[i] = 0;
	for(int i = 1; i <= m; i++) 
		change1(rot[x[i]], 1, n, dfn[y[i]], node({-siz[x[i]] - siz[y[i]], i}));
	for(int i = n; i >= 1; i--) {
		for(auto p : vec[i])
			chkmin(p, query1(rot[i], 1, n, dfn[y[p]], dfn[y[p]] + siz[y[p]] - 1), siz[x[p]] + siz[y[p]]);
		if(f[i]) rot[f[i]] = merge(rot[f[i]], rot[i]);
	}
}
void solve7() {//Case7: sxi-syi-sxj+syj -> xj in xi, yi in yj
	tot = tp = 0;
	for(int i = 1; i <= n; i++) rot[i] = 0;
	for(int i = 1; i <= m; i++) 
		change2(rot[x[i]], 1, n, dfn[y[i]], dfn[y[i]] + siz[y[i]] - 1, node({-siz[x[i]] + siz[y[i]], i}));
	for(int i = n; i >= 1; i--) {
		for(auto p : vec[i])
			res = node(), query2(rot[i], 1, n, dfn[y[p]]), chkmin(p, res, siz[x[p]] - siz[y[p]]);
		if(f[i]) rot[f[i]] = merge(rot[f[i]], rot[i]);
	}
}
void solve8() {//Case8:-sxi+syi+sxj-syj -> xi in xj, yj in yi
	tot = tp = 0;
	for(int i = 1; i <= n; i++) rot[i] = 0;
	for(int i = 1; i <= m; i++) 
		change2(rot[y[i]], 1, n, dfn[x[i]], dfn[x[i]] + siz[x[i]] - 1, node({siz[x[i]] - siz[y[i]], i}));
	for(int i = n; i >= 1; i--) {
		for(auto p : Vec[i])
			res = node(), query2(rot[i], 1, n, dfn[x[p]]), chkmin(p, res, -siz[x[p]] + siz[y[p]]);
		if(f[i]) rot[f[i]] = merge(rot[f[i]], rot[i]);
	}
}
void solve9() {//Case9:-sxi-syi+sxj+syj -> xi in xj, yi in yj
	rt = tot = tp = 0;
	for(int i = 1, u; i <= n + n; i++)
		if(fd[i] < 0) for(; tp > Tp[-fd[i]]; tp--) *st[tp] = Val[tp];
		else {
			Tp[u = fd[i]] = tp;
			for(int j = vec[u].size(), t; j--; ) 
				t = vec[u][j], change2(rt, 1, n, dfn[y[t]], dfn[y[t]] + siz[y[t]] - 1, node(siz[x[t]] + siz[y[t]], t));
			for(int j = vec[u].size(), t; j--; )
				t = vec[u][j], res = node(),query2(rt, 1, n, dfn[y[t]]), chkmin(t, res, -siz[x[t]] - siz[y[t]]);
		}
}
int main() {
#ifdef ddxrS
    freopen("sample.in", "r", stdin);
    freopen("sample.out", "w", stdout);
#endif
    ios::sync_with_stdio(0);							
	cin.tie(0), cout.tie(0);
	cin>>n>>m;
	for(int i = 2; i <= n; i++) cin>>f[i], E[f[i]].push_back(i);
	for(int i = 1; i <= m; i++) cin>>x[i]>>y[i], vec[x[i]].push_back(i), Vec[y[i]].push_back(i), fa[i] = i;
	dfs(1);
	for(int cc = m; cc > 1; ) {
		for(int i = 1; i <= m; i++) pos[i] = 0, w[i] = inf, col[i] = FindSet(i);
		solve1();
		solve2(), solve3(), solve4(), solve5();
		solve6(), solve7(), solve8();                                                                                                                                                         
		solve9(); 
		for(int i = 1; i <= m; i++) if(fa[i] == i && i != FindSet(pos[i])) 
			// cerr<<i<<' '<<pos[i]<<' '<<w[i]<<'\n',
			ans += w[i], fa[FindSet(i)] = FindSet(pos[i]), cc--;
	}
	cout<<ans<<'\n';
	cerr<<clock() * 1.0 / CLOCKS_PER_SEC<<'\n';
	return 0;
}