#include<cstdio>
#include<utility>
#include<vector>
#include<iostream>
using namespace std;
const int o=1e6+10,O=4e6+10;
int n,m,p[o],X[o],Y[o],dfn[o],ed[o],tot,nfd[o*2],h[o],cnt,fa[o],s[o],tp,Tp[o],ptn[o],w[o],col[o];
int rt[o],ls[O],rs[O];long long ans;vector<int> vec[o];
int fr(int x){return fa[x]==x?x:fa[x]=fr(fa[x]);}
struct Edge{int v,p;}e[o];
inline void ad(int U,int V){e[++cnt].v=V;e[cnt].p=h[U];h[U]=cnt;}
struct info{
	pair<int,int> fi,se;
	info(int a=O,int b=0,int c=O,int d=0):fi(make_pair(a,b)),se(make_pair(c,d)){}
	inline void upd(info b){
		if(fi.first>b.fi.first) swap(fi,b.fi);
		if(b.fi.first<=se.first&&col[b.fi.second]-col[fi.second]) se=b.fi;
		if(b.se.first<=se.first&&col[b.se.second]-col[fi.second]) se=b.se;
	}
}tmp,mn[O],Val[O],*st[O];
void dfs(int nw){
	dfn[nw]=++tot;s[nfd[++cnt]=nw]=1;
	for(int i=h[nw];i;i=e[i].p) dfs(e[i].v),s[nw]+=s[e[i].v];
	ed[nw]=tot;nfd[++cnt]=-nw;
}
inline void upd(int x,const info&y,int d){
	x=col[x];
	if(y.fi.first+d<w[x]&&col[y.fi.second]-x) w[x]=y.fi.first+d,ptn[x]=y.fi.second;
	if(y.se.first+d<w[x]&&col[y.se.second]-x) w[x]=y.se.first+d,ptn[x]=y.se.second;
}
void modify1(int&id,int pos,const info&val,int l=1,int r=n){
	if(!id) mn[id=++cnt]=info(),ls[id]=rs[id]=0;
	st[++tp]=mn+id;Val[tp]=*st[tp];mn[id].upd(val);
	if(l==r) return;
	int md=l+r>>1;
	if(pos<=md) modify1(ls[id],pos,val,l,md);
	else modify1(rs[id],pos,val,md+1,r);
}
void query1(int id,int ql,int qr,int l=1,int r=n){
	if(!id) return;
	if(ql<=l&&r<=qr){tmp.upd(mn[id]);return;}
	int md=l+r>>1;
	if(ql<=md) query1(ls[id],ql,qr,l,md);
	if(md<qr) query1(rs[id],ql,qr,md+1,r);
}
void modify2(int&id,int ql,int qr,const info&val,int l=1,int r=n){
	if(!id) mn[id=++cnt]=info(),ls[id]=rs[id]=0;
	if(ql<=l&&r<=qr){st[++tp]=mn+id;Val[tp]=*st[tp];mn[id].upd(val);return;}
	int md=l+r>>1;
	if(ql<=md) modify2(ls[id],ql,qr,val,l,md);
	if(md<qr) modify2(rs[id],ql,qr,val,md+1,r);
}
void query2(int id,int pos,int l=1,int r=n){
	if(!id) return;
	tmp.upd(mn[id]);
	if(l==r) return;
	int md=l+r>>1;
	if(pos<=md) query2(ls[id],pos,l,md);
	else query2(rs[id],pos,md+1,r);
}
int merge(int x,int y){
	if(!x||!y) return x|y;
	mn[x].upd(mn[y]);
	ls[x]=merge(ls[x],ls[y]);rs[x]=merge(rs[x],rs[y]);
	return x;
}
int main(){
	scanf("%d%d",&n,&m);
	for(int i=2;i<=n;++i) scanf("%d",&p[i]),ad(p[i],i);
	for(int i=1;i<=m;++i) scanf("%d%d",&X[i],&Y[i]),vec[X[i]].push_back(i),fa[i]=i;
	cnt=0;dfs(1);
	for(int num=m;num>1;){
		for(int i=1;i<=m;++i) ptn[i]=0,w[i]=O,col[i]=fr(i);
		//case1:sxi+syi+sxj+syj
		tmp=info();
		for(int i=1;i<=m;++i) tmp.upd(info(s[X[i]]+s[Y[i]],i));
		for(int i=1;i<=m;++i) upd(i,tmp,s[X[i]]+s[Y[i]]);
		//case2:sxi+syi+sxj-syj
		rt[0]=cnt=tp=0;
		for(int i=1;i<=m;++i) modify1(rt[0],dfn[Y[i]],info(s[X[i]]-s[Y[i]],i));
		for(int i=1;i<=m;++i) tmp=info(),query1(rt[0],dfn[Y[i]],ed[Y[i]]),upd(i,tmp,s[X[i]]+s[Y[i]]);
		//case3:sxi-syi+sxj+syj
		rt[0]=cnt=tp=0;
		for(int i=1;i<=m;++i) modify2(rt[0],dfn[Y[i]],ed[Y[i]],info(s[X[i]]+s[Y[i]],i));
		for(int i=1;i<=m;++i) tmp=info(),query2(rt[0],dfn[Y[i]]),upd(i,tmp,s[X[i]]-s[Y[i]]);
		//case4:sxi+syi-sxj+syj
		for(int i=1;i<=n;++i) mn[i]=info();
		for(int i=1;i<=m;++i) mn[X[i]].upd(info(-s[X[i]]+s[Y[i]],i));
		for(int i=n;i>1;--i) mn[p[i]].upd(mn[i]);
		for(int i=1;i<=m;++i) upd(i,mn[X[i]],s[X[i]]+s[Y[i]]);
		//case5:-sxi+syi+sxj+syj
		for(int i=1;i<=n;++i) mn[i]=info();
		for(int i=1;i<=m;++i) mn[X[i]].upd(info(s[X[i]]+s[Y[i]],i));
		for(int i=2;i<=n;++i) mn[i].upd(mn[p[i]]);
		for(int i=1;i<=m;++i) upd(i,mn[X[i]],-s[X[i]]+s[Y[i]]);
		//case6:sxi+syi-sxj-syj
		cnt=tp=0;
		for(int i=1;i<=n;++i) rt[i]=0;
		for(int i=1;i<=m;++i) modify1(rt[X[i]],dfn[Y[i]],info(-s[X[i]]-s[Y[i]],i));
		for(int i=n;i;--i){
			for(int j=vec[i].size(),t;j--;)
				t=vec[i][j],tmp=info(),query1(rt[i],dfn[Y[t]],ed[Y[t]]),upd(t,tmp,s[X[t]]+s[Y[t]]);
			if(i>1) rt[p[i]]=merge(rt[p[i]],rt[i]);
		}
		//case7:sxi-syi-sxj+syj
		cnt=tp=0;
		for(int i=1;i<=n;++i) rt[i]=0;
		for(int i=1;i<=m;++i) modify2(rt[X[i]],dfn[Y[i]],ed[Y[i]],info(-s[X[i]]+s[Y[i]],i));
		for(int i=n;i;--i){
			for(int j=vec[i].size(),t;j--;)
				t=vec[i][j],tmp=info(),query2(rt[i],dfn[Y[t]]),upd(t,tmp,s[X[t]]-s[Y[t]]);
			if(i>1) rt[p[i]]=merge(rt[p[i]],rt[i]);
		}
		//case8:-sxi+syi+sxj-syj
		rt[0]=cnt=tp=0;
		for(int i=1,u;i<=n*2;++i)
			if(nfd[i]<0) for(;tp>Tp[-nfd[i]];--tp) *st[tp]=Val[tp];
			else{
				Tp[u=nfd[i]]=tp;
				for(int j=vec[u].size(),t;j--;) t=vec[u][j],modify1(rt[0],dfn[Y[t]],info(s[X[t]]-s[Y[t]],t));
				for(int j=vec[u].size(),t;j--;)
					t=vec[u][j],tmp=info(),query1(rt[0],dfn[Y[t]],ed[Y[t]]),upd(t,tmp,-s[X[t]]+s[Y[t]]);
			}
		//case9:-sxi-syi+sxj+syj
		rt[0]=cnt=tp=0;
		for(int i=1,u;i<=n*2;++i)
			if(nfd[i]<0) for(;tp>Tp[-nfd[i]];--tp) *st[tp]=Val[tp];
			else{
				Tp[u=nfd[i]]=tp;
				for(int j=vec[u].size(),t;j--;) t=vec[u][j],modify2(rt[0],dfn[Y[t]],ed[Y[t]],info(s[X[t]]+s[Y[t]],t));
				for(int j=vec[u].size(),t;j--;)
					t=vec[u][j],tmp=info(),query2(rt[0],dfn[Y[t]]),upd(t,tmp,-s[X[t]]-s[Y[t]]);
			}
		//boruvka
		for(int i=1;i<=m;++i) if(fa[i]==i&&i-fr(ptn[i])) ans+=w[i],fa[fr(i)]=fr(ptn[i]),--num;
	}
	printf("%lld",ans);
	return 0;
}