【JOI 2020 Final】【LOJ3255】奥运公交（最短路）

传送门

---

题解：

直接跑最短路，然后考虑翻转每一条边之后怎么快速算代价。

显然如果这条边不在 $1-n$ 或者 $n-1$ 的最短路上，直接考虑强行走过这条边的最短路就行了，可以通过记录最短路和次短路（要求前驱不同的最短路）来实现。

否则，这样的边只有 $O(n)$ 条，直接重新跑 Dijkstra 即可。

由于图非常稠密，不用堆优化还好一点。

---

代码：

#include<bits/stdc++.h>
#define ll long long
#define re register
#define cs const

namespace IO{
	inline char gc(){
		static cs int Rlen=1<<22|1;static char buf[Rlen],*p1,*p2;
		return (p1==p2)&&(p2=(p1=buf)+fread(buf,1,Rlen,stdin),p1==p2)?EOF:*p1++; 
	}template<typename T>T get_integer(){
		char c;bool f=false;while(!isdigit(c=gc()))f=c=='-';T x=c^48;
		while(isdigit(c=gc()))x=((x+(x<<2))<<1)+(c^48);return f?-x:x;
	}inline int gi(){return get_integer<int>();}
}using namespace IO; 

using std::cerr;
using std::cout;
using pli=std::pair<ll,int>;
#define fi first
#define se second 

cs int N=2e2+7,M=5e4+7;

int n,m;
int fr[M],to[M],c[M],d[M];
bool on_path[M]; 

struct edge{int to,w,id;};
std::vector<edge> G[N];

int ban,pre[N];
int pres[N],pret[N];
ll dis[N],_dis[N];
ll disfrs[N],_disfrs[N];
ll disfrt[N],_disfrt[N];
ll distos[N],_distos[N];
ll distot[N],_distot[N];

bool vs[N];
void Dij(int S,ll *dis,ll *sub,int *pre){
	memset(vs+1,0,sizeof(bool)*n);
	memset(dis+1,0x3f,sizeof(ll)*n);
	memset(sub+1,0x3f,sizeof(ll)*n);
	dis[S]=sub[S]=0;
	for(int re tt=1;tt<n;++tt){
		int u=0;
		for(int re i=1;i<=n;++i)
			if(!vs[i]&&(u==0||dis[i]<dis[u]))
				u=i;
		if(u==0||dis[u]>1e15)break;vs[u]=true;
		for(auto &e:G[u])if(e.id!=ban){
			if(dis[e.to]>dis[u]+e.w){
				sub[e.to]=dis[e.to];
				dis[e.to]=dis[u]+e.w;
				pre[e.to]=e.id;
			}else if(sub[e.to]>dis[u]+e.w)
				sub[e.to]=dis[u]+e.w;
		}
	}
}

void Main(){
	n=gi(),m=gi();
	for(int re i=1;i<=m;++i){
		fr[i]=gi(),to[i]=gi(),c[i]=gi(),d[i]=gi();
		G[fr[i]].push_back({to[i],c[i],i});
	}
	Dij(1,disfrs,_disfrs,pres);
	Dij(n,disfrt,_disfrt,pret);
	if(pres[n])
		for(int re u=n;u!=1;u=fr[pres[u]])
			on_path[pres[u]]=true;
	if(pret[1])
		for(int re u=1;u!=n;u=fr[pret[u]])
			on_path[pret[u]]=true;
	for(int re i=1;i<=n;++i)G[i].clear();
	for(int re i=1;i<=m;++i)
		G[to[i]].push_back({fr[i],c[i],i});
	Dij(1,distos,_distos,pre);
	Dij(n,distot,_distot,pre);
	for(int re i=1;i<=n;++i)G[i].clear();
	for(int re i=1;i<=m;++i)
		G[fr[i]].push_back({to[i],c[i],i});
	ll ans=disfrs[n]+disfrt[1];
	for(int re i=1;i<=m;++i){
		ll res=d[i];
		if(!on_path[i]){
			if(pres[to[i]]==i)
				res+=std::min(disfrs[n],
					_disfrs[to[i]]+c[i]+distot[fr[i]]);
			else 
				res+=std::min(disfrs[n],
					disfrs[to[i]]+c[i]+distot[fr[i]]);
			if(pret[to[i]]==i)
				res+=std::min(disfrt[1],
					_disfrt[to[i]]+c[i]+distos[fr[i]]);
			else 
				res+=std::min(disfrt[1],
					disfrt[to[i]]+c[i]+distos[fr[i]]);
		}else {
			ban=i;G[to[i]].push_back({fr[i],c[i],0});
			Dij(1,dis,_dis,pre),res+=dis[n];
			Dij(n,dis,_dis,pre),res+=dis[1];
			G[to[i]].pop_back();
		}ans=std::min(ans,res);
	}
	if(ans<1e15)cout<<ans<<"\n";
	else cout<<"-1\n";
}

inline void file(){
#ifdef zxyoi
	freopen("bus.in","r",stdin);
#endif
}signed main(){file();Main();return 0;}