畅通工程续
Time Limit: 3000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submission(s): 37269 Accepted Submission(s): 13759
Problem Description
某省自从实行了很多年的畅通工程计划后,终于修建了很多路。不过路多了也不好,每次要从一个城镇到另一个城镇时,都有许多种道路方案可以选择,而某些方案要比另一些方案行走的距离要短很多。这让行人很困扰。
现在,已知起点和终点,请你计算出要从起点到终点,最短需要行走多少距离。
现在,已知起点和终点,请你计算出要从起点到终点,最短需要行走多少距离。
Input
本题目包含多组数据,请处理到文件结束。
每组数据第一行包含两个正整数N和M(0<N<200,0<M<1000),分别代表现有城镇的数目和已修建的道路的数目。城镇分别以0~N-1编号。
接下来是M行道路信息。每一行有三个整数A,B,X(0<=A,B<N,A!=B,0<X<10000),表示城镇A和城镇B之间有一条长度为X的双向道路。
再接下一行有两个整数S,T(0<=S,T<N),分别代表起点和终点。
每组数据第一行包含两个正整数N和M(0<N<200,0<M<1000),分别代表现有城镇的数目和已修建的道路的数目。城镇分别以0~N-1编号。
接下来是M行道路信息。每一行有三个整数A,B,X(0<=A,B<N,A!=B,0<X<10000),表示城镇A和城镇B之间有一条长度为X的双向道路。
再接下一行有两个整数S,T(0<=S,T<N),分别代表起点和终点。
Output
对于每组数据,请在一行里输出最短需要行走的距离。如果不存在从S到T的路线,就输出-1.
Sample Input
3 3 0 1 1 0 2 3 1 2 1 0 2 3 1 0 1 1 1 2
Sample Output
2 -1
思路:
求两点之间的最短路径 ,当然会选择迪杰斯特拉算法
代码[Dijkstra]:
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define size 205
#define INF 0x7fffffff//16进制
int g[size][size];//邻接矩阵 存储图
bool vis[size];//标记函数
int dis[size];//记录最短路径
int n,m;//n表示点的数量 m表示边的数量
void init()
{
for(int i=0;i<n;i++)//初始化整个图
for(int j=0;j<n;j++)
if(i==j)//自身与自身的距离为0
g[i][j]=0;
else
g[i][j]=INF;
}
void Dijkstra(int s,int e)//s起点 e终点 点s 到 点e的最短路径
{
for(int i=0;i<n;i++)//初始化,开始时最短路径dis[]为INF,vis[]为false
{
dis[i]=g[s][i];
vis[i]=false;
}
dis[s]=0;//将起始点s的最短路径的距离记录为0
vis[s]=true;//并且标记访问到的当前点(s);
for(int i=0;i<n;i++)//查找最短路径
{
int minnum=INF;
int temp=s;
for(int j=0;j<n;j++)
if(!vis[j]&&dis[j]<minnum)//寻找一个最短距离,并且记录该点的信息
{
temp=j;
minnum=dis[j];
}
if(minnum==INF)//如果当前的最短距离仍为INF表示此次访问无效 PS这里可有可无 均能AC
break ;
vis[temp]=true;//访问一次,标记一次
for(int j=0;j<n;j++)
if(!vis[j]&&g[temp][j]<INF&&dis[temp]+g[temp][j]<dis[j])//对于每个与点temp相邻接的点j,均访问一次,并且更新最短路径
dis[j]=dis[temp]+g[temp][j];
}
if(dis[e]==INF)//如果到达终点时最短路径仍为INF表示无最短路径,反之则有
printf("-1\n");
else
printf("%d\n",dis[e]);
}
int main()
{
while(~scanf("%d%d",&n,&m))
{
init();//
for(int i=0;i<m;i++)
{
int a,b,x;
scanf("%d%d%d",&a,&b,&x);
if(x<g[a][b])//注意,此处需将原来初始化的(a,b,INF)改为(a,b,c);
g[a][b]=g[b][a]=x;
}
int s,t;
scanf("%d%d",&s,&t);
Dijkstra(s,t);
}
return 0;
}
代码二[Floyd]
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define INF 0xfffffff
#define size 205
int g[size][size];
int n,m;
void Init()
{
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
if(i==j)//自身与自身的距离为0
g[i][j]=0;
else
g[i][j]=INF;
}
void Floyd(){
int i,j,k;
for(k=0;k<n;k++)
for(i=0;i<n;i++)
for(j=0;j<n;j++)
if(g[i][j]>g[i][k]+g[k][j])
g[i][j]=g[i][k]+g[k][j];
}
int main()
{
while(~scanf("%d%d",&n,&m))
{
Init();
for(int i=0;i<m;i++)
{
int a,b,x;
scanf("%d%d%d",&a,&b,&x);
if(x<g[a][b])
g[a][b]=g[b][a]=x;
}
int s,t;
scanf("%d%d",&s,&t);
Floyd();
if(g[s][t]==INF)
printf("-1\n");
else
printf("%d\n",g[s][t]);
}
return 0;
}
代码三[SPFA 思路与Dijkstra 很相似]
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<queue>
#include<algorithm>
using namespace std;
#define INF 0xfffffff
#define size 205
int g[size][size];
int dis[size];
bool vis[size];
int n,m;
void Init()
{
for(int i=0;i<n;i++)
for(int j=0;j<n;j++)
if(i==j)
g[i][j]=0;
else
g[i][j]=INF;
}
int SPFA(int s,int e)
{
for(int i=0;i<n;i++)
{
dis[i]=INF;
vis[i]=false;
}
queue<int> Q;
while(!Q.empty())
Q.pop();
dis[s]=0;
vis[s]=true;
Q.push(s);
int temp=0;
while(!Q.empty())
{
temp=Q.front();
Q.pop();
if(dis[temp]==INF)
break;
vis[temp]=false;
for(int i=0;i<n;i++)
if(dis[i]>dis[temp]+g[temp][i])//注意此处与Dijkstra有点不同
{
dis[i]=dis[temp]+g[temp][i];
if(!vis[i])//在此处进行检验
{
vis[i]=true;
Q.push(i);
}
}
}
if(dis[e]==INF)
printf("-1\n");
else
printf("%d\n",dis[e]);
}
int main()
{
while(~scanf("%d%d",&n,&m))
{
Init();
for(int i=0;i<m;i++)
{
int a,b,x;
scanf("%d%d%d",&a,&b,&x);
if(x<g[a][b])
g[a][b]=g[b][a]=x;
}
int s,t;
scanf("%d%d",&s,&t);
SPFA(s,t);
}
return 0;
}
代码四[你会发现Dijkstra从这上面经过,貌似在严蔚敏的数据结构见过哦]
#include<cstdio>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<algorithm>
using namespace std;
#define INF 0xfffffff
#define size 205
int dis[size];
int n,m,cnt;
struct node
{
int u,v;
int w;
}edge[1005*2];//注意此处
void addedge(int u,int v,int w)
{
edge[cnt].u=u;
edge[cnt].v=v;
edge[cnt].w=w;
cnt++;
edge[cnt].u=v;
edge[cnt].v=u;
edge[cnt].w=w;
cnt++;
}
void Bellman_Ford(int s,int e)
{
for(int i=0;i<n;i++)
dis[i]=INF;
dis[s]=0;
for(int k=0;k<n;k++)
for(int i=0;i<cnt;i++)
if(dis[edge[i].u]!=INF && dis[edge[i].v]>dis[edge[i].u]+edge[i].w)
dis[edge[i].v]=dis[edge[i].u]+edge[i].w;
if(dis[e]==INF)
printf("-1\n");
else
printf("%d\n",dis[e]);
}
int main()
{
while(~scanf("%d%d",&n,&m))
{
cnt=0;
for(int i=0;i<m;i++)
{
int a,b,x;
scanf("%d%d%d",&a,&b,&x);
addedge(a,b,x);
}
int s,t;
scanf("%d%d",&s,&t);
Bellman_Ford(s,t);
}
return 0;
}
Bellmen_Ford算法的运算时间为O(nm)适合n^2>m的稀疏图
当遇到n^2<m的稠密图时最好用Dijkstra算法,其运算时间为O(n^2)。