BZOJ 2407: 探险/4398: 福慧双修

2407: 探险

Description

探险家小T好高兴!X国要举办一次溶洞探险比赛,获奖者将得到丰厚奖品哦!小T虽然对奖品不感兴趣,但是这个大振名声的机会当然不能错过!
比赛即将开始,工作人员说明了这次比赛的规则:每个溶洞和其他某些溶洞有暗道相连。两个溶洞之间可能有多条道路,也有可能没有,但没有一条暗道直接从自己连到自己。参赛者需要统一从一个大溶洞出发,并再次回到这个大溶洞。
如果就这么点限制,那么问题就太简单了,可是举办方又提出了一个条件:不能经过同一条暗道两次。这个条件让大家犯难了。这该怎么办呢?
到了大溶洞口后,小T愉悦地发现这个地方他曾经来过,他还记得有哪些暗道,以及通过每条暗道的时间。小T现在向你求助,你能帮他算出至少要多少时间才能回到大溶洞吗?

Input

第一行两个数n,m表示溶洞的数量以及暗道的数量。

接下来m行,每行4个数s、t、w、v,表示一个暗道连接的两个溶洞s、t,这条暗道正着走(s à t)的所需要的时间w,倒着走(t à s)所需要的时间v。由于溶洞的相对位置不同,wv可能不同。

Output

输出一行一个数t,表示最少所需要的时间。

Sample Input

3 3
1 2 2 1
2 3 4 5
3 1 3 2

Sample Output

8

HINT

N<=10000,M<=200000,1<=W,V<=10000

思路:

  这题是考场上做的,只想出来了四十分做法,同届神犇suika原创本题一种不同于网上大家题解的方法,目前还没更新题解,我的做法是常规做法

>看这里<

我的做法是常规的最短路做法。考虑四十分 : 枚举点1周围的点,以这些点为起点向点1跑最短路,并且屏蔽掉这点到1的边,具体实现有很多种做法,在这里就不赘述了。这样做的话是O(nmlog)的(Dij)可过40分,对于满分,我们考虑这样的情况,当走出第一步之后,剩下的便是从这个点到1的最短路实现,而各个点之间会有很多重复的状态枚举,,考虑新建一张图来抛弃掉这些状态。

我们建立汇点T=n+1 预处理出从1到各个节点的最小代价和第一步分别走了哪些点并用prev数组存下,代价用一次SPFA处理,第一步传递处理就好。

处理好以下之后,枚举每条边进行加图,分成以下几种情况处理:

  1.该边从1到v,若prev[v] == v 不做处理,对答案产生的贡献已经记录到dis里了。

          若prev[v] != v  连接一条从1到v,代价为原代价的边

  2.该边从u到1,若pre[u] != u 直接用dis和该边val更新答案。

         否则 加一条从u到T的,权值为该边val的边

  3.该边从u到v,若u、v的prev相等,那么直接还原这条边

         否则的话,加一条从1到v,代价为f[u]+val的边

  对新图进行SPFA(DIJ)f[t]即为答案

代码如下

#include <cstdio>
#include <cstring>
#include <algorithm>
#include <iostream>
#include <queue>
const int inf=0x3f3f3f3f;
const int N = 41000, M = 400010;
using namespace std;
int f[N],vis[N],val[M],head[N],to[M],next[M],cnt;
int pre[N];
int ra[M],rb[M],rc[M];
void add_edge(int a,int b,int c) {
	to[++cnt] = b;
	next[cnt] = head[a];
	head[a] = cnt;
	val[cnt] = c;
}
queue<int>q;
void spfa() {
	while(!q.empty()) {
		int u=q.front();
		q.pop();
		vis[u] = 0;
		for(int i=head[u];i;i=next[i]) {
			if(f[to[i]]>f[u]+val[i]) {
				pre[to[i]]=pre[u];
				f[to[i]]=f[u]+val[i];
				if(!vis[to[i]]) {
					vis[to[i]]=1;
					q.push(to[i]);
				}
			}
		}
	}
}
int main() {
	int n, m;
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++) {
		int a,b,c,d;
		scanf("%d%d%d%d",&a,&b,&c,&d);
		add_edge(a,b,c);
		ra[cnt]=a;
		rb[cnt]=b;
		rc[cnt]=c;
		add_edge(b,a,d);
		ra[cnt]=b;
		rb[cnt]=a;
		rc[cnt]=d;
	}
	memset(f,0x3f,sizeof(f));
	f[1]=0;
	for(int i=head[1];i;i=next[i]) {
		pre[to[i]]=to[i];
		f[to[i]]=val[i];
		vis[to[i]]=1;
		q.push(to[i]);
	}
	spfa();
	memset(head,0,sizeof(head));
	cnt=0;
	int ans=0x3f3f3f3f;
	for(int i=1;i<=2*m;i++) {
		if(ra[i]==1) {
			if(pre[rb[i]]!=rb[i]) add_edge(1,rb[i],rc[i]);
		}
		else if(rb[i]==1) {
			if(ra[i]!=pre[ra[i]]) ans=min(ans,f[ra[i]]+rc[i]);
			else add_edge(ra[i],n+1,rc[i]);
		}
		else {
			if(pre[ra[i]]==pre[rb[i]]) {
				add_edge(ra[i],rb[i],rc[i]);
			}
			else {
				add_edge(1,rb[i],f[ra[i]]+rc[i]);
			}
		}
	}
	memset(f,0x3f,sizeof f);
	f[1]=0;
	q.push(1);
	spfa();
	ans=min(ans,f[n+1]);
	printf("%d",ans==inf?-1:ans);

}

 

 欢迎来原博客看看 >原文链接<

 

posted @ 2018-05-29 20:27  TOBICHI_ORIGAMI  阅读(105)  评论(0编辑  收藏  举报