【BZOJ 3503】 [Cqoi2014]和谐矩阵|高斯消元|xor方程组

网上题解很多啦~~

很神啊~~

xor也可以gauss啊

我好弱啊~！

#include <cmath>
#include <cstdio>
#include <cstring>
#include <iostream> 
#include <algorithm>
using namespace std;
#define LL long long
const int MAXN=42;
LL init[MAXN][MAXN];
int n,m,f[MAXN][MAXN],val[MAXN][MAXN];
void gauss()
{
	for(int i=1;i<=m;i++)
	{
		bool ok=false;
		for(int j=i;j<=m;j++)
			if(f[j][i]==1)
			{
				ok=true;
				if(j!=i)
					for(int k=i;k<=m;k++)
						swap(f[i][k],f[j][k]);
				break;
			}
		if(ok==false) continue;
		for(int j=i+1;j<=m;j++)
			if(f[j][i]==1)
				for(int k=i;k<=m;k++)
					f[j][k]^=f[i][k];
	}
	for(int i=m;i>=1;i--)
	{
		if(f[i][i]==0)
		{
			val[1][i]=1;
			continue ;	
		}
		for(int j=i+1;j<=m;j++)
			val[1][i]^=f[i][j]&val[1][j];
	}
}
int main()
{
	cin >>n >>m;
	for(int i=0;i<m;i++)
		init[1][i+1]=(1ll<<i);
	for(int i=2;i<=n+1;i++)
		for(int j=1;j<=m;j++)
			init[i][j]=init[i-1][j-1] ^ init[i-1][j] ^ init[i-1][j+1] ^ init[i-2][j];
	for(int i=1;i<=m;i++)
	{
		LL tmp=init[n+1][i];
		for(int j=1;j<=m;j++)
		{
			f[i][j]=tmp&1ll;
			tmp>>=1;	
		}
	}
	gauss();
	for(int i=2;i<=n;i++)
		for(int j=1;j<=m;j++)
			val[i][j]=val[i-1][j-1] ^ val[i-1][j] ^ val[i-1][j+1] ^val[i-2][j];
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
			printf("%d ",(int)val[i][j]);
		printf("\n");
	}
	return 0;
}