poj1222 高斯消元

给了一个01矩阵然后选在一个点1变0或者0变1 然后 与他相邻的 数也相应的变成相反的数，问最后求出一种方案把他们变成全0 将每一个位置上的状态看做一个变元，30个变元，列出30个异或方程

#include <iostream>
#include <algorithm>
#include <string.h>
#include <cstdio>
#include <vector>
using namespace std;

#define eps 1e-9
const int MAXN=;
int a[MAXN][MAXN],x[MAXN];//方程的左边的矩阵和等式右边的值，求解之后x存的就是结果
int equ,var;//方程数和未知数个数
int Gauss()
{
    int i,j,row,col,max_r;
    equ=var=;
    for(row=,col=;row<equ&&col<var;row++,col++)
    {
        max_r=row;
        for(i=row+;i<equ;i++)
          if(a[i][col])
            { max_r=i; break; }
        if(row!=max_r)
        {

            for(j=col;j<var;j++)
              swap(a[row][j],a[max_r][j]);
            swap(x[row],x[max_r]);

        }
        for(i=;i<equ;i++)
          if(i!=row&&a[i][col])
          {
              x[i]^=x[row];
              for(j=col+;j<var;j++)a[i][j]=a[i][j]^a[row][j];
              a[i][col]=;
          }
    }
    return ;
}
int dx[]={,,,-};
int dy[]={,-,,};
int idmark[][];
int main()
{
    for(int i=; i<; i++)
        for(int j=; j<; j++)
           idmark[i][j]=i*+j;
    int cas;
    scanf("%d",&cas);
    for(int cc=; cc<=cas ; cc++)
    {
       for(int i=; i<; i++)
        {
            scanf("%d",&x[i]);
        }
       memset(a,,sizeof(a));
       for(int i=; i<; i++)
       for(int j=; j<; j++){
              a[ idmark[i][j] ][ idmark[i][j] ] = ;

          for(int k=; k<; k++){

              int tx= i+dx[k];

              int ty= j+dy[k];

              if(tx>=&&ty>=&&tx<&&ty<)
              a[ idmark[i][j] ][ idmark[tx][ty] ]=;

          }
       }
       Gauss();
       printf("PUZZLE #%d\n",cc);
       int cur=;
       for(int i=; i<; i++)
       for(int j=; j<; j++)
        printf("%d%c",x[cur++],j==?'\n':' ');
    }

    return ;
}