poj1753 （高斯消元法求异或方程组）

题目链接：http://poj.org/problem?id=1753

题意：经典开关问题，和poj1222一样，进行两次高斯消元即可，只用初始化的时候改一下初始状态。可能存在无解或多解的情况，多解要枚举自由变元的所有状态。

AC代码：

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

const int maxn=;
const int inf=0x3f3f3f3f;
int mp[][],a[maxn][maxn],x[maxn],equ,var,free_x[maxn];
char s[];

void init(int p){
    memset(a,,sizeof(a));
    for(int i=;i<;++i){
        for(int j=;j<;++j){
            int t=i*+j;
            a[t][]=p^mp[i][j];
            a[t][t]=;
            if(i>) a[t][(i-)*+j]=;
            if(i<) a[t][(i+)*+j]=;
            if(j>) a[t][i*+j-]=;
            if(j<) a[t][i*+j+]=;
        }
    }
}

int Gauss(){
    int r=,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[r][i],a[Maxr][i]);
        }
        if(!a[r][c]){
            --r;
            free_x[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;
}

int solve(int t){
    int ret=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_x[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];
        }
        ret=min(ret,cnt);
    }
    return ret;
}

int main(){
    equ=var=;
    for(int i=;i<;++i){
        scanf("%s",s);
        for(int j=;j<;++j)
            if(s[j]=='b') mp[i][j]=;
            else mp[i][j]=;
    }
    init();
    int t1,t2,t3,t4;
    t1=Gauss();
    if(t1!=-) t2=solve(t1);
    init();
    t3=Gauss();
    if(t3!=-) t4=solve(t3);
    if(t1==-&&t3==-)
        printf("Impossible\n");
    else
        printf("%d\n",min(t2,t4));
    return ;
}