2020.10.24--严格次小生成树 --详细讲解

---

2020-1024=996

品，细品

---

10.24祝各位coder编码有成，学业进步，工作顺利，升职加薪！！！

前置知识
1.Prim算法、Kurskal算法求最小生成树
2.LCA求树上两点间的最近公共祖先

概念

如果最小生成树选择的边集是EM，严格次小生成树选择的边集是ES，那么需要满足：
​

​

思路
要求次小生成树 那就是在最小生成树上进行加边和删边的操作求得次小生成树

如果我们求出了次小生成树，我们任意取树上两点u,v 之间连接一条边，那么此时一定生成了一个环 ，我们需要在这个环中删去一条边维护树的结构，那么我删除的这条边要和 添加的这条边的权值尽可能相等（但不能相等），因为这样才可能使新树的权值增加的少。

步骤
1.求出最小生成树
2.求出最小生成树上任意两点间的最大值和次大值
3.枚举多余的边，添加到最小生成树中，然后在产生的环中删除一条边
4.在多次产生的新树的权值中取最小的

具体操作
求任意两点间的最大值和次大值的时候，可以用lca 的思想，树上倍增的方法。
如果添加的一条边连接的是u,v两个点，那么我们可以分别求u到lca的最大值或者次大值，v也同理，然后两者取最大。
因为添加的边的a[i].w的权值一定大于等于树上的最大的边的权值，所以如果求的树上亮点最大值和a[i].w相等，我们取出次小值，只有这样，才能让新树的权值增加的尽可能小

注意

树上任意两点间的次大值要初始化为负无穷！！！

如果不的话，因为如果最小生成树中的边权全部相等，然后添加了一条边的权值和树上的权值也相等的话，求出来的答案是min-0+w;
然后有另一条边的权值是（w+10) 求出来的答案是min-(w)+(w+10);

因为前者并不是答案，但是会被判定成答案，因为他小（在w满足一定情况下）

（自行画一棵树理解一下，比较简单）

Code

#include<iostream>
#include<algorithm>
#include<math.h>
#include<queue>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <stack>
#include <map>
#include <set>
#include <vector>
#include <iomanip>

using namespace std;
typedef long long ll;
typedef pair<int, int> PII;
typedef unsigned long long ull;
const int inf = 0x3f3f3f3f;
const int maxn = 8e5 + 7;
const ll mod = 1000000007;

#define mst(x, a) memset( x,a,sizeof(x) )
#define rep(i, a, b) for(int i=(a);i<=(b);++i)
#define dep(i, a, b) for(int i=(a);i>=(b);--i)

ll read() {
	ll x = 0;
	char ch = getchar();
	while(ch < '0' || ch > '9')ch = getchar();
	while(ch >= '0' && ch <= '9')x = (x << 1) + (x << 3) + (ch ^ 48), ch = getchar();
	return x;
}

void out(ll x) {
	int stackk[40];
	if (x < 0) {
		putchar('-');
		x = -x;
	}
	if (!x) {
		putchar('0');
		return;
	}
	int top = 0;
	while (x) stackk[++top] = x % 10, x /= 10;
	while (top) putchar(stackk[top--] + '0');
}
ll n,m,p[maxn],cnt,head[maxn],fa[maxn][22],max1[maxn][22];
ll max2[maxn][22],dep[maxn],len,vis[maxn],paid;
struct node {
	ll u,v,w;
} a[maxn];
bool cmp(node x,node y) {
	return x.w<y.w;
}
struct nodee {
	ll u,v,w,next;
} e[maxn];
void add(ll u,ll v,ll w) {
	e[cnt].u=u,e[cnt].v=v;
	e[cnt].w=w,e[cnt].next=head[u];
	head[u]=cnt++;
}
void dfs(ll u,ll p) {
	dep[u]=dep[p]+1,fa[u][0]=p;
	for(int i=1 ; i<=len ; i++) fa[u][i]=fa[fa[u][i-1]][i-1];
	for(int i=head[u]; ~i; i=e[i].next) {
		int v=e[i].v;
		if(v==p) continue;
		max1[v][0]=e[i].w;
		dfs(v,u);
	}
}
void cal() {
	for(int i=1 ; i<=len; i++) {
		for(int j=1 ; j<=n; j++)   {
			max1[j][i]=max(max1[j][i-1],max1[fa[j][i-1]][i-1]);
			max2[j][i]=max(max2[j][i-1],max2[fa[j][i-1]][i-1]);
			if(max1[j][i-1]<max1[fa[j][i-1]][i-1]) max2[j][i]=max1[j][i-1];
			else    if(max1[j][i-1]>max1[fa[j][i-1]][i-1]) max2[j][i]=max1[fa[j][i-1]][i-1];
		}
	}
}
ll lca(ll x,ll y) {
	if(dep[x]<dep[y]) swap(x,y);
	for(int i=len ; i>=0 ; i--) if(dep[fa[x][i]]>=dep[y]) x=fa[x][i];
	if(x==y) return x;
	for(int i=len ; i>=0 ; i--)
		if(fa[x][i]!=fa[y][i]) x=fa[x][i],y=fa[y][i];
	return fa[x][0];
}
ll q(ll x,ll y,ll w) {
	ll ans=0;
	for(int i=len ; i>=0 ; i--)
		if(dep[fa[x][i]]>=dep[y]) {
			if(max1[x][i]!=w) ans=max(ans,max1[x][i]);
			else ans=max(ans,max2[x][i]);
			x=fa[x][i];
		}
	return ans;
}
ll find(ll x) {
	if(p[x]==x) return x;
	return p[x]=find(p[x]);
}
int main() {
	mst(head,-1);
	n=read(),m=read();
	len=log2(n)+2;
	for(int i=0 ; i<=n ; i++) p[i]=i;
	for(int i=1 ; i<=m ; i++) {
		a[i].u=read();
		a[i].v=read();
		a[i].w=read();
	}
	sort(a+1,a+1+m,cmp);
	for(int i=1 ; i<=m ; i++) {
		ll u=a[i].u;
		ll v=a[i].v;
		if(find(u)==find(v)) continue;
		p[find(u)]=find(v);
		add(u,v,a[i].w);
		add(v,u,a[i].w);
		vis[i]=1,paid+=a[i].w;
	}
	max2[1][0]=-1e17;
	dfs(1,-1);
	cal();
	ll ans=1e17;
	for(int i=1 ; i<=m ; i++) {
		if(vis[i]) continue;
		ll u=a[i].u;
		ll v=a[i].v;
		ll lc=lca(u,v);
		ll temp1=q(u,lc,a[i].w);
		ll temp2=q(v,lc,a[i].w);
		ans=min(ans,paid-max(temp1,temp2)+a[i].w);
	}
	printf("%lld",ans);
	return 0;
}
/*


*/