POJ 2724 奶酪消毒机 二分匹配 建图 比较难想

#include <cstdio>
#include <algorithm>
#include <cstring>
#include <iostream>
#include <cmath>
#include <vector>
#define MM(a,b)  memset(a,b,sizeof(a))
using namespace std;

vector<int> G[2200];
int match[2200],used[2200];
int g,b,m,cnt,l,r;
int mp[2200],posi[2200];

void add_edge(int u,int v)
{
    G[u].push_back(v);
    G[v].push_back(u);
}

bool dfs(int u)
{
     used[u]=1;
     for(int i=0;i<G[u].size();i++)
     {
       int v=G[u][i];
       int w=match[v];
       if(w<0||!used[w]&&dfs(w))
         {
             match[u]=v;
             match[v]=u;
             return true;
         }
     }
     return false;
}

int  bipartite_match()
{
    MM(match,-1);
    int res=0;
    for(int i=1;i<=1000+r;i++)
      if(match[i]<0)
      {
          memset(used,0,sizeof(used));
          if(dfs(i)) res++;
      }
    return res;
}

bool onlyonedifer(int i,int j)
{
   int c=(i^j);
   return (c&&((c&(c-1))==0));
}//处理i,j,两个数字二进制表示只有一位不同模板

bool jione(int i)
{
   int j=0;
   while(i)
   {
       if(i&1) j++;
       i>>=1;
   }
   return j%2==1;
}

void build()
{
    l=0;r=0;
    for(int i=1;i<=cnt;i++)
       {

          if(jione(mp[i]))
            {
                l++;
                posi[l]=mp[i];
            }
         else
            {
                r++;
                posi[1000+r]=mp[i];
            }
       }

    for(int i=1;i<=1000+r;i++) G[i].clear();

    for(int i=1;i<=l;i++)
        for(int j=1;j<=r;j++)
                  if(onlyonedifer(posi[i],posi[j+1000]))
                      add_edge(i,j+1000);
}

char s[15];
int main()
{
    int n,m;
    while(~scanf("%d %d",&n,&m)&&(n||m))
    {
        cnt=0;
        MM(mp,0);
        for(int i=1;i<=m;i++)
            {
                scanf("%s",s);
                int flag=0;
                ++cnt;
                for(int j=0;j<n;j++)
                {
                        if(s[j]==‘*‘)
                          {
                             flag=j+1;
                             mp[cnt]=(mp[cnt]<<1)+0;
                          }
                        else
                             mp[cnt]=(mp[cnt]<<1)+s[j]-‘0‘;
                }
                if(flag)
                {
                    mp[cnt+1]=(mp[cnt]|(1<<(n-1-(flag-1))));
                    ++cnt;
                }
            }
        sort(mp+1,mp+cnt+1);//unique处理前先排序
        cnt=unique(mp+1,mp+cnt+1)-(mp+1);//去重
        build();
        int w=bipartite_match();
        printf("%d\n",w+cnt-2*w);
    }
    return 0;
}