传送门
题解:
直接跑最短路,然后考虑翻转每一条边之后怎么快速算代价。
显然如果这条边不在 1 − n 1-n 1−n 或者 n − 1 n-1 n−1 的最短路上,直接考虑强行走过这条边的最短路就行了,可以通过记录最短路和次短路(要求前驱不同的最短路)来实现。
否则,这样的边只有 O ( n ) O(n) 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;}