图的邻接矩阵存储
代码如下:
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<malloc.h>
//#include "HeadFile.h"
#define INFINITY 65535
#define MAX_VERTEX_NUM 20
#define MAX_NAME 3
#define MAX_INFO 20
#define ERROR 0
#define OK 1
typedef int ElemType;
typedef int Status;
typedef int VRType;
typedef char InfoType;
typedef enum{DG,DN,AG,AN}GraphKind;
typedef char VertexType[MAX_NAME];
//队列结构
typedef struct node
{
ElemType data;
struct node *next;
}QueueNode;
typedef struct
{
QueueNode *front;
QueueNode *rear;
}LinkQueue;
typedef struct
{
VRType adj;
InfoType *info;
}AdjMatrixNode;
typedef struct
{
VertexType vexs[MAX_VERTEX_NUM];
AdjMatrixNode arcs[MAX_VERTEX_NUM][MAX_VERTEX_NUM];
int vexnum,arcnum;
GraphKind kind;
}MGraph;
//求顶点位置
int LocateVex(MGraph G,VertexType u)
{
int i;
for(i=0;i<G.vexnum;++i)
if(strcmp(u,G.vexs[i])==0)
return i;
return ERROR;
}
//求顶点信息
VertexType *GetVex(MGraph G,int v)
{
if(v>=G.vexnum||v<0)
return ERROR;
return &G.vexs[v];
}
//对顶点赋新值
Status PutVex(MGraph *G,VertexType v,VertexType value)
{
int k;
k=LocateVex(*G,v);
if(k<0)
return ERROR;
strcpy((*G).vexs[k],value);
return 1;
}
//求顶点的第一个邻接顶点的序号
int FirstAdjVex(MGraph G, VertexType v)
{
int i,j=0,k;
k=LocateVex(G,v);
if(G.kind==DN||G.kind==AN)
j=INFINITY;
for(i=0;i<G.vexnum;i++)
if(G.arcs[k][i].adj!=j)
return i;
return ERROR;
}
//求顶点v相对于某顶点w的下一个邻接顶点的序号
int NextAdjVex(MGraph G,VertexType v,VertexType w)
{
int i,j=0,k1,k2;
k1=LocateVex(G,v);
k2=LocateVex(G,w);
if(G.kind==DN||G.kind==AN)
j=INFINITY;
for(i=k2+1;i<G.vexnum;i++)
if(G.arcs[k1][i].adj!=j)
return i;
return ERROR;
}
//插入新顶点
void InsertVex(MGraph *G,VertexType v)
{
int i;
strcpy((*G).vexs[(*G).vexnum],v);
for(i=0;i<=(*G).vexnum;i++)
{
if((*G).kind%2)
{
(*G).arcs[(*G).vexnum][i].adj=INFINITY;
(*G).arcs[i][(*G).vexnum].adj=INFINITY;
}
else
{
(*G).arcs[(*G).vexnum][i].adj=0;
(*G).arcs[i][(*G).vexnum].adj=0;
}
(*G).arcs[(*G).vexnum][i].info=NULL;
(*G).arcs[i][(*G).vexnum].info=NULL;
}
(*G).vexnum+=1;
}
//删除顶点
Status DeleteVex(MGraph *G,VertexType v)
{
int i,j,k;
VRType m=0;
k=LocateVex(*G,v);
if(k<0)
return ERROR;
if((*G).kind==DN||(*G).kind==AN)
m=INFINITY;
for(j=0;j<(*G).vexnum;j++)
if((*G).arcs[j][k].adj!=m)
{
if((*G).arcs[j][k].info)
free((*G).arcs[j][k].info);
(*G).arcnum--;
}
if((*G).kind==DG||(*G).kind==DN)
for(j=0;j<(*G).vexnum;j++)
if((*G).arcs[k][j].adj!=m)
{
if((*G).arcs[k][j].info)
free((*G).arcs[k][j].info);
(*G).arcnum--;
}
for(j=k+1;j<(*G).vexnum;j++)
strcpy((*G).vexs[j-1],(*G).vexs[j]);
for(i=0;i<(*G).vexnum;i++)
for(j=k+1;j<(*G).vexnum;j++)
(*G).arcs[i][j-1]=(*G).arcs[i][j];
for(i=0;i<(*G).vexnum;i++)
for(j=k+1;j<(*G).vexnum;j++)
(*G).arcs[j-1][i]=(*G).arcs[j][i];
(*G).vexnum--;
return OK;
}
//插入新弧
Status InsertArc(MGraph *G,VertexType v,VertexType w)
{
int i,len,v1,w1;
char *info,s[MAX_INFO];
v1=LocateVex(*G,v);
w1=LocateVex(*G,w);
if(v1<0||w1<0)
return ERROR;
(*G).arcnum++;
if((*G).kind%2)
{
printf("请输入此弧或边的权值:");
scanf("%d",&(*G).arcs[v1][w1].adj);
}
else
(*G).arcs[v1][w1].adj=1;
printf("是否有该弧或边的相关信息(0:无 1:有):");
scanf("%d%*c",&i);
if(i)
{
printf("请输入该弧或边的相关信息(<%d个字符):",MAX_INFO);
gets(s);
len=strlen(s);
if(len)
{
info=(char*)malloc((len+1)*sizeof(char));
strcpy(info,s);
(*G).arcs[v1][w1].info=info;
}
}
if((*G).kind>1)
{
(*G).arcs[w1][v1].adj=(*G).arcs[v1][w1].adj;
(*G).arcs[w1][v1].info=(*G).arcs[v1][w1].info;
}
return OK;
}
//删除弧
Status DeleteArc(MGraph *G,VertexType v,VertexType w)
{
int v1,w1;
v1=LocateVex(*G,v);
w1=LocateVex(*G,w);
if(v1<0||w1<0)
return ERROR;
if((*G).kind%2==0)
(*G).arcs[v1][w1].adj=0;
else
(*G).arcs[v1][w1].adj=INFINITY;
if((*G).arcs[v1][w1].info)
{
free((*G).arcs[v1][w1].info);
(*G).arcs[v1][w1].info=NULL;
}
if((*G).kind>=2)
{
(*G).arcs[w1][v1].adj=(*G).arcs[v1][w1].adj;
(*G).arcs[w1][v1].info=NULL;
}
(*G).arcnum--;
return OK;
}
//构造有向图G的邻接矩阵
Status CreateDN(MGraph *G)
{
int i,j,k,w,IncInfo;
char s[MAX_INFO],*info;
VertexType va,vb;
printf("请输入图的类型(有向图:0,有向网:1,无向图:2,无向网:3):");
scanf("%d",&(*G).kind);
printf("请输入的顶点数、弧数、弧是否含其他信息(是:1,否:0):");
scanf("%d %d %d",&(*G).vexnum,&(*G).arcnum,&IncInfo);
printf("请输入%d个顶点的值(<%d个字符):",(*G).vexnum,MAX_NAME);
for(i=0;i<(*G).vexnum;++i)
scanf("%s",(*G).vexs[i]);
for(i=0;i<(*G).vexnum;++i)
for(j=0;j<(*G).vexnum;++j)
{
(*G).arcs[i][j].adj=INFINITY;
(*G).arcs[i][j].info=NULL;
}
if((*G).kind==1||(*G).kind==3)
printf("请顺序输入每条弧(边)的权值、弧尾和弧头(以空格作为间隔):\n");
else
printf("请顺序输入每条弧(边)的弧尾和弧头(以空格作为间隔):\n");
for(k=0;k<(*G).arcnum;++k)
{
printf("第%d条弧",k+1);
if((*G).kind==1||(*G).kind==3)
scanf("%d%s%s",&w,va,vb);
else
scanf("%s%s",va,vb);
i=LocateVex(*G,va);
j=LocateVex(*G,vb);
if((*G).kind==0)
(*G).arcs[i][j].adj=1;
if((*G).kind==1)
(*G).arcs[i][j].adj=w;
if((*G).kind==2)
{
(*G).arcs[i][j].adj=1;
(*G).arcs[j][i].adj=1;
}
if((*G).kind==3)
{
(*G).arcs[i][j].adj=w;
(*G).arcs[j][i].adj=w;
}
if(IncInfo)
{
printf("请输入该弧的相关信息(<%d个字符):",MAX_INFO);
gets(s);
w=strlen(s);
if(w)
{
info=(char*)malloc((w+1)*sizeof(char));
strcpy(info,s);
(*G).arcs[i][j].info=info;
}
}
}
printf("该图的邻接矩阵构造完成!");
return OK;
}
//图的邻接矩阵存储结构的打印输出算法
void Display(MGraph G)
{
int i,j;
char s[7],s1[3];
switch(G.kind)
{
case DG:
strcpy(s,"有向图\0");
strcpy(s1,"弧\0");
break;
case DN:
strcpy(s,"有向网\0");
strcpy(s1,"弧\0");
break;
case AG:
strcpy(s,"无向图\0");
strcpy(s1,"边\0");
break;
case AN:
strcpy(s,"无向网\0");
strcpy(s1,"边\0");
}
printf("%d个顶点%d条%s的%s\n",G.vexnum,G.arcnum,s1,s);
for(i=0;i<G.vexnum;i++)
printf("G.vexs[%d]=%s\n",i,G.vexs[i]);
printf("G.arcs.adj:\n");
for(i=0;i<G.vexnum;i++)
{
for(j=0;j<G.vexnum;j++)
printf("G.arc[%d][%d]=%d\n",i,j,G.arcs[i][j].adj);
printf("\n");
}
printf("G.arcs.info:\n");
for(i=0;i<G.vexnum;i++)
{
for(j=0;j<G.vexnum;j++)
printf("%s\t",G.arcs[i][j].info);
printf("\n");
}
}
int InitQueue(LinkQueue *Q) //链队列的初始化操作
{
Q->front=(QueueNode *)malloc(sizeof(QueueNode));
Q->front->next=NULL;
Q->rear=Q->front;
return OK;
}
int EnQueue(LinkQueue *Q,ElemType e) //链队列的入操作
{
QueueNode *p;
p=(QueueNode *)malloc(sizeof(QueueNode));
p->data=e;
p->next=NULL;
Q->rear->next=p;
Q->rear=p;
return OK;
}
int DeQueue(LinkQueue *Q,ElemType *e) //链队列的出队操作
{
QueueNode *p;
if(Q->front==Q->rear)
return ERROR;
p=Q->front->next;
*e=p->data;
Q->front->next=p->next;
if(Q->rear==p)
Q->rear=Q->front;
free(p);
return OK;
}
int QueueEmpty(LinkQueue Q)
{
if(Q.front==Q.rear)
return 0;
else
return 1;
}
/*Visit(MGraph G,int v)
{
printf("%s",G.vexs[v]);
}*/
//邻接矩阵存储结构上的广度优先遍历算法
void BFSTraverse(MGraph G)
{
int v,u,w;
VertexType w1,u1;
LinkQueue Q;
int visited[MAX_VERTEX_NUM];
for(v=0;v<G.vexnum;v++)
visited[v]=-1;
InitQueue(&Q);
for(v=0;v<G.vexnum;v++)
if(visited[v]==-1)
{
visited[v]=1;
printf("%s ",G.vexs[v]);
EnQueue(&Q,v);
while(!QueueEmpty(Q))
{
DeQueue(&Q,&u);
strcpy(u1,*GetVex(G,u));
for(w=FirstAdjVex(G,u1);w>=0;)
{
if(visited[w]==-1)
{
visited[w]=1;
printf("%s ",G.vexs[w]);
EnQueue(&Q,w);
}
w=NextAdjVex(G,u1,strcpy(w1,*GetVex(G,w)));
}
}
}
printf("\n");
}
//邻接矩阵存储结构上的深度优先遍历的递归算法
/*VisitFunc(MGraph G,int v)
{
printf("%s",G.vexs[v]);
}
int visited[MAX_VERTEX_NUM];
// Status(*VisitFunc)(VertexType);
void DFS(MGraph G,int v)
{
VertexType w1,v1;
int w;
visited[v]=1;
VisitFunc(G,v);
strcpy(v1,*GetVex(G,v));
for(w=FirstAdjVex(G,v1);w>=0;w=NextAdjVex(G,v1,strcpy(w1,*GetVex(G,w))))
if(!visited[w])
DFS(G,w);
}
void DFSTraverse(MGraph G)
{
int visited[MAX_VERTEX_NUM];
//VisitFunc=Visit;
int v;
for(v=0;v<G.vexnum;v++)
visited[v]=0;
for(v=0;v<G.vexnum;v++)
if(!visited[v])
DFS(G,v);
printf("\n");
}*/
void DFSTraverse(MGraph G,int v)
{
int visited[MAX_VERTEX_NUM]={0};
int S[MAX_VERTEX_NUM]={0};
int top=-1;
int j;
printf("%s ",G.vexs[v]);
visited[v]=1;
S[++top]=v;
while(top!=-1)
{
v=S[top];
for(j=0;j<G.vexnum;j++)
{
if(G.arcs[v][j].adj==1 && visited[j]==0)
{
printf("%s ",G.vexs[j]);
visited[j]=1;
S[++top]=j;
break;
}
}
if(j==G.vexnum)
top--;
}
}
int main()
{
int v1=1;
MGraph myG;
CreateDN(&myG);
printf("打印图:\n");
Display(myG);
printf("图的广度优先遍历:");
BFSTraverse(myG);
printf("图的深度优先遍历:");
DFSTraverse(myG,0);
}