POJ1222 EXTENDED LIGHTS OUT 高斯消元 XOR方程组

http://poj.org/problem?id=1222

在学校oj用搜索写了一次，这次写高斯消元，haoi现场裸xor方程消元没写出来，真实zz。

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
#include<map>
using namespace std;
#define LL long long
int a[35]={};
int x[35][35]={};
int id[35]={};
int main(){
    int T;scanf("%d",&T);
    for(int k=1;k<=T;++k){
        int y;memset(x,0,sizeof(x));memset(id,0,sizeof(id));
        for(int i=1;i<=30;i++){
            scanf("%d",&a[i]);y=i%6;if(y==0)y=6;
            x[i][0]=a[i]; x[i][i]=1;
            if(i-6>0)x[i][i-6]=1;
            if(y>1)x[i][i-1]=1;
            if(i+6<=30)x[i][i+6]=1;
            if(y<6)x[i][i+1]=1;
        }
        for(int i=1;i<30;i++){
            if(!x[i][i]){
                for(int j=i+1;j<=30;j++){
                    if(x[j][i]){
                        for(int w=i;w<=30;w++){
                            swap(x[j][w],x[i][w]);
                        }swap(x[j][0],x[i][0]);
                        break;
                    }
                }
            }
            if(!x[i][i])continue;
            for(int j=i+1;j<=30;j++){
                if(!x[j][i])continue;
                for(int w=i;w<=30;w++){
                    x[j][w]^=x[i][w];
                }x[j][0]^=x[i][0];
            }
        }
        for(int i=30;i>=1;i--){
            int z=0;
            for(int j=30;j>i;j--)z^=(x[i][j]*id[j]);
            if(x[i][i]==0)id[i]=0;
            else{
                if(z==x[i][0]){
                    id[i]=0;
                }
                else{
                    id[i]=1;
                }
            }
        }
        printf("PUZZLE #%d
",k);
        for(int i=1;i<=5;i++){
            for(int j=1;j<=6;j++){
                printf("%d ",id[(i-1)*6+j]);
            }printf("
");
        }
    }
    return 0;
}