求多源的最短路径

//无向网和有向网
#include<iostream>
#include<vector>
#include<string>
#define Node_MAX 20
#define INF 0x3f3f3f //无穷

using namespace std;
vector<int> PATH[Node_MAX][Node_MAX];//各个顶点的路径
int D[Node_MAX][Node_MAX];//开始是G的矩阵，最后变成最短路径

typedef struct{
    char vexs[Node_MAX];//顶点名称
    int arc[Node_MAX][Node_MAX];//邻接矩阵，代表长度,在网中没有连接用INF表示
    int vexnum,arcnum;//图的顶点个数，弧的个数
    int kind;//图的类型，0代表无向图，1代表有向图
}MGraph;

int FindV(char v,char *p,int num)//查找顶点名称为v的顶点的位置
{
    for(int t=0;t<num;t++)
        if(v==p[t]) return t;
    cout<<"ERROR:输入空的顶点，结束"<<endl;
    exit(-1);
}

bool CreateUDN(MGraph &G)
{
    cout<<"请输入无向图的顶点个数和弧的数量"<<endl;
    cin>>G.vexnum>>G.arcnum;
    for(int i=0;i<G.vexnum;i++)   //初始化arc邻接矩阵
        for(int j=0;j<G.vexnum;j++)
            if(i==j) G.arc[i][j]=0;
            else G.arc[i][j]=INF;
    cout<<"请依次输入顶点的名称(字符类型)"<<endl;
    for(int t=0;t<G.vexnum;t++)
        cin>>G.vexs[t];
    for(int t=0;t<G.arcnum;t++)
    {
        int i,j,x;
        char v1,v2;
        cout<<"请输入第"<<t+1<<"条弧的两端顶点名称及弧长"<<endl;
        cin>>v1>>v2>>x;
        i=FindV(v1,G.vexs,G.vexnum);
        j=FindV(v2,G.vexs,G.vexnum);
        G.arc[i][j]=G.arc[j][i]=x;
    }
    return true;
}
bool CreateDN(MGraph &G)
{
    cout<<"请输入有向图的顶点个数和弧的数量"<<endl;
    cin>>G.vexnum>>G.arcnum;
    for(int i=0;i<G.vexnum;i++)   //初始化arc邻接矩阵
        for(int j=0;j<G.vexnum;j++)
            if(i==j) G.arc[i][j]=0;
            else G.arc[i][j]=INF;
    cout<<"请依次输入顶点的名称(字符类型)"<<endl;
    for(int t=0;t<G.vexnum;t++)
        cin>>G.vexs[t];
    for(int t=0;t<G.arcnum;t++)
    {
        int i,j,x;
        char v1,v2;
        cout<<"请输入第"<<t+1<<"条弧的弧尾弧头名称及弧长"<<endl;
        cin>>v1>>v2>>x;
        i=FindV(v1,G.vexs,G.vexnum);
        j=FindV(v2,G.vexs,G.vexnum);
        G.arc[i][j]=x;
    }
    return true;
}
bool CreateMGraph(MGraph &G)
{
    cout<<"请输入图的类型：有向图DN，无向图UDN"<<endl;
    string temp;
    cin>>temp;
    if(!temp.compare("DN")) G.kind=1;
    if(!temp.compare("UND")) G.kind=0;
    switch(G.kind)
    {
    case 0:CreateUDN(G);
        break;
    case 1:CreateDN(G);
        break;
    default:
        cout<<"ERROR"<<endl;
        exit(-1);
    }
    return true;
}
void FLOYD(MGraph &G)//最短路径算法
{
    for(int i=0;i<G.vexnum;i++)
        for(int j=0;j<G.vexnum;j++)
        {
            PATH[i][j].push_back(i);
            PATH[i][j].push_back(j);
            D[i][j]=G.arc[i][j];
        }
    for(int t=0;t<G.vexnum;t++)//vivj间序号不小于t的最短路径
        for(int i=0;i<G.vexnum;i++)
            for(int j=0;j<G.vexnum;j++)
            {
                if(D[i][j]>D[i][t]+D[t][j])
                {
                    D[i][j]=D[i][t]+D[t][j];
                    PATH[i][j]=PATH[i][t];
                    PATH[i][j].pop_back();
                    PATH[i][j].insert(PATH[i][j].end(),PATH[t][j].begin(),PATH[t][j].end());
                }
            }
}
int main()
{
    MGraph G;
    CreateMGraph(G);
    FLOYD(G);
    for(int i=0;i<G.vexnum;i++)//输出最短路径及其长度
    {
        for(int j=0;j<G.vexnum;j++)
        {
            cout<<G.vexs[i]<<"->"<<G.vexs[j]<<endl;
            for(vector<int>::iterator it=PATH[i][j].begin();it!=PATH[i][j].end();it++)
            {
                cout<<G.vexs[*it]<<" ";
            }
            cout<<endl<<D[i][j]<<endl;
        }

    }
    return 0;
}