非朴素最短路

非朴素最短路

次短路

struct node{
    int id,len,flag;
    bool operator<(const node &a)const{
        return len>a.len;
    }
};
void dijkstra(int s){
    //dis[i][0]表示从s到i的最短路径长度
    //dis[i][1]表示从s到i的次短路径长度
    for(int i=0;i<N;i++)dis[i][0]=1e9,dis[i][1]=1e9;
    priority_queue<node>q;
    q.push(node{s,0,0}),dis[s][0]=0;
    while(!q.empty()){
        int u=q.top().id,t=q.top().flag;
        q.pop();
        if(used[u])continue;
        used[u]=true;
        for(int i=Last[u];i;i=Next[i]){
            int v=End[i];
            if(dis[v][0]>dis[u][t]+Len[i]){
                dis[v][1]=dis[v][0];
                dis[v][0]=dis[u][t]+Len[i];
                q.push(node{v,dis[v][1],1});
                q.push(node{v,dis[v][0],0});
            }else if(dis[v][1]>dis[u][t]+Len[i]){
                dis[v][1]=dis[u][t]+Len[i];
                q.push(node{v,dis[v][1],1});
            }
        }
    }
}

多维最短路

分层图问题。

例一

话说puga来到了神秘岛。他发现，这是一个由n个小岛和一个中心岛组成的群岛，群岛之间有m座桥。令他感到惊讶的是，这些桥并不是固定不变的，经较长时间的观察，发现它们会随时间作周期性的变化（即桥的两端会不断更换）。

puga用望远镜看到到远远的那个中心岛上有一间小屋，架在一棵好大好大的树上。于是他决定前往中心岛上的那间空中楼阁。puga当然希望越早到越好，那么，你能帮帮他们吗？

为方便计算，puga把小岛按1…n编号，0表示中心岛。puga一开始在编号为1的小岛上。在岛上行走的时间忽略不计，过桥的时间为1个单位。岛上的桥变化的周期为T，在n*T+i(n=0,1,2,…;i=1,2,…,T)时刻岛上的桥为第i种状态，一开始的时刻为1。两个小岛间可能有多条桥相连。在任一时刻，puga可以选择过桥，也可以原地不动。当然，如果无桥可过，puga只能在原地等待。

这道题相当于每一个时间都有一张图，如果把每一张图按照时间叠起来，那么相邻两层图之间的路就是对应的时间点（比如第1秒的p1和第2秒的p1相连，第1秒的p2和第2秒的p2相连…）

这道题比朴素的最短路多的那一维就是时间。其实代码并没有什么大的差别。

#include<iostream>
#include<cstdio>
#include<queue>
#define N 100
using namespace std;
int n,m,T;
int dis[N][N];
bool used[N][N],Map[N][N][N];
struct node{int id,ti;};
void SPFA(){
    for(int i=0;i<N;i++)
        for(int j=0;j<N;j++)dis[i][j]=1e9;
    queue<node>q;
    q.push(node{1,0}),dis[1][0]=0,used[1][0]=true;
    while(!q.empty()){
        int u=q.front().id,t=q.front().ti;
        q.pop(),used[u][t]=false;
        for(int v=0;v<=n;v++){
            int tt=(t+1)%T;
            if(!Map[u][v][t])continue;
            if(dis[u][t]+1<dis[v][tt]){
                dis[v][tt]=dis[u][t]+1;
                if(!used[v][tt]){
                    used[v][tt]=true;
                    q.push(node{v,tt});
                }
            }
        }
    }
}
int main(){
    int Ans=1e9;
    scanf("%d%d%d",&n,&m,&T);
    for(int i=0;i<T;i++)
        for(int j=1;j<=m;j++){
            int x,y;scanf("%d%d",&x,&y);
            Map[x][y][i]=1,Map[y][x][i]=1;
        }
    for(int i=0;i<=n;i++)
        for(int j=0;j<T;j++)Map[i][i][j]=1;
    SPFA();
    for(int i=0;i<T;i++)Ans=min(Ans,dis[0][i]);
    if(Ans==1e9)puts("Poor puga!");
    else printf("%d",Ans);
    return 0;
} 

例二

每天,农夫John需要经过一些道路去检查牛棚N里面的牛.

农场上有M(1<=M<=50,000)条双向泥土道路,编号为1..M. 道路i连接牛棚P1_i和P2_i (1 <= P1_i <= N; 1 <= P2_i<= N).

John需要T_i (1 <= T_i <= 1,000,000)时间单位用道路i从P1_i走到P2_i或者从P2_i
走到P１_i 他想更新一些路经来减少每天花在路上的时间.具体地说,他想更新K (1 <= K <=
20)条路经，将它们所须时间减为０．帮助FJ选择哪些路经需要更新使得从１到N的时间尽量少.

最短路+动规。

f[u][j]表示从起点s到u号点更新j条路径所需的最短时间。f[i][j]=min{ f[k][j] + Len[i][k] , f[k][j-1] }(i k之间有路).

#include<iostream>
#include<cstdio>
#include<queue>
#include<cstring>
#define N 100050
using namespace std;
//dis[i][j]表示从起点到i点，修改j条路径的最短用时
//used[i][j]表示从起点到i点，修改j条路径的情况是否已经被用过
int dis[N][50],used[N][50],n,m,k;
int Len[N],Next[N],Last[N],End[N],cnt;
struct node{
    int id,v,k;
    node(int a,int b,int c){id=a,v=b,k=c;}
    bool operator<(const node &a)const{
        return v>a.v;
    }
};
void Ins(int x,int y,int l){
    End[++cnt]=y;
    Len[cnt]=l;
    Next[cnt]=Last[x];
    Last[x]=cnt;
}
void dijkstra(int s){
    memset(dis,0x3f,sizeof(dis));
    priority_queue<node>q;
    q.push(node(s,0,0));dis[s][0]=0;
    while(!q.empty()){
        int u=q.top().id,t=q.top().k;
        q.pop();
        if(used[u][t])continue;
        used[u][t]=true;
        for(int i=Last[u];i;i=Next[i]){
            int v=End[i];
            if(dis[v][t]>dis[u][t]+Len[i]){
                dis[v][t]=dis[u][t]+Len[i];
                q.push(node(v,dis[v][t],t));
            }
            if(t+1<=k && dis[v][t+1]>dis[u][t]){
                dis[v][t+1]=dis[u][t];
                q.push(node(v,dis[v][t+1],t+1));
            }
        }
    }
}
int main(){
    scanf("%d%d%d",&n,&m,&k);
    for(int i=1;i<=m;i++){
        int x,y,z;
        scanf("%d%d%d",&x,&y,&z);
        Ins(x,y,z),Ins(y,x,z);
    }
    dijkstra(1);
    printf("%d",dis[n][k]);
    return 0;
}