题解：P2680 [NOIP 2015 提高组] 运输计划

同步发布至 洛谷

题解——关于这篇树上差分和树上倍增的经典题目

根据题目，我们不难看出这道题我们其实就是要找一条最优的边，使其被改造为虫洞之后，所有的运输计划的时间最小。

原本的话，其实可以打贪心的，也就是说选择经过最多的点，但是嘛，这道题是一道和手的练习树上问题的题目，所以说还是打差分吧。

因为任意一条边都有可能，所以说我们要枚举每条边，然后来判断。但是呢，这样的时间复杂度为 $O(n)$，如果里面再套一层循环啥的，就废了，铁打的超时。

所以呢，我们不难可以想到，找到最大的一条边不就行？然后根据最大的这条边，来二分答案，也就是直接枚举答案，看当前的这个值能不能达到，这样的时间复杂度就是 $O(\log n)$，那么现在的问题就是怎么检查。

我们先令所有边中最大的一条边为 $totmaxn$，当前二分所得到的答案为 $mid$。

首先我们一定要找出来所有的大于当前答案的边，因为这些边是一定要被减去时间的，然后我们可以根据树上差分，来标记。
标记完所有的答案之后呢，就对树上差分求子树和，然后找到所有点所经过的一条或多条边（易迷惑点 1），然后在这些边中找到一个最大的（假设为 $maxn$），然后用最大的一条边 $totmaxn$ 减去这条最大的可以节省时间的边 $maxn$（易迷惑点 2），看看是不是小于或等于我当前二分的答案（即 $maxn\le mid$），如果是，说明当前的答案是可以的，如果不是，代表当前答案就是不是可以的。

这样检查的时间复杂度为 $O(n+m)$，不算太高，配上外围的 $O(\log n)$ 刚刚好。

解释一下上面的两个可能会让人迷惑的点：

易迷惑点 1：
为什么是要找到所有点都要经过的一条边？请注意，这里的所有点指的是大于当前答案的点，即大于 $mid$ 的点，如果说一个要被开成虫洞的点没有经过当前所有大于 $mid$ 的点的话，是一定不行的，这是很显然的吧。

易迷惑点 2：
为什么要用最大的一条边来减当前的答案？请注意，如果说最大的一条边都可以成功的话，为什么比他要小的点成功不了呢？换言之，如果最大的一条边都不能成功，说明，当前的答案也是不可以的啊。

其他要注意的事项：

1. 我们要先用树上倍增处理出来每个运输计划的时间，以方便我们下面使用。
2. 注意在写树上差分的时候要先求出根节点。
3. 除了这些的话，还是开一个读写优化吧，要不然要废掉啊。

我相信一定是没有人喜欢看上面的讲解吧，所以我们直接上代码吧。

ACcode

#include<bits/stdc++.h>
using namespace std;
const int INF=3e5+10;

struct Node{
	int point,num;
};
vector<Node> mp[INF];
int a[INF],b[INF],dp[INF][30],w[INF][30],p[INF],edge[INF];
int deep[INF],dis[INF];
int maxlen=INT_MIN,maxtot=INT_MIN,ans=INT_MAX;
void prepare(int x,int fa){
	for (int j=1;(1<<j)<=deep[x]-1;j++){
		dp[x][j]=dp[dp[x][j-1]][j-1];
		w[x][j]=w[dp[x][j-1]][j-1]+w[x][j-1];
	}
	int len=mp[x].size();
	for (int i=0;i<len;i++){
		if (mp[x][i].point==fa)continue;
		int t=mp[x][i].point;
		deep[t]=deep[x]+1,dp[t][0]=x,w[t][0]=mp[x][i].num;
		prepare(t,x); 
	}
}

int lca_length(int x,int y){
	int tot=0;
	if (deep[x]<deep[y])swap(x,y);
	int index=__lg(deep[x]-deep[y]);
	for (int i=index;i>=0;i--){
		if (deep[dp[x][i]]>=deep[y])tot+=w[x][i],x=dp[x][i];
		if (deep[x]==deep[y])break;
	}
	if (x==y)return tot;
	for (int i=20;i>=0;i--){
		if (dp[x][i]!=dp[y][i])tot+=w[x][i]+w[y][i],x=dp[x][i],y=dp[y][i];
	}
	tot+=w[x][0]+w[y][0];
	return tot;
}

int lca_root(int x,int y){
	if (deep[x]<deep[y])swap(x,y);
	int index=__lg(deep[x]-deep[y]);
	for (int i=index;i>=0;i--){
		if (deep[dp[x][i]]>=deep[y])x=dp[x][i];
		if (deep[x]==deep[y])break;
	}
	if (x==y)return x;
	for (int i=20;i>=0;i--){
		if (dp[x][i]!=dp[y][i])x=dp[x][i],y=dp[y][i];
	}
	return dp[x][0];
}

void get(int x,int fa){
	int len=mp[x].size();
	for (int i=0;i<len;i++){
		if (mp[x][i].point==fa)continue;
		int t=mp[x][i].point;
		get(t,x);
		p[x]+=p[t];
	}
} 

int main(){
	ios::sync_with_stdio(0);
	cin.tie(0),cout.tie(0);
	int n,m;
	cin>>n>>m;
	for (int i=1;i<n;i++){
		int u,v,w;
		cin>>u>>v>>w;
		mp[u].push_back({v,w});
		mp[v].push_back({u,w});
		maxlen=max(maxlen,w);
	}
	deep[1]=1;
	prepare(1,-1);
	for (int i=1;i<=m;i++){
		cin>>a[i]>>b[i];
		dis[i]=lca_length(a[i],b[i]);
		maxtot=max(dis[i],maxtot);
	}
//	cout<<maxtot<<endl;
	int l=maxtot-maxlen,r=300000000;
	while (l<=r){
		int maxn=0,tim=0;
		memset(p,0,sizeof(p));
		int mid=(l+r)>>1;
		for (int i=1;i<=m;i++){
			if (dis[i]>mid){
				int root=lca_root(a[i],b[i]);
				p[a[i]]++,p[b[i]]++,p[root]-=2;
				tim++;
			}
		}
		if (tim==0){
			ans=mid;
			r=mid-1;
			continue; 
		}
		get(1,-1);
		for (int i=1;i<=n;i++){
			int len=mp[i].size();
			for (int j=0;j<len;j++){
				if (deep[mp[i][j].point]>deep[i]&&p[mp[i][j].point]==tim)maxn=max(maxn,mp[i][j].num);
			} 
		}
		if (maxtot-maxn<=mid){
			ans=mid;
			r=mid-1;
		}else l=mid+1;
	}
	cout<<ans;
	return 0;
}