#include<iostream>
#include<vector>
#include<set>
using namespace std;
#define N 200010
int n,to[N],m,a[N],b[N],qq[N],ans,rr[N];
int qwq,cid[N],td[N],cnt,dep[N],dfn[N],ed[N],ys[N];
bool vis[N],ins[N],inc[N];
vector<int> g[N],pi[N],pe[N],qi[N],qe[N];
set<int> s;
multiset<int> s2;
struct node
{
	int len,sum,ma;
	friend node operator+(node a,node b)
	{
		return (node){a.len+b.len,a.sum+b.sum,max(a.ma-b.len+b.sum,b.ma)};
	}
};
class sgtree
{
	public:
		#define sN 30000010
		int cnt,rt[N],ls[sN],rs[sN],qwq=600000;
		node a[sN];
		void pushdown(int o,int l,int r)
		{
			if(l==r||ls[o]||rs[o])
			 return;
			int mid=l+r>>1;
			ls[o]=++cnt;
			a[ls[o]]=(node){mid-l+1,0,0};
			rs[o]=++cnt;
			a[rs[o]]=(node){r-mid,0,0};
		}
		void upd(int o,int l,int r)
		{
			int mid=l+r>>1;
			node w1=(node){mid-l+1,0,0};
			node w2=(node){r-mid,0,0};
			if(ls[o])
			 w1=a[ls[o]];
			if(rs[o])
			 w2=a[rs[o]];
			a[o]=w1+w2;
		}
		int update(int o,int l,int r,int x,int k)
		{
			if(l>x||r<x||!o)
			 return o;
			int d=++cnt;
			a[d]=a[o];
			if(l==r)
			{
				a[d].sum+=k;
				a[d].ma=max(a[d].sum,0);
				return d;
			}
			pushdown(o,l,r);
			int mid=l+r>>1;
			ls[d]=update(ls[o],l,mid,x,k);
			rs[d]=update(rs[o],mid+1,r,x,k);
			upd(d,l,r);
			return d;
		}
		int add(int o,int x,int k)
		{
			int tmp=update(o,1,qwq,x,k);
			return tmp;
		}
		int merge(int o1,int o2,int l,int r)
		{
			if(!o1||!o2)
			 return o1+o2;
			if(l==r)
			{
				a[o1].sum+=a[o2].sum;
				a[o1].ma=max(a[o1].sum,0);
				return o1;
			}
			int mid=l+r>>1;
			ls[o1]=merge(ls[o1],ls[o2],l,mid);
			rs[o1]=merge(rs[o1],rs[o2],mid+1,r);
			upd(o1,l,r);
			return o1;
		}
		int Merge(int x,int y)
		{
			return merge(x,y,1,qwq);
		}
		node query(int o,int l,int r,int x,int y)
		{
			if(l>y||r<x)
			 return (node){0,0,0};
			if(!o)
			 return (node){min(y,r)-max(l,x)+1,0,0};
			if(l>=x&&r<=y)
			 return a[o];
			int mid=l+r>>1;
			node q1=query(ls[o],l,mid,x,y);
			node q2=query(rs[o],mid+1,r,x,y);
			return q1+q2;
		}
		int ask(int rt)
		{
//			cerr<<"UUU:"<<a[rt].sum<<endl;
			int l=0,r=qwq,ans=qwq;
			while(l<=r)
			{
				int mid=l+r>>1;
				node res=query(rt,1,qwq,1,mid);
//				cerr<<"umnik"<<mid<<"::"<<res.len<<" "<<res.sum<<" "<<res.ma<<endl;
				if(res.sum==a[rt].sum&&res.ma<=0)
				{
					ans=mid;
					r=mid-1;
				}
				else
				{
					l=mid+1;
				}
			}
			return ans;
		}
		int init()
		{
			int u=++cnt;
			ls[u]=0;
			rs[u]=0;
			a[u]=(node){qwq,0,0};
			return u;
		}
}tr;
void dfs0(int u)
{
	vis[u]=1;
	ins[u]=1;
	if(ins[to[u]])
	 inc[to[u]]=1;
	if(vis[to[u]])
	{
		ins[u]=0;
		return;
	}
	dfs0(to[u]);
	ins[u]=0;
}
void dfs1(int u)
{
	dfn[u]=++cnt;
	for(int v:g[u])
	{
		cid[v]=cid[u];
		dep[v]=dep[u]+1;
		td[v]=td[u];
		dfs1(v);
	}
	ed[u]=cnt;
}
void init()
{
	for(int i=1;i<=n;i++)
	if(!vis[i])
	 dfs0(i);
	for(int i=1;i<=n;i++)
	if(inc[i]&&!inc[to[i]])
	{
		cid[i]=++qwq;
		int now=to[i];
		while(!inc[now])
		{
			inc[now]=1;
			cid[now]=cid[i];
			now=to[now];
		}
	}
	for(int i=1;i<=n;i++)
	if(!inc[i])
	 g[to[i]].push_back(i);
	for(int i=1;i<=n;i++)
	if(inc[i])
	{
		td[i]=i;
		dfs1(i);
	}
}
bool check(int x,int y)
{
	if(cid[x]!=cid[y])
	 return 0;
	if(!inc[y]&&!(dfn[x]>=dfn[y]&&dfn[x]<=ed[y]))
	 return 0;
	return 1;
}
void dfs2(int u)
{
	if(!g[u].size())
	 tr.rt[u]=tr.init();
	for(int v:g[u])
	{
		dfs2(v);
		tr.rt[u]=tr.Merge(tr.rt[u],tr.rt[v]);
	}
	for(int x:pi[u])
	 tr.rt[u]=tr.add(tr.rt[u],dep[a[x]],1);
	for(int x:pe[u])
	 tr.rt[u]=tr.add(tr.rt[u],dep[a[x]],-1);
	int w=tr.ask(tr.rt[u]);
	ans=max(ans,w-dep[u]);
}
void solvec(int bbb)
{
	if(to[bbb]==bbb)
	 return;
	vector<int> d;
	while(!vis[bbb])
	{
		vis[bbb]=1;
		d.push_back(bbb);
		bbb=to[bbb];
	}
	s.clear();
//	tr.cnt=0;
	tr.rt[0]=tr.init();
	for(int i=0;i<d.size();i++)
	 ys[d[i]]=d.size()-i;
	for(int u:d)
	{
		for(int x:qi[u])
		{
//			cerr<<"insA:"<<x<<"|"<<a[x]<<" "<<b[x]<<"||"<<ys[td[a[x]]]<<"+"<<qq[x]<<"+"<<d.size()<<endl;
			s.insert(x);
			rr[x]=1;
			tr.rt[0]=tr.add(tr.rt[0],ys[td[a[x]]]+qq[x]+d.size(),1);
		}
		for(int x:qe[u])
		{
			if(s.find(x)!=s.end())
			{
				s.erase(x);
//				cerr<<"ersA:"<<x<<endl;
				tr.rt[0]=tr.add(tr.rt[0],ys[td[a[x]]]+qq[x]+d.size(),-1);
			}
		}
	}
	for(int u:d)
	{
		for(int x:qi[u])
		{
//			cerr<<"insB:"<<x<<"|"<<a[x]<<" "<<b[x]<<"||"<<ys[td[a[x]]]<<"+"<<qq[x]<<endl;
			s.insert(x);
			rr[x]=0;
			tr.rt[0]=tr.add(tr.rt[0],ys[td[a[x]]]+qq[x],1);
		}
		for(int x:qe[u])
		{
			if(s.find(x)!=s.end())
			{
				s.erase(x);
//				cerr<<"ersB:"<<x<<"||"<<rr[x]<<endl;
				if(rr[x])
				 tr.rt[0]=tr.add(tr.rt[0],ys[td[a[x]]]+qq[x]+d.size(),-1);
				else
				 tr.rt[0]=tr.add(tr.rt[0],ys[td[a[x]]]+qq[x],-1);
			}
		}
//		cerr<<"set:-----"<<endl;
//		for(int ww:s)
//		 cerr<<ww<<" ";
//		cerr<<endl;
//		if(u==d.back())
//		 continue;
//		cerr<<"#2 "<<u<<":"<<tr.ask(tr.rt[0])<<" - "<<ys[u]<<endl;
		ans=max(ans,tr.ask(tr.rt[0])-ys[u]);
	}
	return;
	s2.clear();
	for(int x:s)
	 s2.insert(ys[td[a[x]]]-1+qq[x]);
	int las=0,sum=0;
	for(int x:s2)
	{
//		cerr<<"!!!"<<x<<endl;
		sum=max(0,sum-(x-las)+1);
		las=x;
	}
//	cerr<<"#3 >>"<<d.back()<<":"<<las<<"+"<<sum<<endl;
	ans=max(ans,las+sum);
}
int main()
{
	ios::sync_with_stdio(false);
	cin.tie(0);
	cout.tie(0);
	cin>>n;
	for(int i=1;i<=n;i++)
	 cin>>to[i];
	init();
	cin>>m;
	for(int i=1;i<=m;i++)
	{
		cin>>a[i]>>b[i];
		if(!check(a[i],b[i]))
		{
			cout<<"-1"<<endl;
			return 0;
		}
		if(!inc[a[i]])
		 pi[a[i]].push_back(i);
		if(!inc[b[i]])
		 pe[b[i]].push_back(i);
		else if(!inc[a[i]])
		 pe[td[a[i]]].push_back(i);
		if(inc[b[i]])
		{
//			cerr<<i<<"-->"<<td[a[i]]<<endl;
			qq[i]=dep[a[i]];
			qi[td[a[i]]].push_back(i);
			qe[b[i]].push_back(i);
		}
	}
	for(int i=1;i<=n;i++)
	if(inc[i])
	 dfs2(i);
	for(int i=1;i<=n;i++)
	 vis[i]=0;
	for(int i=1;i<=n;i++)
	if(inc[i]&&!vis[i])
	{
		solvec(i);
	}
	cout<<ans<<endl;
}