第十三周项目—算法验证(2)
来源:互联网 发布:软件售后服务流程图 编辑:程序博客网 时间:2024/05/18 00:12
- Copyright (c) 2017, 烟台大学计算机学院
- *All rights reserved.
- *作 者:张行
- *完成日期:2017年12月30日
- *版 本 号:v1.0
- *问题描述:1:Dijkstra算法的验证
- 2:Floyd算法验证
- 3:拓扑排序算法验证
- *结果显示:如图所示
main.cpp
#include <stdio.h>#include <malloc.h>#include "graph.h"#define MaxSize 100//Dijkstra算法的验证void Ppath(int path[],int i,int v) //前向递归查找路径上的顶点{ int k; k=path[i]; if (k==v) return; //找到了起点则返回 Ppath(path,k,v); //找顶点k的前一个顶点 printf("%d,",k); //输出顶点k}void Dispath(int dist[],int path[],int s[],int n,int v){ int i; for (i=0; i<n; i++) if (s[i]==1) { printf(" 从%d到%d的最短路径长度为:%d\t路径为:",v,i,dist[i]); printf("%d,",v); //输出路径上的起点 Ppath(path,i,v); //输出路径上的中间点 printf("%d\n",i); //输出路径上的终点 } else printf("从%d到%d不存在路径\n",v,i);}void Dijkstra(MGraph g,int v){ int dist[MAXV],path[MAXV]; int s[MAXV]; int mindis,i,j,u; for (i=0; i<g.n; i++) { dist[i]=g.edges[v][i]; //距离初始化 s[i]=0; //s[]置空 if (g.edges[v][i]<INF) //路径初始化 path[i]=v; else path[i]=-1; } s[v]=1; path[v]=0; //源点编号v放入s中 for (i=0; i<g.n; i++) //循环直到所有顶点的最短路径都求出 { mindis=INF; //mindis置最小长度初值 for (j=0; j<g.n; j++) //选取不在s中且具有最小距离的顶点u if (s[j]==0 && dist[j]<mindis) { u=j; mindis=dist[j]; } s[u]=1; //顶点u加入s中 for (j=0; j<g.n; j++) //修改不在s中的顶点的距离 if (s[j]==0) if (g.edges[u][j]<INF && dist[u]+g.edges[u][j]<dist[j]) { dist[j]=dist[u]+g.edges[u][j]; path[j]=u; } } Dispath(dist,path,s,g.n,v); //输出最短路径}//Floyd算法验证void Ppath(int path[][MAXV],int i,int j) //前向递归查找路径上的顶点{ int k; k=path[i][j]; if (k==-1) return; //找到了起点则返回 Ppath(path,i,k); //找顶点i的前一个顶点k printf("%d,",k); Ppath(path,k,j); //找顶点k的前一个顶点j}void Dispath(int A[][MAXV],int path[][MAXV],int n){ int i,j; for (i=0; i<n; i++) for (j=0; j<n; j++) { if (A[i][j]==INF) { if (i!=j) printf("从%d到%d没有路径\n",i,j); } else { printf(" 从%d到%d=>路径长度:%d 路径:",i,j,A[i][j]); printf("%d,",i); //输出路径上的起点 Ppath(path,i,j); //输出路径上的中间点 printf("%d\n",j); //输出路径上的终点 } }}void Floyd(MGraph g){ int A[MAXV][MAXV],path[MAXV][MAXV]; int i,j,k; for (i=0; i<g.n; i++) for (j=0; j<g.n; j++) { A[i][j]=g.edges[i][j]; path[i][j]=-1; } for (k=0; k<g.n; k++) { for (i=0; i<g.n; i++) for (j=0; j<g.n; j++) if (A[i][j]>A[i][k]+A[k][j]) { A[i][j]=A[i][k]+A[k][j]; path[i][j]=k; } } Dispath(A,path,g.n); //输出最短路径}int main(){ MGraph g; int A[7][7]= { {0,4,6,6,INF,INF,INF}, {INF,0,1,INF,7,INF,INF}, {INF,INF,0,INF,6,4,INF}, {INF,INF,2,0,INF,5,INF}, {INF,INF,INF,INF,0,INF,6}, {INF,INF,INF,INF,1,0,8}, {INF,INF,INF,INF,INF,INF,0} }; ArrayToMat(A[0], 7, g); printf("Dijkstra序列:\n"); Dijkstra(g,0); printf(" Floyd序列:\n"); Floyd(g); return 0;}
graph.cpp
#include <stdio.h>#include <malloc.h>#include "graph.h"//功能:由一个反映图中顶点邻接关系的二维数组,构造出用邻接矩阵存储的图//参数:Arr - 数组名,由于形式参数为二维数组时必须给出每行的元素个数,在此将参数Arr声明为一维数组名(指向int的指针)// n - 矩阵的阶数// g - 要构造出来的邻接矩阵数据结构void ArrayToMat(int *Arr, int n, MGraph &g){ int i,j,count=0; //count用于统计边数,即矩阵中非0元素个数 g.n=n; for (i=0; i<g.n; i++) for (j=0; j<g.n; j++) { g.edges[i][j]=Arr[i*n+j]; //将Arr看作n×n的二维数组,Arr[i*n+j]即是Arr[i][j],计算存储位置的功夫在此应用 if(g.edges[i][j]!=0 && g.edges[i][j]!=INF) count++; } g.e=count;}void ArrayToList(int *Arr, int n, ALGraph *&G){ int i,j,count=0; //count用于统计边数,即矩阵中非0元素个数 ArcNode *p; G=(ALGraph *)malloc(sizeof(ALGraph)); G->n=n; for (i=0; i<n; i++) //给邻接表中所有头节点的指针域置初值 G->adjlist[i].firstarc=NULL; for (i=0; i<n; i++) //检查邻接矩阵中每个元素 for (j=n-1; j>=0; j--) if (Arr[i*n+j]!=0) //存在一条边,将Arr看作n×n的二维数组,Arr[i*n+j]即是Arr[i][j] { p=(ArcNode *)malloc(sizeof(ArcNode)); //创建一个节点*p p->adjvex=j; p->info=Arr[i*n+j]; p->nextarc=G->adjlist[i].firstarc; //采用头插法插入*p G->adjlist[i].firstarc=p; } G->e=count;}void MatToList(MGraph g, ALGraph *&G)//将邻接矩阵g转换成邻接表G{ int i,j; ArcNode *p; G=(ALGraph *)malloc(sizeof(ALGraph)); for (i=0; i<g.n; i++) //给邻接表中所有头节点的指针域置初值 G->adjlist[i].firstarc=NULL; for (i=0; i<g.n; i++) //检查邻接矩阵中每个元素 for (j=g.n-1; j>=0; j--) if (g.edges[i][j]!=0) //存在一条边 { p=(ArcNode *)malloc(sizeof(ArcNode)); //创建一个节点*p p->adjvex=j; p->info=g.edges[i][j]; p->nextarc=G->adjlist[i].firstarc; //采用头插法插入*p G->adjlist[i].firstarc=p; } G->n=g.n; G->e=g.e;}void ListToMat(ALGraph *G,MGraph &g)//将邻接表G转换成邻接矩阵g{ int i,j; ArcNode *p; g.n=G->n; //根据一楼同学“举报”改的。g.n未赋值,下面的初始化不起作用 g.e=G->e; for (i=0; i<g.n; i++) //先初始化邻接矩阵 for (j=0; j<g.n; j++) g.edges[i][j]=0; for (i=0; i<G->n; i++) //根据邻接表,为邻接矩阵赋值 { p=G->adjlist[i].firstarc; while (p!=NULL) { g.edges[i][p->adjvex]=p->info; p=p->nextarc; } }}void DispMat(MGraph g)//输出邻接矩阵g{ int i,j; for (i=0; i<g.n; i++) { for (j=0; j<g.n; j++) if (g.edges[i][j]==INF) printf("%3s","∞"); else printf("%3d",g.edges[i][j]); printf("\n"); }}void DispAdj(ALGraph *G)//输出邻接表G{ int i; ArcNode *p; for (i=0; i<G->n; i++) { p=G->adjlist[i].firstarc; printf("%3d: ",i); while (p!=NULL) { printf("-->%d/%d ",p->adjvex,p->info); p=p->nextarc; } printf("\n"); }}graph.h
#ifndef GRAPH_H_INCLUDED#define GRAPH_H_INCLUDED#define MAXV 100 //最大顶点个数#define INF 32767 //INF表示∞typedef int InfoType;//以下定义邻接矩阵类型typedef struct{ int no; //顶点编号 InfoType info; //顶点其他信息,在此存放带权图权值} VertexType; //顶点类型typedef struct //图的定义{ int edges[MAXV][MAXV]; //邻接矩阵 int n,e; //顶点数,弧数 VertexType vexs[MAXV]; //存放顶点信息} MGraph; //图的邻接矩阵类型//以下定义邻接表类型typedef struct ANode //弧的结点结构类型{ int adjvex; //该弧的终点位置 struct ANode *nextarc; //指向下一条弧的指针 InfoType info; //该弧的相关信息,这里用于存放权值} ArcNode;typedef int Vertex;typedef struct Vnode //邻接表头结点的类型{ Vertex data; //顶点信息 int count; //存放顶点入度,只在拓扑排序中用 ArcNode *firstarc; //指向第一条弧} VNode;typedef VNode AdjList[MAXV]; //AdjList是邻接表类型typedef struct{ AdjList adjlist; //邻接表 int n,e; //图中顶点数n和边数e} ALGraph; //图的邻接表类型//功能:由一个反映图中顶点邻接关系的二维数组,构造出用邻接矩阵存储的图//参数:Arr - 数组名,由于形式参数为二维数组时必须给出每行的元素个数,在此将参数Arr声明为一维数组名(指向int的指针)// n - 矩阵的阶数// g - 要构造出来的邻接矩阵数据结构void ArrayToMat(int *Arr, int n, MGraph &g); //用普通数组构造图的邻接矩阵void ArrayToList(int *Arr, int n, ALGraph *&); //用普通数组构造图的邻接表void MatToList(MGraph g,ALGraph *&G);//将邻接矩阵g转换成邻接表Gvoid ListToMat(ALGraph *G,MGraph &g);//将邻接表G转换成邻接矩阵gvoid DispMat(MGraph g);//输出邻接矩阵gvoid DispAdj(ALGraph *G);//输出邻接表G#endif // GRAPH_H_INCLUDED
阅读全文
0 0
- 第十三周--项目--验证算法(2)
- 第十三周项目-验证算法(2)
- 第十三周项目—算法验证(2)
- 第十三周实践项目—算法验证(2)Kruskal算法的验证
- 第十三周项目-验证算法(2-Kruskal算法)
- 第十三周——【项目 - 验证算法】
- 第十三周——【项目 - 验证算法】
- 第十三周项目2Kruskal算法验证
- 第十三周项目2 Kruskal算法验证
- 第十三周项目1-验证算法2
- 第十三周--项目--验证算法(1)
- 第十三周--项目--验证算法(3)
- 第十三周项目-验证算法(1)
- 第十三周项目-验证算法(3)
- 第十三周项目-验证算法(4)
- 第十三周项目-验证算法(5)
- 第十三周项目2—Kruskal算法的验证
- 第十三周项目2—Dijkstra算法的验证
- 什么叫身心灵?
- 从Ubuntu 14.04 LTS版升级到Ubuntu 16.04 LTS
- C++标志粘贴##怎么使用
- android中输入一行字符,分别统计出其中英文字母、空格、数字和其它字符的个数
- 第九周项目三 利用二叉树遍历解决思想问题
- 第十三周项目—算法验证(2)
- 明德扬FPGA-培训班公开课-第05期-《边缘检测工程讲解》
- 状态机生成工具——SMC
- [文摘]信息熵
- javaWeb应用的生命周期
- IntelliJ Idea 常用快捷键列表
- 第十一周-项目4-1-应用图的深度优先遍历思路求解问题
- vue 快捷轮播组件代码实现
- LeetCode--122. Best Time to Buy and Sell Stock II