poj 1224

题意:有一个5 * 6的矩阵，每个位置表示灯，1表示灯亮，0表示灯灭。 
然后如果选定位置i，j点击，则位置i，j和其上下左右的灯的状态都会反转。 
现在要你求出一个5 * 6的矩阵，1表示这个灯被点击过，0表示没有。 
要求这个矩阵能够使得原矩阵的灯全灭。

题解：构造一个30 X 31增广矩阵（一个有30个目标结果，所以方程组的个数为30，然后可用的变量个数为30 ，增广一下就是30 X 31了）然后异或高斯就可以了。

ac代码：

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std;
int dir[][]={{,},{,-},{,},{-,}};
int mp[][];
int x[];
int a[][];
void init()
{
    memset(x,,sizeof(x));
    memset(a,,sizeof(a));
}
int check(int x,int y)
{
    if(x<= || x>= || y<= || y>=) return -;
    return ;
}
void Debug(int n,int m)
{
    for(int i=;i<=n;i++)
    {
        for(int j=;j<=m;j++) cout<<a[i][j]<<' ';
        cout<<endl;
    }
}
void gauss()
{
    int max_r;
    int col,k;
    int i,j;
    for(int z=;z<=;z++)
    {
        for(int zz=;zz<=;zz++)
        {
            i=(z-)*+zz;// row
            j=(z-)*+zz;
            a[i][]=mp[z][zz];
            a[i][j]=;
            for(int low=;low<;low++)// 里面的一一对应 求集中注意力写/。。。
            {
                int xx=z+dir[low][];
                int yy=zz+dir[low][];
                if(check(xx,yy)==) //attention !! 这个条件都写错了，， 里面不成立写的是-1 所以外面要注意一点。。
                {
                    j=(xx-)*+yy;
                  //  cout<<low<<' '<<xx<<' '<<yy<<endl;
                    a[i][j]=;
                }
            }
        }
    }

   // Debug(30,31);
    col=k=;
    while( k<= && col<=) // 阶梯化
    {
        max_r=k;
        for(int u=k;u<=;u++) if(a[u][col])
        {
            max_r=u; // find max
            break;
        }
        if(a[max_r][col]!=)
        {
            if(max_r!=k) // not the frist ont swap it
            {
                for(int u=;u<=;u++) swap(a[max_r][u],a[k][u]);
            }
            for(int u=k+;u<=;u++)
            {
                if(a[u][col])
                {
                    for(int z=;z<=;z++) a[u][z]^=a[k][z];
                }
            }
            k++;
        }
        col++;
    }
    //Debug(30,31);
    int temp=;
    for(int z=;z>=;z--)
    {
       x[z]=a[z][];
       for(int zz=z+;zz<=;zz++) if(a[z][zz])
       {
           x[z]^=x[zz];
       }
    }
}
int main()
{
    int t;
    scanf("%d",&t);
    int Case=;
    while(t--)
    {
        for(int i=;i<=;i++)
        {
            for(int j=;j<=;j++) scanf("%d",&mp[i][j]);
        }
        init();//a x
        gauss(); // 还有一个点 就是方程的构造问题
        printf("PUZZLE #%d\n",++Case);
        for(int i=;i<=;i++)
        {
            if(i%!=) printf("%d ",x[i]);
            else printf("%d\n",x[i]);
        }
    }
    return ;
}