UVA 12035 War Map
让我先把微笑送给出题人
这个题最基础的一个想法:先找出一个度数和为总度数和的1/2的点集，然后判断这个点集和这个点集的补集能否形成二分图。但是就算我们把判断的复杂度看成O(1)，这个算法的复杂度也是 O(T * 2^n)，根本过不了。
下面我就来列举一下我这两个月做的优化:
1.把度数序列排序。
这个看起来没啥用，但是后面好多地方都要用到这个。
2.meet in the middle。
我们可以先搜出前一半的所有点集，把它们存在以它们的度数和为下标的链表里；
然后我们再搜后一半，再从 以 总度数/2 - 搜出的点集度数和 为下标的链表里找点集 ，然后这两个点集加起来就是二分图的左边.
虽然不知道这个效果咋样，但肯定是比直接搜点集要优很多的。
3.网络流判断。
这个其实不能算什么优化，因为判断能不能形成二分图的时候，我实在想不到什么其他办法23333。
也算是开了一下脑洞把，，，，搜索+网络流也是玄学。。。
4.询问判重构。
直接把排完序的度数序列表示成23进制数，这样就可以直接hash(可以用map记，因为这里的log 是并行的不会记到复杂度里)判重构了，如果一组询问以前出现过那么直接输出以前计算出来的答案。
5.选出的点集判重构。
这个可不能用map记了，这里如果多个log是要记在复杂度里的，，，很可能就被卡了2333.
但稍微想一想就能发现，这个是可以用二进制表示的。我们只需要让每种度数对应一个二进制位，并且不产生冲突即可。
比如说 1,2,2,3,5这个序列吧，我们只需要把它们的对应二进制位设成 2^0,2^1,2^1,2^3,2^4，就可以了。
6.二分图左右点集的最大度数和点数关系。
比如说二分图左边的最大度数是6 ，而右边只有5个点的话，是肯定不能形成合法二分图的，所以我们就提前判断而不要跑一遍网络流。
而且感觉上述6种优化缺一不可，，，因为中间14次UNACCEPT大部分就是因为少了某种优化。。。
后记：某位同校大佬还想到了一个剪枝，那就是当 总度数/2＞最多可能的边数的时候，直接判不可行。
#include<bits/stdc++.h>
#define ll long long
#define pb push_back
using namespace std;
vector<int> g[35];
map<int,int> mmp;
struct lines{
int to,flow,cap;
}l[2333];
int S,T,t,d[35],q[35],tp,tl,cur[35];
bool v[35],used[1100005];
inline void add(int from,int to,int cap){
l[++t]=(lines){to,0,cap},g[from].pb(t);
l[++t]=(lines){from,0,0},g[to].pb(t);
}
inline bool BFS(){
memset(v,0,sizeof(v));
q[tp=tl=1]=S,v[S]=1,d[S]=0;
int x; lines e;
while(tp<=tl){
x=q[tp++];
for(int i=g[x].size()-1;i>=0;i--){
e=l[g[x][i]];
if(e.flow<e.cap&&!v[e.to]){
v[e.to]=1,d[e.to]=d[x]+1;
q[++tl]=e.to;
}
}
}
return v[T];
}
int dfs(int x,int A){
if(x==T||!A) return A;
int flow=0,f,sz=g[x].size();
for(int &i=cur[x];i<sz;i++){
lines &e=l[g[x][i]];
if(d[e.to]==d[x]+1&&(f=dfs(e.to,min(A,e.cap-e.flow)))){
A-=f,flow+=f;
e.flow+=f,l[g[x][i]^1].flow-=f;
if(!A) break;
}
}
return flow;
}
inline int max_flow(){
int an=0;
while(BFS()){
memset(cur,0,sizeof(cur));
an+=dfs(S,1<<30);
}
return an;
}
int G[2333],sum[2333],num,zt[1100005];
int Q,n,de[25],cnt,ci[33],HD,val[25];
int hd[445],ne[2333],HF,now,win;
inline bool can(int s){
int tt=0;
for(int i=0;i<n;i++) if(ci[i]&s) tt+=ci[val[i]];
if(used[tt]) return 0;
used[tt]=1,zt[++num]=tt;
int OT=0,OM=0,ZT=0,ZM=0;
for(int i=0;i<n;i++)
if(ci[i]&s) OT++,OM=max(OM,de[i]);
else ZT++,ZM=max(ZM,de[i]);
if(OT<ZM||ZT<OM) return 0;
S=0,T=n+1,t=-1;
for(int i=0;i<=T;i++) g[i].clear();
for(int i=0;i<n;i++)
if(ci[i]&s) add(S,i+1,de[i]);
else add(i+1,T,de[i]);
for(int i=0;i<n;i++) if(ci[i]&s)
for(int j=0;j<n;j++) if(!(ci[j]&s)) add(i+1,j+1,1);
return max_flow()==HD;
}
inline void front_search(){
for(int i=0;i<HF;i++) sum[ci[i]]=de[i];
for(int i=0,N,L;i<ci[HF];i++){
if(i){
N=i&-i,L=N^i;
sum[i]=sum[N]+sum[L];
}
G[++cnt]=i,ne[cnt]=hd[sum[i]],hd[sum[i]]=cnt;
}
}
inline void second_search(){
for(int i=0;i<n-HF;i++) sum[ci[i]]=de[HF+i];
for(int i=0,N,L;i<ci[n-HF];i++){
if(i){
N=i&-i,L=N^i;
sum[i]=sum[N]+sum[L];
}
if(HD-sum[i]>=0){
N=HD-sum[i];
int x;
for(int j=hd[N];j;j=ne[j]){
x=G[j];
if(can(x+i*ci[HF])){
win=1;
return;
}
}
}
}
}
inline void solve(int CA){
for(int i=1;i<=num;i++) used[zt[i]]=0;
num=cnt=now=win=0,memset(hd,0,sizeof(hd));
scanf("%d",&n),HF=n>>1,HD=0;
for(int i=0;i<n;i++) scanf("%d",de+i),HD+=de[i];
sort(de,de+n),printf("Case %d: ",CA);
val[0]=0;
for(int i=1;i<n;i++)
if(de[i]==de[i-1]) val[i]=val[i-1];
else val[i]=i;
for(int i=0;i<n;i++) now=now*23+de[i];
if(mmp.count(now)){
if(mmp[now]==1) puts("YES");
else puts("NO");
return;
}
if(HD&1){
puts("NO");
return;
}
HD>>=1;
if(HD>HF*(n-HF)){
puts("NO");
return;
}
front_search();
second_search();
mmp[now]=win;
if(win) puts("YES");
else puts("NO");
}
int main(){
// freopen("warmap.in","r",stdin);
// freopen("warmap.out","w",stdout);
ci[0]=1;
for(int i=1;i<=20;i++) ci[i]=ci[i-1]<<1;
scanf("%d",&Q);
for(int i=1;i<=Q;i++) solve(i);
return 0;
}
