NOIP2015DAY2T3【运输计划】

Description

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

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

小 P 掌管一家物流公司， 该公司有很多个运输计划，每个运输计划形如：有一艘物流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然，飞船驶过一条航道是需要时间的，对于航道 j，任意飞船驶过它所花费的时间为 tj，并且任意两艘飞船之间不会产生任何干扰。

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

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

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

Input Format

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

接下来 n−1 行描述航道的建设情况，其中第 i 行包含三个整数 ai,bi 和 ti，表示第 i 条双向航道修建在 ai 与 bi 两个星球之间，任意飞船驶过它所花费的时间为 ti。数据保证 1≤ai,bi≤n 且 0≤ti≤1000。

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

Output Format

输出文件只包含一个整数，表示小 P 的物流公司完成阶段性工作所需要的最短时间。

Sample Input

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

Sample Output

11

Hint

将第 1 条航道改造成虫洞： 则三个计划耗时分别为：11,12,11，故需要花费的时间为 12。

将第 2 条航道改造成虫洞： 则三个计划耗时分别为：7,15,11，故需要花费的时间为 15。

将第 3 条航道改造成虫洞： 则三个计划耗时分别为：4,8,11，故需要花费的时间为 11。

将第 4 条航道改造成虫洞： 则三个计划耗时分别为：11,15,5，故需要花费的时间为 15。

将第 5 条航道改造成虫洞： 则三个计划耗时分别为：11,10,6，故需要花费的时间为 11。

故将第 3 条或第 5 条航道改造成虫洞均可使得完成阶段性工作的耗时最短，需要花费的时间为 11。 

【题解】

 二分答案+在线lca

 由于答案很难直接算出来运用二分答案将问题转化为可行性

 若经改造之后m条路径的长度都小于ans则ans可行

接下来考虑如何算每条路径长度及删除哪条边

求路径可以在树上倍增求lca时实现dist=dist[u]+dist[v]-2*dist[lca(u,v)]（不会lca欢迎观看本博客其他lca文章中的详细介绍）（我这题是用离线lca tarjan大家可以自行百度）

假设未改造之前有k条路径长度大于ans那么肯定是删除被k条路径都经过的最长边最优

首先考虑如何找到被k条路径都经过的边  可以用差分打标记的方法 在u,v节点各打上1，在lca(u,v)打上-1 用dfs遍历一遍（树型动规）

然后把所有经过k次的边挑出来，排序一下找最长边，这条边就是删除的边，若最大路径长-删除边边长<=ans，ans就可行

详见代码（好像有点冗余部分）

#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <string>
#include <cstdio>
#include <queue>
#include <ctime>
#include <cmath>
using namespace std;
int firstr[300005],firstq[300005],lca[300005],rec[300005],vis[300005],fa[300005],numr,numq;
int dis[300005],f[600005],a[300005];
int i,j,k,l,m,n,r,mid,ans,x,y,z;
struct info
  {
  	int fr,ar,l,next;
  }road[600005],ques[600005];
void add(int x,int y,int z)
  {
  	numr++;road[numr]=(info){x,y,z,firstr[x]};firstr[x]=numr;
  }
void add1(int x,int y,int z)
  {
  	numq++;ques[numq]=(info){x,y,z,firstq[x]};firstq[x]=numq;
  }

int getfa(int x)
  {
  	if (fa[x]==x) return x;else return fa[x]=getfa(fa[x]);
  }
void tarjan(int u)
  {
  	int i,v;
  	vis[u]=1;fa[u]=u;
  	for (i=firstq[u];i;i=ques[i].next)
  	  {
		v=ques[i].ar;
		if (vis[v]) 
		  {
		  	lca[ques[i].l]=getfa(v);
		  	rec[ques[i].l]=dis[u]+dis[v]-2*dis[lca[ques[i].l]];
		  }
	  }
	for (i=firstr[u];i;i=road[i].next)
  	  {
		v=road[i].ar;
		if (!vis[v]) 
		  {
		  	dis[v]=dis[u]+road[i].l;
		  	tarjan(v);
		  	fa[getfa(v)]=u;
		  }
	  }
  }
void dfs(int u,int fa)
  {
    int i;
    for (i=firstr[u];i;i=road[i].next)
      if (road[i].ar!=fa)
        {
           dfs(road[i].ar,u);
           a[u]+=a[road[i].ar];
           f[i]=a[road[i].ar];
		}
  }
bool pd(int l)
  {
  	int i,j,k,mx;
  	memset(a,0,sizeof(a));
  	memset(f,0,sizeof(f));
  	for (i=1,k=0;i<=m;i++)
	   if (rec[i]>l) 
	    {
	      k++;a[ques[2*i].fr]++;a[ques[2*i].ar]++;a[lca[i]]-=2;
		}
  	dfs(1,-1);
  	for (i=1,mx=0;i<=numr;i++)
  	  if (f[i]==k) mx=max(mx,road[i].l);
  	for (i=1;i<=m;i++) if (rec[i]-mx>l) return 0;
  	return 1;
  }
int main()
  {
  	scanf("%d%d",&n,&m);
  	for (i=1;i<n;i++) 
	  scanf("%d%d%d",&x,&y,&z),add(x,y,z),add(y,x,z);
	for (i=1;i<=m;i++) 
	  scanf("%d%d",&x,&y),add1(x,y,i),add1(y,x,i);
	tarjan(1);
    for (i=1;i<=m;i++) r=max(r,rec[i]);
    for (;l<=r;)
      {
      	mid=(l+r)/2;
      	if (pd(mid)) r=mid-1,ans=mid;else l=mid+1;
	  }
	printf("%d",ans);
  }