[CQOI2014] 和谐矩阵

题目描述：

QAQ…

题目分析：

a[i][j]^a[i][j-1]^a[i-1][j]^a[i+1][j]^a[i][j+1]=0
暴力构造 n*m 个异或方程，bitset优化高斯消元即可…
高消和bitset常数真小…

题目链接：

BZOJ 3503
Luogu 3164

Ac 代码：

#include <cstring>
#include <algorithm>
#include <bitset>
#include <cstdio>
#include <iostream>
const int maxm=1610;
std::bitset <maxm+100> a[maxm];
int ans[maxm];
int dx[]={0,1,-1,0,0,0},dy[]={0,0,0,1,-1,0};
inline void Gauss(int n)
{
    for(int i=1;i<=n;i++)
    {
        int j=i;
        while(j<=n&&!a[j][i]) j++;
        if(j==n+1) continue;
        std::swap(a[j],a[i]);
        for(int k=i+1;k<=n;k++)
         if(a[k][i]) a[k]^=a[i];
    }
    for(int i=n;i>=1;i--)
    if(!a[i][i]) ans[i]=1;
    else for(int j=i+1;j<=n;j++) if(a[i][j]) ans[i]^=ans[j];
}
int n,m;
int id[maxm][maxm],tot;
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      id[i][j]=++tot;
    for(int i=1;i<=n;i++)
     for(int j=1;j<=m;j++)
      for(int k=1;k<=5;k++)
      {
        int x=i+dx[k],y=j+dy[k];
        if(x<1||x>n||y<1||y>m) continue;
        a[id[i][j]][id[x][y]]=1;
      }
    Gauss(n*m);
    for(int i=1;i<=n;i++)
    {
        for(int j=1;j<=m;j++)
         printf("%d ",ans[id[i][j]]);
        puts("");
    }

    return 0;
}