前置知识:
关于图的几个概念定义:
- 连通图:在无向图中,若任意两个顶点vi与vj都有路径相通,则称该无向图为连通图。
- 强连通图:在有向图中,若任意两个顶点vi与vj都有路径相通,则称该有向图为强连通图。
- 连通网:在连通图中,若图的边具有一定的意义,每一条边都对应着一个数,称为权;权代表着连接连个顶点的代价,称这种连通图叫做连通网。
- 生成树:一个连通图的生成树是指一个连通子图,它含有图中全部n个顶点,但只有足以构成一棵树的n-1条边。一颗有n个顶点的生成树有且仅有n-1条边,如果生成树中再添加一条边,则必定成环。
- 最小生成树:在连通网的所有生成树中,所有边的代价和最小的生成树,称为最小生成树。
一、定义
(看板子题就好)
给出一个无向图,求出最小生成树,如果该图不连通,则输出orz
输入格式:
第一行包含两个整数N、M,表示该图共有N个结点和M条无向边。(N<=5000,M<=200000)
接下来M行每行包含三个整数Xi、Yi、Zi,表示有一条长度为Zi的无向边连接结点Xi、Yi
输出格式:
输出包含一个数,即最小生成树的各边的长度之和;如果该图不连通则输出orz
二、有两种算法
kruskal 和 prim
两者区别:Prim在稠密图中比Kruskal优,在稀疏图中比Kruskal劣。
Prim是以更新过的节点的连边找最小值,Kruskal是直接将边排序。
kruskal
此算法又可称为“加边法”
初始最小生成树边数为0
每迭代一次就选择一起满足条件的最小代价边,加入到最小生成树的边集合中
1.把图中的所有边按代价从小到大排序
2.把图中的n个顶点看成独立的n棵树
3.按权值从小到大选择边,所选的边连接的两个顶点ui,viui,vi,应属于两颗不同的树,则成为最小生成树的一条边,并将这两颗树合并作为一颗树。
4. 重复(3),直到所有顶点都在一颗树内或者有n-1条边为止。
#include#include #include using namespace std;inline int read(){ int sum = 0,p = 1; char ch = getchar(); while(ch < '0' || ch > '9') { if(ch == '-') p = -1; ch = getchar(); } while(ch >= '0' && ch <= '9') { (sum *= 10)+= ch - '0'; ch = getchar(); } return sum * p;}const int maxm = 200005,maxn = 5005;struct Edge{ int from,to,wei;}edge[maxm];int fa[maxn],n,m,ans,u,v,cnt;inline bool cmp(Edge a,Edge b){ return a.wei < b.wei;}int find(int o){ if(o == fa[o]) return o; else return fa[o] = find(fa[o]);}inline void kruskal(){ sort(edge+1,edge+m+1,cmp); for(int i = 1;i <= m;i++) { u = find(edge[i].from); v = find(edge[i].to); if(u == v) continue; ans += edge[i].wei; fa[v] = u; if(++cnt == n - 1) break; }}int main(){ n = read(),m = read(); for(int i = 1;i <= n;i++) fa[i] = i; for(int i = 1;i <= m;i++) { edge[i].from = read(); edge[i].to = read(); edge[i].wei = read(); } kruskal(); printf("%d",ans); return 0;}
Prim
此算法也可以成为“加点法”
每次迭代选择代价最小的边对应的点加入到最小生成树中
算法从某一个顶点s开始逐渐长大覆盖整个连通网的所有顶点
1.图的所有顶点集合为v;初始令集合u={s},v=V-u;
2.在两个集合u,v能组成的边中,选择我一条代价最小的边(u0,v0),加入到最小生成树中,并把v0并入到集合u中
3.重复上述步骤,知道最小生成树有n-1条边或者n个顶点为止
由于不断向集合u中加点,所以最小代价边必须同步更新;需要建立一个辅助数组closedge,用来维护集合v中每个顶点与集合u中最小代价边信息
#include#include #include using namespace std;inline int read(){ int sum = 0,p = 1; char ch = getchar(); while(ch < '0' || ch > '9') { if(ch == '-') p = -1; ch = getchar(); } while(ch >= '0' && ch <= '9') { (sum *= 10)+= ch - '0'; ch = getchar(); } return sum * p;}const int maxn = 5005,maxm = 200005,inf = 123456789;struct edge{ int v,w,next;}e[maxm<<1];int head[maxn],dis[maxn],cnt,n,m,tot,now = 1,ans;bool vis[maxn];inline void add(int u,int v,int w){ e[++cnt].v = v; e[cnt].w = w; e[cnt].next = head[u]; head[u] = cnt;}int prim(){ for(int i = 2;i <= n;i++) dis[i] = inf; for(int i = head[1];i;i = e[i].next) dis[e[i].v] = min(dis[e[i].v],e[i].w); while(++tot <= n;i++) if(!vis[i] && minn > dis[i]) { minn = dis[i]; now = i; } ans += minn; for(int i = head[now];i;i = e[i].next) { int v = e[i].v; if(dis[v] > e[i].w && !vis[v]) dis[v] = e[i].w; } } return ans;}int main(){ n = read(),m = read(); int u,v,w; for(int i = 1;i <= m;i++) { u = read(),v = read(),w = read(); add(u,v,w); add(v,u,w); } printf("%d",prim()); return 0;}