2019牛客暑期多校训练营（第五场）F maximum clique 1

题目链接：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;
}