（模板）poj1681 高斯消元法求异或方程组（无解、唯一解、多解）

题目链接：https://vjudge.net/problem/POJ-1681

题意：类似于poj1222，有n×n的01矩阵，翻转一个点会翻转其上下左右包括自己的点，求最少翻转多少点能使得矩阵全0。

思路：

　　同样的可以枚举第一行的状态，这里不说了。

　　用高斯消元法来解这道题，每个点的状态表示一个变量，那么有n*n个方程，n*n个变量的方程组，用高斯消元法来解，可能存在无解，唯一解，多解的情况。多解的时候要枚举自由变元的状态。

AC代码：

/*
    poj1681
    开关问题，高斯消元法解异或方程组
    求最少要翻转的开关使得矩阵全0
    存在无解，唯一解，多解的情况
    多解时要枚举自由变元的状态
   */
#include<cstdio>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<cstring>
using namespace std;

const int maxn=;
const int inf=0x3f3f3f3f;
int T,n,equ,var,a[maxn][maxn],x[maxn],free_xx[maxn];
int ans;
char s[];

void init(){    //初始化
    memset(a,,sizeof(a));
    for(int i=;i<n;++i){
        for(int j=;j<n;++j){
            int t=i*n+j;
            a[t][t]=;
            if(i>) a[t][(i-)*n+j]=;
            if(i<n-) a[t][(i+)*n+j]=;
            if(j>) a[t][i*n+j-]=;
            if(j<n-) a[t][i*n+j+]=;
        }
    }
}

int Gauss(){
    int r=,cnt=;  //cnt表示自由变元个数
    for(int c=;r<equ&&c<var;++r,++c){
        int Maxr=r;
        for(int i=r+;i<equ;++i)
            if(abs(a[i][c])>abs(a[Maxr][c]))
                Maxr=i;
        if(Maxr!=r){
            for(int i=c;i<var+;++i)
                swap(a[Maxr][i],a[r][i]);
        }
        if(!a[r][c]){
            --r;
            free_xx[cnt++]=c;
            continue;
        }
        for(int i=r+;i<equ;++i){
            if(!a[i][c]) continue;
            for(int j=c;j<var+;++j)
                a[i][j]^=a[r][j];
        }
    }
    for(int i=r;i<equ;++i)
        if(a[i][var])
            return -;    //无解
    return var-r;         //返回自由变元的个数,cnt=var-r
}

int solve(int t){
    ans=inf;
    for(int i=;i<(<<t);++i){  //枚举自由变元的状态
        int cnt=;        //要翻转的个数
        memset(x,,sizeof(x));
        for(int j=;j<t;++j){
            if((i>>j)&){
                ++cnt;
                x[free_xx[j]]=;
            }
        }
        for(int j=var-t-;j>=;--j){
            int tmp=a[j][var],tp,ok=;
            for(int k=j;k<var;++k){
                if(!a[j][k]) continue;
                if(ok){   //找主元
                    ok=;
                    tp=k;
                }
                else{
                    tmp^=x[k];
                }
            }
            x[tp]=tmp;
            cnt+=x[tp];
        }
        ans=min(ans,cnt);  //取最小
    }
    return ans;
}

int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%d",&n);
        init();
        equ=var=n*n;
        for(int i=;i<n;++i){
            scanf("%s",s);
            for(int j=;j<n;++j)
                if(s[j]=='y') a[i*n+j][n*n]=;
                else a[i*n+j][n*n]=;
        }
        int t=Gauss();
        if(t==-)
            printf("inf\n");
        else
            printf("%d\n",solve(t));
    }
    return ;
}