[BZOJ1923][Sdoi2010]外星千足虫（高斯消元解异或方程组）

Description

Input

第一行是两个正整数 $N$ ， $M$ 。 接下来 $M$ 行，按顺序给出 Charles 这 $M$ 次使用“点足机”的统计结果。每行包含一个“01”串和一个数字，用一个空格隔开。“01”串按位依次表示每只虫子是否被放入机器：如果第 $i$ 个字符是“0”则代表编号为 $i$ 的虫子未被放入，“1”则代表已被放入。后面跟的数字是统计的昆虫足数 $\bmod 2$ 的结果。 由于 NASA 的实验机器精确无误，保证前后数据不会自相矛盾。即给定数据一定有解。

Output

在给定数据存在唯一解时有 $N+1$ 行，第一行输出一个不超过 $M$ 的正整数 $K$ ，表明在第 $K$ 次统计结束后就可以确定唯一解；接下来 $N$ 行依次回答每只千足虫的身份，若是奇数条足则输出“?y7M#”（火星文），偶数条足输出“Earth”。如果输入数据存在多解，输出“Cannot Determine”。 所有输出均不含引号，输出时请注意大小写。

Sample Input

3 5
011 1
110 1
101 0
111 1
010 1

Sample Output

4
Earth
?y7M#
Earth

HINT

对于 20% 的数据，满足 $N=M\le20$ ；
对于 40% 的数据，满足 $N=M\le500$ ；
对于 70% 的数据，满足 $N\le500,M\le1000$ ；
对于 100% 的数据，满足 $N\le1000,M\le2000$ 。

Solution

把 $N$ 只虫子看作 $N$ 个 0/1 变量。
那么问题就转化成， $M$ 个方程，每个方程都是一些变量异或得到的结果。
套用高斯消元的模板解异或方程组即可，如果有自由元则有多解。
但题目还有一个要求：
求：至少知道前多少个方程就能得出唯一解。
由于高斯消元的过程中需要不停地交换两行，因此还需要记一个 $ans$ ，假设现在对第 $i$ 个变量进行消元，找到第 $i$ 行到第 $M$ 行中，满足第 $i$ 个数为 $1$ 的第一个行号 $j$ ，那么就把 $ans$ 与 $j$ 取 $\max$ 。
此外，高斯消元 $O(N^2M)$ 的复杂度是过不去的，因此需要用 bitset 压位。

Source

#include <cmath>
#include <bitset>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
#define For(i, a, b) for (i = a; i <= b; i++)
using namespace std;
const int N = 1005, M = 2005;
int n, m; bitset<N> a[M]; char s[N];
int solve() {
    int i, j; int res = -1; For (i, 1, n) {
        int pos = -1; For (j, i, m) if (a[j][i]) {pos = j; break;}
        if (pos == -1) return -1; res = max(res, pos);
        if (pos != i) swap(a[i], a[pos]);
        For (j, 1, m) if (j != i && a[j][i]) a[j] ^= a[i];
    }
    return res;
}
int main() {
    int i, j; cin >> n >> m; For (i, 1, m) {
        scanf("%s", s + 1); For (j, 1, n) a[i][j] = s[j] - '0';
        scanf("%s", s + 1); a[i][n + 1] = s[1] - '0';
    }
    int ans = solve(); if (ans == -1) return puts("Cannot Determine"), 0;
    cout << ans << endl; For (i, 1, n) puts(a[i][n + 1] ? "?y7M#" : "Earth");
    return 0;
}