Floyd算法

intvexnum,arcnum;//图的当前顶点数和弧数
GraphKindkind;//图的种类标志
}MGraph;
void find(int P[MAX_VERTEX_NUM][MAX_VERTEX_NUM][MAX_VERTEX_NUM],MGraph G,int a,int b);
void main(){
printf("%c到%c最短路径为",65+a,65+b);
printf("%c\t",65+a);
find(P,G,a,b);
printf("%c\t",65+b);
printf("长度为%d",D[a][b]);
printf("\n");
}
}
void find(int P[MAX_VERTEX_NUM][MAX_VERTEX_NUM][MAX_VERTEX_NUM],MGraphG,int a,int b){
G.kind=DG;
printf("请输入邻接矩阵\n");
for (v = 0; v < G.vexnum; v++)
for (w = 0; w < G.vexnum; w++)
scanf("%d",&G.arcs[v][w]); //读入邻接矩阵
// P[v][w][k]为TRUE，则从v到w的最短路径中含有k节点
D[v][w] = D[v][k] + D[k][w];
for (i = 0; i < G.vexnum; i++)
P[v][w][i] = P[v][k][i] || P[k][w][i];
}
for(a=0; a<G.vexnum; a++)
for(b=0; b<G.vexnum; b++)
if(D[a][b] < INFINITY && a!=b){
问题描述
实现Floyd算法，并求所示有向图中各顶点之间的最短路径及其长度。
算法思想
采用图的邻接矩阵存储，实现Floyd算法~，数组P[][][]存储是否存在中间点使长度缩短。
设计描述
数据存储结构类型的定义：
typedef struct MGraph{
char vexs[MAX_VERTEX_NUM];
#define TRUE 1
#define FALSE 0
typedef enum{DG, DN, UDG, UDN}
GraphKind; //四种图类型
typedef struct MGraph{
charvexs[MAX_VERTEX_NUM];//顶点向量
intarcs[MAX_VERTEX_NUM][MAX_VERTEX_NUM];//邻接矩阵
MGraph G;
intD[MAX_VERTEX_NUM][MAX_VERTEX_NUM],P[MAX_VERTEX_NUM][MAX_VERTEX_NUM][MAX_VERTEX_NUM];
int v,w,k,a,b,i;
printf("请输入顶点数和弧数");
scanf("%d %d",&G.vexnum,&G.arcnum);
if (D[v][w] < INFINITY)
P[v][w][v] = P[v][w][w] = TRUE;
}
for (k = 0; k < G.vexnum; k++)
for (v = 0; v < G.vexnum; v++)
for (w = 0; w < G.vexnum; w++)
if (D[v][k] + D[k][w] < D[v][w]){
问题主要出在输出打印方面
Voidfind(int P[MAX_VERTEX_NUM][MAX_VERTEX_NUM][MAX_VERTEX_NUM],MGraph G,int a,int b){
int k;
for(k = 0; k < G.vexnum; k++)
if(P[a][b][k]==TRUE && k!=a && k!=b){百度文库
find(P,G,a,k);
printf("\t%c",65+k);
find(P,G,k,b);
}
}
采用递归的思想一次寻找是否存在中间点然后打印出来。
问题出在红色部分，判定是否采用递归，未考虑k是否与a b相同，结果导致无限递归。从而发现stack overflow的错误提示有可能出于递归无法跳出，导致栈的溢出问题。
int arcs[MAX_VERTEX_NUM][MAX_VERTEX_NUM];
int vexnum,arcnum;
GraphKind kind;
}MGraph;
源程序
#include<stdlib.h>
#include<stdio.h>
#define INFINITY 1000 //最大值
#define MAX_VERTEX_NUM 20 //最大顶点个数
// D[v][w]从v到w的最短路径的长度
for (v = 0; v < G.vexnum; v++)
for (w = 0; w < G.vexnum; w++){
D[v][w] = G.arcs[v][w];
for (k = 0; k < G.vexnum; k++)
P[v][w][k] = FALSE;
int k;
for(k = 0; k < G.vexnum; k++)
if(P[a][b][k]==TRUE && k!=a && k!=b){
find(P,G,a,k);
printf("%c\t",65+k);
find(P,G,k,b);
}
}
测试结果
输入
（1000为无穷！~
输出
心得体会~。！！
懂得了floyd算法的思想，用邻接矩阵存储带权值的图。