码迷,mamicode.com
首页 > 编程语言 > 详细

次短路算法

时间:2020-06-14 16:38:23      阅读:68      评论:0      收藏:0      [点我收藏+]

标签:return   else   const   min   isp   span   入队   cout   ref   

技术图片

次短路

(一)
  • \(u\)(父节点)到\(v\)(子节点)次短路直接更新(通常在最短路已经确定的情况下才进行直接更新次短路)
  • \(u\)(父节点)到\(v\)(子节点)最短路不更新,但是距离比次短路距离小,更新次短路
  • \(u\)(父节点)到\(v\)(子节点)最短路更新,原来的最短路就成了次短路

数组\(dis1\)存放的是到各点的最短路,\(dis2\)存放的是到各点的次短路

\[dis_2[v]=\begin{cases} dis_2[u]+w \dis_1[u]+w \ \ \ \ \ \ \ \ \ dis_1[u]+w<dis_1[v]\dis_1[v] \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ dis_1[v]>dis_1[u]+w \end{cases}\]

(二)换边法
  • 1到n的次短路长度必然产生于:从\(1\)走到\(x\)的最短路 + \(edge[x][y]\) + \(y\)\(n\)的最短路(与次小生成树有些相似)
    1.首先预处理好1到每一个节点的最短路,和n到每一个节点的最短路
    2.然后枚举每一条边作为中间边\((x,y)\)或者\((y,x)\),如果加起来长度等于最短路长度则跳过,否则更新。
    \(1\)走到\(x\)的最短路 + \(edge[x][y]\) + \(y\)\(n\)的最短路,只比\(dist[n]\)大,比其他的都小的那个值。

Roadblocks

Spfa

  • 一直更新记录,找出的是各点的最短路和次短路
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
int dis[5500];
int dis2[5500];
int vis[5500];
int head[200100],cnt;
struct node
{
    int v,w;
    int nxt;
} edge[200100];
void add(int x,int y,int w)
{
    edge[++cnt].nxt=head[x];
    edge[cnt].v=y;
    edge[cnt].w=w;
    head[x]=cnt;
}
void Dijkstra(int x)
{
    queue<int>qu;
    memset(vis,0,sizeof(vis));
    memset(dis,INF,sizeof(dis));
    memset(dis2,INF,sizeof(dis2));
    dis[x]=0;
    vis[x]=1;
    qu.push(x);
    while(!qu.empty( ))
    {
        int u=qu.front( );
        qu.pop( );
        vis[u]=0;
        for(int i=head[u]; i!=-1; i=edge[i].nxt)
        {
            int v=edge[i].v;
            int w=edge[i].w;
            if(dis[v]>dis[u]+w)//最短路更新,原来的最短路变成次短路
            {
                dis2[v]=dis[v];
                dis[v]=dis[u]+w;
                if(!vis[v])
                {
                    qu.push(v);
                    vis[v]=1;
                }
            }
            if(dis2[v]>dis2[u]+w)//次短路直接更新(通常是最短路确定后才进行)
            {
                dis2[v]=dis2[u]+w;
                if(!vis[v])
                {
                    qu.push(v);
                    vis[v]=1;
                }
            }
            if(dis2[v]>dis[u]+w&&dis[v]<dis[u]+w)//最短路不更新,但距离比次短路小
            {
                dis2[v]=dis[u]+w;
                if(!vis[v])
                {
                    qu.push(v);
                    vis[v]=1;
                }
            }
        }
    }
}
int main( )
{
    int n,r,x,y,w;
    memset(head,-1,sizeof(head));
    cnt=0;
    scanf("%d%d",&n,&r);
    for(int i=1; i<=r; i++)
    {
        scanf("%d%d%d",&x,&y,&w);
        add(x,y,w);
        add(y,x,w);
    }
    Dijkstra(1);
    for(int i=1;i<=n;i++)
    {
        printf("%d ",dis[i]);
    }
    cout<<endl;
    for(int i=1;i<=n;i++)
    {
        printf("%d ",dis2[i]);
    }
    return 0;
}

利用dijkstra求出\(1\)到各点的最小距离,\(n\)到各点的距离,然后枚举每条边,找出第二小的值

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
const int M=200010;
int head[M],cnt;
struct node
{
    int v,w,nxt;
}edge[M];
int dis1[5050];
int disn[5050];
int vis[5050];
void add(int x,int y,int w)
{
    edge[++cnt].nxt=head[x];
    edge[cnt].v=y;
    edge[cnt].w=w;
    head[x]=cnt;
}
void Dijkstra(int x,int dis[ ])
{
    memset(vis,0,sizeof(vis));
    priority_queue<pair<int,int>,vector<pair<int,int> >,greater<pair<int,int> > >qu;
    dis[x]=0;
    qu.push(make_pair(0,x));
    while(!qu.empty( ))
    {
        int u=qu.top( ).second;
        qu.pop( );
        if(!vis[u])
        {
            vis[u]=1;
            for(int i=head[u];i!=-1;i=edge[i].nxt)
            {
                int v=edge[i].v;
                if(dis[v]>dis[u]+edge[i].w)
                {
                    dis[v]=dis[u]+edge[i].w;
                    qu.push(make_pair(dis[v],v));
                }
            }
        }
    }
}
int main( )
{
    int n,r;
    scanf("%d%d",&n,&r);
    memset(head,-1,sizeof(head));
    cnt=0;
    int x,y,w;
    for(int i=1;i<=r;i++)
    {
        scanf("%d%d%d",&x,&y,&w);
        add(x,y,w);
        add(y,x,w);
    }
    memset(dis1,INF,sizeof(dis1));
    memset(disn,INF,sizeof(disn));
    Dijkstra(1,dis1);
    Dijkstra(n,disn);
    // for(int i=1;i<=n;i++)
    // {
    //     cout<<dis1[i]<<" ";
    // }
    // cout<<endl;
    // for(int i=1;i<=n;i++)
    // {
    //     cout<<disn[i]<<" ";
    // }
    // cout<<endl;
    int ans=dis1[n];
    int temp=INF;
    for(int i=1;i<=n;i++)
    {
        for(int j=head[i];j!=-1;j=edge[j].nxt)
        {
            int v=edge[j].v;
            int w=dis1[i]+edge[j].w+disn[v];
            if(w==ans)
            {
                continue;
            }
            else
            {
                temp=min(temp,w);
            }
            
        }
    }
    cout<<temp<<endl;
    return 0;
}

一边更新最短路,一边更新次短路

#include<iostream>
#include<cstdio>
#include<cstring>
#include<queue>
#define INF 0x3f3f3f3f
using namespace std;
const int M=200000+10;
int head[M],cnt;
struct node
{
    int v,w,nxt;
}edge[M];
int dis1[5005];
int dis2[5005];
void add(int x,int y,int w)
{
    edge[++cnt].nxt=head[x];
    edge[cnt].v=y;
    edge[cnt].w=w;
    head[x]=cnt;
}
void Dijkstra(int x)
{
    memset(dis1,INF,sizeof(dis1));
    memset(dis2,INF,sizeof(dis2));
    priority_queue<pair<int,int>,vector<pair<int,int> >,greater<pair<int,int> > >qu;
    dis1[x]=0;//最短路初始值为0,次短路无穷大
    qu.push(make_pair(0,x));
    while(!qu.empty( ))
    {
        int w=qu.top( ).first;//弹出最小值,或许是最短路,或许是次短路
        int u=qu.top( ).second;
        qu.pop( );
        if(dis2[u]<w)//弹出来的值比当前的次短路大,就可以跳过这个
        {
            continue;
        }
        for(int i=head[u];i!=-1;i=edge[i].nxt)
        {
            int v=edge[i].v;
            int cost=w+edge[i].w;//u到v的花费
            if(dis1[v]>cost)//花费大于原来的最小值,更新最短路
            {
                swap(dis1[v],cost);//交换值
                qu.push(make_pair(dis1[v],v));//压入队列
            }
            if(dis2[v]>cost&&cost>dis1[v])//交换次短路
            {
                swap(dis2[v],cost);
                qu.push(make_pair(dis2[v],v));//压入队列,之所以次短路要压入队列是因为后面更新需要。
                                             //例子:dis[2] = 10, dis2[2] = 20 有一条边 2 到 6 的边权值为 5                 
                                            //如果不把 dis2 入队,那么之后的算法中 dis[6] = 15, dis2[6] = INF              
                                           //只有当队列里有 20 这个值,才能 20+5 得出 25,然后更新 dis2[6] = 25 
            }
        }
    }
}
int main( )
{
    int n,r;
    scanf("%d%d",&n,&r);
    memset(head,-1,sizeof(head));
    cnt=0;
    int x,y,w;
    for(int i=1;i<=r;i++)
    {
        scanf("%d%d%d",&x,&y,&w);
        add(x,y,w);
        add(y,x,w);
    }
    Dijkstra(1);
    printf("%d\n",dis2[n]);
    return 0;
}

次短路算法

标签:return   else   const   min   isp   span   入队   cout   ref   

原文地址:https://www.cnblogs.com/lcbwwy/p/13125183.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!