【CF839E】Mother of Dragons 折半状压

【CF839E】Mother of Dragons

题意:给你一张n个点,m条边的无向图。你有k点能量,你可以把能量分配到任意一些点上,每个点分到的能量可以是一个非负实数。定义总能量为:对于所有边<a,b>,a的能量*b的能量 的和。让你最大化总能量。

$n\le 40,k\le 1000$

题解:容易发现,最后的分配方案一定是给一个大小为cnt的完全子图中的每个点都分配$k\over cnt$点能量。

那么本题就变成了一般图最大团问题,可以用随机化搞定,这里给出一种meet in the middle算法。

对于前一半的点,用状压预处理出f[S]表示S中的最大团是什么。对于后一半的点,我们可以枚举其中的所有最大团,然后把与这些点都相邻的前一半的点拿出来,用前一半点的f值+后一半点的团大小更新答案即可。

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
typedef long long ll;
int n,m,n1,n2;
double k;
const int N=(1<<20)+4;
ll v[50];
int f1[N],Log[N],g[N],f2[N],cnt[N];
int main()
{
	scanf("%d%lf",&n,&k);
	int i,j,a;
	for(i=0;i<n;i++)	for(j=0;j<n;j++)
	{
		scanf("%d",&a);
		if(a)	v[i]|=1ll<<j,v[j]|=1ll<<i;
	}
	n1=n/2,n2=n-n/2;
	for(i=0;i<n2;i++)	Log[1<<i]=i;
	for(i=1;i<(1<<n1);i++)
	{
		a=Log[i&-i];
		f1[i]=max(f1[i^(1<<a)],f1[i&v[a]]+1);
	}
	m=f1[(1<<n1)-1];
	g[0]=1,f2[0]=(1<<n1)-1;
	for(i=1;i<(1<<n2);i++)
	{
		a=Log[i&-i];
		f2[i]=f2[i^(1<<a)]&v[a+n1],cnt[i]=cnt[i^(1<<a)]+1;
		if(g[i^(1<<a)]&&(((v[a+n1]>>n1)&i)==(i^(1<<a))))	g[i]=1,m=max(m,f1[f2[i]]+cnt[i]);
	}
	printf("%.10f",k*k*(m-1.0)/m/2.0);
	return 0;
}

原文地址:https://www.cnblogs.com/CQzhangyu/p/8283333.html

时间: 2024-10-18 15:04:08

【CF839E】Mother of Dragons 折半状压的相关文章

CF839E Mother of Dragons 最大团 Bron-Kerbosch算法

题意简述 给你一个\(n\)个节点的无向图\(G=\{V,E\}\)的邻接矩阵\(g\)和每个点的点权为\(s_i\),且\(\sum_{i=1}^n s_i = K\),要你求出\(\mathrm{max} \{ \sum_{u,v \in E} s_u \times s_v\}\) 做法 设两个不相邻的点\(u\),\(v\)的点权为\(s_u\)和\(s_v\),令\(a_u = \sum_{g[u][i]=1} s_i, a_v=\sum_{g[v][i]=1} s_i\),此时这对点\

v&#183;y「状压dp」

一直对状压dp怀有一种恐惧感 不会打,不会调,关键是不会调 做了这两道题,虽然还是不会状压dp,但总比之前好了一些 y 普通状压应该很好打 复杂度$O(2^d*n*(n+m))$ for(ll i=2;i<=d;i++){ for(ll state=0;state<=maxn;state++){ for(ll x=1;x<=n;x++){ for(ll j=head[x];j;j=nxt[j]){ ll y=ver[j]; if(!f[i-1][x][state]) continue ;

ZOJ3305Get Sauce 状压DP,

状压DP的题目留个纪念,首先题意一开始读错了,搞了好久,然后弄好了,觉得DFS可以,最后超时,修改了很久还是超时,没办法看了一下n的范围,然后觉得状压可以,但是没有直接推出来,就记忆化搜索了一下,可是一直错,莫名奇妙,然后没办法看了一下题解,发现了下面这个比较好的方法,然后按照这个方程去推,然后敲,也是WA了好多把,写的太搓了,没人家的清楚明了,唉~也算是给自己留个纪念,状压一直做的都不太好~唉~还好理解了, 参考了  http://blog.csdn.net/nash142857/articl

poj 2411 Mondriaan&#39;s Dream(状压DP)

Mondriaan's Dream Time Limit: 3000MS   Memory Limit: 65536K Total Submissions: 12232   Accepted: 7142 Description Squares and rectangles fascinated the famous Dutch painter Piet Mondriaan. One night, after producing the drawings in his 'toilet series

HDU 5025 Saving Tang Monk(状压搜索)

钥匙是必须有序,蛇是不要求有序的.所以一个需要状压一个不用 因为时间计算和步数计算不同.所以要遍历整个空间,或者使用优先队列.优先时间短的. 风格就这样了. #include <iostream> #include <cstring> #include <cstdio> #include <algorithm> #include <queue> #define inf 0x3f3f3f3f #define maxn 110 using names

URAL 1326. Bottle Taps(简单的状压dp)

题目不太好读懂,就是先给你一个n代表要从n个物品中买东西,然后告诉你这n个东西的单价,在给你m个集合的情况,就是每个结合中有x件物品,他们合起来买的价格是k.这x件物品依次是:p1--px.之后给你一个kk,表示你要买的物品的编号.让你求出来如何花费最少的钱买到要求的序列. 20,可以状压啊,注意一开始的时候先把单价的状态处理出来...之后就是水题了啊. 1326. Bottle Taps Time limit: 3.0 second Memory limit: 64 MB Programmer

奇怪的道路(状压)

[Jxoi2012]奇怪的道路 时间限制: 1 Sec  内存限制: 128 MB提交: 55  解决: 23[提交][状态][讨论版] 题目描述 小宇从历史书上了解到一个古老的文明.这个文明在各个方面高度发达,交通方面也不例外.考古学家已经知道,这个文明在全盛时期有n座城市,编号为1..n.m条道路连接在这些城市之间,每条道路将两个城市连接起来,使得两地的居民可以方便地来往.一对城市之间可能存在多条道路.据史料记载,这个文明的交通网络满足两个奇怪的特征.首先,这个文明崇拜数字K,所以对于任何一

(状压dp)uva 10817 Headmaster&#39;s Headache

题目地址 1 #include <bits/stdc++.h> 2 typedef long long ll; 3 using namespace std; 4 const int MAX=1e5+5; 5 const int INF=1e9; 6 int s,m,n; 7 int cost[125]; 8 //char sta[MAX]; 9 string sta; 10 int able[125]; 11 int dp[125][1<<8][1<<8]; 12 in

【计算几何】【状压dp】Codeforces Round #226 (Div. 2) D. Bear and Floodlight

读懂题意发现是傻逼状压. 只要会向量旋转,以及直线求交点坐标就行了.(验证了我这俩板子都没毛病) 细节蛮多. #include<cstdio> #include<algorithm> #include<cstring> #include<cmath> using namespace std; const double PI=acos(-1.0); #define EPS 0.00000001 struct Point { double x,y; Point(