#include<bits/stdc++.h>
using namespace std;
#define max_n 400020
int n;
int nn;
int q;
struct Q{int u,v;}que[max_n+2];
int eto[max_n+2];
int der[max_n+2];
bool cir[max_n+2];
vector<int>g[max_n+2];
int fa[max_n+2];
int de[max_n+2];
int dfn[max_n+2],tim=0;
int edfn[max_n+2];
vector<pair<int,int>>op[max_n+2];
struct S{int sum,ma,lma,rma,sma;};
S operator+(S a,S b){
	S res;
	res.sum=a.sum+b.sum; res.ma=max(a.ma,b.ma);
	res.lma=max(a.lma,a.sum+b.lma); res.rma=max(b.rma,b.sum+a.rma);
	res.sma=max(max(a.sma,b.sma),a.rma+b.lma);
	return res; 
}
struct T{
	int ls,rs;
	S s;
}tr[max_n*52+20];
int ti=0,rt[max_n+2];
int ans=0;
inline int read(){
	int res=0;
	char c=getchar();
	while(c<'0'||c>'9')c=getchar();
	while(c>='0'&&c<='9'){res=res*10+(c^48); c=getchar();}
	return res;
}
inline void topo(){
	memset(cir,1,sizeof(cir));
	queue<int>q;
	for(int i=1;i<=n;++i)if(!der[i])q.push(i);
	while(!q.empty()){
		int u=q.front(); q.pop(); cir[u]=false;
		if(!--der[eto[u]])q.push(eto[u]);
	}
}
inline void dfs1(int u,int f){
	fa[u]=f; de[u]=de[f]+1; dfn[u]=++tim;
	for(auto v:g[u])dfs1(v,u);
	edfn[u]=tim;
}
inline bool checkin(int u,int v){return dfn[v]<=dfn[u]&&dfn[u]<=edfn[v];}//u in v
inline void pu(int o,int l,int r){
	int mid=(l+r)>>1;
	S sl,sr;
	if(tr[o].ls)sl=tr[tr[o].ls].s; else sl=(S){-(mid-l+1),-1,0,0,0};
	if(tr[o].rs)sr=tr[tr[o].rs].s; else sr=(S){-(r-mid),-1,0,0,0};
	tr[o].s=sl+sr;
}
inline void upd(int &o,int l,int r,int x,int w){
	if(!o){o=++ti; tr[o].s=(S){-(r-l+1),-1,0,0,0};}
	if(l>=r){
		tr[o].s.sum+=w; tr[o].s.ma+=w;
		tr[o].s.lma=tr[o].s.rma=tr[o].s.sma=max(tr[o].s.ma,0);
		return;
	}
	int mid=(l+r)>>1;
	if(x<=mid)upd(tr[o].ls,l,mid,x,w);
	else upd(tr[o].rs,mid+1,r,x,w);
	pu(o,l,r);
}
inline int merge(int o1,int o2,int l,int r){
	if(!o1||!o2)return o1+o2;
	int o=++ti;
	if(l>=r){
		tr[o].s.sum=tr[o1].s.sum+tr[o2].s.sum+1;
		tr[o].s.ma=tr[o1].s.ma+tr[o2].s.ma+1;
		tr[o].s.lma=tr[o].s.rma=tr[o].s.sma=max(tr[o].s.ma,0);
		return o;
	}
	int mid=(l+r)>>1;
	tr[o].ls=merge(tr[o1].ls,tr[o2].ls,l,mid);
	tr[o].rs=merge(tr[o1].rs,tr[o2].rs,mid+1,r);
	pu(o,l,r); return o;
}
S query(int o,int l,int r,int ql,int qr){
	if(ql<=l&&r<=qr)return o?tr[o].s:((S){-(r-l+1),-1,0,0,0});
	int mid=(l+r)>>1;
	if(qr<=mid)return query(tr[o].ls,l,mid,ql,qr);
	if(ql>mid)return query(tr[o].rs,mid+1,r,ql,qr);
	return query(tr[o].ls,l,mid,ql,qr)+query(tr[o].rs,mid+1,r,ql,qr);
}
int queryr(int o,int l,int r){
	if(!o||tr[o].s.ma<=-1)return 0;
	if(l>=r)return l;
	int mid=(l+r)>>1;
	int res=queryr(tr[o].rs,mid+1,r);
	if(!res)res=queryr(tr[o].ls,l,mid);
	return res;
}
inline void dfs2(int u){
	for(auto v:g[u]){dfs2(v); rt[u]=merge(rt[u],rt[v],1,n);}
	for(auto p:op[u])upd(rt[u],1,n,p.first,p.second);
	if(query(rt[u],1,n,1,n).ma<=-1)return;
	int pos=queryr(rt[u],1,n);
	S s=query(rt[u],1,n,1,pos);
	ans=max(ans,s.rma+(pos-de[u]+1));
}
int main(){
	#ifdef dzy
	freopen("package.in","r",stdin);
	freopen("package.out","w",stdout);
	#endif
	n=read();
	for(int i=1;i<=n;++i){eto[i]=read(); ++der[eto[i]];}
	q=read();
	for(int i=1;i<=q;++i){que[i].u=read(); que[i].v=read();}
	topo();
	for(int i=1;i<=n;++i)eto[i+n]=eto[i]+n;
	for(int i=1;i<=n;++i){
		if(!cir[i])continue;
		for(int u=i;;u=eto[u]){
			if(!cir[u])break;
			cir[u]=false;
		}
		eto[i]+=n; eto[i+n]=n*2+1;
	}
	nn=n; n=n*2+1;
	for(int i=1;i<n;++i)g[eto[i]].push_back(i);
	dfs1(n,0);
	for(int i=1,laq=q;i<=laq;++i){
		if(!checkin(que[i].u,que[i].v)){
			if(!checkin(que[i].u,que[i].v+nn)){printf("-1\n"); return 0;}
			que[i].v+=nn;
			que[++q]=(Q){que[i].u+nn,n};
		}
		else que[++q]=(Q){que[i].u+nn,que[i].v+nn};
	}
	for(int i=1;i<=q;++i){
		op[que[i].u].push_back({de[que[i].u],1});
		op[que[i].v].push_back({de[que[i].u],-1});
	}
	tr[0].s=(S){-1,0,0}; dfs2(n);
	printf("%d\n",ans);
	return 0;
}