#include<iostream>
#include<vector>
#include<algorithm>
#define BUFSIZE 10000000
struct read{
	char buf[BUFSIZE],*p1,*p2,c;
	read():p1(buf),p2(buf){}
	char gc(void){
	    return p1==p2?(std::cin.read(p1=buf,BUFSIZE),p2=p1+std::cin.gcount(),p1==p2?EOF:*p1++):*p1++;
	}
	read& operator >>(int& x){
		c=gc(),x=0;
		for(;c<'0'||c>'9';c=gc());
		for(;c>='0'&&c<='9';c=gc())x=x*10+(c-'0');
		return *this;
	}
}cin;
using std::cout;
std::vector<int> v[1000010],x,y;
std::vector<std::pair<int,int>> t[1000010],h[1000010];
std::vector<std::pair<int,int>>*g;
int n,m,a[100010],b[100010],sz[1000010],xl[1000010],xr[1000010],yl[1000010],yr[1000010];
int f[100010],pcnt,col[100010],cs,*L,*R;
int find(int x){
	while(x!=f[x]) x=f[x]=f[f[x]];
	return x;
}
int rk[1000010];
int p[2000010],cct;
void dfs(int x){
	static int cnt=0;
	rk[x]=++cnt,sz[x]=1,p[++cct]=x;
	for(auto u:v[x]) dfs(u),sz[x]+=sz[u];
	p[++cct]=-x;
}
int const inf=1e9;
struct pnt{
	int mn,mnt;
	pnt():mn(inf),mnt(){}
	pnt(int x,int y):mn(x),mnt(y){}
	friend bool operator <(pnt const &a,pnt const &b){
		return a.mn<b.mn;
	}
	pnt operator -(int x)const{
		return pnt(mn-x,mnt);
	}
	pnt operator +(int x)const{
		return pnt(mn+x,mnt);
	}
}res[100010];
struct node{
	pnt m,xm;
	node():m(),xm(){}
	node(pnt const &x):m(x),xm(){}
	node(int x,int y):m(x,y),xm(){}
	node(pnt const &x,pnt const &y):m(x),xm(y){}
	void operator +=(node const &y){
		if(y.m<m) xm=std::min(y.xm,col[m.mnt]==col[y.m.mnt]?xm:m),m=y.m;
		else xm=std::min(xm,col[y.m.mnt]==col[m.mnt]?y.xm:y.m);
	} 
	pnt query(int x)const{
		if(col[m.mnt]==col[x]) return xm;
		else return m;
	}
};
namespace sub1{
	node f[1000010],r[1000010];
	void work(){
		f[1]=node();
		for(int i=1;i<=cct;i++){
			int x=p[i];
			if(x<0){
				x=-x;
				for(auto u:v[x]) r[x]+=r[u];
				for(auto [i,u]:g[x]) res[i]=std::min(res[i],std::min(r[x].query(i)+sz[x]+sz[u],f[x].query(i)-sz[x]+sz[u]));
			}else{
				r[x]=node();
				for(auto [i,u]:g[x]) f[x]+=node(sz[x]+sz[u],i),r[x]+=node(-sz[x]+sz[u],i);
				for(auto u:v[x]) f[u]=f[x];
			}
		}
	}
}
node tr[1800010];
namespace sub2{
	int cnt,ls[1800010],rs[1800010],rt[1000010];
	void update(int p,node const &v,int &x,int l=1,int r=cs){
		if(!x) x=++cnt,ls[x]=rs[x]=0;
		tr[x]+=v;
		if(l==r) return;
		int mid=(l+r)>>1;
		if(p<=mid) update(p,v,ls[x],l,mid);
		else update(p,v,rs[x],mid+1,r);
	}
	void merge(int &x,int y){
		if(!x||!y) return x|=y,void();
		tr[x]+=tr[y];
		merge(ls[x],ls[y]),merge(rs[x],rs[y]);
	}
	pnt query(int pl,int pr,int i,int x,int l=1,int r=cs){
		if(l==pl&&r==pr) return tr[x].query(i);
		int mid=(l+r)>>1;
		if(pr<=mid) return query(pl,pr,i,ls[x],l,mid);
		else if(pl>mid) return query(pl,pr,i,rs[x],mid+1,r);
		else return std::min(query(pl,mid,i,ls[x],l,mid),query(mid+1,pr,i,rs[x],mid+1,r));
	}
	void work(){
		for(int i=1;i<=cct;i++){
			int x=p[i];
			if(x<0){
				x=-x;
				rt[x]=0;	
				for(auto u:v[x]) merge(rt[x],rt[u]);
				for(auto [i,u]:g[x]) update(L[u],node(-sz[x]-sz[u],i),rt[x]);
				for(auto [i,u]:g[x]) res[i]=std::min(res[i],query(L[u],R[u],i,rt[x])+sz[x]+sz[u]);
			}
		}
		for(int i=1;i<=cnt;i++) tr[i]=node();
		cnt=0;
	}
}
int M;
int tp,s[3600010],top[1000010];
node val[3600010];
namespace sub3{
	void update(int p,node const &v){
		for(p+=M;p;p>>=1) ++tp,s[tp]=p,val[tp]=tr[p],tr[p]+=v;
	}
	pnt query(int x,int y,int i){
		pnt ans;
		for(x+=M-1,y+=M+1;x^y^1;x>>=1,y>>=1){
			if(~x&1) ans=std::min(ans,tr[x^1].query(i));
			if(y&1) ans=std::min(ans,tr[y^1].query(i));
		}
		return ans;
	}
	void work(){
		M=1;
		while(M<cs) M<<=1;
		for(int i=1;i<=cct;i++){
			int x=p[i];
			if(x<0){
				x=-x;
				while(tp!=top[x]) tr[s[tp]]=val[tp],--tp;
			}else{
				top[x]=tp;
				for(auto [i,u]:g[x]) update(L[u],node(sz[x]-sz[u],i));
				for(auto [i,u]:g[x]) res[i]=std::min(res[i],query(L[u],R[u],i)-sz[x]+sz[u]);
			}
		}
	}
}
namespace sub4{
	void update(int x,int y,node const &v){
		for(x+=M-1,y+=M+1;x^y^1;x>>=1,y>>=1){
			if(~x&1) ++tp,s[tp]=x^1,val[tp]=tr[x^1],tr[x^1]+=v;
			if(y&1) ++tp,s[tp]=y^1,val[tp]=tr[y^1],tr[y^1]+=v;
		}
	}
	pnt query(int p,int i){
		pnt ans;
		for(p+=M;p;p>>=1) ans=std::min(ans,tr[p].query(i));
		return ans;
	}
	void work(){
		M=1;
		while(M<cs) M<<=1;
		for(int i=1;i<=cct;i++){
			int x=p[i];
			if(x<0){
				x=-x;
				while(tp!=top[x]) tr[s[tp]]=val[tp],--tp;
			}else{
				top[x]=tp;
				for(auto [i,u]:g[x]) update(L[u],R[u],node(sz[x]+sz[u],i));
				for(auto [i,u]:g[x]) res[i]=std::min(res[i],query(L[u],i)-sz[x]-sz[u]);
			}
		}
	}
}
int main(){
	std::ios::sync_with_stdio(false),std::cin.tie(nullptr);
	cin>>n>>m;
	for(int i=1,x;i<n;i++) cin>>x,v[x].push_back(i+1);
	dfs(1);
	for(int i=1;i<=m;i++){
		cin>>a[i]>>b[i],t[a[i]].emplace_back(i,b[i]),h[b[i]].emplace_back(i,a[i]),f[i]=i;
		x.push_back(rk[b[i]]),y.push_back(rk[a[i]]);
	}
	std::sort(x.begin(),x.end());
	std::sort(y.begin(),y.end());
	x.resize(std::unique(x.begin(),x.end())-x.begin());
	y.resize(std::unique(y.begin(),y.end())-y.begin());
	for(int i=1;i<=m;i++)
		xl[b[i]]=std::lower_bound(x.begin(),x.end(),rk[b[i]])-x.begin()+1,xr[b[i]]=std::upper_bound(x.begin(),x.end(),rk[b[i]]+sz[b[i]]-1)-x.begin(),
		yl[a[i]]=std::lower_bound(y.begin(),y.end(),rk[a[i]])-y.begin()+1,yr[a[i]]=std::upper_bound(y.begin(),y.end(),rk[a[i]]+sz[a[i]]-1)-y.begin();
	pcnt=m;
	long long ans=0;
	while(pcnt>1){
		for(int i=1;i<=m;i++) col[i]=find(i);
		node tt;
		for(int i=1;i<=m;i++) tt+=node(sz[a[i]]+sz[b[i]],i);
		for(int i=1;i<=m;i++) res[i]=tt.query(i)+sz[a[i]]+sz[b[i]];
		g=t,cs=x.size(),L=xl,R=xr;
		sub1::work();
		sub2::work();
		sub3::work();
		sub4::work();
		g=h,cs=y.size(),L=yl,R=yr;
		sub1::work();
		sub3::work();
		static pnt link[100010];
		static int c[100010];
		for(int i=1;i<=m;i++) link[i]=pnt(),c[i]=0;
		for(int i=1;i<=m;i++) if(res[i]<link[col[i]]) link[col[i]]=res[i],c[col[i]]=i;
		for(int i=1;i<=m;i++) if(col[i]==i){
			int u=find(link[i].mnt),v=find(c[i]);
			if(u!=v) f[u]=v,ans+=link[i].mn,--pcnt;
		}
	}
	cout<<ans<<'\n';
	return 0;
}