poj 1222 EXTENDED LIGHTS OUT高斯消元。

把每个格子的状态看成一个变元。

列出30个状态方程。 高斯消元 求解异或方程即可。

#include <vector>
#include <list>
#include <map>
#include <set>
#include <deque>
#include <stack>
#include <cstring>
#include <bitset>
#include <algorithm>
#include <functional>
#include <numeric>
#include <utility>
#include <sstream>
#include <iostream>
#include <iomanip>
#include <cstdio>
 
#include <cstdlib>
#include <ctime>
#include <assert.h>
#include <queue>
#define REP(i,n) for(int i=0;i<n;i++)
#define TR(i,x) for(typeof(x.begin()) i=x.begin();i!=x.end();i++)
#define ALLL(x) x.begin(),x.end()
#define SORT(x) sort(ALLL(x))
#define CLEAR(x) memset(x,0,sizeof(x))
#define FILLL(x,c) memset(x,c,sizeof(x))
using namespace std;
 
const double EPS = 1e-8;

#define LL long long 
#define pb push_back

//const int maxn = 30;
int a[31][31];
int ans[31];
bool left(int i){
    if(i%6 ==0)return false;
    return true;
}
bool right(int i){
    if(i%6 ==5) return false;
    return true;
}
bool up(int i){
    if(i<=5)return false;
    return true;
}
bool down(int i){
    if(i>=24)return false;
    return true;
}
void init(){
    for(int i=0;i<30;i++){
        a[i][i] = 1;
        if(left(i)){
            a[i-1][i] =1;
        }
        if(right(i)){
            a[i+1][i] =1;
        }
        if(up(i)){
            a[i-6][i] =1;
        }
        if(down(i)){
            a[i+6][i] =1;
        }
    }
}
void guass(){
    for(int i=0;i<30;i++){
        int k = i;
        for(;k<30;k++){
            if(a[k][i]>0)break;
        }
        for(int j=0;j<=30;j++){
            swap(a[i][j],a[k][j]);
        }
        for(int j=0;j<30;j++){
            if(i!=j&&a[j][i]){
                for(int k =0;k<=30;k++){
                    a[j][k] = a[j][k]^a[i][k];
                }
            }
        }
    }
    
}
int main(){
    int t ;
    cin >>t;
    int cas =0 ;
    while(t--){
        cas ++ ;
        printf("PUZZLE #%d\n",cas);
            init();
            for(int i=0;i<30;i++){
                    scanf("%d",&a[i][30]);
                    ans[i] = 0;
            }
            guass();
            for(int i=0;i<30;i++){
                    ans[i] = a[i][30];
                    printf("%d",ans[i]);
                    if(i%6 ==5){
                        printf("\n");
                    }else{
                        printf( " ");
                    }
            }
        
    }
    return 0;
}