许多实际问题在解决时都可以抽象成一张图,以图的形式来对数据进行处理(如洛谷P1135),在图上有一类操作十分常见,即求从一点到另一点的最短路。解决这个问题通常有dijkstra和spfa两种算法,在这里我们使用dijkstra.
dijkstra是基于贪心思想的最短路算法,每次都选择最短的路径对可达的点进行更新,因此可以保证每一个被到达的点都能得到最优解。
dijkstra在建立路径时采取“松弛”操作,例如:我们已经知道起点到某一个节点的距离为 n n n,那么起点到这个节点的下一个节点的距离就是 n n n+边权,这是十分自然的。
分析时间复杂度:我们需要对所有 n n n个点进行遍历,每次查询最短路径的耗时为 O ( l o g n ) O(logn) O(logn)(使用堆优化),所以总复杂度是 O ( n l o g n ) O(nlogn) O(nlogn).
代码如下(原题:洛谷P4779):
#include<cstdio>
#include<queue>
#define ma 1000000000
using namespace std;
priority_queue<pair<int,int> > a;//pair<点权值,点编号>
int n,m,s,x,y,z,dis[100005],head[100005],ecnt=1,cnt;
bool vis[100005];
struct ljb//与head[]一起构成邻接表
{//r 终点 w 权值 n 下一条边
int r,w,n;
}ljb[500005];
void add(int e,int f,int g)
{
ljb[ecnt].r=f;
ljb[ecnt].w=g;
ljb[ecnt].n=head[e];
head[e]=ecnt;
ecnt++;
//head[e]:下一次加入从该点出发的边时,该边为第e条被添加到图中的边
}
void dijkstra(int p)
{
int k=p;//在下面的程序里可以直接使用p代替k,没有影响
while(a.size())
{
k=a.top().second;//取出点编号
a.pop();//已使用的信息出堆
if(vis[k]) continue;//如果已经扫过,则该点一定有最优解,此时可以跳过这一点
vis[k]=1;//标记此点已被到达
for(int i=head[k];ljb[i].r!=0;i=ljb[i].n)//遍历每个出边
if(dis[k]+ljb[i].w<dis[ljb[i].r])//判断新到达的点是否有更优的解
{
dis[ljb[i].r]=dis[k]+ljb[i].w;//更新
a.push(make_pair(-dis[ljb[i].r],ljb[i].r));
//priority_queue默认为大根堆,用相反数法保证其顺序
}
}
}
int main()
{
scanf("%d%d%d",&n,&m,&s);
a.push(make_pair(0,1));
for(int i=1;i<=n;i++) dis[i]=ma;//除起点外的其他点距离为无穷大
dis[1]=0;//起点的距离为0
for(int i=1;i<=m;i++)
{
scanf("%d%d%d",&x,&y,&z);
add(x,y,z);
}
dijkstra(s);
for(int i=1;i<=n;i++)
printf("%d ",dis[i]);
return 0;
}