Bubble Cup 8 finals C. Party (575C)

时间:2023-03-09 09:36:36
Bubble Cup 8 finals C. Party (575C)

题意:

给定n个人,分两天晚上去夜总会开派对,要求每天恰好有n/2个人去,且每人去的夜总会各不相同。

每个人对不同的晚上不同的夜总会有不同的满意度,求一个方案使得所有人的满意度之和最大。

夜总会数量=人的数量=n,2<=n<=20,且n是偶数。

0<=每一项满意度<=10^6。

时间限制2s,空间限制4MB。

题解:

在这题上卡了很久…

初看这题觉得是费用流…写完发现图建错了…

然后改成暴力枚举哪些人在第一天晚上去再跑费用流…

每个人只和对应晚上的夜总会连边,然后两天晚上的夜总会再连到一个辅助点,再连到汇…

就避免了跑出来的东西有两人在不同的晚上去同一个夜总会…

理所当然地超时了…

后来换成比较擅长这种图的zkw费用流…快了很多,但还是严重超时…

然后开始怀疑自己的做法…试图写DP…很快放弃了…

然后突然发现自己之前建图时傻掉了,枚举夜总会在哪个晚上有人就可以把图变成二分图…

于是改改改…改成KM…还是严重超时…大数据跑100s…

于是弃疗了…去膜拜Tourist的代码…结果发现就是KM…只不过是非递归的KM…

都是KM…差距就是这么大…

主要原因就是点是一个个加进二分图的,如果在最后跑KM,就完全没有用到子图的信息…

用非递归KM的话…一层层下去,冗余运算就减少了很多。

简单来说这题的做法就是:

枚举夜总会,一边向图中加点一边跑KM,然后就可以了…(CF评测机快,2s可过…

 1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
#include <cstdio>
#include <cstring>
#define mv(a,b) memcpy(a,b,(n<<2)+4)
inline int read()
{
int s = ; char c; while((c=getchar())<'0'||c>'9');
do{s=s*+c-'0';}while((c=getchar())>='0'&&c<='9');
return s;
}
const int N = , INF = 0x7f7f7f7f;
int n,n2,w1[N][N],w2[N][N],w[N][N],ans,tans,lx[N],ly[N],sl[N],lk[N],pat[N];
bool by[N];
void dfs(int nn,int d1,int d2)
{
if(nn>n){ if(ly[]>ans) ans = ly[]; return; }
int _lx[N],_ly[N],_lk[N],_a,i,cy,cx,ny;
for(int ch=;ch<;ch++)
{
if(d1==n2&&ch==) continue;
if(d2==n2&&ch==) continue;
mv(_lx,lx), mv(_ly,ly), mv(_lk,lk);
if(ch) mv(w[nn],w1[nn]); else mv(w[nn],w2[nn]);
memset(sl,0x7f,(n<<)+), memset(by,,n+);
cy = ; lk[] = nn;
while(lk[cy])
{
by[cy] = ; cx = lk[cy]; int tmp = INF;
for(i=;i<=n;i++)
if(!by[i])
{
int t = -lx[cx]-ly[i]-w[cx][i];
if(t<sl[i]) sl[i] = t, pat[i] = cy;
if(sl[i]<tmp) tmp = sl[i], ny = i;
}
for(i=;i<=n;i++)
{
if(by[i]) lx[lk[i]] += tmp, ly[i] -= tmp;
else sl[i] -= tmp;
}
cy = ny;
}
while(cy) lk[cy] = lk[pat[cy]], cy = pat[cy];
dfs(nn+,d1+(ch),d2+(ch^));
mv(lx,_lx), mv(ly,_ly), mv(lk,_lk);
}
}
int main()
{
int i,j;
n2 = (n = read())>>;
for(i=;i<=n;i++) for(j=;j<=n;j++) w1[i][j] = read();
for(i=;i<=n;i++) for(j=;j<=n;j++) w2[i][j] = read();
dfs(,,);
printf("%d\n",ans);
return ;
}