-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathShortestPath.c
263 lines (216 loc) · 5.59 KB
/
ShortestPath.c
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
/***********************************************************************
************************************************************************
author:wangkang qq:417301568 e-mail: wangkang806@163.com
begin time:2010-5-30
end time:
************************************************************************
************************************************************************/
#include <stdio.h>
#define VEXNUM_MAX 20 //定义最大顶点数。
#define QUANZHI_NOT_REACHE 10000 //定义权值不可达时的值。
#define CityNameLength 20
typedef int Quanzhi;
typedef struct{
int value;
char info[CityNameLength];//城市名字
}Vex,*Vex_point;
typedef struct{
Vex v[VEXNUM_MAX]; //顶点数组
Quanzhi quan[VEXNUM_MAX][VEXNUM_MAX]; //权值
int vexnum;//顶点数
int quannum;//权的数量 ,即弧数
} *Graph_point, Graph;
/***********************************************************************
************************************************************************
创建邻接矩阵
************************************************************************
************************************************************************/
void createJuZhen(Graph_point m)
{
int i,j,k,citynum,x,y;
Quanzhi z;
char c[CityNameLength];
printf("请你输入顶点数目: ");
scanf("%d",&i);
m->vexnum=i;
//给矩阵初值。
for(i=0;i<m->vexnum;i++)
{
for(j=0;j<m->vexnum;j++)
{
if(i==j) m->quan[i][j]=0;
else m->quan[i][j]=QUANZHI_NOT_REACHE;
}
}
printf("请你输入权的数目: ");
scanf("%d",&j);
m->quannum=j;
for(citynum=0;citynum<m->vexnum;citynum++)
{
printf("请输入第%d点顶点代表的城市:",citynum);
scanf("%s",c);
for(k=0;k<CityNameLength;k++)
{
m->v[citynum].info[k]=c[k];
}
m->v[citynum].value=citynum;
}
for(i=0;i<m->quannum;i++)
{
printf("请输入两个顶点:(中间加逗号)");
scanf("%d,%d",&x,&y);
printf("输入两个顶点的权值(城市之间坐车所用费用):");
scanf("%d",&z);
m->quan[x][y]=z;
m->quan[y][x]=z;
}
//printf("%s\n",m->v[0].info);
//printf("%d",j);
}
/***********************************************************************
************************************************************************
显示邻接矩阵的
************************************************************************
************************************************************************/
void display(Graph_point m)
{
int i,j,k;
printf("\t");
for(i=0;i<m->vexnum;i++) //输出顶点
{
//printf("%8d",m->v[i].value); //输出顶点
printf("%8s",m->v[i].info); //输出城市名字
}
printf("\n");
for(i=0;i<m->vexnum;i++)
{
//printf("%8d",m->v[i].value);
printf("%8s",m->v[i].info);
for(j=0;j<m->vexnum;j++)
{
printf("%8d",m->quan[i][j]);
}
printf("\n");
}
}
/***********************************************************************
************************************************************************
最短路径 用到了Dijstra算法。
************************************************************************
************************************************************************/
typedef struct{
Quanzhi value; //权值
int pre;//前一个顶点
}dijstra,*dijstra_point;
void shortertpath_Dijstra(Graph_point m,int vi,int vj)
{
Quanzhi min; //最短路径
dijstra dist[VEXNUM_MAX];
int s[VEXNUM_MAX]; //某个顶点是否访问过,访问过为1,没访问过为0
int i,j,k,n;
for(i=0;i<m->vexnum;i++)
{
//
s[i]=0;
if(i!=vi)
{
dist[i].value=m->quan[vi][i];
if(m->quan[vi][i]==QUANZHI_NOT_REACHE) dist[i].pre=-1;
else dist[i].pre=vi;
}
}
s[vi]=1;
//主循环
for(i=0;i<m->vexnum;i++)
{
min=QUANZHI_NOT_REACHE;
for(n=0;n<m->vexnum;n++)
{
if(s[n]==0&& dist[n].value<min)
{
min=dist[n].value;
k=n;
}
}
if(k==vj) break;//如果计算到vj顶点,就跳出循环。
//
s[k]=1;
for(j=0;j<m->vexnum;j++)
{
if(s[j]==0 && m->quan[k][j]<QUANZHI_NOT_REACHE && (dist[k].value+m->quan[k][j])<dist[j].value )
{
dist[j].value=dist[k].value+m->quan[k][j];
dist[j].pre=k;
}
}
}
//打印输出
printf("最小花费为:%d\n",dist[vj]);
printf("%d%s<---",vj,m->v[vj].info);
j=vj;
while(j!=vi &&j!=-1)
{
j=dist[j].pre;
printf("%d%s<---",j,m->v[j].info);
}
printf("\b\b\b\b");
printf(" \n");
}
void main()
{
int i,j,select,k=0;
FILE *fp;
Graph g;
Graph_point s;
s=&g;
printf("1. 创建矩阵。2.求最短路径。 3.退出。请选择你要进行的操作,输入(1,2,3): ");
scanf("%d",&select);
while(k==0)
{
switch(select)
{
case 1:
{
createJuZhen(s);
fp=fopen("short.bin","wb");
fwrite(s,sizeof(Graph),1,fp);
fclose(fp);
printf("矩阵创建成功\n");
printf("矩阵显示如下--------\n");
display(s);
printf("\n\n\n");
printf("1. 创建矩阵。2.求最短路径。 3.退出。请继续选择你要进行的操作,输入(1,2,3): ");
scanf("%d",&select);
break;
}
case 2:
{
printf("请输入起始顶点:");
scanf("%d",&i);
printf("请输入末顶点:");
scanf("%d",&j);
fp=fopen("short.bin","rb");
fread(s,sizeof(Graph),1,fp);
fclose(fp);
printf("矩阵显示如下--------\n");
display(s);
printf("\n\n\n");
printf("最短路径如下-------\n");
shortertpath_Dijstra(s,i,j);
printf("1. 创建矩阵。2.求最短路径。 3.退出。请继续选择你要进行的操作,输入(1,2,3): ");
scanf("%d",&select);
break;
}
case 3:
{
k=1;
break;
}
default:
{
printf("输入有误,请重新输入:");
scanf("%d",&select);
}
}
}
}