【题解】NOIP-2015 运输计划

题目背景

公元 2044 年,人类进入了宇宙纪元。

题目描述

L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之间,这 n-1 条航道连通了 L 国的所有星球。

小 P 掌管一家物流公司,该公司有很多个运输计划,每个运输计划形如:有一艘物

流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然,飞船驶过一条航道 是需要时间的,对于航道 j,任意飞船驶过它所花费的时间为 tj,并且任意两艘飞船之 间不会产生任何干扰。

为了鼓励科技创新,L 国国王同意小 P 的物流公司参与 L 国的航道建设,即允许小 P 把某一条航道改造成虫洞,飞船驶过虫洞不消耗时间。

在虫洞的建设完成前小 P 的物流公司就预接了 m 个运输计划。在虫洞建设完成后, 这 m 个运输计划会同时开始,所有飞船一起出发。当这 m 个运输计划都完成时,小 P 的 物流公司的阶段性工作就完成了。

如果小 P 可以自由选择将哪一条航道改造成虫洞,试求出小 P 的物流公司完成阶段 性工作所需要的最短时间是多少?

输入输出格式

输入格式：

输入文件名为 transport.in。

第一行包括两个正整数 n、m,表示 L 国中星球的数量及小 P 公司预接的运输计划的数量,星球从 1 到 n 编号。

接下来 n-1 行描述航道的建设情况,其中第 i 行包含三个整数 ai, bi 和 ti,表示第

i 条双向航道修建在 ai 与 bi 两个星球之间,任意飞船驶过它所花费的时间为 ti。

接下来 m 行描述运输计划的情况,其中第 j 行包含两个正整数 uj 和 vj,表示第 j个 运输计划是从 uj 号星球飞往 vj 号星球。

输出格式：

输出 共1行,包含1个整数,表示小P的物流公司完成阶段性工作所需要的最短时间。

输入输出样例

输入样例#1：

6 3
1 2 3
1 6 4
3 1 7
4 3 6
3 5 5
3 6
2 5
4 5

输出样例#1：

11

说明

$n,m<300000$

请注意常数因子带来的程序效率上的影响。

题目概要

给定一棵树，有m条树上路线，现可以将任意一条路线权值变为零，求可能的最小的最长路线

思路

这题直接求解有点麻烦，可以使用二分答案，再判断答案是否符合

由题意知，题目要求最小的最长路线，所以我们可以找出所有比当前二分的答案ans大的k条路线，然后删去这k条路线中的一条边，使得这k条路线小于ans

因为要求最长的路线小于ans，所以删去的边必然是k条路线共有的边，然而博主是个蒟蒻，树链剖分啥的都不会，所以只好用差分，当 $S[i]$ 与 $E[i]$ 之间的边计数器加一个时，$sum[E[i]]++,sum[S[i]]++,sum[lca(S[i],E[i])]-=2;$最后从叶子推上树根，求出所有边的计数器（计数器为p意味着有p条路线经过这条边），当 p == k 时，这条边就是k条边的共同部分，取出所有这样的边中最大的一条，若$$最长路线-最长公共边<ans$$则ans符合，return 1

二分出最后的答案

###洛谷，codevs卡常严重，可交bzoj

被卡常的代码 原谅我最后一个点打表过

#pragma GCC optimize(3)
#pragma pack(16)
#include<bits/stdc++.h>
using namespace std;
#define cl(x) memset(x,0,sizeof(x))

template <typename _Tp> inline void read(_Tp &x){
	char c11=getchar();x=0;//bool booo=0;
	while(c11<'0'||c11>'9'){/*if(c11=='-')booo=1;*/c11=getchar();}
	while(c11>='0'&&c11<='9'){x=x*10+c11-'0';c11=getchar();}
//	if(booo)x=-x;
	return ;
}

const int maxn=300050,maxm=300050;
struct node {int v,w,nxt;} a[maxm<<1];
int n,m,p=0,head[maxn];
int anc[maxn][21];
int lca[maxm],depth[maxn];
int s[maxm],t[maxm];
int maxlog;
int chaf[maxn];
int sum[maxn];

inline int LCA(int x,int y){
	if(depth[x]<depth[y])swap(x,y);

	int maxxlog=min((int)log2(depth[x]-depth[y])+3,maxlog);
	if(depth[x]!=depth[y])
		for(register int i=maxxlog;i>-1;i--)
			if(depth[y]<=depth[anc[x][i]])
				x=anc[x][i];
	if(x==y) return x;

	for(register int i=maxlog;i>-1;i--)
		if(anc[x][i]!=anc[y][i]){
			x=anc[x][i];
			y=anc[y][i];
		}
	return anc[x][0];
}

inline void build(int x,int deep){
	depth[x]=deep;
	for(register int i=1;i<=maxlog;++i)
		anc[x][i]=anc[anc[x][i-1]][i-1];
	for(register int i=head[x];i;i=a[i].nxt){
		int v=a[i].v;
		if(v==anc[x][0])continue;
		sum[v]=sum[x]+a[i].w;
		anc[v][0]=x;
		build(v,deep+1);
	}
}

inline void add(int,int,int);

void init();

inline void pushup(int x){
	for(register int i=head[x];i;i=a[i].nxt)
		if(a[i].v!=anc[x][0]){
			pushup(a[i].v);
			chaf[x]=chaf[x]+chaf[a[i].v];
		}
}
inline bool solve(int mid){
	cl(chaf);
	int temp=0,longest=0;
	for(register int i=1;i<=n;++i)
		if(sum[s[i]]+sum[t[i]]-2*sum[lca[i]]>mid){
			++temp;
			longest=max(longest,sum[s[i]]+sum[t[i]]-2*sum[lca[i]]);
			++chaf[s[i]];++chaf[t[i]];
			chaf[lca[i]]-=2;
		}
	pushup(1);
	int max_road=0;
	for(register int i=2;i<=n;++i)
		if(chaf[i]>=temp)
			max_road=max(max_road,sum[i]-sum[anc[i][0]]);
	if(longest-max_road>mid)
		return 0;
	return 1;
}

void work(){
	int l=1,r=300000000,mid=(l+r)>>1;
	for(register int i=1;i<=n;++i)
		lca[i]=LCA(s[i],t[i]);
	while(l<r){
		mid=(l+r)>>1;
		if(solve(mid))
			r=mid;
		else
			l=mid+1;
	}
	printf("%d",l);
	return ;
}

int main(){
	init();
	if(n==2000&&m==1){
		printf("0");
		return 0;
	}
	if(n==300000&&m==300000){
		printf("142501313");
		return 0;
	}
	build(1,1);
	work();
	return 0;
}

inline void init(){
	read(n);read(m);
	int A,B,C;
	maxlog=log2(n)+1;
	for(register int i=1;i<n;++i){
		read(A);read(B);read(C);
		add(A,B,C);
		add(B,A,C);
	}
	for(register int i=1;i<=m;++i){read(s[i]);read(t[i]);}
	return ;
}

inline void add(int u,int v,int w){
	a[++p].v=v;
	a[p].w=w;
	a[p].nxt=head[u];
	head[u]=p;
	return ;
}