#include <cstdio>
#include <cstring>
#include <iostream>
#include <queue>
#include <vector>
using namespace std;
const int maxn=10010;
int n,m,tot,cnt,S,T,ans;
int fa[maxn],to[maxn*240],next[maxn*240],val[maxn*240],head[maxn*80],d[maxn*80],v[maxn],rt[maxn],siz[maxn],son[maxn];
vector<int> ch[maxn];
struct sag
{
	int ls,rs;
}s[maxn*80];
queue<int> q;
inline int rd()
{
	int ret=0,f=1;	char gc=getchar();
	while(gc<‘0‘||gc>‘9‘)	{if(gc==‘-‘)	f=-f;	gc=getchar();}
	while(gc>=‘0‘&&gc<=‘9‘)	ret=ret*10+(gc^‘0‘),gc=getchar();
	return ret*f;
}
inline void add(int a,int b,int c)
{
	to[cnt]=b,val[cnt]=c,next[cnt]=head[a],head[a]=cnt++;
	to[cnt]=a,val[cnt]=0,next[cnt]=head[b],head[b]=cnt++;
}
void insert(int &x,int y,int l,int r,int a,int b)
{
	x=++tot;
	if(l==r)
	{
		add(x+n,b,1<<30);
		if(y)	add(x+n,y+n,1<<30);
		return ;
	}
	int mid=(l+r)>>1;
	if(a<=mid)	s[x].rs=s[y].rs,insert(s[x].ls,s[y].ls,l,mid,a,b);
	else	s[x].ls=s[y].ls,insert(s[x].rs,s[y].rs,mid+1,r,a,b);
}
void getall(int x,int y)
{
	insert(rt[x],rt[x],1,n,v[y],y);
	for(int i=0;i<(int)ch[y].size();i++)	getall(x,ch[y][i]);
}
void DFS(int x)
{
	siz[x]=1;
	for(int i=0;i<(int)ch[x].size();i++)
	{
		DFS(ch[x][i]),siz[x]+=siz[ch[x][i]];
		if(siz[ch[x][i]]>siz[son[x]])	son[x]=ch[x][i];
	}
	insert(rt[x],rt[son[x]],1,n,v[x],x);
	for(int i=0;i<(int)ch[x].size();i++)	if(ch[x][i]!=son[x])	getall(x,ch[x][i]);
}
void query(int x,int l,int r,int a,int b,int y)
{
	if(!x)	return ;
	if(a<=l&&r<=b)
	{
		add(y+n,x+n,1<<30);
		return ;
	}
	int mid=(l+r)>>1;
	if(a<=mid)	query(s[x].ls,l,mid,a,b,y);
	if(b>mid)	query(s[x].rs,mid+1,r,a,b,y);
}
inline int dfs(int x,int mf)
{
	if(x==T)	return mf;
	int i,k,temp=mf;
	for(i=head[x];i!=-1;i=next[i])	if(d[to[i]]==d[x]+1&&val[i])
	{
		k=dfs(to[i],min(temp,val[i]));
		if(!k)	d[to[i]]=0;
		val[i]-=k,val[i^1]+=k,temp-=k;
		if(!temp)	break;
	}
	return mf-temp;
}
int bfs()
{
	while(!q.empty())	q.pop();
	memset(d,0,sizeof(d));
	int i,u;
	q.push(S),d[S]=1;
	while(!q.empty())
	{
		u=q.front(),q.pop();
		for(i=head[u];i!=-1;i=next[i])	if(!d[to[i]]&&val[i])
		{
			d[to[i]]=d[u]+1;
			if(to[i]==T)	return 1;
			q.push(to[i]);
		}
	}
	return 0;
}
int main()
{
	//freopen("bz3681.out","w",stdout);
	n=rd(),m=rd();
	int i,a,b,x,y;
	memset(head,-1,sizeof(head));
	for(i=2;i<=n;i++)	fa[i]=rd(),ch[fa[i]].push_back(i);
	for(i=1;i<=n;i++)	v[i]=rd();
	DFS(1);
	S=0;
	for(i=1;i<=tot;i++)
	{
		if(s[i].ls)	add(i+n,s[i].ls+n,1<<30);
		if(s[i].rs)	add(i+n,s[i].rs+n,1<<30);
	}
	for(i=1;i<=m;i++)
	{
		a=rd(),b=rd(),x=rd(),y=rd(),add(S,(++tot)+n,y);
		query(rt[x],1,n,a,b,tot);
	}
	T=tot+n+1;
	for(i=1;i<=n;i++)	add(i,T,1);
	while(bfs())	ans+=dfs(S,1<<30);
	printf("%d",ans);
	return 0;
}