2019暑假牛客多校五 F maximum clique 1

https://ac.nowcoder.com/acm/contest/885/F
题意：给出n个不同的数字，求最大的一个子集，满足子集内任意两个数的二进制形式至少有两位不同。
思路：n个数字均不同，所以任意两个数二进制形式不可能相同，那么至少2位不同的对立面就是恰好1位不同。我们把每个数字看作一个点，恰好一位不同就连一条边，题目转化为无向图求最多点，使得任意两点不相邻。但这个问题是难以解决的，需要转化为二分图的最大独立集，事实上，与网格模型转化（依据（行+列）%2分成两列）类似，这个问题依据二进制形式中1的个数%2为0或1分成两列。因为%2同为1或0一定不相邻，只有分别为0和1才可能相邻。
因此，子集大小就是二分图最大独立集。
输出方案的问题，参考https://www.cnblogs.com/jianglangcaijin/p/6035945.html。
从Y列的非匹配点出发，依次走非匹配边，匹配边，非匹配边，匹配边…最后以匹配边结尾，并把路径上的点均标记。这样，{Y列未标记+X列标记}标记就是最小覆盖，{Y列标记+X列未标记}就是最大独立集。

#include <bits/stdc++.h>
using namespace std;
const int maxn=5500;
const int INF=0x3f3f3f3f;

struct Edge{
	int from,to,cap,flow;
};

struct Dinic{
	int n,n1,n2,m,s,t;
	vector<Edge> edges;
	vector<int> G[maxn];
	bool vis[maxn];
	int d[maxn];
	int cur[maxn];
	int val[maxn];
	vector<int> L,R;
	bool book[maxn];

	void init()
	{
		int x,cnt=0;
		cin>>n;
		s=0,t=n+1;
		for(int i=1;i<=n;i++)
		{
			scanf("%d",&x);
			if(__builtin_popcount(x)&1)n1=++cnt,AddEdge(s,n1,1),L.push_back(n1),val[n1]=x;
			else n2=++cnt,AddEdge(n2,t,1),R.push_back(n2),val[n2]=x;
		}
	 	for(int i=0;i<L.size();i++)for(int j=0;j<R.size();j++)
		{
			int l=L[i],r=R[j];
			if(__builtin_popcount(val[l]^val[r])==1)AddEdge(l,r,1);
		}

	}

	void AddEdge(int f,int t,int c)
	{
		edges.push_back((Edge){f,t,c,0});
		edges.push_back((Edge){t,f,0,0});
		m=edges.size();
		G[f].push_back(m-2);
		G[t].push_back(m-1);
	}

	bool bfs()
	{	
		memset(vis,0,sizeof(vis));
		queue<int> Q;
		Q.push(s);
		d[s]=0;
		vis[s]=1;
		while(!Q.empty())
		{
			int x=Q.front();Q.pop();
			for(int i=0;i<G[x].size();i++)
			{
				Edge& e=edges[G[x][i]];
				if(!vis[e.to] && e.cap>e.flow)
				{
					vis[e.to]=1;
					d[e.to]=d[x]+1;
					Q.push(e.to);
				}
			}
		}
		return vis[t];
	}

	int dfs(int x,int a)
	{
		if(x==t || a==0)return a;
		int flow=0,f;
		for(int& i=cur[x];i<G[x].size();i++)
		{
			Edge& e=edges[G[x][i]];
			if(d[x]+1==d[e.to] && (f=dfs(e.to,min(a,e.cap-e.flow)))>0)
			{
				e.flow+=f;
				edges[G[x][i]^1].flow-=f;
				flow+=f;
				a-=f;
				if(!a)break;
			}
		}
		return flow;
	}

	int MaxFlow()
	{
		int flow=0;	
		while(bfs())
		{
			memset(cur,0,sizeof(cur));
			flow+=dfs(s,INF);
		}
		return flow;
	}

	void dfs2(int u,bool flag)
	{
		book[u]=1;
		for(int i=0;i<G[u].size();i++)
		{
			Edge& e=edges[G[u][i]];
			if(flag==1&&e.flow==0&&e.to!=t&&book[e.to]==0)dfs2(e.to,0);
			else if(flag==0&&e.flow==1&&e.to!=s&&book[e.to]==0)dfs2(e.to,1);
		}
	}

	void print()
	{
		for(int i=0;i<m;i+=2)
		{
			Edge& e=edges[i];
			if(e.to==t&&e.flow==0)dfs2(e.from,1);
		}
		for(int i=0;i<L.size();i++)if(book[L[i]]==0)printf("%d ",val[L[i]]);
		for(int i=0;i<R.size();i++)if(book[R[i]]==1)printf("%d ",val[R[i]]);
	}
}ans;

int main()
{
	//freopen("input.in","r",stdin);
	ans.init();	
	cout<<ans.n-ans.MaxFlow()<<"\n";
	ans.print();
	return 0;
}