#include "tree.h"
#include <set>
#include <cstdio>
#include <cassert>
#include <algorithm>
using namespace std;
typedef long long ll;
typedef vector<int> vi;
const int N=200003;
const int INF=0x3f3f3f3f;
typedef pair<int,ll> pil;
int n,L,R;
vi vec[N];
multiset<pil> st[N];
vi w;ll res;
ll dfs(int u){
	st[u].clear();
	if(vec[u].empty()){
		res+=(ll)L*w[u];
		st[u].emplace(INF,L);
		return L;
	}
	else{
		ll sm=0,all=0;
		for(int v:vec[u]){
			sm+=dfs(v);
			if(st[v].size()>st[u].size()) st[u].swap(st[v]);
			for(auto [x,c]:st[v]) st[u].emplace(x,c);
		}
		while(!st[u].empty()){
			auto it=prev(st[u].end());
			if((it->first)>w[u]){
				all+=it->second;
				st[u].erase(it);
			}
			else break;
		}
		st[u].emplace(w[u],all);
		while(!st[u].empty()&&sm>R){
			auto it=st[u].begin();
			auto [x,c]=*it;
			st[u].erase(it);
			if(c>sm-R){
				c-=sm-R;sm=R;
				res+=(ll)(sm-R)*x;
				st[u].emplace(x,c);
			}
			else sm-=c,res+=(ll)c*x;
		}
		/*
		all=L;
		while(!st[u].empty()&&all){
			auto it=prev(st[u].end());
			auto [x,c]=*it;
			st[u].erase(it);
			if(all<c){
				c-=all;all=0;
				st[u].emplace(x,c);
				break;
			}
			else all-=c;
		}
		st[u].emplace(INF,L);
		*/
		return sm;
	}
}
void init(vi P,vi W){
	n=P.size();w=W;
	for(int i=1;i<n;++i) vec[P[i]].emplace_back(i);
}
ll query(int _L,int _R){
	L=_L;R=_R;res=0;dfs(0);
	return res;
}