bzoj 3445: [Usaco2014 Feb] Roadblock （堆优化dijkstra）

最新推荐文章于 2023-05-19 01:00:35 发布

clover_hxy

最新推荐文章于 2023-05-19 01:00:35 发布

阅读量2.3k

点赞数

CC 4.0 BY-SA版权

分类专栏：图论

本文链接：https://blog.youkuaiyun.com/clover_hxy/article/details/53186282

图论专栏收录该内容

102 篇文章

订阅专栏

本文介绍了一个USACO竞赛题目“Roadblock”的解决思路，该问题要求找到延长特定路径的方法。通过两次Dijkstra算法计算最短路径并对比结果，最终确定最佳策略。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

3445: [Usaco2014 Feb] Roadblock

Time Limit: 10 Sec Memory Limit: 128 MB
Submit: 123 Solved: 91
[ Submit][ Status][ Discuss]

Description

有一个无向图，共N个节点，编号1至N，共M条边。FJ在节点1，它想到达节点N。FJ总是会选择最短路径到达节点N。作为捣蛋的奶牛Bessie，它想尽量延迟FJ到达节点N的时间，于是Bessie决定从M条边之中选择某一条边，使得改边的长度变成原来的两倍，由于智商的问题，Bessie不知道选择加倍哪条边的长度才能使得FJ到达N号节点的时间最迟。注意：不管Bessie选择加倍哪条边的长度，FJ总是会从1号节点开始走最短路径到达N号点。

Input

第一行，两个整数N和M. 1 <=N<=250, 1<=M<=250000。

接下来有M行，每行三个整数：A，B，L，表示节点A和节点B之间有一条长度为L的无向边。1<=L<=1000000。

Output

一个整数。Bessie选择了加倍某一条边的长度后，奶牛FJ从节点1到达节点N的最短路径是多少。但是输出的格式有变化，假设Bessie没有加倍某一条边的长度之前，FJ从1号节点到达N号节点的最短路径是X；在Bessie加倍某一条边的长度之后，FJ从1号节点到达N号节点的最短路径是Y，那么你输出的结果是Y-X。

Sample Input

5 7
2 1 5
1 3 1
3 2 8
3 5 7
3 4 3
2 4 7
4 5 2
INPUT DETAILS: There are 5 fields and 7 pathways. Currently, the shortest path from the house (field 1) to the barn (field 5) is 1-3-4-5 of total length 1+3+2=6.

Sample Output

2
（把节点3到节点4的边从原来的长度3变成长度6）

HINT

Source

Gold

[ Submit][ Status][ Discuss]

题解：堆优化dijkstra

一遍扫出最短路，并记录下路径。只有路径上的点发生变化才会影响最短路，然后之间枚举路径上的边，然后计算出此时的最短路，更新答案即可。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<queue>
#define N 300
#define M 500003
#define LL long long
#define pa pair<LL,int>
using namespace std;
int n,a[N],tot,c[N],cnt,vis[M],m;
int pre[M],next[M],v[M],point[N];
LL dis[N],len[M];
void add(int x,int y,LL z)
{
	tot++; next[tot]=point[x]; point[x]=tot; v[tot]=y; len[tot]=z;
	tot++; next[tot]=point[y]; point[y]=tot; v[tot]=x; len[tot]=z;
}
void dijkstra()
{
	priority_queue<pa,vector<pa>,greater<pa> > p;
	memset(dis,127/3,sizeof(dis));
	dis[1]=0; p.push(make_pair(dis[1],1));
	while (!p.empty()){
		int now=p.top().second; p.pop();
		for (int i=point[now];i!=-1;i=next[i])
		{
		 LL t=0;
		 if (vis[i]) t=len[i]*2LL;
		 else t=len[i];
		 if (dis[v[i]]>dis[now]+t){
		 	dis[v[i]]=dis[now]+t; pre[v[i]]=i;
		 	p.push(make_pair(dis[v[i]],v[i]));
		 }
	    }
	}
}
void solve()
{
	int now=n;
	while (now!=1) {
		c[++cnt]=pre[now];
		now=v[pre[now]^1];
	}
}
int main()
{
	freopen("rblock.in","r",stdin);
	freopen("rblock.out","w",stdout);
	scanf("%d%d",&n,&m);
	tot=-1;
	memset(point,-1,sizeof(point));
	for (int i=1;i<=m;i++){
		int x,y; LL z; scanf("%d%d%I64d",&x,&y,&z);
		add(x,y,z);
	}
	dijkstra(); LL ans=dis[n]; LL maxn=0;
	solve();
	for (int i=1;i<=cnt;i++) {
		vis[c[i]]=1; vis[c[i]^1]=1;
		dijkstra();  maxn=max(maxn,dis[n]-ans);
		vis[c[i]]=0; vis[c[i]^1]=0;
	}
	printf("%I64d\n",maxn);
}