UVA1663 净化器 Purifying Machine(二分图匹配)

Problem

Solution

• 得到了二分图匹配的提示才想到怎么做，还是太笨了 。
• 一个模板可能匹配一个串或两个串，要让匹配两个串的模板尽可能的多，建立二分图，能用一个模板匹配的串连边，设最大匹配为ans，给定的模板集合所能匹配的串数为ele，答案即为ele-ans+ans/2。
• 位运算小技巧：x第y位置0： x|=(1<<y)；x第y位置1：x&=~(1<<y)。

Code

#include<iostream>
#include<cstring>
#include<queue>
#include<cstdio>
using namespace std;
const int maxn=4005,maxm=2000000,inf=1e7;
int cnt=1,n,last[maxn],dep[maxn],ansn,m,s,t,cur[maxn],ele;
bool inq[maxn];
struct edge{
	int v,c,next;
}e[maxm*2];
inline void add(int u,int v,int c)
{
	e[++cnt].c=c;
	e[cnt].next=last[u];
	e[cnt].v=v;
	last[u]=cnt;
}
bool bfs()
{
	queue<int>q;
	memset(dep,0,sizeof(dep));
	dep[s]=1;
	q.push(s);
	while(!q.empty())
	{
		int u=q.front();
		q.pop();
		for(int i=last[u];i;i=e[i].next)
		{
			int v=e[i].v,c=e[i].c;
			if(!dep[v]&&c>0)
			{
				q.push(v);
				dep[v]=dep[u]+1;
			}
		}
	}
	return dep[t]!=0;
}
int dfs(int u,int dis)
{
	if(u==t)
	return dis;
	for(int i=cur[u];i;i=e[i].next)
	{
		cur[u]=e[i].next;
		int v=e[i].v,c=e[i].c;
		if(dep[v]==dep[u]+1&&c)
		{
			int d=dfs(v,min(dis,c));
			if(d>0)
			{
				e[i].c-=d;
				e[i^1].c+=d;
				return d;
			}
		}
	}
	return 0;
}
int dinic()
{
	int ans = 0;
	while(bfs())
	{
		for(int i=0;i<=t;i++)
		cur[i] = last[i];
		while(int d=dfs(s,inf))
		ans+=d;
	}
	return ans;
}
void init()
{
	memset(last,0,sizeof(last));
	cnt=1;
	ele=0;
	memset(inq,0,sizeof(inq));
}
int main()
{
//	freopen("1.txt","w",stdout);
	while(cin>>n>>m&&n&&m)
	{
		for(int i=1;i<=m;i++)
		{
			getchar();
			int ans1=0,ans2=0;
			for(int i=1;i<=n;i++)
			{
				ans1<<=1;ans2<<=1;
				char c=getchar();
				if(c=='1')
				{
					ans1++;
					ans2++;
				}
				if(c=='*')
				ans1++;
			}
			inq[ans1]=inq[ans2]=1;
		}
		ansn=(1<<n);
		s=2*ansn+1;t=2*ansn+2;
		for(int i=0;i<(1<<n);i++)
		{
			if(!inq[i]) continue;
			ele++;
			int tmp=1,tmp0;
			for(int j=1;j<=n;j++)
			{
				tmp0=((i|tmp)==i)?(i&(~tmp)):(i|tmp);
				tmp<<=1;
				if(!inq[tmp0]) continue;
				add(i,tmp0+ansn,1);
				add(tmp0+ansn,i,0);
				
			}
			add(s,i,1);
			add(i,s,0);
			add(i+ansn,t,1);
			add(t,i+ansn,0);
		}
		int ans=dinic();
		cout<<ele-ans+ans/2<<endl;
		init();
	}
	return 0;
}