POJ 1222 EXTENDED LIGHTS OUT （异或消元）

分析：一共有$30$个开关，$30$个灯，设$a[i][j]$代表第$j$个开关对第$i$个灯的贡献，那么$a[i][j]=1$或者$0$，其中$a[i][i]=1$，周围$4$个灯的贡献为$1$，其他都是$0$，用$x[i]$表示第$i$个开关有无按下，$c[i]$表示第$i$个灯的初始状态，那么对于第$i$个灯，得到异或方程：

$$c[i]\ xor\ (a[i][0]*x[0])\ xor\ (a[i][1]*x[1])\ xor\ ...\ xor\ (a[i][29]*x[29])=0\\ \implies{(a[i][0]*x[0])\ xor\ (a[i][1]*x[1])\ xor\ ...\ xor\ (a[i][29]*x[29])=0\ xor\ c[i]}$$

得到一个异或方程组，这个时候异或消元，找出方程组的解即可。

这里也不需要按照高斯消元来写（上一次我就是按照高斯消元写的。。。），详见代码：

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#define LL long long
#define FOR(i,x,y)  for(int i = x;i < y;++ i)
#define IFOR(i,x,y) for(int i = x;i > y;-- i)

using namespace std;

const int maxn = 30;

int a[maxn];
int n,mat[5][6],c[30];

void init(){
    FOR(i,0,5){
        FOR(j,0,6)  scanf("%d",&mat[i][j]),c[6*i+j] = 0^mat[i][j];
    }
    FOR(i,0,5){
        FOR(j,0,6){
            int u = i*6+j;
            a[u] = (1<<(30-u))^c[u];
            if(i > 0)   a[u] ^= (1<<(36-u));
            if(i < 4)   a[u] ^= (1<<(24-u));
            if(j > 0)   a[u] ^= (1<<(31-u));
            if(j < 5)   a[u] ^= (1<<(29-u));
        }
    }
}

int xorguass(int n){
    int row = 0;
    for(int i = 30;i >= 0;-- i){
        int j;
        for(j = row;j < n;++ j){
            if(a[j] & (1<<i))   break;
        }
        if(j != n){
            swap(a[j],a[row]);
            for(j = 0;j < n;++ j){
                if(j == row)    continue;
                if(a[j] & (1<<i))   a[j] ^= a[row];
            }
        }
        ++ row;
    }
    return row;
}

int ans[5][6];

void work(){
    int n = xorguass(30);
    FOR(i,n,30) a[i] = 0;
    FOR(i,0,30){
        ans[i/6][i%6] = a[i]%2;
    }
    FOR(i,0,5){
        FOR(j,0,5)  printf("%d ",ans[i][j]);
        printf("%d\n",ans[i][5]);
    }
}

int main()
{
    //freopen("test.in","r",stdin);
    int T,tCase = 0;    scanf("%d",&T);
    while(T--){
        printf("PUZZLE #%d\n",++tCase);
        init();
        work();
    }
    return 0;
}