#include<bits/stdc++.h>
using namespace std;
typedef long long LL;
const int N=1e3+3,M=5e3;
int n,x,a[N],h[N],idx,f[N],s[2][M+3];
struct node{
	int v,ne;
}e[N];
void add(int x,int y){
	e[++idx]={y,h[x]};
	h[x]=idx;
}
void dfs(int u){
	int son=0,sum=0;
	for(int i=h[u];i;i=e[i].ne){
		int v=e[i].v;
		dfs(v);
		sum+=a[v]+f[v];
	}
	for(int i=0;i<=a[u];i++) s[0][i]=s[1][i]=0;
	s[0][0]=1;
	for(int i=h[u];i;i=e[i].ne){
		int v=e[i].v;son++;
		for(int j=a[u];j>=0;j--) 
		    s[son&1][j]=0;
		for(int j=a[u];j>=f[v];j--)
		    s[son&1][j]|=s[(son-1)&1][j-f[v]];
		for(int j=a[u];j>=a[v];j--)
		    s[son&1][j]|=s[(son-1)&1][j-a[v]];
	}
	//cout<<u<<' '<<sum<<endl;
	//for(int i=0;i<=a[u];i++) cout<<s[son&1][i]<<' ';cout<<endl;
	for(int i=a[u];i>=0;i--)
    	if(s[son&1][i]){
    		f[u]=sum-i;
    		//cout<<f[u]<<endl;
    		return;
		}
	puts("IMPOSSIBLE");
	exit(0);	
}
int main()
{
    scanf("%d",&n);
    for(int i=2;i<=n;i++){
    	scanf("%d",&x);
    	add(x,i);
	}
	for(int i=1;i<=n;i++) scanf("%d",&a[i]);
	dfs(1);
	puts("POSSIBLE");
 	return 0;
}