POJ1222 高斯消元法解抑或方程

第一次学怎么用高斯消元法解抑或方程组，思想其实很简单，方法可以看下面的链接：http://blog.csdn.net/zhuichao001/article/details/5440843

有了这种思想之后，一些简单的翻牌问题也算是有了头绪，还记得之前做一到翻一次牌影响曼哈顿距离为k的点的题，现在看来是有思路，但那个貌似是900个点，不好搞呀，自己回头再想想吧。。先贴一记水题的代码

#include<iostream>
#include<cstring>
#include<cstdio>
#include<string>
#include<algorithm>
using namespace std;

int m[][];
int ans[][];
int eq[][];

void gauss(int a[][])
{
    for(int i=;i<;++i){
        int k=i;
        for(;k<;++k){
            if(a[k][i]!=){
                break;
            }
        }
        for(int j=;j<=;++j){
            swap(a[i][j],a[k][j]);
        }
        for(int j=;j<;++j){
            if(i!=j&&a[j][i]){
                for(int k=;k<=;++k){
                    a[j][k]=a[i][k]^a[j][k];
                }
            }
        }
    }
    for(int i=;i<;++i){
        ans[i/][i%]=eq[i][];
    }
}

int main()
{
    int T;cin>>T;int ca=;
    while(T--)
    {
        memset(eq,,sizeof(eq));
        for(int i=;i<;++i){
            for(int j=;j<;++j){
                scanf("%d",&m[i][j]);
                eq[i*+j][]=m[i][j];
                eq[i*+j][i*+j]=;
                if(j->=) eq[i*+j][i*+j-]=;
                if(j+<)  eq[i*+j][i*+j+]=;
                if(i*+j->=) eq[i*+j][i*+j-]=;
                if(i*+j+<=) eq[i*+j][i*+j+]=;
            }
        }
        gauss(eq);
        printf("PUZZLE #%d\n",++ca);
        for(int i=;i<;++i){
            printf("%d",ans[i][]);
            for(int j=;j<;++j){
                printf(" %d",ans[i][j]);
            }
            puts("");
        }
    }
    return ;
}