#include<iostream>
#include<algorithm>
#include<vector>
#include<string>
#include<set>
#include<queue>
#include<stack>
#include<bitset>
#include<functional>
using namespace std;
int p = 1000000007;
#define int long long
#define vel vector<long long>
#define vvel vector<vel>
#define rep(i,n) for(long long i=0;i<n;i++)
#define sor(v) sort(v.begin(),v.end())
#define mmax(a,b) a=max(a,b)
#define mmin(a,b) a=min(a,b)
#define mkp make_pair
#define pin pair<int,int>
#define V vector
#define Endl endl
#define veb vector<bool>
#define sq(a) (a)*(a)
#define rev(s) reverse(s.begin(),s.end())
#define end_pr(s) cout << s <<endl;return 0
#define itn int
#define bs bitset<5001>
int kai_size = 1000001;
vel kai(kai_size, 1);
vel ink(kai_size, 1);
vel dist;
int RE() {
	vel v(3, 2);
	return v.at(4); 
}
int ru(int a, int r) {
	if (r == 0) { return 1; }
	int ans = ru(a, r / 2);
	ans *= ans; ans %= p;
	if (r % 2 == 1) { ans *= a; }
	return ans % p;
}
int inv(int a) {
	return ru(a, p - 2);
}
void make_kai() {
	rep(i, kai_size-1) { kai[i + 1] = (kai[i] * (i + 1)) % p; }
	rep(i, kai_size) { ink[i] = inv(kai[i]); }
}
int com(int n, int r) {
	if (r < 0 || n < r) { return 0; }
	int ans = kai[n] * ink[r];
	ans %= p; ans *= ink[n - r]; ans %= p;
	return ans;
}
vel dis(int mid1, vvel &way) {
	int n = way.size();
	vel dist(n, -1); dist[mid1] = 0;
	queue<int> q;
	q.push(mid1);
	while (!q.empty()) {
		int st = q.front(); q.pop();
		rep(i, way[st].size()) {
			int to = way[st][i];
			if (dist[to] == -1) {
				dist[to] = dist[st] + 1;
				q.push(to);
			}
		}
	}
	return dist;
}
pin most_far(int now, int n, vvel &way) {
	vel dist1 = dis(now, way);
	pin ans = mkp(-1, 0);
	rep(i, n) {
		if (dist1[i] > ans.first) { ans = mkp(dist1[i], i); }
	}
	return ans;
}

int per(int a, int b) {
	int ans = a % b;
	if (ans < 0) { ans += b; }
	return ans;
}
vel uni(vel v) {
	sor(v);
	vel ans(1, v[0]);
	for (int i = 1; i < v.size(); i++) {
		if (v[i] != v[i-1]) { ans.push_back(v[i]); }
	}
	return ans;
}
int s_gcd(int a, int b) {
	if (b == 0) { return a; }
	return s_gcd(b, a%b);
}
int gcd(int a, int b) {
	if (a < b) { swap(a, b); }
	return s_gcd(a, b);
}
bool is_prime(int i) {
	for (int j = 2; j*j <= i; j++) {
		if (i%j == 0) { return false; }
	}
	return true;
}
int get(int i, int j, vvel &ans) {
	int n = ans.size();
	if (i < 0 || n <= i || j < 0 || n <= j) { return 1; }
	return ans[i][j];
}
int lcm(int a, int b) {
	return a * b / gcd(a, b);
}
void del(vel &v, int a) {
	vel ans;
	for (auto x : v) {
		if (x != a) { ans.push_back(x); }
	}
	v = ans;
}
void sz_tree(int ce0, V<V<pin>> &chi,vel &sz) {
	for (auto x : chi[ce0]) {
		sz_tree(x.first, chi, sz);
		sz[ce0] += sz[x.first];
	}
}
int ord(int n) {
	if (n % 2 == 1) { return 0; }
	return 1 + ord(n / 2);
}
int max1(vel &a, int lim) {
	bs smt;
	smt[0] = 1;
	for (auto x : a) {
		smt |= (smt << x);
	}
	int ans = 0;
	for (int i = 0; i <= lim; i++) {
		if (smt[i]) { ans = i; }
	}
	return ans;
}
bool solve(int st, vvel &chi, vel &x, vel &y) {
	int sum_min = 0;
	vel diff(chi[st].size()); int co = 0;
	int all_sum = 0;
	for (auto to : chi[st]) {
		if (!solve(to, chi, x, y)) { return false; }
		if (x[to] < y[to]) { swap(x[to], y[to]); }
		sum_min += y[to];
		diff[co] = x[to] - y[to];
		co++;
		all_sum += x[to] + y[to];
	}
	if (x[st] < sum_min) { return false; }
	y[st] = all_sum - (max1(diff, x[st] - sum_min)+sum_min);
	return true;
}
signed main() {
	int n; cin >> n;
	vvel chi(n);
	rep(i, n - 1) {
		int p; cin >> p;
		chi[p - 1].push_back(i + 1);
	}
	vel x(n);
	rep(i, n) { cin >> x[i]; }
	vel y(n);
	if (solve(0, chi, x, y)) { end_pr("POSSIBLE"); }
	cout << "IMPOSSIBLE" << Endl;
	return 0;
}

Measured memory_bytes: 16000