// 采用邻接矩阵表示图

// 简单起见，邻接矩阵（图中各结点间的距离权值）在main函数中直接输入了

#include <stdio.h>
#include <stdlib.h>

#define MAXVEX 9
#define INFINITY 65535

typedef int Patharc[MAXVEX];  //存放最短路径上的结点编号
typedef int ShortPathTable[MAXVEX]; //各节点到源点的路径(经过中间结点)长度

struct MGraph{
	int numVertexes;
	int *vex;
	int arc[MAXVEX][MAXVEX];
};

void ShortestPath_Dijkstra(MGraph *G,int v0,int vv,Patharc *P,ShortPathTable *D)
{
	int v,w,k,min;
	int result[MAXVEX];

	int n=0;

	for(v=0;v<G->numVertexes;++v)
	{
		result[v]=0;
		(*D)[v]=G->arc[v0][v];  //初始化，其余各节点到源点的距离
		(*P)[v]=-1;
	}
	(*D)[v0]=0;
	result[v0]=1; //源点加入
	(*P)[n++]=v0;

	for(v=0;v<G->numVertexes;++v)
	{
		min=INFINITY;
		for(w=0;w<G->numVertexes;++w) //寻找未加入的所有结点中到源点距离最短的结点k
		{
			if(!result[w] && (*D)[w]<min)
			{
				k=w;
				min=(*D)[w];
			}
		}
		result[k]=1; //结点k加入
		(*P)[n++]=k;

		//判断是否到达终点
		if(k==vv)
			break;

		for(w=0;w<G->numVertexes;++w) //更新其余各节点到源点的距离的距离数组
		{
			//if(!result[w] && (min+G->arc[k][w]<(*D)[w])) //新节点w到源点的距离=k到源点的距离+k到w的距离
			if(!result[w])
			{
				(*D)[w]=min+G->arc[k][w];
				//(*P)[w]=k;
			}
		}
	}
	//输出结果
	printf("\n");
	for(v=0;v<G->numVertexes;++v)
	{
		printf("%d-->",(*P)[v]);
	}
}

void main()  
{  
    MGraph *my_g=(struct MGraph*)malloc(sizeof(struct MGraph));  
    int i,j;  
    int t=0;

	int v0=0;
	int vv=8;

    my_g->numVertexes=9;  
    my_g->vex=(int*)malloc(sizeof(char)*my_g->numVertexes);  
    if(!my_g->vex) return;  
    for(i=0;i<my_g->numVertexes;++i)  //一维数组(图中各结点)初始化{0,1,2,3,4,5,6,7,8}  
        my_g->vex[i]=i++;

    for(i=0;i<my_g->numVertexes;++i)  
        for(j=0;j<my_g->numVertexes;++j)  
            my_g->arc[i][j]=INFINITY;

	// 无向图的权值二维数组为对称矩阵
    my_g->arc[0][1]=1;  my_g->arc[0][2]=5;
    my_g->arc[1][2]=3;  my_g->arc[1][3]=7;  my_g->arc[1][4]=5;  
    my_g->arc[2][4]=1;  my_g->arc[2][5]=7;
    my_g->arc[3][4]=2;  my_g->arc[3][6]=3;
	my_g->arc[4][5]=3;  my_g->arc[4][6]=6;  my_g->arc[4][7]=9;
	my_g->arc[5][7]=5;
	my_g->arc[6][7]=2;  my_g->arc[6][8]=7;
	my_g->arc[7][8]=4;
    for(i=0;i<my_g->numVertexes;++i)  
        for(j=0;j<=i;++j)  
        {  
            if(i==j)  
            {  
                my_g->arc[i][j]=0;  
                continue;  
            }  
            my_g->arc[i][j]=my_g->arc[j][i];  
        }  
    for(i=0;i<my_g->numVertexes;++i)  //二维数组表示图中各结点间连接边的weight  
    {  
        for(j=0;j<my_g->numVertexes;++j)  
            printf("%5d  ",my_g->arc[i][j]);  
        printf("\n");  
    }  
    printf("\n\n");


	Patharc P;
	ShortPathTable D;
	ShortestPath_Dijkstra(my_g,v0,vv,&P,&D);

	free(my_g->vex);
	//free(my_g->arc);
}

运行结果表明从节点0到结点8的最短路径所经过的结点编号（在本例中输出了初始化时的-1，可优化不输出-1）