[BZOJ2407]探险——最短路、贪心

[BZOJ2407]探险

终于做到有链接的题了…

题解

题目要求进过点1的最小环，并且这个环边不重复。

考虑一个简化的问题，即边的两个方向边权相等，那么我们可以考虑建出以1为根的最短路树。此时的最小环必定是树上的两个lca等于1的点到1的路径长度和+两点间的最小边权（非树边）。如果我们通过枚举非树边来统计答案的话，容易发现这样虽然会统计漏一些环，但是一定不会漏掉最小环。

然后考虑更复杂的问题，即边的两个方向边权不一定相等。这个时候显然不能只建一棵最短路树，需要正反方向各建一棵。此时我们仍然可以枚举 $(u,v)$，计算第一棵树上 $u$ 到1的距离与第二棵树上 $v$ 到1的距离和，加上边权。这时只需要满足 $u,v$ 的最短路的第一条边不相同且不等于枚举的边即可，因为即使满足条件后两条路径实际上有重复边，那也一定存在一个更优的合法解。

实现精细一点就可以涵盖掉走两条共端点但不同的边的情况。

代码

空间别开太大，不然Checker会T

#include<bits/stdc++.h>//JZM yyds!!
#define ll long long
#define uns unsigned
#define IF (it->first)
#define IS (it->second)
#define END putchar('\n')
using namespace std;
const int MAXN=10005;
const ll INF=1e18;
inline ll read(){
	ll x=0;bool f=1;char s=getchar();
	while((s<'0'||s>'9')&&s>0){if(s=='-')f^=1;s=getchar();}
	while(s>='0'&&s<='9')x=(x<<1)+(x<<3)+(s^48),s=getchar();
	return f?x:-x;
}
int ptf[50],lpt;
inline void print(ll x,char c='\n'){
	if(x<0)putchar('-'),x=-x;
	ptf[lpt=1]=x%10;
	while(x>9)x/=10,ptf[++lpt]=x%10;
	while(lpt)putchar(ptf[lpt--]^48);
	if(c>0)putchar(c);
}
inline ll lowbit(ll x){return x&-x;}

struct edge{
	int v,to,w;edge(){}
	edge(int V,int T,int W){v=V,to=T,w=W;}
}e[400005];
int EN=1,G[MAXN];
inline void addedge(int u,int v,int w1,int w2){
	e[++EN]=edge(v,G[u],w1),G[u]=EN;
	e[++EN]=edge(u,G[v],w2),G[v]=EN;
}
struct itn{
	int s,id;itn(){}
	itn(int S,int D){s=S,id=D;}
	inline bool operator<(const itn&b)const{return s>b.s;}
}f[MAXN],g[MAXN];
#define pti pair<itn,int>
#define fi first
#define se second
priority_queue<pti>q;
int n,m;
signed main()
{
	n=read(),m=read();
	for(int i=1;i<=m;i++){
		int u=read(),v=read(),w1=read(),w2=read();
		addedge(u,v,w1,w2);
	}
	for(int i=1;i<=n;i++)f[i]=g[i]=itn(1e9,0);
	f[1]=g[1]=itn(0,0);
	while(!q.empty())q.pop();
	q.push(pti(f[1],1));
	while(!q.empty()){
		int u=q.top().se;
		itn a=q.top().fi;
		q.pop();
		if(a<f[u])continue;
		for(int i=G[u];i;i=e[i].to){
			int v=e[i].v,w=e[i].w;
			itn to=itn(a.s+w,a.id);
			if(!a.id)to.id=i>>1;
			if(f[v]<to)f[v]=to,q.push(pti(f[v],v));
		}
	}
	while(!q.empty())q.pop();
	q.push(pti(g[1],1));
	while(!q.empty()){
		int u=q.top().se;
		itn a=q.top().fi;
		q.pop();
		if(a<g[u])continue;
		for(int i=G[u];i;i=e[i].to){
			int v=e[i].v,w=e[i^1].w;
			itn to=itn(a.s+w,a.id);
			if(!a.id)to.id=i>>1;
			if(g[v]<to)g[v]=to,q.push(pti(g[v],v));
		}
	}
	int ans=1e9;
	for(int x=1;x<=n;x++)
		for(int i=G[x];i;i=e[i].to){
			int v=e[i].v;
			if((i>>1)==f[x].id||(i>>1)==g[v].id||f[x].id==g[v].id)
				continue;
			ans=min(ans,f[x].s+g[v].s+e[i].w);
		}
	print(ans);
	return 0;
}