题目链接:https://ac.nowcoder.com/acm/contest/885/F
题意:
给你5000个1e9范围内的数,现在要你在这些数中取一个子集,使得这个子集中任意两个元素用二进制表示有至少两个不同的位,要你输出这个子集的大小以及值。
做法:
图论选手居然被卡在计算几何上...这样的二分匹配题还是强悍的队友做的= =..老了老了。。
我们会发现,两个数不能同时存在的数字只能是有一位不一样的,那么一个数不能同时兼容的数最多只会存在30个(但事实上5000个数不会都是这种极限数据的),那么假设我们对两个不能兼容的数进行连边,那么题目就转化成了求最小独立集,最小独立集=点数-最小点覆盖,最小点覆盖=最大匹配数,那么个数我们就能轻而易举的算出。
在这个基础上,我们会发现点集并不好求(后来还是问了队友,心塞塞),结果是,我们会发现,我们对于一个没有连边的点,让它去找一个非匹配的边上的点,再用那个点去找匹配点,找到的那个匹配点再去找非匹配的边上的点....如此循环往复(前提是我们要先对点进行分类,即按照1的奇偶性进行是否开始匹配,否则后来答案不好找,因为二分图本来就是对两个不同的集合进行匹配),最后我们会发现,答案就是另一边vis到的有匹配的点和开始匹配的这一侧的没匹配到的点。
#include<bits/stdc++.h>
#define rep(i,a,b) for(int i=(int)a;i<=(int)b;i++)
using namespace std;
const int maxn=5005;
const int maxm=2*maxn*30;
int n,a[maxn],head[maxn],obj[maxn];
int nex[maxm],to[maxm],cnt,vis[maxn];
void add(int u,int v){
to[cnt]=v,nex[cnt]=head[u];
head[u]=cnt++;
}
int match(int x){
for(int i=head[x];~i;i=nex[i]){
int v=to[i];
if(!vis[v]){
vis[v]=1;
if(obj[v]==0||match(obj[v])){
obj[x]=v;
obj[v]=x;
return 1;
}
}
}
return 0;
}
int ck(int x){
return (x&(-x))==x;
}
void dfs(int u,int f){
for(int i=head[u];~i;i=nex[i]){
int v=to[i];
if(!vis[v]&&f==0&&v!=obj[u]){
vis[v]=1;
dfs(v,1);
}
if(!vis[v]&&f&&v==obj[u]){
vis[v]=1;
dfs(v,0);
}
}
}
int main(){
memset(head,-1,sizeof(head));
scanf("%d",&n);
rep(i,1,n){
scanf("%d",&a[i]);
rep(j,1,i-1){
if(ck(a[i]^a[j])){
add(i,j); add(j,i);
}
}
}
int sum=0,ans;
rep(i,1,n){
if(__builtin_popcount(a[i])%2==0) continue;
memset(vis,0,sizeof(vis));
if(match(i)) sum++;
}
ans=n-sum;
printf("%d\n",ans);
memset(vis,0,sizeof(vis));
for(int i=1;i<=n;i++){
if(__builtin_popcount(a[i])%2==0&&obj[i]==0){
vis[i]=1;
dfs(i,0);
}
}
int ff=0;
rep(i,1,n){
if(__builtin_popcount(a[i])%2==0&&vis[i]){
if(!ff) ff=1;
else printf(" ");
printf("%d",a[i]);
}
else if(__builtin_popcount(a[i])%2==1&&!vis[i]){
if(!ff) ff=1;
else printf(" ");
printf("%d",a[i]);
}
}
printf("\n");
return 0;
}