【ybtoj 4.5 倍增问题课堂过关】【洛谷P2680】E.运输计划【树上差分】

题目链接：gmoj，洛谷

分析

最大值最小，可以二分。
二分答案，然后判断的时候把所有大于$mid$值的路径记录下来，找出被所有这样路径覆盖的最长的道路。
如果没有这样的道路 false；如果这样的道路被减去之后依然大于mid false ，找出被所有路径覆盖的道路。
在树中将所有路径起、始权值+1，LCA权值-2，从所有叶节点往上累加，最终权值为路径数的点到其父亲的边为所求边 $dis[i]$表示i到根的距离；$t[i]$表示$i$这个点通往父亲的边，目的是记录这条边被遍历的次数 ；$num[i]$的作用是找到叶节点向上累加

上代码

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
using namespace std;

struct node
{
	int to,next,w;
}e[600001];

struct lwx
{
	int u,v,lcaa,diss;
}a[600001];

int n,m,dis[300001],fa[300001][24],dep[300001];//fa数组大小要注意看循环多少 
int tot,hd[300001],vis[300001];
int sum,cnt;
int t[300001],num[300001];

void add(int x,int y,int z)
{
	e[++tot]=(node){y,hd[x],z};
	hd[x]=tot;
}

void dfs(int x,int fx)
{
	num[++cnt]=x;
	dep[x]=dep[fx]+1;
	fa[x][0]=fx;
	for(int i=1;i<=23;i++)
	{
		fa[x][i]=fa[fa[x][i-1]][i-1];
	}
	for(int i=hd[x];i>0;i=e[i].next)
	{
		int k=e[i].to;
		if(k!=fx)
		{
			dis[k]=dis[x]+e[i].w;
			dfs(k,x);
		}
	}
}

int lca(int x,int y)	 
{
	if(dep[x]<dep[y]) swap(x,y);
	for(int i=23;i>=0;i--)
	{
		if(dep[fa[x][i]]>=dep[y])
		{
			x=fa[x][i];
		}
	}
	if(x==y) return x;
	for(int i=23;i>=0;i--)
	{
		if(fa[x][i]!=fa[y][i])
		{
			x=fa[x][i];
			y=fa[y][i];
		}
	}
	return fa[x][0];
}

int pd(int mid)
{
	int ans=0,cnt=0;
	memset(t,0,sizeof(t));
	for(int i=1;i<=m;i++)
	{
		if(a[i].diss>mid)
		{
			t[a[i].u]++;
			t[a[i].v]++;
			t[a[i].lcaa]-=2;
			ans=max(ans,a[i].diss-mid);
			cnt++; 
		}
	}
	if(cnt==0) return 1;
	for(int i=n;i>=1;i--)
	{
		t[fa[num[i]][0]]+=t[num[i]];
	}
	for(int i=2;i<=n;i++)
	{
		if(t[i]==cnt&&dis[i]-dis[fa[i][0]]>=ans) return 1;
	}
	return 0;
}

int main()
{
	cin>>n>>m;
	for(int i=1;i<=n-1;i++)
	{
		int x,y,w;
		scanf("%d%d%d",&x,&y,&w);
		add(x,y,w);
		add(y,x,w);
		sum+=w;//权值和 
	}
	dis[1]=0;
	dfs(1,0);
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&a[i].u,&a[i].v);
		a[i].lcaa=lca(a[i].u,a[i].v);
		a[i].diss=dis[a[i].u]+dis[a[i].v]-2*dis[a[i].lcaa];//点距离 
	}
	int l=0,r=sum;
	while(l<r)
    {
    	int mid=(l+r)/2;
    	if(pd(mid)) r=mid;
    	else l=mid+1;
	}
	cout<<l;
	return 0;
}