[ACNOI2021]《普林斯普的呼唤》卷三-CSDN博客

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

题目

题目背景
在破除了一些封印之后，仍然有很多的封印。看来普通的方法已经行不通了，还有最后一个办法：血祭。以吾之血，浇灌足下之土地，盛放死亡之花。

题目描述
有 $n$ 朵血之花。其中有 $m$ 对 $a_i,b_i)$ ，能够提供血祭的能量。即，如果在第 $i$ 朵血之花处，浇灌 $\omega_i$ 单位的鲜血，那么血祭的能量就会增加 $\omega_{a_i}\times \omega_{b_i}$ 。

你只能浇灌最多 $x$ 单位的鲜血，需要最大化血祭的能量。注意 $\omega_i$ 可以为实数。

数据范围与提示
$n\leqslant 40$ 且 $m\leqslant\frac{n(n-1)}{2}$ ，而 $x\leqslant 100$ 。

思路

用调整法可知，可以让两个 $\omega\ne 0$ 的值满足周围点的 $\sum\omega$ 相同。

然而我没想到的是：如果二者不相邻，那么即使周围点的 $\sum\omega$ 相同，也可以让其中一个为 $0$ ，另一个为 $\omega_1+\omega_2$ 。也就是说，所有非零点都相邻。

相当于就是个团呗！那么大小为 $k$ 的团的最优方案肯定是 $\omega={x\over k}$ ，因为任意点的贡献是 $\omega_i(x-\omega_i)$ ，求和即 $x^2-\sum\omega_i^2$ ，显然是均分 $\omega$ 最优。同时我们也看到，答案就是 $x^2(1-{1\over k})$ ，这是关于 $k$ 单调的，所以问题转为 求最大团。

很显然， $n\leqslant 40$ 就是要我们折半搜索。时间复杂度 $\mathcal O(n2^{n/2})$ 。

代码

注意上面那个式子中，每一对被算了两次，输出时除以二即可。

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

const int MAXN = 20;
int ga[1<<MAXN], gb[1<<MAXN], gm[1<<MAXN];
int dpa[1<<MAXN], dpb[1<<MAXN];

inline void fuckZxy(int g[],int n){
	rep(S,1,(1<<n)-1) g[S] = g[S&-S] bitand g[S&(S-1)];
}
void getBool(int g[],int dp[],int n){
	rep(S,(dp[0]=1),(1<<n)-1)
		dp[S] = dp[S&(S-1)] and (S&g[S&-S]) == (S&(S-1));
}
void getInt(int dp[],int n){
	rep(S,0,(1<<n)-1){
		dp[S] = (-dp[S])&(__builtin_popcount(S));
		rep(i,0,n-1) getMax(dp[S],(S>>i&1)*dp[S^(1<<i)]);
	}
}

int main(){
	freopen("nanami.in","r",stdin);
	freopen("nanami.out","w",stdout);
	int n = readint(), m = readint();
	int x = readint(), ans = 0;
	const int HALF = n>>1;
	for(int a,b; m; --m){
		a = readint()-1, b = readint()-1;
		if(a > b) std::swap(a,b);
		if(b < HALF) // in left part
			ga[1<<a] |= (1<<b), ga[1<<b] |= (1<<a);
		else if(HALF <= a){ // in right part
			gb[1<<(a-HALF)] |= (1<<(b-HALF));
			gb[1<<(b-HALF)] |= (1<<(a-HALF));
		}
		else gm[1<<(b-HALF)] |= (1<<a);
	}
	getBool(ga,dpa,HALF); getInt(dpa,HALF);
	gm[0] = (1<<HALF)-1; fuckZxy(gm,n-HALF);
	getBool(gb,dpb,n-HALF);
	for(int S=0; S<(1<<(n-HALF)); ++S) if(dpb[S])
		getMax(ans,__builtin_popcount(S)+dpa[gm[S]]);
	printf("%.6f\n",x*x*(1-double(1)/ans)/2);
	return 0;
}