xor高斯消元 poj2965 The Pilots Brothers' refrigerator

用的是kuangbin的模板，，也是66666

数组A是增广矩阵，下标从0开始，到n下标，其中n下标代表的是右边的常数

最后求得的答案在第n列中（下标从0开始）

对于这道题，先列出增广矩阵，也就是一个点 会有哪些点是否翻转而影响它，可以建立起一个这样的矩阵

1 1 1 1 1 0 0 0 1 0 0 0 1 0 0 0
1 1 1 1 0 1 0 0 0 1 0 0 0 1 0 0
1 1 1 1 0 0 1 0 0 0 1 0 0 0 1 0
1 1 1 1 0 0 0 1 0 0 0 1 0 0 0 1
1 0 0 0 1 1 1 1 1 0 0 0 1 0 0 0
0 1 0 0 1 1 1 1 0 1 0 0 0 1 0 0
0 0 1 0 1 1 1 1 0 0 1 0 0 0 1 0
0 0 0 1 1 1 1 1 0 0 0 1 0 0 0 1
1 0 0 0 1 0 0 0 1 1 1 1 1 0 0 0
0 1 0 0 0 1 0 0 1 1 1 1 0 1 0 0
0 0 1 0 0 0 1 0 1 1 1 1 0 0 1 0
0 0 0 1 0 0 0 1 1 1 1 1 0 0 0 1
1 0 0 0 1 0 0 0 1 0 0 0 1 1 1 1
0 1 0 0 0 1 0 0 0 1 0 0 1 1 1 1
0 0 1 0 0 0 1 0 0 0 1 0 1 1 1 1
0 0 0 1 0 0 0 1 0 0 0 1 1 1 1 1

其中，每一行表示对应的某个点，每行的1表示哪些点对这个点有影响，也就是在同一行或者同一列

然后呢，对于不同的样例，带上不同的常数列，如测试样例是这样的

-+--
----
----
-+--

那么我们得到的增广矩阵就应该是这样的

1 1 1 1 1 0 0 0 1 0 0 0 1 0 0 0 0
1 1 1 1 0 1 0 0 0 1 0 0 0 1 0 0 1
1 1 1 1 0 0 1 0 0 0 1 0 0 0 1 0 0
1 1 1 1 0 0 0 1 0 0 0 1 0 0 0 1 0
1 0 0 0 1 1 1 1 1 0 0 0 1 0 0 0 0
0 1 0 0 1 1 1 1 0 1 0 0 0 1 0 0 0
0 0 1 0 1 1 1 1 0 0 1 0 0 0 1 0 0
0 0 0 1 1 1 1 1 0 0 0 1 0 0 0 1 0
1 0 0 0 1 0 0 0 1 1 1 1 1 0 0 0 0
0 1 0 0 0 1 0 0 1 1 1 1 0 1 0 0 0
0 0 1 0 0 0 1 0 1 1 1 1 0 0 1 0 0
0 0 0 1 0 0 0 1 1 1 1 1 0 0 0 1 0
1 0 0 0 1 0 0 0 1 0 0 0 1 1 1 1 0
0 1 0 0 0 1 0 0 0 1 0 0 1 1 1 1 1
0 0 1 0 0 0 1 0 0 0 1 0 1 1 1 1 0
0 0 0 1 0 0 0 1 0 0 0 1 1 1 1 1 0

可以发现，在第n列中，只有第1列(代表的是点(0,1) )和第13列(代表是点(3,1) )是1，因为是+号，其他都是0

这样，我们就列出了方程组了

然后直接丢给模板，就能解出方程，然后整理下就能输出了

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;

const int MX = 20 + 5;
const int INF = 0x3f3f3f3f;

typedef bool Matrix[MX][MX];

Matrix A;
char S[MX][MX];
int ax[MX], ay[MX];

int gauss(int equ, int var) {
    int max_r, col, k;
    for(k = 0, col = 0; k < equ && col < var; k++, col++) {
        max_r = k;
        for(int i = k + 1; i < equ; i++) {
            if(A[i][col] > A[max_r][col]) {
                max_r = i;
            }
        }
        if(A[max_r][col] == 0) {
            k--;
            continue;
        }
        if(max_r != k) {
            for(int j = col; j < var + 1; j++) {
                swap(A[k][j], A[max_r][j]);
            }
        }
        for(int i = k + 1; i < equ; i++) {
            if(A[i][col] != 0) {
                for(int j = col; j < var + 1; j++) {
                    A[i][j] ^= A[k][j];
                }
            }
        }
    }
    for(int i = k; i < equ; i++) {
        if(A[i][col] != 0) return -1;
    }
    if(k < var) return var - k;

    for(int i = var - 1; i >= 0; i--) {
        for(int j = i + 1; j < var; j++) {
            A[i][var] ^= (A[i][j] && A[j][var]);
        }
    }
    return 0;
}

int main() {
    //freopen("input.txt", "r", stdin);
    memset(A, 0, sizeof(A));
    for(int i = 0; i < 4; i++) {
        for(int j = 0; j < 4; j++) {
            A[i * 4 + j][i * 4 + j] ^= 1;
            for(int k = 0; k < 4; k++) {
                A[i * 4 + j][i * 4 + k] ^= 1;
                A[i * 4 + j][k * 4 + j] ^= 1;
            }
        }
    }


    while(~scanf("%s%s%s%s", S[0], S[1], S[2], S[3])) {
        for(int i = 0; i < 4; i++) {
            for(int j = 0; j < 4; j++) {
                A[i * 4 + j][16] = S[i][j] == '+' ? 1 : 0;
            }
        }

        gauss(16, 16);
        int tot = 0;
        for(int i = 0; i < 4; i++) {
            for(int j = 0; j < 4; j++) {

                if(A[i * 4 + j][16]) {
                    tot++;
                    ax[tot] = i;
                    ay[tot] = j;
                }
            }
        }

        printf("%d\n", tot);
        for(int i = 1; i <= tot; i++) {
            printf("%d %d\n", ax[i] + 1, ay[i] + 1);
        }
    }
    return 0;
}