UVA 1663(p381)----Purifying Machine

最新推荐文章于 2022-02-23 14:12:40 发布

wang2147483647

最新推荐文章于 2022-02-23 14:12:40 发布

阅读量560

点赞数

本文链接：https://blog.csdn.net/wang2147483647/article/details/50749264

版权

算法竞赛入门经典（第二版）----图论模型与算法专栏收录该内容

9 篇文章 0 订阅

订阅专栏

该博客主要介绍了如何解决UVA 1663题目的过程，涉及图论和位运算的知识。通过读取输入的字符串并进行位操作，构建一个邻接矩阵来表示二分图，并使用Kuhn-Munkres算法找到最大匹配，最终输出匹配的数量。博客中包含了完整的C++代码实现。

摘要由CSDN通过智能技术生成

# include<iostream>
# include<cstdio>
# include<vector>
# include<cstring>
# include<algorithm>
using namespace std;
const int maxn=3000+50;
vector<int>b;
int vis[maxn];
int n,m,a[maxn][maxn],t[maxn];
int A[11]= {1,2,4,8,16,32,64,128,256,512,1024};
void f(string p)
{
    int r=0;
    for(int i=0; i<m; ++i)
        r=r*2+p[i]-'0';
    b.push_back(r);
}
bool ok(int x,int y)
{
    int z=x^y;
    for(int i=0; i<11; ++i)
        if(z==A[i]) return true;
    return false;
}
void make()
{
    memset(a,0,sizeof(a));
    for(int i=0; i<n; ++i)
        for(int j=i+1; j<n; ++j)
            if(ok(b[i],b[j]))
            {
                a[i][j]=1;
                a[j][i]=1;
            }
}
bool match(int u)
{
    for(int v=0; v<n; ++v)
    {
        if(a[u][v]&&!vis[v])
        {
            vis[v]=true;
            if(t[v]==-1||match(t[v]))
            {
                t[v]=u;
                return true;
            }
        }
    }
    return false;
}
void solve()
{
    int ans=0;
    memset(t,-1,sizeof(t));
    for(int i=0; i<n; ++i)
    {
        memset(vis,0,sizeof(vis));
        if(match(i)) ans++;
    }
    ans/=2;
    for(int i=0; i<n; i++)
        if(t[i]==-1) ans++;
    printf("%d\n",ans);
}
int main()
{
    string p;
    while(scanf("%d%d",&m,&n)&&(n+m))
    {
        b.clear();
        for(int i=0; i<n; ++i)
        {
            cin>>p;
            int j;
            for(j=0; j<m; ++j) if(p[j]=='*') break;
            if(j==m) f(p);
            else
            {
                p[j]='0';
                f(p);
                p[j]='1';
                f(p);
            }
        }
        sort(b.begin(),b.end());
        n=unique(b.begin(),b.end())-b.begin();
        make();
        /*for(int i=0;i<2*n-1;i++)
            for(int j=0;j<2*n-1;j++)
              if(a[i][j]) cout<<i<<" "<<j<<endl;*/
        solve();
    }
    return 0;
}