数据结构实验之图论七：驴友计划

/*

测试数据 教科书 P189 G6 的邻接矩阵 其中 数字 1000000 代表无穷大

6

1000000 1000000 10 100000 30 100

1000000 1000000 5 1000000 1000000 1000000

1000000 1000000 1000000 50 1000000 1000000

1000000 1000000 1000000 1000000 1000000 10

1000000 1000000 1000000 20 1000000 60

1000000 1000000 1000000 1000000 1000000 1000000

结果：

D[0]   D[1]   D[2]   D[3]   D[4]   D[5]

  0   1000000   10     50     30     60

*/

#include <iostream>

#include <cstdio>

#define MAX 1000000

using  namespace  std;

int  arcs[10][10]; //邻接矩阵

int  D[10]; //保存最短路径长度

int  p[10][10]; //路径

int  final[10]; //若final[i] = 1则说明 顶点vi已在集合S中

int  n = 0; //顶点个数

int  v0 = 0; //源点

int  v,w;

void  ShortestPath_DIJ()

{

      for  (v = 0; v < n; v++)  //循环 初始化

      {

           final[v] = 0; D[v] = arcs[v0][v];

           for  (w = 0; w < n; w++) p[v][w] = 0; //设空路径

           if  (D[v] < MAX) {p[v][v0] = 1; p[v][v] = 1;}

      }

      D[v0] = 0; final[v0]=0;  //初始化 v0顶点属于集合S

      //开始主循环 每次求得v0到某个顶点v的最短路径 并加v到集合S中

      for  ( int  i = 1; i < n; i++)

      {

           int  min = MAX;

           for  (w = 0; w < n; w++)

           {

                //我认为的核心过程--选点

                if  (!final[w])  //如果w顶点在V-S中

                {

                     //这个过程最终选出的点 应该是选出当前V-S中与S有关联边

                     //且权值最小的顶点 书上描述为 当前离V0最近的点

                     if  (D[w] < min) {v = w; min = D[w];}

                }

           }

           final[v] = 1;  //选出该点后加入到合集S中

           for  (w = 0; w < n; w++) //更新当前最短路径和距离

           {

                /*在此循环中 v为当前刚选入集合S中的点

                则以点V为中间点 考察 d0v+dvw 是否小于 D[w] 如果小于 则更新

                比如加进点 3 则若要考察 D[5] 是否要更新 就 判断 d(v0-v3) + d(v3-v5) 的和是否小于D[5]

                */

                if  (!final[w] && (min+arcs[v][w]<D[w]))

                {

                     D[w] = min + arcs[v][w];

                    // p[w] = p[v];

                     p[w][w] = 1;  //p[w] = p[v] +　[w]

                }

           }

      }

}

 

 

int  main()

{

     cin >> n;

     for  ( int  i = 0; i < n; i++)

     {

          for  ( int  j = 0; j < n; j++)

          {

               cin >> arcs[i][j];

          }

     }

     ShortestPath_DIJ();

     for  ( int  i = 0; i < n; i++)  printf ( "D[%d] = %d\n" ,i,D[i]);

     return  0;

}

 

/*用迪杰斯特拉算法求有向网G的V0顶点到其他顶点的最短路径P，以及其带权长度D。其中P是二维数组，行号表示终点，列号表示经过的路径。P[v][w]为TRUE的意思就是从v0到v,要经过w点)。D是一维数组，表示某顶点到v0点的路径长(D[v] == 10表示从v0到v要经过的路径长度为10。final存放已经求得的路径结果(比如final[v]为TRUE表示已经找到v0到v的最短路径)。*/
void  ShorttestPath_DIJ(  MGraph  G,  int  v0,  PathMatrix  &P,  ShortPathTable  &D)
{
        for(  v  =  0;  v  <  G.vexnum;  ++v  )  
       {
                final[v]  =  FALSE; 
                D[v]  =  G.arcs[v0][v];
                for(  w  =  0;  w  <  G.vexnum;  ++w  )  
                {
                        P[v][w]  =  FALSE;
                }
                if(  D[v]  <  INFINITY  )  
                {  //如果有直接互通的两个顶点，直接将这个路径赋值到数组P[v]。
                        P[v][v0]  =  TRUE;
                        P[v][v]  =  TRUE;
                }
        }
        D[v0]  =  0;  final[v]  =  TRUE;

       /*下面开始主循环，每次求得v0到某个v顶点的最短路径，同时刷新之前的最短路径。*/
        for(  i  =  1;  i  <  G.vexnum;  ++i  )  
        {  //  对于除了v0之外的顶点(这个循环仅仅限制次数，i的值不用).
                min  =  INFINITY;  //  假定初始的“最小值”为无穷大。
                for(  w  =  0;  w  <  G.vexnum;  ++w  )  
                {
                        if(  !final[w]  )  //  w顶点在V - S中，即还未确定的顶点。  
                                if(  D[w]  <  min  )  
                                {
                                        v  =  w;
                                        min  =  D[w];  //  随着循环进行，依与v0的距离大小，从小到大取得顶点v，并标记进final。
                                }
                 }  
                final[v]  =  TRUE;  //  标记已经找到
                for(  w  =  0;  w  <  G.vexnum;  w++  )  
                {  //  更新路径
                        if(  !final[w]  &&  (min  +  G.arcs[v][w]  <  D[w]) )  
                        {
                                D[w]  =  min  +  G.arcs[v][w];
                                P[w]  =  P[v];  //  把一行都给赋值了
                                P[w][w]  =  TRUE;
                         }
                }
        }

/---------------------------------------------------------------------------------------------------------------------------------------------------------/

数据结构实验之图论七：驴友计划

Time Limit: 1000MS  Memory Limit: 65536KB

Submit Statistic

Problem Description

做为一个资深驴友，小新有一张珍藏的自驾游线路图，图上详细的标注了全国各个城市之间的高速公路距离和公路收费情况，现在请你编写一个程序，找出一条出发地到目的地之间的最短路径，如果有多条路径最短，则输出过路费最少的一条路径。

Input

连续T组数据输入，每组输入数据的第一行给出四个正整数N,M，s，d，其中N(2 <= N <= 500)是城市数目，城市编号从0～N-1，M是城市间高速公路的条数，s是出发地的城市编号，d是目的地的城市编号；随后M行，每行给出一条高速公路的信息，表示城市1、城市2、高速公路长度、收费额，中间以空格间隔，数字均为整数且不超过500，输入数据均保证有解。 

Output

在同一行中输出路径长度和收费总额，数据间用空格间隔。 

Example Input

1
4 5 0 3
0 1 1 20
1 3 2 30
0 3 4 10
0 2 2 20
2 3 1 20

Example Output

3 40

Hint

  

[plain]  view plain  copy

 

1. #include<stdio.h>  
2. #include<math.h>  
3. #define max 9999999  
4. int n,m,s,d;  
5.   
6. struct node  
7. {  
8.     int length;  
9.     int money;  
10. }road[500][500];  
11.   
12. //初始化邻接矩阵图  
13. void InitG()  
14. {  
15.     for(int i=0;i<n;i++)  
16.     {  
17.         for(int j=0;j<n;j++)  
18.         {  
19.              
20.             road[i][j].length=max;//两点之间的长度和路费全部初始化为无穷大  
21.             road[i][j].money=max;  
22.              
23.         }  
24.     }  
25.   
26. }  
27.   
28. //寻找最短路  
29. //此函数运行完成后，任意两个点之间的路线长度均为最短  
30. void findroad()  
31. {  
32.     for(int i=0;i<n;i++)  
33.     {  
34.         for(int j=0;j<n;j++)  
35.         {  
36.             for(int k=0;k<n;k++)  
37.             {  
38.                 if(road[k][j].length>road[k][i].length+road[i][j].length)  
39.                 {  
40.                     road[k][j].length=road[k][i].length+road[i][j].length;  
41.                     road[k][j].money=road[k][i].money+road[i][j].money;  
42.                 }  
43.                 else if(road[k][j].length==road[k][i].length+road[i][j].length)//路线长度相同还要比较价格  
44.                 {  
45.                     if(road[k][j].money>road[k][i].money+road[i][j].money)  
46.                     {  
47.                         road[k][j].money=road[k][i].money+road[i][j].money;  
48.                     }  
49.                 }  
50.             }  
51.         }  
52.     }  
53.     printf("%d %d\n",road[s][d].length,road[s][d].money);  
54. }  
55.   
56. int main()  
57. {  
58.     int t;  
59.     scanf("%d",&t);  
60.     while(t--)  
61.     {  
62.         scanf("%d %d %d %d",&n,&m,&s,&d);  
63.         InitG();  
64.         for(int i=0;i<m;i++)  
65.         {  
66.             int a,b,l,m;  
67.             scanf("%d %d %d %d",&a,&b,&l,&m);  
68.             road[a][b].length=road[b][a].length=l;  
69.             road[a][b].money=road[b][a].money=m;  
70.         }  
71.         findroad();  
72.     }  
73.     return 0;  
74. }  

[plain]  view plain  copy

 

1. //迪杰斯特拉算法11  
2. #include <stdio.h>  
3. #include <algorithm>  
4. #include <string.h>  
5. using namespace std;  
6. const int inf=0x3f3f3f3f;  
7. int mp[600][600][2],n;//0 luchag 1lufei  
8. int dist[600][2],vis[600];  
9. void dij(int s)  
10. {  
11.     int i,j;  
12.     memset(vis,0,sizeof(vis));  
13.     vis[s]=1;  
14.     for(i=0;i<n;i++)  
15.     {  
16.         dist[i][0]=mp[s][i][0];  
17.         dist[i][1]=mp[s][i][1];  
18.     }  
19.     int max1,max2,u=s;  
20.     for(i=1;i<n;i++)  
21.     {  
22.         max1=max2=inf;  
23.         for(j=0;j<n;j++)  
24.         {  
25.             if(!vis[j] && ((dist[j][0]<max1) || (dist[j][0]==max1 && dist[j][1]<max2)) )  
26.             {  
27.                 max1=dist[j][0];max2=dist[j][1];  
28.                 u=j;  
29.             }  
30.         }  
31.         vis[u]=1;  
32.         for(j=0;j<n;j++)  
33.         {  
34.             if(!vis[j] && ((dist[j][0]>mp[u][j][0]+dist[u][0]) || (dist[j][0]==mp[u][j][0]+dist[u][0] && dist[j][1]>mp[u][j][1]+dist[u][1])) )  
35.             {  
36.                 dist[j][0]=mp[u][j][0]+dist[u][0];dist[j][1]=mp[u][j][1]+dist[u][1];  
37.             }  
38.         }  
39.     }  
40. }  
41.   
42. int main()  
43. {  
44.     int m,u,v,w,h;  
45.     int begin1,end1;  
46.     int t;  
47.     scanf("%d",&t);  
48.     while(t--)  
49.     {  
50.         scanf("%d%d%d%d",&n,&m,&begin1,&end1);  
51.         for(int i=0;i<n;i++)  
52.         {  
53.             for(int j=0;j<n;j++)  
54.             {  
55.                 mp[i][j][0]=mp[i][j][1]=inf;  
56.             }  
57.         }  
58.         while(m--)  
59.         {  
60.             scanf("%d%d%d%d",&u,&v,&w,&h);  
61.             mp[u][v][0]=w;  
62.             mp[v][u][0]=w;  
63.             mp[u][v][1]=h;  
64.             mp[v][u][1]=h;  
65.         }  
66.         dij(begin1);  
67.         printf("%d %d\n",dist[end1][0],dist[end1][1]);  
68.     }  
69.     return 0;  
70. }