【AT4162 [ARC099C]】 Independence 题解

Link:(洛谷)AT4162 [ARC099C] Independence

首先在此感谢 @kymru 大佬提供的学术支持。

Solution

题意

给一张图, N N N 个点, M M M 条边,要求将图分成两个团,使得两个顶点都在同一个团中的边最少,求最小值。(团是一个两两之间有边的顶点集合。)

1

发现把所有节点分层两个团的本质就是对这些点进行黑白染色。设最后染成黑色的节点数为 a a a,白色的节点数为 b b b,则最终输出的答案即为 a × ( a − 1 ) + b × ( b − 1 ) 2 \frac{a\times(a-1)+b\times(b-1)}{2} 2a×(a1)+b×(b1)

发现直接对原图上手太复杂了,那些边也很难处理。所以考虑对原图建它的补图。根据补图的定义,知道在补图中有一条边连接的两个节点在原图中不直接联通。所以,这样一来就把题目中“团是一个两两之间有边的顶点集合”转化为“集合内的节点两两不直接联通”。说道这里,显然,就变成将补图转化为一个二分图了。

在一个二分图内,试将其某一部点都染成白色,另一部点都染成黑色。而这两部点就是所求的两个团。但是,原图的补图会有多个二分图,即不止一个连通块。那如何确定左右部点各自所染的颜色呢?

2

回过头,答案的最终形式是 a × ( a − 1 ) + b × ( b − 1 ) 2 \frac{a\times(a-1)+b\times(b-1)}{2} 2a×(a1)+b×(b1),运用简单小学数学知识,可以将它转化为 ( a + b ) × ( a + b − 1 ) − 2 a b (a+b)\times(a+b-1)-2ab (a+b)×(a+b1)2ab,即 n × ( n − 1 ) − 2 a b n\times(n-1)-2ab n×(n1)2ab。再运用小学数学知识,知道“和一定,差小积大”,所以目标就变成让 a a a b b b 的差尽可能小。

因此,在遍历一个连通块将其转化为一二分图的过程中,记录下每个二分图它的左右部点数量的差值,并全部相加,记为 s u m sum sum。接着要做的就是适当分配这些差值(分成两部分 x x x y y y),使这两部分的和,它们两的差尽可能小。绕这么多,其实就是为了 a a a b b b 的差尽可能小。

显然,我们想要 x x x y y y 都尽可能接近 s u m 2 \frac{sum}{2} 2sum。这里运用到最优性转可行性的方法——背包 dp。

3

通过简单背包,可以求得 x x x y y y 的最小差值,即 a a a b b b 的最小差值。同时,已知 a a a b b b 的和(即 n n n),组合起来,就又是一个小学数学知识,就成为了一个简单的和差问题。往下的解和差问题就不多说了。

至此,此题被完美解决,复杂度可过。

Code

#include
using namespace std;

#define rep(i, a, b) for(register int i = a; i <= b; ++i)
#define per(i, a, b) for(register int i = a; i >= b; --i)
const int maxn = 705, maxm = 500005;
int n, m, sum;
bool flg, g[maxn][maxn], f[maxn];
int cnt, hd[maxn];
struct edge{
	int to, nxt;
}e[maxm];
int del[maxn], col[maxn], cur;
int u, v, tmp;

inline void add(int u, int v){
	e[++cnt] = (edge){v, hd[u]}, hd[u] = cnt;
}

inline void dfs(int u, int t){
	if(flg) return;
	del[t] += (col[u] == 1 ? 1 : -1);
	for(int v, i = hd[u]; i; i = e[i].nxt)
		if(!col[v = e[i].to])
			col[v] = 3 - col[u], dfs(v, t);
		else if(col[v] == col[u]){
			flg = 1; return;
		}
}

int main(){
	scanf("%d%d", &n, &m);
	rep(i, 1, m)
		scanf("%d%d", &u, &v), g[u][v] = g[v][u] = 1;
	rep(i, 1, n) rep(j, 1, n)
		if(!g[i][j] and i != j) add(i, j); 
	rep(i, 1, n){
		if(!col[i])
			col[i] = 1, dfs(i, ++cur),
			del[cur] = abs(del[cur]), sum += del[cur];
		if(flg){
			printf("-1\n"); return 0;
		}
	}
	f[0] = 1;
	rep(i, 1, cur) per(j, sum / 2, del[i])
		f[j] |= f[j - del[i]];
	per(i, sum / 2, 0) if(f[i]){
		tmp = i; break;
	}
	tmp = abs(sum - tmp * 2);
	int a = (n + tmp) >> 1, b = (n - tmp) >> 1;
	printf("%d\n", a * (a - 1) / 2 + b * (b - 1) / 2);
	return 0;
} 

感谢阅读。

你可能感兴趣的:(题解,算法,图论)