【uva 1663】Purifying Machine 二分图匹配

来源：互联网发布：九州邮币卡软件下载编辑：程序博客网时间：2024/05/16 14:56

一开始给你m个n长度的原串，*可以代替1或者0 于是可以匹配出一个新串，现在要你找到一个新的模板串使模板最少。

好吧，就是一个基本上裸的二分图匹配可是由于自己还是太脑残，在debug的时候，改掉一个地方以后就会有连锁反应，导致其他的地方也会错，所以下一次debug的时候一定要关注一下上下面有没有牵连例如if 和 else 语句，注意，注意‘

其次，在找到原来串的时候我们用到状态压缩1的位置在哪里出现过来判重也是一个很不错的技巧

代码比较冗杂，又不想优化了

#include<cstdio>#include<cstring>#include<iostream>#define maxn 2100using namespace std;int n,m,link[maxn],tot;bool vis[maxn];char s[1500];bool used[maxn*maxn],can[maxn][maxn];struct node{//-1代表品* int ch[15];void clear(){memset(ch,0,sizeof(ch));}bool operator ==(const node& b){int _cnt=0;for(int i=1;i<=n;i++){if(b.ch[i]!=ch[i]){_cnt++;if(_cnt>1)return false; }}return true;}}nod[maxn];int dfs(int u){for(int i=1;i<=tot;i++){if(i!=u&&!vis[i]&&can[u][i]){vis[i]=1;if(!link[i]||dfs(link[i])){link[i]=u;return 1;}}}return 0;}int main(){while(scanf("%d%d",&n,&m)!=EOF&&(n||m)){memset(link,0,sizeof(link));memset(used,0,sizeof(used));memset(can,0,sizeof(can));tot=0;for(int i=1;i<=m;i++){scanf("%s",&s);int pos=0,y=-1;for(int j=0;j<n;j++){if(s[j]=='1')pos|=1<<j;else if(s[j]=='*')y=j;}if(y==-1){if(!used[pos]){used[pos]=true;tot++;for(int j=0;j<n;j++){if(s[j]=='1')nod[tot].ch[j+1]=1;else nod[tot].ch[j+1]=0;}}}else{if(!used[pos]){used[pos]=true;tot++;for(int j=0;j<n;j++){if(s[j]=='1')nod[tot].ch[j+1]=1;if(s[j]=='0')nod[tot].ch[j+1]=0;if(s[j]=='*')nod[tot].ch[j+1]=0;}}pos|=1<<y;if(!used[pos]){used[pos]=true;tot++;for(int j=0;j<n;j++){if(s[j]=='1')nod[tot].ch[j+1]=1;if(s[j]=='0')nod[tot].ch[j+1]=0;if(s[j]=='*')nod[tot].ch[j+1]=1;}}}}/*cout<<endl;for(int i=1;i<=tot;i++){for(int k=1;k<=n;k++)printf("%d",nod[i].ch[k]);cout<<endl;}*/for(int i=1;i<=tot;i++){for(int j=1;j<=tot;j++){if(i!=j&&nod[i]==nod[j])can[i][j]=true;}}int ans=0;for(int i=1;i<=tot;i++){memset(vis,0,sizeof(vis));ans+=dfs(i);}//for(int i=1;i<=tot;i++)printf("%d : %d\n",i,link[i]);printf("%d\n",tot-ans/2);}return 0;}/*5 61011101010*10001*01010*0111*110 0*/

1 1