BZOJ 3503 高斯消元

来源：互联网发布：护肤品淘宝店产品下架编辑：程序博客网时间：2024/06/06 16:27

思路：
高斯消元就好啦
注意每个格子最多只能和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');    }}

0 0