Submission #1070289

#	Time	Username	Problem	Language	Result	Execution time	Memory
1070289		pcc	Closing Time (IOI23_closing)	C++17	9 / 100	42 ms	10844 KiB

closing

#include "closing.h"
#include <vector>
#include <bits/stdc++.h>
using namespace std;

#define pll pair<ll,ll>
#define pii pair<int,int>
#define fs first
#define sc second
#define ll long long
const int mxn = 3020;
const ll inf = 2e18;

vector<pii> tree[mxn];
pii eul[mxn];
int ptr = 0;
ll dp[mxn][4][mxn*2];
ll arr[mxn],brr[mxn],crr[mxn];
ll dist[mxn];
int N;
int X,Y;
ll K;
int sz[mxn];
ll tmp[4][mxn*2];

void init(){
	ptr = 0;
	for(int i = 0;i<=N;i++)tree[i].clear();
}

void dfs(int now,int par){
	if(now == par)dist[now] = 0;
	for(auto [nxt,w]:tree[now]){
		if(nxt == par)continue;
		dist[nxt] = dist[now]+w;
		dfs(nxt,now);
	}
	return;
}

void dfs2(int now,int par){
	eul[now].fs = ++ptr;
	for(auto [nxt,w]:tree[now]){
		if(nxt == par)continue;
		dfs2(nxt,now);
	}
	eul[now].sc = ptr;
	return;
}

void calc(ll dp1[],ll dp2[],int s1,int s2,ll rdp[]){
	assert(max(s1,s2)<=N*2);
	for(int i = 0;i<=s1-s2;i++){
		for(int j = 0;j<=s2;j++){
			rdp[i+j] = min({rdp[i+j],dp1[i]+dp2[j]});
		}
	}
	return;
}

void dfs3(int now,int par){
	sz[now] = 2;
	for(int i = 0;i<4;i++)fill(dp[now][i],dp[now][i]+mxn*2,inf);
	dp[now][0][0] = 0;
	dp[now][1][1] = arr[now];
	dp[now][2][1] = brr[now];
	dp[now][3][2] = crr[now];
	//cerr<<now<<":"<<arr[now]<<' '<<brr[now]<<' '<<crr[now]<<' '<<chain<<endl;
	for(auto [nxt,w]:tree[now]){
		if(nxt == par)continue;
		dfs3(nxt,now);
		sz[now] += sz[nxt];

		for(int i = 0;i<4;i++){
			fill(tmp[i],tmp[i]+mxn*2,inf);
			for(int j = 0;j<=sz[now];j++)tmp[i][j] = dp[now][i][j];
		}

		bool chain = false;
		if(eul[Y].fs>=eul[nxt].fs&&eul[Y].sc<=eul[nxt].sc)chain = true;

		if(!chain){

			calc(dp[now][1],dp[nxt][1],sz[now],sz[nxt],tmp[1]);
			calc(dp[now][3],dp[nxt][1],sz[now],sz[nxt],tmp[3]);
			calc(dp[now][2],dp[nxt][2],sz[now],sz[nxt],tmp[2]);
			calc(dp[now][3],dp[nxt][2],sz[now],sz[nxt],tmp[3]);
			calc(dp[now][3],dp[nxt][3],sz[now],sz[nxt],tmp[3]);

		}
		else{

			calc(dp[now][0],dp[nxt][0],sz[now],sz[nxt],tmp[0]);
			calc(dp[now][1],dp[nxt][0],sz[now],sz[nxt],tmp[1]);
			calc(dp[now][1],dp[nxt][1],sz[now],sz[nxt],tmp[1]);
			calc(dp[now][0],dp[nxt][2],sz[now],sz[nxt],tmp[0]);
			calc(dp[now][1],dp[nxt][2],sz[now],sz[nxt],tmp[1]);
			calc(dp[now][2],dp[nxt][2],sz[now],sz[nxt],tmp[2]);
			calc(dp[now][3],dp[nxt][2],sz[now],sz[nxt],tmp[3]);
			calc(dp[now][1],dp[nxt][3],sz[now],sz[nxt],tmp[1]);
			calc(dp[now][3],dp[nxt][3],sz[now],sz[nxt],tmp[3]);

		}

		for(int i = 0;i<4;i++){
			for(int j = 0;j<=sz[now];j++)dp[now][i][j] = min(dp[now][i][j],tmp[i][j]);
		}
	}


	/*
	for(int i =0;i<4;i++){
		for(int j = 0;j<=sz[now];j++){
			cerr<<now<<' '<<i<<' '<<j<<":"<<dp[now][i][j]<<endl;
		}
	}
	*/

	//for(int i = 0;i<4;i++)dp[now][i][0] = 0;

	return;
}

int max_score(int NN, int XX, int YY, long long KK,
              std::vector<int> U, std::vector<int> V, std::vector<int> W){

	X = XX,Y = YY;
	K = KK;
	init();
	N = NN;
	for(int i = 0;i<U.size();i++){
		int a = U[i],b = V[i],w = W[i];
		tree[a].push_back(pii(b,w));
		tree[b].push_back(pii(a,w));
	}
	dfs(X,X);
	for(int i = 0;i<N;i++)arr[i] = dist[i];
	dfs(Y,Y);
	for(int i = 0;i<N;i++)brr[i] = dist[i];
	for(int i = 0;i<N;i++)crr[i] = max(arr[i],brr[i]);

	dfs2(X,X);
	dfs3(X,X);

	int ans = 0;

	for(int i =0;i<4;i++){
		for(int j = 0;j<=N*2;j++){
			if(dp[X][i][j]<=K)ans = max(ans,j);
		}
	}

	return ans;
}

Compilation message (stderr)

closing.cpp: In function 'int max_score(int, int, int, long long int, std::vector<int>, std::vector<int>, std::vector<int>)':
closing.cpp:131:17: warning: comparison of integer expressions of different signedness: 'int' and 'std::vector<int>::size_type' {aka 'long unsigned int'} [-Wsign-compare]
  131 |  for(int i = 0;i<U.size();i++){
      |                ~^~~~~~~~~

• Subtask 0: 0 / 0
• Subtask 1: 0 / 8
• Subtask 2: 0 / 9
• Subtask 3: 0 / 12
• Subtask 4: 0 / 14
• Subtask 5: 9 / 9
• Subtask 6: 0 / 11
• Subtask 7: 0 / 10
• Subtask 8: 0 / 10
• Subtask 9: 0 / 17