結果

問題 No.1524 Upward Mobility
コンテスト
ユーザー 沙耶花
提出日時 2021-05-05 13:41:44
言語 C++17
(gcc 13.3.0 + boost 1.89.0)
結果
MLE  
実行時間 -
コード長 1,930 bytes
記録
記録タグの例:
初AC ショートコード 純ショートコード 純主流ショートコード 最速実行時間
コンパイル時間 4,769 ms
コンパイル使用メモリ 262,152 KB
最終ジャッジ日時 2025-01-21 07:19:27
ジャッジサーバーID
（参考情報） 		judge4 / judge5
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 24 TLE * 3 MLE * 4
権限があれば一括ダウンロードができます
コンパイルメッセージ
main.cpp: In function ‘int main()’:
main.cpp:86:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   86 |                 scanf("%d",&p);
      |                 ~~~~~^~~~~~~~~
main.cpp:95:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   95 |                 scanf("%d",&a[i]);
      |                 ~~~~~^~~~~~~~~~~~
main.cpp:98:22: warning: ignoring return value of ‘int scanf(const char*, ...)’ declared with attribute ‘warn_unused_result’ [-Wunused-result]
   98 |         rep(i,N)scanf("%lld",&b[i]);
      |                 ~~~~~^~~~~~~~~~~~~~

ソースコード

diff #
raw source code 

#include <stdio.h>
#include <bits/stdc++.h>
#include <atcoder/all>
using namespace atcoder;
using mint = modint998244353;
using namespace std;
#define rep(i,n) for (int i = 0; i < (n); ++i)
#define Inf 1000000000

int N;
vector<vector<int>> E;
vector<int> a;
vector<long long> b;
vector<vector<int>> ps;

vector<vector<long long>> v;

void dfs(int cur){
	if(E[cur].size()==0){
		ps[cur] = {a[cur],N};
		v[cur] = {b[cur],0};
		return;
	}
	
	vector<int> pp(1,a[cur]);
	
	rep(i,E[cur].size()){
		int to = E[cur][i];
		dfs(to);
		rep(j,ps[to].size())pp.push_back(ps[to][j]);
	}
	sort(pp.begin(),pp.end());
	pp.erase(unique(pp.begin(),pp.end()),pp.end());
	
	v[cur].resize(pp.size());
	
	rep(i,E[cur].size()){
		int to = E[cur][i];
		rep(j,ps[to].size()){
			int d0 = -1;
			if(j!=0)d0 = distance(pp.begin(),lower_bound(pp.begin(),pp.end(),ps[to][j-1]));
			int d1 = distance(pp.begin(),lower_bound(pp.begin(),pp.end(),ps[to][j]));
			v[cur][d1] += v[to][j];
			if(d0!=-1)v[cur][d0] -= v[to][j];
		}
	}
	
	for(int i=v[cur].size()-1;i>=1;i--)v[cur][i-1] += v[cur][i];
	for(int i=v[cur].size()-1;i>=1;i--)v[cur][i-1] = max(v[cur][i-1],v[cur][i]);
	
	int d = distance(pp.begin(),lower_bound(pp.begin(),pp.end(),a[cur]));
	v[cur][d] += b[cur];
	for(int i=v[cur].size()-1;i>=1;i--)v[cur][i-1] = max(v[cur][i-1],v[cur][i]);
		
	
	rep(i,E[cur].size()){
		int to = E[cur][i];
		ps[to].clear();
		v[to].clear();
	}
	
	{
		vector<long long> vv;
		rep(i,pp.size()){
			if(vv.size()>0 && vv.back()==v[cur][i]){
				ps[cur].back() = pp[i];
			}
			else{
				vv.push_back(v[cur][i]);
				ps[cur].push_back(pp[i]);
			}
		}
		v[cur] = vv;
	}
	
}

int main(){
	
	cin>>N;
	
	E.resize(N);
	
	rep(i,N-1){
		int p;
		scanf("%d",&p);
		p--;
		E[p].push_back(i+1);
	}
	
	a.resize(N);
	b.resize(N);
	
	rep(i,N){
		scanf("%d",&a[i]);
		a[i]--;
	}
	rep(i,N)scanf("%lld",&b[i]);
	
	ps.resize(N);
	v.resize(N);
	dfs(0);
	
	cout<<v[0][0]<<endl;
	
    return 0;
}
0