【uva 1663】Purifying Machine 二分图匹配

最新推荐文章于 2021-11-10 20:36:45 发布

原创最新推荐文章于 2021-11-10 20:36:45 发布 · 326 阅读

0 ·

CC 4.0 BY-SA版权

ac之路专栏收录该内容

303 篇文章

订阅专栏

这篇博客讨论了如何解决UVA 1663编程题，即寻找最小子串（模板串）以实现原串的匹配。博主指出问题本质为二分图匹配，并分享了在调试过程中遇到的挑战，强调了调试时要注意连锁反应，特别是if和else语句的影响。此外，文章提到了使用状态压缩来避免重复的一种技巧，但提到代码较为冗长且未进行优化。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

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

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

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

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

#include<cstdio>
#include<cstring>
#include<iostream>
#define maxn 2100
using 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 6
10111
01010
*1000
1*010
10*01
11*11
0 0

*/