问题引入:

  我们先来回想一下生成树是如何定义的,生成树就是用n - 1条边将图中的所有n个顶点都连通为一个连通分量,这样的边连成子树称为生成树。

  最小生成树很明显就是生成树中权值最小的生成树,那么我们即将要学的次小生成树或者K小生成树是怎么定义的呢,很明显就是生成树中权值第k小的生成树。

  下面给出刘老师书中对次小生成树的定义,我是用自己的话描述的。

  对于一个无向图G(V, E),其定义了边权为W(u, v),若T为他的一颗最小生成树,那么我们假设存在一颗生成树T1,不存在任意一颗G的生成树T2满足W(T) <= W(T2) < W(T1),那么我们就称T1为

G的次小生成树。

非严格次小生成树的求解:

  我们知道最小生成树我们是通过Prim和Kruskal这样的贪心算法求得的,那么次小生成树我们只是对这两种算法进行我们需要的修改就可以进行次小生成树的求解。

  我们很容易可以想到最小生成树和次小生成树应该是有联系的,那么是如何联系的呢?次小生成树就是图G的所有生成树中权值第二小的生成树,也就是说我们只需要替换最小生成树的一条边(u, v)

就可以得到次小生成树,显然这条边肯定不可以属于原最小生成树,如果我们将一条不属于原最小生成树的边(u, v)加入T,那么此时T中就形成了一个环,我们在环中选一个除(u, v)权值最大的边进行删除

(想一下为什么是选择权值最大的那条边),得到的树依然是一颗图G的生成树,我们将所有边逐个加入原最小生成树T,得出并记录所有的生成树的权值T1,那么最后T1中最小的那个值即是次小生成树的

权值。

下面只对与求解最小生成树中不同的部分进行说明:

Prim:

  我们知道Prim算法是以给定的任意点作为起始点运用一定的方法对所有点进行贪心处理,缩点从而生成一颗最小生成树,那我们只需要用数组用来描述最小生成树中每条边的访问情况以及最小

生成树中每两个顶点之间的最大边权还需要保存最小生成树中每个顶点的父亲顶点,从而就可以方便用于计算次小生成树。

  具体操作:

    初始化:初始化所有点( i )距离最小生成树子树的距离为cost[source][ i ],所有边初始化为未访问,所有顶点之间的最大边权初始化为0。

    加边:每次加入一条安全边(这里不对安全边进行解释,有不了解的可以查阅博主的上一篇博客),并将最小生成子树中顶点之间的最大边权进行更新,接着更新lowc即可。

    求解次小生成树:我们逐一枚举出所有不属于最小生成树的边(u, v),并且用w(u, v)来替代最大边权和Max(u, v),怎么个替代法?

      SecondMST = min(MST + w(u, v) - Max(u, v))    ((u, v) not belong to MST)。

      OK?

参考代码:

  

 #include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std; const int maxn = + , INF = 0x3f3f3f3f;
int n, m, lowc[maxn], pre[maxn], Max[maxn][maxn], cost[maxn][maxn];
bool vis[maxn], used[maxn][maxn]; int Prim() {
int ans = ;
memset(vis, false, sizeof vis);
memset(Max, , sizeof Max);
memset(used, false, sizeof used);
vis[] = true;
pre[] = -;
for(int i = ; i <= n; i ++) {
lowc[i] = cost[][i];
pre[i] = ;
}
lowc[] = ;
for(int i = ; i <= n; i ++) {
int MIN = INF, p = -;
for(int j = ; j <= n; j ++) {
if(!vis[j] && MIN > lowc[j]) {
MIN = lowc[j];
p = j;
}
}
if(MIN == INF) return -;
ans += MIN;
vis[p] = true;
used[p][pre[p]] = used[pre[p]][p] = true;
for(int j = ; j <= n; j ++) {
if(vis[j] && j != p) Max[j][p] = Max[p][j] = max(Max[j][pre[p]], lowc[p]);
if(!vis[j] && lowc[j] > cost[p][j]) {
lowc[j] = cost[p][j];
pre[j] = p;
}
}
}
return ans;
} int Second_Prim(int MST) {
int ans = INF;
for(int i = ; i <= n; i ++)
for(int j = i + ; j <= n; j ++)
if(!used[i][j] && cost[i][j] != INF) ans = min(ans, MST - Max[i][j] + cost[i][j]);
return ans;
} int main() {
int t, a, b, c;
scanf("%d", &t);
while(t --) {
scanf("%d %d", &n, &m);
for(int i = ; i <= n; i ++)
for(int j = ; j <= n; j ++)
cost[i][j] = INF;
for(int i = ; i <= m; i ++) {
scanf("%d %d %d", &a, &b, &c);
cost[a][b] = cost[b][a] = c;
}
int MST = Prim();
int Second_MST = Second_Prim(MST);
printf("%d\n", Second_MST);
}
return ;
}

Kruskal:

  Kruskal算法是将图G的所有边进行排序,从小到大满足边的两个顶点有一个不在subMST中就将其加入MST,在求解次小生成树问题时我们也需要记录MST中结点的连接情况,以及MST中两个顶点

间的最大边权。

  具体操作:

    初始化:初始化并查集,初始化在subMST中每个结点 i 直接或者间接相连的边为i。

    加边:每次加入一条边时,我们更新subMST中所有与u, v相连的结点间的最大边权,接着将所有与结点v相连的边都与结点u也连起来就行了(前提是在合并时head[ head[ v ] ] = head[ u ])。

    求解次小生成树:

      SecondMST = min(MST + w(u, v) - Max(u, v))    ((u, v) not belong to MST)。滑稽.jpg

参考代码:

 #include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std; const int maxn = + , maxe = * / + , INF = 0x3f3f3f3f;
int n, m, pre[maxn], head[maxn], Max[maxn][maxn];
struct Edge {
int u, v, w;
bool vis;
}edge[maxe];
vector<int> G[maxn]; bool cmp(const Edge &a, const Edge &b) {
return a.w < b.w;
} void Init() {
for(int i = ; i <= n; i ++) {
G[i].clear();
G[i].push_back(i);
head[i] = i;
}
} int Find(int x) {
if(head[x] == x) return x;
return head[x] = Find(head[x]);
} int Kruskal() {
sort(edge + , edge + + m, cmp);
Init();
int ans = , cnt = ;
for(int i = ; i <= m; i ++) {
if(cnt == n - ) break;
int fx = Find(edge[i].u), fy = Find(edge[i].v);
if(fx != fy) {
cnt ++;
edge[i].vis = true;
ans += edge[i].w;
int len_fx = G[fx].size(), len_fy = G[fy].size();
for(int j = ; j < len_fx; j ++)
for(int k = ; k < len_fy; k ++)
Max[G[fx][j]][G[fy][k]] = Max[G[fy][k]][G[fx][j]] = edge[i].w;
head[fx] = fy;
for(int j = ; j < len_fx; j ++)
G[fy].push_back(G[fx][j]);
}
}
return ans;
} int Second_Kruskal(int MST) {
int ans = INF;
for(int i = ; i <= m; i ++)
if(!edge[i].vis)
ans = min(ans, MST + edge[i].w - Max[edge[i].u][edge[i].v]);
return ans;
} int main() {
int t;
scanf("%d", &t);
while(t --) {
scanf("%d %d", &n, &m);
for(int i = ; i <= m; i ++) {
scanf("%d %d %d", &edge[i].u, &edge[i].v, &edge[i].w);
edge[i].vis = false;
}
int MST = Kruskal();
int Second_MST = Second_Kruskal(MST);
printf("%d\n", Second_MST );
}
return ;
}

次小生成树(Prim + Kruaskal)的更多相关文章

  1. UVA 10462 Is There A Second Way Left?(次小生成树&Prim&Kruskal)题解

    思路: Prim: 这道题目中有重边 Prim可以先加一个sec数组来保存重边的次小边,这样不会影响到最小生成树,在算次小生成树时要同时判断次小边(不需判断是否在MST中) Kruskal: Krus ...

  2. The Unique MST POJ - 1679 次小生成树prim

    求次小生成树思路: 先把最小生成树求出来  用一个Max[i][j] 数组把  i点到j 点的道路中 权值最大的那个记录下来 used数组记录该条边有没有被最小生成树使用过   把没有使用过的一条边加 ...

  3. hdu 4081 Qin Shi Huang's National Road System(次小生成树prim)

    题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4081 题意:有n个城市,秦始皇要修用n-1条路把它们连起来,要求从任一点出发,都可以到达其它的任意点. ...

  4. HDOJ-4081(次小生成树+Prim算法)

    Qin Shi Huang's National Road System HDOJ-4081 本题考查的是次小生成树的问题,这里的解决方法就是先使用Prim算法求解最小生成树. 在求解最小生成树的时候 ...

  5. poj1789--最小生成树(prim)

    水题... 题目大意: 用一个7位的字符串代表一个编号,两个编号之间的distance代表这两个编号之间不同字母的个数.一个编号只能由另一个编号“衍生”出来,代价是这两个编号之间相应的distance ...

  6. ZOJ-1586 QS Network---最小生成树Prim

    题目链接: https://vjudge.net/problem/ZOJ-1586 题目大意: 首先给一个t,代表t个测试样例,再给一个n,表示有n个QS装置,接下来一行是n个QS装置的成本.接下来是 ...

  7. POJ-1789 Truck History---最小生成树Prim算法

    题目链接: https://vjudge.net/problem/POJ-1789 题目大意: 用一个7位的string代表一个编号,两个编号之间的distance代表这两个编号之间不同字母的个数.一 ...

  8. 10-最小生成树-Prim算法

    #include <iostream> #include <cstring> #include <cstdio> using namespace std; #def ...

  9. 【次小生成树】【Kruskal】【prim】【转】

    原博客出处:https://blog.csdn.net/yasola/article/details/74276255 通常次小生成树是使用Prim算法进行实现的,因为可以在Prim算法松弛的同时求得 ...

随机推荐

  1. 实现单选框点击label标记中的文字也能选中

    实例: <label for="man"> <input type="radio" value="男" name=&quo ...

  2. atomic,nonatomic

    atomic和nonatomic用来决定编译器生成的getter和setter是否为原子操作.         atomic 设置成员变量的@property属性时,默认为atomic,提供多线程安全 ...

  3. ulimit 管理系统资源

    具体的 options 含义以及简单示例可以参考以下表格. 选项 含义 例子 -H 设置硬资源限制,一旦设置不能增加. ulimit – Hs 64:限制硬资源,线程栈大小为 64K. -S 设置软资 ...

  4. $_POST 和 php://input 的区别

    手册中摘取的几句话: 当 HTTP POST 请求的 Content-Type 是 application/x-www-form-urlencoded 或 multipart/form-data 时, ...

  5. JavaWeb(四):JDBC

    数据持久化(persistence) 把数据保存到可掉电式存储设备中以供之后使用. 大多数情况下,特别是企业级应用,数据持久化意味着将内存中的数据保存到硬盘上加以”固化”,而持久化的实现过程大多通过各 ...

  6. [BZOJ] 地精部落

    问题描述 传说很久以前,大地上居住着一种神秘的生物:地精. 地精喜欢住在连绵不绝的山脉中.具体地说,一座长度为 N 的山脉 H 可分 为从左到右的 N 段,每段有一个独一无二的高度 Hi,其中 Hi ...

  7. 【leetcode】1034. Coloring A Border

    题目如下: Given a 2-dimensional grid of integers, each value in the grid represents the color of the gri ...

  8. chrles设置断点

    1.选择你要断点的接口,右键Breakpoints 2.配置断点接口proxy>Breakpoint settings query设置为* 3.开始断点,重新抓取接口 修改入参.请求头 修改出参 ...

  9. NIO编程模式示例

    1. 服务端 import java.net.InetSocketAddress; import java.nio.ByteBuffer; import java.nio.channels.*; im ...

  10. php stripos()函数 语法

    php stripos()函数 语法 作用:寻找字符串中某字符最先出现的位置,不区分大小写.直线电参数 语法:stripos(string,find,start) 参数: 参数 描述 string  ...