洛谷 P2680 运输计划

最新推荐文章于 2022-07-23 12:46:10 发布
最新推荐文章于 2022-07-23 12:46:10 发布
阅读量432 收藏 1
点赞数 1
分类专栏： 倍增 搜索 图论 文章标签： 搜索 lca 图论
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Rlt1296/article/details/53022618
版权
图论 同时被 3 个专栏收录
89 篇文章 0 订阅
订阅专栏
搜索
45 篇文章 0 订阅
订阅专栏
倍增
5 篇文章 0 订阅
订阅专栏

题目背景

公元 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 号星球。

输出格式：

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

说明

所有测试数据的范围和特点如下表所示


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


预处理出询问的lca，枚举在那条路建立虫洞，dfs暴力更新到根节点的距离，可以得50分。


#include<algorithm>
#include<iostream>
#include<cstdio>
using namespace std;
const int N=300005;
int n,m,cnt,ans=1e9+7,hd[N],dep[N],f[N][21],dis[N],ban;
struct edge
{
	int to,nxt,val,num;
}v[2*N];
struct node
{
	int x,y,lcaa;
}e[N];
void addedge(int x,int y,int z,int id)
{
	++cnt;
	v[cnt].to=y;
	v[cnt].nxt=hd[x];
	v[cnt].val=z;
	v[cnt].num=id;
	hd[x]=cnt;
}
void init(int u,int fa)
{
	f[u][0]=fa;
	dep[u]=dep[fa]+1;
	for(int i=hd[u];i;i=v[i].nxt)
		if(v[i].to!=fa)
			init(v[i].to,u);
}
int lca(int x,int y)
{
	if(dep[x]<dep[y])
		swap(x,y);
	for(int i=18;i>=0;i--)
		if(dep[f[x][i]]>=dep[y])
			x=f[x][i];
	if(x==y)
		return x;
	for(int i=18;i>=0;i--)
		if(f[x][i]!=f[y][i])
			x=f[x][i],y=f[y][i];
	return f[x][0];
}
void dfs(int u,int fa)
{
	for(int i=hd[u];i;i=v[i].nxt)
		if(v[i].to!=fa)
		{
			if(v[i].num!=ban)
				dis[v[i].to]=dis[u]+v[i].val;
			else
				dis[v[i].to]=dis[u];
			dfs(v[i].to,u);
		}
}
int main()
{
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n-1;i++)
	{
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		addedge(x,y,z,i);
		addedge(y,x,z,i);
	}
	init(1,0);
	for(int i=1;i<=18;i++)
		for(int j=1;j<=n;j++)
			f[j][i]=f[f[j][i-1]][i-1];
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&e[i].x,&e[i].y);
		e[i].lcaa=lca(e[i].x,e[i].y);
	}
	for(int i=1;i<=n-1;i++)//删边 
	{
		ban=i;
		dfs(1,0);
		int mx=0;
		for(int j=1;j<=m;j++)
			mx=max(mx,dis[e[j].x]+dis[e[j].y]-2*dis[e[j].lcaa]);
		ans=min(ans,mx);
	}
	printf("%d\n",ans);
	return 0;
}

迷之二分答案，转变为可行性问题。
如果所有的边都比lim小于等于，一定可行，否则从两端点搜索到lca，将经过的边++，如果有一条边所有的计划都经过了，且mx-边权<=lim，就可行。80分。

 
#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int N=300005;
int n,m,cnt,hd[N],dep[N],f[N][21],dis[N],ac[N],val[N];
struct edge
{
	int to,nxt,num;
}v[2*N];
struct node
{
	int x,y,lca,dist;
}e[N];
void addedge(int x,int y,int id)
{
	++cnt;
	v[cnt].to=y;
	v[cnt].nxt=hd[x];
	v[cnt].num=id;
	hd[x]=cnt;
}
void init(int u,int fa)
{
	f[u][0]=fa;
	dep[u]=dep[fa]+1;
	for(int i=hd[u];i;i=v[i].nxt)
		if(v[i].to!=fa)
		{
			dis[v[i].to]=dis[u]+val[v[i].num];
			init(v[i].to,u);
		}
}
int getlca(int x,int y)
{
	if(dep[x]<dep[y])
		swap(x,y);
	for(int i=17;i>=0;i--)
		if(dep[f[x][i]]>=dep[y])
			x=f[x][i];
	if(x==y)
		return x;
	for(int i=17;i>=0;i--)
		if(f[x][i]!=f[y][i])
			x=f[x][i],y=f[y][i];
	return f[x][0];
}
void dfs(int u,int ban)
{
	if(u==ban)
		return ;
	for(int i=hd[u];i;i=v[i].nxt)
		if(dep[v[i].to]<dep[u])
		{
			ac[v[i].num]++;
			dfs(v[i].to,ban);
		}
}
bool pd(int lim)
{
	memset(ac,0,sizeof(ac));
	int sum=0,mx=0;
	for(int i=1;i<=m;i++)
		if(e[i].dist>lim)
		{
			sum++;
			mx=max(mx,e[i].dist);
			dfs(e[i].x,e[i].lca);
			dfs(e[i].y,e[i].lca);
		}
	if(sum==0)
		return 1;
	for(int i=1;i<=n-1;i++)
		if(ac[i]==sum&&mx-val[i]<=lim)
			return 1;
	return 0;
}
int main()
{
	int l=0,r=0,mid;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=n-1;i++)
	{
		int x,y,z;
		scanf("%d%d%d",&x,&y,&z);
		addedge(x,y,i);
		addedge(y,x,i);
		val[i]=z;
	}
	init(1,0);
	for(int i=1;i<=17;i++)
		for(int j=1;j<=n;j++)
			f[j][i]=f[f[j][i-1]][i-1];
	for(int i=1;i<=m;i++)
	{
		scanf("%d%d",&e[i].x,&e[i].y);
		e[i].lca=getlca(e[i].x,e[i].y);
		e[i].dist=dis[e[i].x]+dis[e[i].y]-2*dis[e[i].lca];
		r=max(r,e[i].dist);
	}
	while(l<r)
	{
		mid=(l+r)/2;
		if(!pd(mid))
			l=mid+1;
		else
			r=mid;
	}
	printf("%d\n",l);
	return 0;
}

想一想，每次都直接dfs太zz了，可以把它们存起来，从根节点只dfs一遍，就是差分（不懂），详见代码。（终于ac了） 

#include<algorithm>
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
const int N=300005;
int n,m,cnt,hd[N],dep[N],f[N][21],dis[N],ac[N],val[N],hh[N];
struct edge
{
    int to,nxt,num;
}v[2*N];
struct node
{
    int x,y,lca,dist;
}e[N];
void addedge(int x,int y,int id)
{
    ++cnt;
    v[cnt].to=y;
    v[cnt].nxt=hd[x];
    v[cnt].num=id;
    hd[x]=cnt;
}
void init(int u,int fa)
{
    f[u][0]=fa;
    dep[u]=dep[fa]+1;
    for(int i=hd[u];i;i=v[i].nxt)
        if(v[i].to!=fa)
        {
            dis[v[i].to]=dis[u]+val[v[i].num];
            init(v[i].to,u);
        }
}
int getlca(int x,int y)
{
    if(dep[x]<dep[y])
        swap(x,y);
    for(int i=17;i>=0;i--)
        if(dep[f[x][i]]>=dep[y])
            x=f[x][i];
    if(x==y)
        return x;
    for(int i=17;i>=0;i--)
        if(f[x][i]!=f[y][i])
            x=f[x][i],y=f[y][i];
    return f[x][0];
}
int dfs(int u,int fa)
{
    int t=0;
    for(int i=hd[u];i;i=v[i].nxt)
        if(v[i].to!=fa)
        {
            int x=dfs(v[i].to,u);
            ac[v[i].num]+=x;
            t+=x;
        }
    return t+hh[u];
}
bool pd(int lim)
{
    memset(ac,0,sizeof(ac));
    memset(hh,0,sizeof(hh));
    int sum=0,mx=0;
    for(int i=1;i<=m;i++)
        if(e[i].dist>lim)
        {
            sum++;
            mx=max(mx,e[i].dist);
            hh[e[i].x]++;
            hh[e[i].y]++;
            hh[e[i].lca]-=2;
        }
    if(sum==0)
        return 1;
    dfs(1,0);
    for(int i=1;i<=n-1;i++)
        if(ac[i]==sum&&mx-val[i]<=lim)
            return 1;
    return 0;
}
int main()
{
    int l=0,r=0,mid;
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n-1;i++)
    {
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        addedge(x,y,i);
        addedge(y,x,i);
        val[i]=z;
    }
    init(1,0);
    for(int i=1;i<=17;i++)
        for(int j=1;j<=n;j++)
            f[j][i]=f[f[j][i-1]][i-1];
    for(int i=1;i<=m;i++)
    {
        scanf("%d%d",&e[i].x,&e[i].y);
        e[i].lca=getlca(e[i].x,e[i].y);
        e[i].dist=dis[e[i].x]+dis[e[i].y]-2*dis[e[i].lca];
        r=max(r,e[i].dist);
    }
    while(l<r)
    {
        mid=(l+r)/2;
        if(!pd(mid))
            l=mid+1;
        else
            r=mid;
    }
    printf("%d\n",l);
    return 0;
}

「已注销」
关注
  • 1
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
洛谷P1422解答（C++）
08-12
洛谷P1422解答（C++）请勿在未经授权的情况下上传任何涉及著作权侵权的资源，除非该资源完全由您个人创作
洛谷P1563玩具谜题
08-26
洛谷P1563玩具谜题
[luogu]P2680 运输计划
Myka1991的博客
10-28 131
原题链接 ：P2680 运输计划 分析 题意很简单，给定一张连通图，n个点，n-1条边，很显然是一棵树。 现在给定m条链\((u,v)\)。 现在可以cut掉一条边(边权置为0)。 现在求最长链的最小值。 55pts 蒟蒻的我肯定是拿不到满分的。。我们直接考虑部分分。 对于m=1的点，我们只需要枚举链上的全部边，然后cut掉最长的那一个就可以了。 裸搜一遍dfs，然后求出链长，最后减去最...
P2680 [NOIP2015 提高组] 运输计划
You_are_hanson的博客
07-23 388
P2680 [NOIP2015 提高组] 运输计划
洛谷 2680 运输计划
ToheartZhang的博客
08-24 350
//二分+倍增LCA+树上前缀和 //尽力缩小二分范围，能多拿分就多拿分 //掌握树上前缀和的初始化和区间修改 //读入优化 //二分答案，然后找出超出mid的计划，记录他们的公共重边，这里用树上前缀和找前缀和exc==超出计划数 的最大的边，删掉后再判断是否<=mid。 #include<bits/stdc++.h> using namespace std; const int inf = 0
【NOIP2015】洛谷2680 运输计划【解法二】
YHaX的博客
10-28 957
二分答案+LCA+差分
洛谷 P1164 小A点菜讲解PPT
03-12
洛谷 P1164 小A点菜讲解PPT 本资源是关于洛谷 P1164 题目的讲解PPT，题目名称为小A点菜，通过对题目的分析和解决方法的讲解，帮助学习者更好地理解和掌握该题目的解决思路和方法。 知识点一：01背包问题 在本题目...
洛谷P1028.cpp
03-21
洛谷P1028.cpp
洛谷P4035参考源码
01-11
思路1（洛谷题解） 设n维球体为α，其半径为r（注意，这是一个设而不求。），其球心X的坐标为(x_1,x_2,…,x_n )。∀A_1,A_2,…,A_(n+1)∈α，点A_i (1≤i≤n+1)坐标为(a_((i,1) ),a_((i,2) ),…,a_((i,n) ) )。由n维...
【NOIP2015】洛谷2680 运输计划【解法一】
YHaX的博客
10-28 1235
二分答案+LCA
洛谷 P2680 运输计划 树链剖分+最近公共祖先
liangzihao1的博客
08-03 295
题目背景 公元2044年，人类进入了宇宙纪元。 题目描述 公元2044年，人类进入了宇宙纪元。 LLL国有nnn个星球，还有n−1n−1n−1条双向航道，每条航道建立在两个球之间，这n−1n−1n−1条航道连通了LLL国的所有星球。 小PPP掌管一家物流公司，该公司有很多个运输计划，每个运输计划形如：有一艘物流飞船需要从uiuiu_i号星球沿最快的宇航路径飞行到viviv_i号星球去。显...
洛谷P2680 运输计划（树上差分+二分）
weixin_33816300的博客
09-19 97
传送门 考虑树上乱搞 首先这是满足二分性质的，如果在某个时间可以完成工作那么比他更长的时间肯定也能完成工作 然后考虑二分，设当前答案为$mid$，如果有一条链的长度大于$mid$，那么这条链上必须得删去一条边。我们可以贪心的删去所有可以删去的边中最长的，然后看看最长边减去删去的边是否小于等于$mid$，如果成立说明可行 然后考虑怎么solve，我们可以用树上差分，就是对于每个点把它看...
P2680 运输计划
weixin_30873847的博客
09-09 81
题目背景 公元 2044 年,人类进入了宇宙纪元。 题目描述 L 国有 n 个星球,还有 n-1 条双向航道,每条航道建立在两个星球之间,这 n-1 条航道连通了 L 国的所有星球。 小 P 掌管一家物流公司,该公司有很多个运输计划,每个运输计划形如:有一艘物 流飞船需要从 ui 号星球沿最快的宇航路径飞行到 vi 号星球去。显然,飞船驶过一条航道 是需要时间的,对于航道 j,任意飞船驶...
P2680 【运输计划】
weixin_42404190的博客
08-22 142
绝对是一个全新的暴力思路。~只不过要开O2才能过。 接下来进入正题。 看到有一条边可以变成0，很自然地就想到要删边。 那么怎么判断一条边是不是在一条路径上呢？ 设一条边的两个端点为u,v; 那条路径的两个端点为x,y; 我们把这个问题拆成两个小问题来做， 1先判定一个点是不是在那条路径上； 2再判断另一个点是不是也在那条路径上； 判定一个点是不是在一条路径上满足如下等式即可 设dis(i,j)为两...
洛谷P2680运输计划
weixin_30487317的博客
11-08 70
传送门啦 要求的就是，把树上的一条边的权值设为0之后，所有路径中的最大值的最小值。 首先二分最大值，假设某次二分的最大值为x，我们首先找出所有大于x的路径（也就是我们需要通过改权缩短的路径），并把路径上的所有边都标记一下。 在标记完成后，枚举所有边，如果存在一条边位于所有长度大于于x的路径上，并且删除之后能使所有路径都满足 $ length<=x $ ，则返回 $ true $ ，否则 $ ...
P2680 运输计划（二分+树上差分）
weixin_30872499的博客
05-27 96
P2680 运输计划 链接 分析： 　　二分+树上差分。 　　首先可以二分一个答案，那么所有比这个答案大的路径，都需要减去些东西才可以满足这个答案。 　　那么减去的这条边一定在所有的路径的交集上。 　　那么如果求快速的求出这个交集并判断呢，树剖可以，把所有大于的路径都标记一下，然后判断，复杂度太大了。 　　于是用到了树上差分，get新技能。 　　在两个端点出标记+1，在lca...
luogu P2680 运输计划
zsyz_ZZY的博客
03-21 280
背景： 好久之前的坑… 题目传送门： https://www.luogu.org/problemnew/show/P2680 题意： 有一棵树，现在让你选取一条边的距离改为000，求树上mmm个点对到达的用时（走单位111的距离用时为111）。 思路： 谨慎阅读，非正解+卡常。 显然我们不能暴力。于是我们考虑将哪一条边的权值改为000的问题求最小值就变成了二分最后的时间（有单调性）。 ...
洛谷p5721c语言
最新发布
11-13
这里提供两个C语言程序的例子，一个是数数的手指，另一个是洛谷p5721题目的解答。 1. 数数的手指 ```c #include int main() { int i; for (i = 1; i ; i++) { printf("%d ", i); } return 0; } ``` 这个...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值