算法:最小生成树(prime,kruskal,floyd)

本文详细介绍了三种求解最小生成树的算法:Floyd算法、Kruskal算法和Prime算法。针对有向图和无向图,分别展示了如何处理负权边和避免负权回路。提供了每种算法的思路、伪代码及C++实现,适用于处理含有重边和自环的图。同时,针对每种算法的数据复杂度进行了分析,包括O(n^3)的Floyd算法和O(mlogm)的Kruskal算法等。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

最小生成树

Floyd算法O(n^3)- 动态规划

思路

  • f[i, j, k] 表示从i走到j的路径上除ij点外只经过1k的点的所有路径的最短距离。那么f[i, j, k] = min(f[i, j, k - 1], f[i, k, k - 1] + f[k, j, k - 1])
  • 读入邻接矩阵,每次通过动态规划转换成从i到j的最短距离矩阵
  • 在下面代码中,判断从a到b是否是无穷大距离时,需要进行if(t > INF/2)判断,而非是if(t == INF)判断,原因是INF是一个确定的值,并非真正的无穷大,会随着其他数值而受到影响,t大于某个与INF相同数量级的数即可。
  • 注意邻接矩阵的初始化。

题目描述

给定一个n个点m条边的有向图,图中可能存在重边和自环,边权可能为负数。

再给定k个询问,每个询问包含两个整数x和y,表示查询从点x到点y的最短距离,如果路径不存在,则输出“impossible”。

数据保证图中不存在负权回路。

输入格式

第一行包含三个整数n,m,k

接下来m行,每行包含三个整数x,y,z,表示存在一条从点x到点y的有向边,边长为z。

接下来k行,每行包含两个整数x,y,表示询问点x到点y的最短距离。

输出格式

共k行,每行输出一个整数,表示询问的结果,若询问两点间不存在路径,则输出“impossible”。

数据范围

1≤n≤200,
1≤k≤n2
1≤m≤20000,
图中涉及边长绝对值均不超过10000。

输入样例:
3 3 2
1 2 1
2 3 2
1 3 1
2 1
1 3
输出样例:
impossible
1

代码
#include <iostream>
#include <cstring>
using namespace std;
const int N = 210, INF = 1e9;
int d[N][N];

int main()
{
    int n, m, q;
    cin >> n >> m >> q;
    for (int i = 1; i <= n; ++i)
    {
        for (int j = 1; j <= n; ++j)
        {
            if (i == j) d[i][j] = 0;
            else d[i][j] = INF;
        }
    }
    while (m --)
    {
        int x, y, k;
        cin >> x >> y >> k;
        d[x][y] = min(d[x][y], k);
    }
    for (int k = 1; k <= n; ++k)
        for (int i = 1; i <= n; ++i)
            for (int j = 1; j <= n; ++j)
                d[i][j] = min(d[i][j], d[i][k] + d[k][j]);
    while(q --)
    {
        int x, y;
        cin >> x >> y;
        if (d[x][y] > INF / 2) cout << "impossible" << endl;
        else cout << d[x][y] << endl;
    }
    return 0;
}

原题链接

Kruskal算法O(mlogm)

思路

  • 将所有边按权重从小到大排序
  • 枚举每条边a,b权重c
    • if a, b两点不连通
      将a, b边加入集合中

注意:枚举每条边的操作是并查集操作。
需要使用变量cnt来记录加进集合的边数,若cnt < n - 1表示不能遍历所有点

问题描述

给定一个n个点m条边的无向图,图中可能存在重边和自环,边权可能为负数。

求最小生成树的树边权重之和,如果最小生成树不存在则输出impossible。

给定一张边带权的无向图G=(V, E),其中V表示图中点的集合,E表示图中边的集合,n=|V|,m=|E|。

由V中的全部n个顶点和E中n-1条边构成的无向连通子图被称为G的一棵生成树,其中边的权值之和最小的生成树被称为无向图G的最小生成树。

输入格式

第一行包含两个整数n和m。

接下来m行,每行包含三个整数u,v,w,表示点u和点v之间存在一条权值为w的边。

输出格式

共一行,若存在最小生成树,则输出一个整数,表示最小生成树的树边权重之和,如果最小生成树不存在则输出impossible。

数据范围
1≤n≤105,
1≤m≤2∗105,
图中涉及边的边权的绝对值均不超过1000。

输入样例:
4 5
1 2 1
1 3 2
1 4 3
2 3 2
3 4 4
输出样例:
6

代码

#include <iostream>
#include <algorithm>
using namespace std;
const int N = 1000010;

int n, m;
int p[N];
struct Edge
{
    int a, b, w;
    
    bool operator< (const Edge &W)const
    {
        return w < W.w;
    }
}edges[N];

int find(int x)
{
    if (x != p[x]) p[x] = find(p[x]);
    return p[x];
}


int main()
{
    int n, m;
    cin >> n >> m;
    for (int i = 1; i <= n; ++i)
        p[i] = i;
    for (int i = 0; i < m; ++i)
    {
        int a, b, w;
        cin >> a >> b >> w;
        edges[i] = {a, b, w};
    }
    sort(edges, edges+m);
    int ans = 0, cnt = 0;
    for (int i = 0; i < m; ++i)
    {
        int a = find(edges[i].a);
        int b = find(edges[i].b);
        if (a != b){
            p[a] = b;
            ans += edges[i].w;
            cnt ++;
        }
    }
    if (cnt < n - 1) cout << "impossible";
    else cout << ans << endl;
    return 0;
}

原题链接

Prime算法O(n^2)

思路

dist[i]距离设置为无穷大
s:当前已经在连通块中的所有点
for i 1…n
(1)t 找到集合外距离最近的点 总共O(n^2)次
(2)st 对该店进行标记
(3)用t更新其他点到集合的距离
dist[1] = 0

问题描述

给定一个n个点m条边的无向图,图中可能存在重边和自环,边权可能为负数。

求最小生成树的树边权重之和,如果最小生成树不存在则输出impossible。

给定一张边带权的无向图G=(V, E),其中V表示图中点的集合,E表示图中边的集合,n=|V|,m=|E|。

由V中的全部n个顶点和E中n-1条边构成的无向连通子图被称为G的一棵生成树,其中边的权值之和最小的生成树被称为无向图G的最小生成树。

输入格式

第一行包含两个整数n和m。

接下来m行,每行包含三个整数u,v,w,表示点u和点v之间存在一条权值为w的边。

输出格式

共一行,若存在最小生成树,则输出一个整数,表示最小生成树的树边权重之和,如果最小生成树不存在则输出impossible。

数据范围
1≤n≤500,
1≤m≤105,
图中涉及边的边权的绝对值均不超过10000。

输入样例:
4 5
1 2 1
1 3 2
1 4 3
2 3 2
3 4 4
输出样例:
6

代码

#include <iostream>
#include <cstring>
using namespace std;

int n, m;
const int N = 510, INF = 0x3f3f3f3f;
int g[N][N];
int st[N];
int dist[N];

int prime()
{
    int ans = 0;
    memset (dist, 0x3f, sizeof dist);
    for (int i = 0; i < n; ++i)
    {
        int t = -1;
        for (int j = 1; j <= n; ++j)
            if (!st[j] && (t == -1 || dist[t] > dist[j]))
                t = j;
        if (i && dist[t] == INF) return INF;
        if (i) ans += dist[t];   
        st[t] = true;
        for (int j = 1; j <= n; ++j)
        {
            dist[j] = min(dist[j], g[t][j]);
        }
    }
    return ans;
}


int main()
{
    cin >> n >> m;
    memset(g, 0x3f, sizeof g);
    while (m -- )
    {
        int x, y, c;
        cin >> x >> y >> c;
        g[x][y] = g[y][x] = min(g[x][y], c);
    }
    int t = prime();
    if (t == INF) cout << "impossible" << endl;
    else cout << t << endl;
    return 0;
}

原题链接

好的,下面是使用C语言实现的Kruskal算法、Prim算法、Dijkstra算法Floyd算法的完整代码。这些代码将输出最小生成树和最短路径长度矩阵及其路径。 ### Kruskal算法 ```c #include <stdio.h> #include <stdlib.h> #define MAX 100 typedef struct { int u, v, w; } Edge; int parent[MAX]; int find(int i) { while (parent[i] != i) i = parent[i]; return i; } void union1(int i, int j) { int a = find(i); int b = find(j); parent[a] = b; } void kruskal(Edge edges[], int n, int e) { int i, j, u, v, min; Edge temp; int sum = 0; for (i = 0; i < n; i++) parent[i] = i; for (i = 0; i < e - 1; i++) { min = i; for (j = i + 1; j < e; j++) { if (edges[j].w < edges[min].w) min = j; } temp = edges[i]; edges[i] = edges[min]; edges[min] = temp; u = find(edges[i].u); v = find(edges[i].v); if (u != v) { printf("Edge %d: %d - %d = %d\n", i + 1, edges[i].u, edges[i].v, edges[i].w); union1(u, v); sum += edges[i].w; } } printf("Total weight of the minimum spanning tree: %d\n", sum); } ``` ### Prim算法 ```c #include <stdio.h> #include <stdbool.h> #include <limits.h> #define V 5 int minKey(int key[], bool mstSet[]) { int min = INT_MAX, min_index; for (int v = 0; v < V; v++) if (mstSet[v] == false && key[v] < min) min = key[v], min_index = v; return min_index; } void printMST(int parent[], int graph[V][V]) { printf("Edge \tWeight\n"); for (int i = 1; i < V; i++) printf("%d - %d \t%d \n", parent[i], i, graph[i][parent[i]]); } void primMST(int graph[V][V]) { int parent[V]; int key[V]; bool mstSet[V]; for (int i = 0; i < V; i++) key[i] = INT_MAX, mstSet[i] = false; key[0] = 0; parent[0] = -1; for (int count = 0; count < V - 1; count++) { int u = minKey(key, mstSet); mstSet[u] = true; for (int v = 0; v < V; v++) if (graph[u][v] && mstSet[v] == false && graph[u][v] < key[v]) parent[v] = u, key[v] = graph[u][v]; } printMST(parent, graph); } ``` ### Dijkstra算法 ```c #include <stdio.h> #include <stdbool.h> #include <limits.h> #define V 5 int minDistance(int dist[], bool sptSet[]) { int min = INT_MAX, min_index; for (int v = 0; v < V; v++) if (sptSet[v] == false && dist[v] <= min) min = dist[v], min_index = v; return min_index; } void printSolution(int dist[], int n) { printf("Vertex \tDistance from Source\n"); for (int i = 0; i < V; i++) printf("%d \t\t %d\n", i, dist[i]); } void dijkstra(int graph[V][V], int src) { int dist[V]; bool sptSet[V]; for (int i = 0; i < V; i++) dist[i] = INT_MAX, sptSet[i] = false; dist[src] = 0; for (int count = 0; count < V - 1; count++) { int u = minDistance(dist, sptSet); sptSet[u] = true; for (int v = 0; v < V; v++) if (!sptSet[v] && graph[u][v] && dist[u] != INT_MAX && dist[u] + graph[u][v] < dist[v]) dist[v] = dist[u] + graph[u][v]; } printSolution(dist, V); } ``` ### Floyd算法 ```c #include <stdio.h> #define V 5 void printSolution(int dist[][V]) { printf("The matrix of shortest path lengths is:\n"); for (int i = 0; i < V; i++) { for (int j = 0; j < V; j++) { if (dist[i][j] == INT_MAX) printf("INF\t"); else printf("%d\t", dist[i][j]); } printf("\n"); } } void floydWarshall(int graph[][V]) { int dist[V][V]; for (int i = 0; i < V; i++) for (int j = 0; j < V; j++) dist[i][j] = graph[i][j]; for (int k = 0; k < V; k++) { for (int i = 0; i < V; i++) { for (int j = 0; j < V; j++) { if (dist[i][k] != INT_MAX && dist[k][j] != INT_MAX && dist[i][k] + dist[k][j] < dist[i][j]) dist[i][j] = dist[i][k] + dist[k][j]; } } } printSolution(dist); } ``` 这些代码实现了Kruskal算法、Prim算法、Dijkstra算法Floyd算法,并输出相应的结果。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

evil心安

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值