#include<bits/stdc++.h>
using namespace std;
const int N=1010,M=5050;
int i,j,k,n,m,En;
int h[N],fa[N],f[N],X[N];
bitset<M>g;
struct edge
{
    int s,n;
} E[N];
void E_add(int x,int y)
{
    E[++En].s=y;
    E[En].n=h[x];
    h[x]=En;
}
bool dfs(int x)
{
    for(int k=h[x]; k; k=E[k].n)
        if(!dfs(E[k].s))
            return 0;
    g.reset();
    g[0]=1;
    int sum=0;
    for(int k=h[x]; k; k=E[k].n)
    {
        g=g<<X[E[k].s]|g<<f[E[k].s];
        sum+=X[E[k].s]+f[E[k].s];
    }
    for(int i=X[x]; i>=0; i--)
        if(g[i])
        {
            f[x]=sum-i;
            return 1;
        }
    return 0;
}
int main()
{
    scanf("%d",&n);
    for(i=2; i<=n; i++)
        scanf("%d",&fa[i]),E_add(fa[i],i);
    for(i=1; i<=n; i++)
        scanf("%d",&X[i]);
    puts(dfs(1)?"POSSIBLE":"IMPOSSIBLE");
    return 0;
}
