#include <iostream>
#include <vector>
#include <bitset>
using namespace std;

int X[5010], Y[5010];
vector<vector<int>> G;
vector<vector<bitset<5010>>> dp;

bool dfs(int key){
    int sum=0;
    for(int i=0; i<G[key].size(); ++i){
        if(!dfs(G[key][i])) return false;
        sum += X[G[key][i]]+Y[G[key][i]];
    }
    dp[key].resize(G[key].size()+1);
    dp[key][0][0]=1;
    for(int i=1; i<5010; ++i) dp[key][0][i]=0;
    for(int i=0; i<G[key].size(); ++i){
        dp[key][i+1]=(dp[key][i]<<(X[G[key][i]]))|(dp[key][i]<<(Y[G[key][i]]));
    }
    for(int i=X[key]; i>=0; --i){
        if(dp[key][G[key].size()][i]){
            Y[key]=sum-i;
            return true;
        }
    }
    return false;
}

int main(void){
    int N;
    cin >> N;
    G.resize(N);
    dp.resize(N);
    for(int i=1; i<N; ++i){
        int P;
        cin >> P;
        G[P-1].emplace_back(i);
    }
    for(int i=0; i<N; ++i) cin >> X[i];
    cout << (dfs(0) ? "" : "IM") << "POSSIBLE" << endl;
    return 0;
}
