思路:
小路大路分开存储。小路跑一遍Floyd算法,并将连续小路平方后等价于带小路标记的直达路,并入大路中。但之后用Dijkstra算法时将不可叠加计算,即小路->小路为不允许情况,直接跳过即可。
本算法有漏洞过不了自己编的例子,之后再做修改,不过测评已经100分了。
#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll INF=1e18;
const int MAXV=501;
int n,m;
int pre[MAXV];
ll dis[MAXV][MAXV];
ll G[MAXV][MAXV];
int tag[MAXV][MAXV];
ll d[MAXV];
bool vis[MAXV];
inline void Floyd(){
for(register int i=1;i<=n;i++)dis[i][i]=0;
for(register int k=1;k<=n;k++){
for(register int i=1;i<=n;i++){
for(register int j=1;j<=n;j++){
if(dis[i][k]!=INF&&dis[k][j]!=INF&&dis[i][k]+dis[k][j]<dis[i][j]){
dis[i][j]=(dis[i][k]+dis[k][j]);
}
}
}
}
for(register int i=1;i<=n;i++){ //将连续小路等价于带小路标记的直达路,之后用Dijkstra算法时将不可叠加计算
for(register int j=1;j<=n;j++){
if(dis[i][j]!=INF){
dis[i][j]*=dis[i][j];
if(dis[i][j]<G[i][j]){
G[i][j]=dis[i][j];
tag[i][j]=tag[j][i]=1;
}
}
}
}
}
inline void Dijkstra(int s){
fill(d,d+MAXV,INF);
for(int i=1;i<=n;i++){
pre[i]=i;
}
d[s]=0;
for(int i=1;i<=n;i++){
ll u=-1,MIN=INF;
for(int j=0;j<=n;j++){
if(vis[j]==false&&d[j]<MIN){
u=j;
MIN=d[j];
}
}
if(u==-1)return;
vis[u]=true;
for(int v=1;v<=n;v++){
if(pre[u]!=u&&tag[pre[u]][u]==1&&tag[u][v]==1)continue; //小路->小路 不允许跳过
if(vis[v]==false&&G[u][v]!=INF&&d[u]+G[u][v]<d[v]){ // 路 + 路 优化
d[v]=d[u]+G[u][v];
pre[v]=u;
}
}
}
}
int main(){
//输入
ios::sync_with_stdio(false);
fill(dis[0],dis[0]+MAXV*MAXV,INF);
fill(G[0],G[0]+MAXV*MAXV,INF);
cin>>n>>m;
int t,a,b,c;
for(int i=0;i<m;i++){
cin>>t>>a>>b>>c;
if(!t){
if(c<G[a][b]){
G[a][b]=G[b][a]=c;
tag[a][b]=tag[b][a]=t;
}
}
else {
if(c<dis[a][b]){
dis[a][b]=dis[b][a]=c;
}
}
}
//处理
Floyd();
Dijkstra(1);
//输出
cout<<d[n];
return 0;
}
样例:
/*
样例输入
6 7
1 1 2 3
1 2 3 2
0 1 3 30
0 3 4 20
0 4 5 30
1 3 5 6
1 5 6 1
4 4
0 1 3 10
1 1 2 2
1 2 3 1
1 3 4 1
*/