数据结构——图

一、Dijkstra算法

1. 算法过程

首先定义一个辅助集合$S$，$S$中包括了已经加入到最短路径中的顶点，例如使用一个数组来实现s[]，s[i] = 1代表第$i$个顶点加入到了最短路径。然后定义两个辅助数组：

• dist[]：记录了从源点$v_0$开始，到剩余顶点的最短路径长度，初值为dist[i] = arcs[0][i]
• path[]：path[i]记录了从源点$v_0$到顶点$v_i$的最短路径的前驱结点，再每次算法结束时，可以回溯$v_0$到$v_i$的最短路径。

算法过程：

• 初始化：集合$S$初始化为{0}，初始化dist[]
• 从顶点集合$V-S$中选出$v_j$，满足 d i s t [ j ] = M i n { d i s t [ i ] ∣ v i ∈ V − S } dist[j] = Min\{dist[i] | v_i \in V-S\} dist[j]=Min{dist[i]∣vi​∈V−S}，这里就得到了一条从$v_0$到$v_j$的最短路径。$S=S\cup j$
• 修改从$v_0$出发到$V-S$任意一点$v_k$的最短路径长度：

  if (dist[j] + arcs[j][k] < dist[k])
  	dist[k] = dist[j] + arcs[j][k];
  

• 重复上述步骤，直到所有的顶点都在$S$中。

2. 代码

template <typename T>
void Graph<T>::_dijkstra(int v, int *path, int *shortpath) const {
    bool *visited = new bool [vexnum];
    for (int i = 0; i < vexnum; ++i)
        visited[i] = false;

    visited[v] = true;
    path[v] = -1;
    shortpath[v] = 0;

    int k = v;
    int nextk = -1;
    EdgeNode *p;
    int count = 0;

    for (int i = 1; i < vexnum; ++i) {
        for (p = vexlist[k].first; p; p = p->next) {
            int j = p->adjvex;
            if (!visited[j] && shortpath[j] > shortpath[k] + p->weight)
                shortpath[j] = shortpath[k] + p->weight;
        }
        nextk = _findMinLen(shortpath, visited);    // 查找源点邻点中权值最小的
        if (nextk == -1)
            break;
        else {
            visited[nextk] = true;                  // 将找到的点加入到最短路径中
            path[nextk] = k;
            k = nextk;
            count++;
        }
    }

    std::cout << "shortest path: " << std::endl;
    int lastVer;
    std::stack<int> tmp;
    for (int i = 1; i <= count; ++i) {
        std::cout << shortpath[k] << " : ";

        tmp.push(k);
        lastVer = path[k];
        while (lastVer != v) {
            tmp.push(lastVer);
            lastVer = path[lastVer];
        }
        tmp.push(v);

        std::cout << tmp.top();
        tmp.pop();
        while (!tmp.empty()) {
            std::cout << "->" << tmp.top();
            tmp.pop();
        }
        std::cout << "\n";
        k = path[k];
    }
}

二、Prim算法