#include<bits/stdc++.h>
using namespace std;
typedef long long int uli;
const int mx=1e3+10;
const uli oo=1e18;
vector<int>g[mx];
uli x[mx];
uli f[2][5000+10];
uli dp[mx];
bool ok=true;
void dfs(int u){
  uli sm=0;
  for(int v:g[u]){
    dfs(v);
    sm+=dp[v]+x[v];
  }
  if(g[u].empty()){
    dp[u]=0;
    return;
  }
  int b=0;
  for(int s=0;s<=x[u];s++)f[b][s]=sm-s;
  for(int v:g[u]){
    b^=1;
    for(int s=0;s<=x[u];s++){
      uli p=dp[v],q=x[v];
      f[b][s]=oo;
      for(int it=0;it<2;it++){
        if(p+s<=x[u]){
          f[b][s]=min(f[b][s],f[b^1][s+p]);
        }
        swap(p,q);
      }
    }
  }
  dp[u]=f[b][0];
  if(dp[u]>=oo)ok=false;
}
int main(){
  int n;
  scanf("%d",&n);
  for(int i=1;i<n;i++){
    int pi;
    scanf("%d",&pi);
    --pi;
    g[pi].push_back(i);
  }
  for(int i=0;i<n;i++)scanf("%lld",x+i);
  dfs(0);
  if(ok)puts("POSSIBLE");
  else puts("IMPOSSIBLE");
  return 0;
}
