bzoj4326 noip 2015 运输计划

4326: NOIP2015 运输计划

Time Limit: 30 Sec   Memory Limit: 128 MB
Submit: 1841   Solved: 1173
[ Submit][ Status][ Discuss]

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

第一行包括两个正整数 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

输出文件只包含一个整数，表示小 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。

Source

[ Submit][ Status][ Discuss] 

HOME Back

刚开始看到题，最大值最小，树上路径问题又不带修改，直接树上差分+二分。然后还要在树上做下路径以及lca，尽管知识点很多，很杂，但是不难想，算是noip中比较宽容的一个d2t3吧。。。

先二分一个答案mid，然后把路径长度大于mid的运输计划取出来做树上差分，判断所有路径都有交集的点（因为只有这样才能把所有路径长度减到mid以下），if(f[j]==i&&mid>=p[1].dis-dist[j][0])如果最长路径可行说明下面都可行（事先按照降序排好）（dist[i][j]表示i起点向上2^1长度的路径总长）。如果这个答案可行，更新ans，然后r=mid，否则l=mid，一直更新答案即可。

具体细节见代码

#include <cstdio>
#include <vector>
#include <algorithm>
#include <cstring>
using namespace std;
vector<pair<int,int> >e[300005];
int n,m,ai,bi,ti,l=0,r,f[300005],anc[524289][20],dist[524289][20],deep[300005]={0,1},cnt,ans=2147483647;
struct trans
{
	int u,v,dis,Lca;
}p[300005];
int read()
{
    int ans=0,f=1;
    char c=getchar();
    while(c<'0'||c>'9')
	{
		if(c=='-') f=-1;
		c=getchar();
	}
    while(c>='0'&&c<='9')
    {
    	ans=ans*10+c-'0';
    	c=getchar();
    }
    return ans*f;
}
bool cmp(trans a,trans b)
{
	return a.dis>b.dis;
}
void dfs1(int now)
{
	for(int i=1;i<=19;i++)
	{
		anc[now][i]=anc[anc[now][i-1]][i-1];
		dist[now][i]=dist[anc[now][i-1]][i-1]+dist[now][i-1];
	}
	for(int i=0;i<e[now].size();i++)
		if(e[now][i].first!=anc[now][0])
		{
			deep[e[now][i].first]=deep[now]+1;
			dist[e[now][i].first][0]=e[now][i].second;
			anc[e[now][i].first][0]=now;
			dfs1(e[now][i].first);
		}
}
void lca(int cur,int u,int v)
{
	int ret=0;
	if(deep[u]<deep[v]) swap(u,v);
    for(int i=19;i>=0;i--)
        if(deep[v]<=deep[anc[u][i]])
        {
        	ret+=dist[u][i];
            u=anc[u][i];
        }
    if(u==v)
    {
		p[cur].Lca=u;
		p[cur].dis=ret;
		return;
    }
    for(int i=19;i>=0;i--)
        if(anc[u][i]!=anc[v][i])
        {
        	ret+=dist[u][i]+dist[v][i];
            u=anc[u][i];
            v=anc[v][i];
        }
    ret+=dist[u][0]+dist[v][0];
    p[cur].Lca=anc[u][0];
    p[cur].dis=ret;
}
void dfs2(int now)
{
	for(int i=0;i<e[now].size();i++)
		if(e[now][i].first!=anc[now][0])
		{
			dfs2(e[now][i].first);
			f[now]+=f[e[now][i].first];
		}
}
int main()
{
	n=read();m=read();
	for(int i=1;i<n;i++)
	{
		ai=read();bi=read();ti=read();
		e[ai].push_back(make_pair(bi,ti));
		e[bi].push_back(make_pair(ai,ti));
	}
	dfs1(1);
	for(int i=1;i<=m;i++)
	{
		p[i].u=read();p[i].v=read();
		lca(i,p[i].u,p[i].v);//把lca和dis预处理好，避免在线计算 
		r=max(r,p[i].dis);
	}
	sort(p+1,p+m+1,cmp);
	while(l<r)
	{
		int mid=cnt?r:(l+r)>>1,i=1,ok=0;//防止死循环以及漏掉答案
		memset(f,0,sizeof(f));
		for(;i<=m;i++)
		{
			if(p[i].dis<=mid) break;
			f[p[i].Lca]-=2;
			f[p[i].u]++;
			f[p[i].v]++;//树上差分 
		}
		i--;
		if(!i)
		{
			ans=min(ans,r=mid);
			continue;
		}
		dfs2(1);
		for(int j=2;j<=n;j++)
			if(f[j]==i&&mid>=p[1].dis-dist[j][0]) ok=1;
		if(ok&&!cnt) ans=min(ans,r=mid);
		else if(!ok&&!cnt) l=mid;
		if(r-l==1) cnt++;//防止死循环以及漏掉答案
		if(cnt==2) break;//防止死循环以及漏掉答案
	}
	printf("%d\n",ans==2147483647?0:ans);//细节要注意 
}