Dijkstra 查找最短路径
Dijkstra是一种解决 单源最短路径 的算法,其可扩展性较强,扩展后适用于大多数最短路径问题,但不能解决带负权边的问题
问题:
求0号城市到其他各城市的最短路径,结果包存在dis[ ]数组里面。
输入:
第一行为 n(城市数量)和 m(道路数量)
接下来 m 行,每行为 城市 a,城市 b,它们之间的距离
输入样例:
6 9
0 1 1
0 2 12
1 2 9
1 3 3
2 4 5
3 2 4
3 4 13
3 5 15
4 5 4
输出样例:
dis[]= 0 1 8 4 13 17
本问题我使用了 邻接矩阵 和 邻接表(vector实现)两种数据结构储存图,vector确实好用,嘿嘿
邻接矩阵(C):
#include <stdio.h>
#define INF 99999999
int main()
{
int map[50][50],dis[50],book[50];
int n,m,i,j,lenth,index,index1,index2,min;
scanf ("%d%d",&n,&m);
for (i=0;i<n;i++)//初始化矩阵和标记数组
{
for (j=0;j<n;j++)
{
if (i==j)
map[i][j]=0;
else
map[i][j]=INF;
}
book[i]=0;
}
for (i=0;i<m;i++)//输入数据
{
scanf ("%d%d%d",&index1,&index2,&lenth);
map[index1][index2]=lenth;
// map[index2][index1]=lenth;//无向图
}
for (i=0;i<n;i++)//初始化dis[]
dis[i]=map[0][i];//初始时dis[]中为顶点0到其他顶点的估计值
book[0]=1;//标记源点到源点的距离为确定值
for (i=0;i<n;i++)
{
min=INF;//找出估计值中最小值的下标
for (j=0;j<n;j++)
{
if (book[j]==0&&dis[j]<min)
{
min=dis[j];
index=j;//index为最小值的下标
}
}
book[index]=1;//标记为确定值
for (j=0;j<n;j++)
{//找出此点的所有出边
if (map[index][j]!=INF)
{//如果源点到 j点的距离大于源点到index点与index点到 j点的和,则缩小估计值
if (dis[j]>dis[index]+map[index][j])
dis[j]=dis[index]+map[index][j];
}
}
}
printf ("dis[]=");
for (i=0;i<n;i++)
printf (" %d",dis[i]);
return 0;
}
邻接表(C++):
#include <iostream>
#include <vector>
#define INF 99999999
using namespace std;
struct Node{
int to;//另一个顶点
int weight;//权值
};
int main()
{
vector <struct Node> map[1000];
int dis[1000],book[1000]={0};
int n,m,i,j,min,index,index1,index2,weight;
cin>>n>>m;//n个顶点,m条边
struct Node temp;
for (i=0;i<m;i++)
{
cin>>index1>>index2>>weight;
temp.weight=weight;
temp.to=index2;
map[index1].push_back(temp);
// temp.to=index1;//无向图
// map[index2].push_back(temp);
}
for (i=0;i<n;i++)//初始化dis[]数组
dis[i]=INF;
dis[0]=0;
for (i=0;i<map[0].size();i++)
dis[map[0][i].to]=map[0][i].weight;
book[0]=1;
for (i=0;i<n;i++)
{
min=INF;
index=-1;
for (j=0;j<n;j++)
{//遍历dis[]找估计值中的最小值
if (book[j]==0&&dis[j]<min)
{
min=dis[j];
index=j;
}
}
if (index==-1)
break;
book[index]=1;
for (j=0;j<map[index].size();j++)
{
if (dis[map[index][j].to]>dis[index]+map[index][j].weight)
dis[map[index][j].to]=dis[index]+map[index][j].weight;
}
}
cout<<"dis[]=";
for (i=0;i<n;i++)
cout<<" "<<dis[i];
return 0;
}