bzoj 3503 和谐矩阵

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<bitset>
#define md
#define ll long long
#define inf (int) 1e9
#define eps 1e-8
#define N 50
using namespace std;
int lx[4]={-1,0,0,0},ly[4]={0,-1,0,1};
ll ju

;
bool ans

;
bitset<N> a
;
int m,n;
void gaosi()
{
	for (int i=1;i<=m;i++)
	{
		int k=i;
		while (k<=m&&!a[k][i]) k++;
		if (k>m) continue;
		swap(a[i],a[k]);
		for (int j=i+1;j<=m;j++)
		  if (a[j][i]) a[j]^=a[i];
	}
	for (int i=m;i;i--)
	{
		if (a[i][i]==0) ans[1][i]=a[i][0]=1; else ans[1][i]=a[i][0];
		for (int j=1;j<i;j++) if (a[j][i]) a[j][0]=a[i][0]^a[j][0];
	}
}

int main()
{
	scanf("%d%d",&n,&m);
	for (int i=1;i<=m;i++) ju[1][i]=(ll)1<<i;
	for (int i=2;i<=n+1;i++)
	  for (int j=1;j<=m;j++)
	  {
	  	ll now=0; 
	  	for (int k=0;k<=3;k++)
	  	{
	  		int x=i-1+lx[k],y=j+ly[k];
	  		if (x<1||y<1||x>n||y>m) continue;
	  		now^=ju[x][y];
	  	}
	  	ju[i][j]=now;
	  }
	for (int i=1;i<=m;i++)
	  for (int j=1;j<=m;j++)
	    a[i][j]=(ju[n+1][i]>>j)&1;
	gaosi();
	for (int i=2;i<=n;i++)
	  for (int j=1;j<=m;j++)
	  {
	  	ans[i][j]=0;
	  	for (int k=0;k<=3;k++)
	  	{
	  		int x=i-1+lx[k],y=j+ly[k];
	  		if (x<1||y<1||x>n||y>m) continue;
	  		ans[i][j]^=ans[x][y];
	  	}
	  }
	for (int i=1;i<=n;i++)
	{
		printf("%d",ans[i][1]);
		for (int j=2;j<=m;j++)
		  printf(" %d",ans[i][j]);
		printf("\n");
	}
	return 0;
}