UVA1663/ POJ - 2724 Purifying Machine(二分图)

题目比较巧妙，给了一堆模板串，*可以被替换成0/1，经过替换可以得到一堆01串，问你能不能减少模板串，比如：
*01
100
011
这三个模板串可以得到四个串 001 101 100 011 ，但是这四个串实际上可以由 10–和0–1两个模板串得到。

这题有两个点，一个是二分图建图问题和二分图匹配问题
首先不容易想到是二分图(多做一些可能就能容易点)，那么如何二分图建图呢，我们可以将同一个模板得到的两个字符串进行连线，这样进行匹配的时候，一条连线就是一个匹配，就是一个模板串。
建图的一个巧妙方法使用位运算转化成数字，具体方法如下(另一种方法也可以写个函数直接将01串变成数字，记录*的位置，存两次，但不如这个方法直接)：

for (int i = 0; i < m; i++) {
            string s;
            getchar();
            cin >> s;
            int res = 0;
            int pos = -1;
            for (int j = 0; j< s.length(); j++) { //从前往后或者从后往前都行，因为数字是不重复的，虽然实际的数字应该从后往前遍历
                if(s[j] == '1') res |= 1<<j;
                if(s[j] == '*') pos = j;
            }
            num[res] = 1;
            if(pos != -1) num[res|(1<<pos)] = 1;
            //cout << res <<  " " << (res|(1<<pos)) << endl;
        }
        for (int i = 0; i < MAXN; i++) {  
            if(num[i]) {
                cnt ++;
                for (int j = 0; j < n; j++) {//遍历每一位，看是否有对应的串，1000看是否存在0000 1100 1010 1001 
                    if(num[i^(1<<j)]) {//异或，将这一位反过来
                        map_[i][i^(1<<j)] = 1;
                    }
                }
            }
        }

建图完成后使用匈牙利算法就可以解决匹配问题了，得到匹配结果后，由于我们建图不是严格的二分图建图，因此我们得到的匹配数目实际上是两倍，所以我们要剪掉一半，比如

得到的是双倍的匹配，所以匹配数要除2，最后的答案就是所有的串的数量-匹配数/2,因为没匹配到的串本身就是一个模板串
完整代码：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>

using namespace std;

const int MAXN = 2500+10;
int map_[MAXN][MAXN];
int num[MAXN];
int n, m, e;
int T[MAXN], Left[MAXN];
int ans = 0;
int cnt = 0;//一共几个串

void init() {
    cnt = 0;ans =0;
    memset(num, 0, sizeof(num));
    memset(T, 0, sizeof(T));
    memset(Left, -1, sizeof(Left));
    memset(map_, 0, sizeof(map_));
//    cout << Left[100] << endl;
}

bool DFS(int u) {
    for (int i = 0; i < MAXN; i++) {
        if(map_[u][i] && !T[i]) {
            T[i] = 1;
            if(Left[i] == -1 || DFS(Left[i])) {
                Left[i] = u;
                return true;
            }
        }
    }
    return false;
}

int main() {
    while(cin >> n >> m) {
        if((m+n) == 0) return 0;
        init();
        for (int i = 0; i < m; i++) {
            string s;
            getchar();
            cin >> s;
            int res = 0;
            int pos = -1;
            for (int j = 0; j< s.length(); j++) { //从前往后或者从后往前都行，因为数字是不重复的，虽然实际的数字应该从后往前遍历
                if(s[j] == '1') res |= 1<<j;
                if(s[j] == '*') pos = j;
            }
            num[res] = 1;
            if(pos != -1) num[res|(1<<pos)] = 1;
            //cout << res <<  " " << (res|(1<<pos)) << endl;
        }
        for (int i = 0; i < MAXN; i++) {
            if(num[i]) {
                cnt ++;
                for (int j = 0; j < n; j++) {
                    if(num[i^(1<<j)]) {//异或，将这一位反过来
                        map_[i][i^(1<<j)] = 1;
                    }
                }
            }
        }
        int cns = 0;// 匹配数目
        for (int i = 0; i < MAXN; i++) {
            memset(T, 0, sizeof(T));
            if(DFS(i)) cns ++;
        }
        ans = cnt-cns/2;
        cout << ans << endl;
    }
}