BZOJ 3503 高斯消元

思路：

高斯消元就好啦

注意每个格子最多只能和4个相邻

所以是 n*m*n*m*5 的 并不会TLE

//By SiriusRen
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int n,m,a[44][44],xx[]={0,0,1,-1,0},yy[]={1,-1,0,0,0},eli[1666][1666],b[1666],ans[1666];
int main(){
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            a[i][j]=1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            for(int k=0;k<=4;k++){
                int dx=i+xx[k],dy=j+yy[k];
                if(a[dx][dy]){
                    int tempa=(i-1)*m+j,tempb=(dx-1)*m+dy;
                    eli[tempa][tempb]=1;
                }
            }
    for(int i=1;i<=n*m;i++){
        bool flag=0;
        for(int j=i;j<=n*m;j++){
            if(eli[j][i]){
                for(int k=1;k<=n*m;k++)
                    swap(eli[i][k],eli[j][k]);
                flag=1;break;
            }
        }
        if(!flag)continue;
        for(int j=i+1;j<=n*m;j++){
            if(eli[j][i])
                for(int k=1;k<=n*m;k++)
                    eli[j][k]^=eli[i][k];
        }
    }
    for(int i=n*m;i;i--){
        ans[i]=eli[i][i]?b[i]:1;
        if(ans[i])for(int j=1;j<=i-1;j++)if(eli[j][i])b[j]^=1;
    }
    for(int i=1;i<=n*m;i++){
        printf("%d ",ans[i]);
        if(i%m==0)putchar('\n');
    }
}