图论最短路

图论最短路算法是一种家喻户晓的算法，并且一直受到各种出题人的青睐，在此简单介绍下几种常见算法
先看看个奇怪的概念(十分钟重要)：
松弛
以上图为例，黑点到红点有条路，你可以把它看成黑点和红点的连线，也可以看成黑红间一些点的连线，总之，已知求出的最短路长度即为7
这时候，我们引入一点(绿点)，发现比起原来求出的，通过从黑点用绿点中转到达红点，结果更优，这个操作，就叫做松弛

多源最短路：Floyd

多源最短路，即需要求图中任意两个点的最短路径，这种最短路目前貌似只有Floyd比较好用
这个算法的实现和思想都比较简单，是一种基于动规的思想
首先令点 i 到点 j 的最短路径为dis[ i ][ j ]，若 i 和 j 不连通，可以先设
dis[ i ][ j ]=INF(一个很大的值即可)
初始化后进入主题部分，我们每次枚举三个数:k,i,j，看看是原来求出的
dis[ i ][ j ]更优，还是通过k点中转，即dis[ i ][ k ]+dis[ k ][ j ]更优，这样就可以得出我们的转移方程，即dis[ i ][ j ]=min{ dis[ i ][ j ],dis[ i ][ k ]+dis[ k ][ j ] }
于是乎，核心代码也是相当的简单

#include<bits/stdc++.h>
using namespace std;
const int N=110;
int n,m;
struct node{
	int x,y;
}a[N];
int dis[N][N];//邻接矩阵存图
void Floyed()//核心代码
{
	for(int k=1;k<=n;k++)
		for(int i=1;i<=n;i++)
			for(int j=1;j<=n;j++)
				if(i!=j && j!=k)
					if(dis[i][k]+dis[k][j]<dis[i][j])
						dis[i][j]=dis[i][k]+dis[k][j];
}
int main()
{
	cin>>n;
	for(int i=1;i<=n;i++)
	{
		int x,y;
		scanf("%d%d",&x,&y);
		a[i].x=x;
		a[i].y=y;
	}
	memset(dis,1<<30,sizeof(dis));//初始化，赋一个极大值
	cin>>m;
	for(int i=1;i<=m;i++)
	{
		int x,y,u;
		cin>>x>>y>>u;
		dis[x][y]=dis[y][x]=u;
	}
	int s,t;
	cin>>s>>t;//求点s到点t的最短路
	Floyed();
	cout<<dis[s][t]<<endl;
	return 0;	
} 

时间复杂度显而易见为O(n^3)
Floyd虽然粗暴了点，但是用起来很方便，预处理后直接O(1)拿出来用，十分方便

单源最短路

单源最短路，即求指定的一点到其它点的最短路

法一：Bellman-Ford算法

由于本人对这个算法运用不是很熟练，所以就暂时跳过了，有兴趣的话可以去看看别的大佬的讲解

法二：Dijkstra

迪杰在OIer是最受欢迎的最短路算法，理解起来可能有点难，但是它的速度很快(优化后)，且较之SPFA，它很难被卡掉
初始化：起点s，dis[ s ]=0，其它赋一个很大的值
算法流程：
①在没被访问的节点中，找一个u，使得dis[ u ]最小
②将它标记为已知的最短路
③对所有和u相连的点做松弛操作
上面的文字略显枯燥，我们可以动手做一做

比如我们有这么一个图，要求单源最短路，起点为1
初始化：dis={0,∞,∞,∞,∞}

我们首先找到dis最小的，即起点，将它标记
然后发现，它可以对2，3，4进行松弛
修改：dis[2]=2,dis[3]=4,dis[4]=7
dis={0,2,4,7,∞}

接着，在未标记的点，找到dis最小的点2并标记，发现它可以对3，5进行松弛
修改：dis[3]=3,dis[5]=4
dis={0,2,3,7,4}

然后找到3并标记，发现它可以对4，5进行松弛，但由于数值比原来大，所以对5的松弛失败
改变：dis[4]=4
dis={0,2,3,4,4}
简单来说，就是每一次找到距离最小的点，然后对它相连的点进行松弛

#include<bits/stdc++.h>
using namespace std;
const int N=1e4+10;
int n,m,s,dis[N],vis[N];
struct node{
    int to,u;
};
vector<node> v[N];//邻接表存图
void Dijkstra()//算法主体
{
    for(int i=1;i<=n;i++)
    {
        int k=0;
        for(int j=1;j<=n;j++)//寻找距离最小点
            if(vis[j]==0 && dis[j]<dis[k])
                k=j;
	    vis[k]=1;//标记
        for(int j=0;j<v[k].size();j++)//寻找相连的点
        {
            int to=v[k][j].to,u=v[k][j].u;
            if(dis[k]+u<dis[to])//松弛
                dis[to]=dis[k]+u;
        }
    }
}
int main()
{
    cin>>n>>m>>s;
    v[s].push_back((node){s,0});
    for(int i=0;i<=n;i++)
    	dis[i]=INT_MAX;
    dis[s]=0;
    for(int i=1;i<=m;i++)
    {
        int x,y,u;
        scanf("%d%d%d",&x,&y,&u);
        v[x].push_back((node){y,u});
        v[y].push_back((node){x,u});
    }
    Dijkstra();
    for(int i=1;i<=n;i++)
        cout<<dis[i]<<' ';
}

算法复杂度为O(n^2),不能解决有边权为负数的图

Dijkstra的优化

其实，每次都循环找最小值十分浪费，c++提供给我们一种数据结构，堆，当我们用这个数据结构在线查询时，复杂度会小很多
改完后变成这样的了

#include<bits/stdc++.h>
using namespace std;
const int N=1e5+10;
int n,m,s,dis[N],vis[N];
struct node{
    int to,u;
};
vector<node> v[N];
priority_queue<pair<int,int> > q;//优先队列(堆)，pair的第一关键字存dis的值的相反数(以此保证小的在堆顶，实现小根堆，第二关键字是节点的编号
void Dijkstra()
{
    q.push(make_pair(0,s));//压入初始节点
    while(q.size())
    {
        int x=q.top().second;//查询最小的，并存下它的节点编号
        q.pop();
        if(vis[x]==1) continue;
        vis[x]=1;
        for(int i=0;i<v[x].size();i++)
        {
            int to=v[x][i].to,u=v[x][i].u;
            if(dis[to]>dis[x]+u)//松弛
            {
                dis[to]=dis[x]+u;
                q.push(make_pair(-dis[to],to));//将松弛后的信息压入堆
            }
        }
    }
    
}
int main()
{
    cin>>n>>m>>s;
    v[s].push_back((node){s,0});
    for(int i=0;i<=n;i++)
    	dis[i]=INT_MAX;
    dis[s]=0;
    for(int i=1;i<=m;i++)
    {
        int x,y,u;
        scanf("%d%d%d",&x,&y,&u);
        v[x].push_back((node){y,u});
        v[y].push_back((node){x,u});
    }
    Dijkstra();
    for(int i=1;i<=n;i++)
        cout<<dis[i]<<' ';
}

创作不易，求赞