[UOJ176]新年的繁荣

最新推荐文章于 2023-12-20 20:27:56 发布

OneInDark

最新推荐文章于 2023-12-20 20:27:56 发布

阅读量335

点赞数 1

分类专栏：图论 # 树数学文章标签：算法 c语言数据结构

本文链接：https://blog.csdn.net/qq_42101694/article/details/104065299

版权

数学同时被 3 个专栏收录

280 篇文章 1 订阅

订阅专栏

图论

122 篇文章 0 订阅

订阅专栏

树

98 篇文章 0 订阅

订阅专栏

题目

传送门 to UOJ

题目概要
$n$ 个点的完全图，给定整数序列 $a_1,a_2,a_3,\dots,a_n$ 满足边 $\lang x,y\rang$ 的边权为 $(a_x\operatorname{and}a_y)$ ，其中 $\text{and}$ 是按位与运算符。求最大生成树。

数据范围与约定
将会给定一个 $m$ ，满足 $\forall x\in[1,n],\;a_x\in[0,2^m)$ 。

对于所有数据， $n\le 10^5$ 且 $m\le 18$ 。

思路壹

感谢@CzxingcHen 的博客提供了思路与代码。

如果 $a$ 值有重复，譬如 $a_x=a_y$ ，那么所有与 $a_y$ 相连的边，都可以嫁接到 $x$ 身上（而边权不变）。毕竟 $a_x$ 和 $a_y$ 是一定要连边的（根据相关定理，与 $a_x$ 相连的所有边中，边权最大的一个一定属于最大生成树）。所以，可以认为 $a$ 的值互不相同。

为了方便叙述，认为 $a_x$ 就是点的编号。也就是说点的编号并不是连续的。

考虑使用 $\tt kruskal$ 算法求解。边权的大小 $v\in[0,2^m)$ ，枚举一个边权，问题转化为了：找到有哪些边的边权为 $v$ 。

首先边权为 $v$ 的边一定是点集 $\{x\;|\;v\subseteqq x\}$ 之间的边。其次，这个点集内的边的权不小于 $v$ ，所以连接之后它们会成为一个连通块。类似的，点集 $\{x\;|\;(v\operatorname{or}2^r)\subseteqq x\}$ 在考虑边权为 $v$ 的边之前，就已经成为了一个连通块。一个连通块只需要选取一个点。

同时注意到 $r\in[0,m)$ 就已经把点集都划分完了。所以我们只需要找这 $m$ 个连通块，每个选出一个代表，然后尝试和 $v$ 连边。反正两两连接，和直接与 $v$ 连接，效果是一样的，本质上是计算了连通块个数。

我们可以用 $f (v)$ 表示，点集 $\{x\;|\;v\subseteqq x\}$ 这个连通块中，任取的一个代表。事实上，只要随意找到一个点即可，因为形成连通块是必然的。那么可以递推，在 $f(v\operatorname{or}2^r)$ 中找一个就行。

时间复杂度 $\mathcal O(m2^m+n)$ ，相当快。

代码

#include <cstdio>
#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long int_;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(int x){
	if(x > 9) writeint(x/10);
	putchar((x-x/10*10)^48);
}

const int MaxN = 1<<18;
namespace UFS{
	int fa[MaxN], rnk[MaxN];
	void init(int n){
		rep(i,1,n)
			fa[i] = i, rnk[i] = 1;
	}
	inline int find(int a){
		if(fa[a] != a)
			fa[a] = find(fa[a]);
		return fa[a];
	}
	bool merge(int a,int b){
		a = find(a), b = find(b);
		if(a == b) return false;
		if(rnk[a] > rnk[b]) swap(a,b);
		fa[a] = b, rnk[b] += rnk[a];
		return true;
	}
}

int f[MaxN];
int main(){
	int n = readint(), m = readint();
	int_ ans = 0;
	for(int a; n; --n){
		a = readint();
		if(f[a]) ans += a;
		else f[a] = a;
	}
	UFS::init((1<<m)-1);
	for(int i=(1<<m)-1; ~i; --i){
		for(int j=0; j<m&&(!f[i]); ++j)
			f[i] = f[i|(1<<j)];
		if(!f[i]) continue; // nothing
		for(int j=0; j<m; ++j)
			if(!(i>>j&1) && f[i^(1<<j)] &&
			UFS::merge(f[i],f[i^(1<<j)]))
				ans += i; // one more edge
	}
	printf("%lld\n",ans);
	return 0;
}

思路贰（正确性未知）

这是我的第一思路。似乎和思路一的思考方向差不多。但是一点也不显然。

鉴于 正确性未知，只简单介绍一下。对于每一个二进制位，试图将所有该位上为 $1$ 的数连接成一个连通块。利用 $\tt LCT$ ，当新加入的边形成环的时候，就删去最短的一条边。

思路叁

使用 ~~B开头的算法（我们老师这样称呼它）~~ $\tt{Boruvka}$ 算法。简单来说就是，一个点的邻边中，边权最大的一定在最大生成树上，对于所有点都连接这条边，每次至少会连 $n\over 2$ 条边，然后缩点。最多连接 $\mathcal O(\log n)$ 轮。

要求按位与的最大值，就建 $0/1\;\tt trie$ 呗。建完之后，自底向上，把 $1$ 儿子往 $0$ 儿子里复制一份，这是套路了，因为 $\operatorname{and}$ 操作中，当另一个数字是 $0$ 时，该数字是 $0$ 或 $1$ 等效。时间复杂度是 $\mathcal O(m2^m)$ 的。