基于堆优化的Prim,另附kruskal解法

最小生成树的Prim算法也是贪心算法的一大经典应用。Prim算法的特点是时刻维护一棵树,算法不断加边,加的过程始终是一棵树。

Prim算法过程:

一条边一条边地加, 维护一棵树。

初始 E = {}空集合, V = {任意节点}

循环(n – 1)次,每次选择一条边(v1,v2), 满足:v1属于V , v2不属于V。且(v1,v2)权值最小。

E = E + (v1,v2)
V = V + v2

最终E中的边是一棵最小生成树, V包含了全部节点。
最后,我们来提供输入输出数据,由你来写一段程序,实现这个算法,只有写出了正确的程序,才能继续后面的课程。

输入

第1行:2个数N,M中间用空格分隔,N为点的数量,M为边的数量。(2 <= N <= 1000, 1 <= M <= 50000)
第2 - M + 1行:每行3个数S E W,分别表示M条边的2个顶点及权值。(1 <= S, E <= N,1 <= W <= 10000)

输出

输出最小生成树的所有边的权值之和。

输入示例

9 14
1 2 4
2 3 8
3 4 7
4 5 9
5 6 10
6 7 2
7 8 1
8 9 7
2 8 11
3 9 2
7 9 6
3 6 4
4 6 14
1 8 8

输出示例

37

对于prim, 如果每次都遍历包含在已构造生成树的点集的minCost[v],需要O(V2)时间。
但是如果用堆来维护minCost[v],则可以把时间降为O(ElogV)

Prim实现:

#include <cstdio>
#include <set>
#include <queue>
#include <vector>
#include <utility>
#include <cstring>
using namespace std; 

typedef long long ll;
typedef pair<int, int> pii;
const int EMAX = 5000 + 10;
const int VMAX = 1000 + 10;

struct Node {
  int a, b;
  ll w;
  Node(int x, int y, ll z) {
    a = x, b = y, w = z;
  }
  bool operator < (const Node &n) const {
    return w > n.w;
  }
};

int vNum, eNum;
bool look[VMAX];
vector<pii> e[EMAX];
set<int> v;

void Init() {
  memset(look, 0, sizeof(look));
  for (int i = 0; i < EMAX; ++i) {
    e[i].clear();
  }
  v.clear();
}

ll Prim(int s) {
  priority_queue<Node> que;
  que.push(Node(s, s, 0));
  ll ans = 0;
  while (!que.empty()) {
    Node head = que.top(); que.pop();
    if (look[head.b]) continue;
    look[head.b] = true;
    ans += head.w;
    for (int i = 0; i < e[head.b].size(); ++i) {
      pii &edge = e[head.b][i];
      if (!look[edge.first]) {
        que.push(Node(head.b, edge.first, edge.second));
      }
    }
  }
  return ans;
}

int main() {
#ifndef ONLINE_JUDGE
  freopen("in.txt", "r", stdin);
#endif
  while (scanf("%d%d", &vNum, &eNum) != EOF) {
    Init();
    int a, b, w;
    for (int i = 0; i < eNum; ++i) {
      scanf("%d%d%d", &a, &b, &w);
      e[a].push_back(pii(b, w));
      e[b].push_back(pii(a, w));
      v.insert(a);
      v.insert(b);
    }
    printf("%lld\n", Prim(a));
  }
  return 0;
}

Kruskal实现:

#include <cstdio>
#include <set>
#include <queue>
#include <vector>
#include <utility>
#include <cstring>
#include <algorithm>
using namespace std; 

typedef long long ll;
const int EMAX = 5000 + 10;
const int VMAX = 1000 + 10;

struct Edge {
  int a, b, w;
  Edge(int x, int y, int z) {
    a = x, b = y, w = z;
  }
  bool operator < (const Edge &arg) const {
    return w < arg.w;
  }
};

int vNum, eNum, father[VMAX];
vector<Edge> e;

void Init() {
  e.clear();
  for (int i = 0; i < VMAX; ++i) {
    father[i] = i;
  }
}

int Find(int x) {
  return father[x] == x ? x : father[x] = Find(father[x]);
}

void Unite(int a, int b) {
  a = Find(a), b = Find(b);
  if (a != b) {
    father[b] = a;
  }
}

ll Kruskal() {
  sort(e.begin(), e.end());
  ll ans = 0;
  for (int i = 0; i < eNum; ++i) {
    if (Find(e[i].a) != Find(e[i].b)) {
      ans += e[i].w;
      Unite(e[i].a, e[i].b);
    }
  }
  return ans;
}

int main() {
#ifndef ONLINE_JUDGE
  //freopen("in.txt", "r", stdin);
#endif
  while (scanf("%d%d", &vNum, &eNum) != EOF) {
    Init();
    int a, b, w;
    for (int i = 0; i < eNum; ++i) {
      scanf("%d%d%d", &a, &b, &w);
      e.push_back(Edge(a, b, w));
    }
    printf("%lld\n", Kruskal());
  }
  return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值