Bzoj3503--Cqoi2014和谐矩阵

#include<bits/stdc++.h>
#define INF 1000000000
#define LNF 100000000000000ll
#define eps 1e-9
#define LL long long
#define MOD 1000000007
inline int _max(int a,int b) {return a>b?a:b;}
inline double _fabs(double a) {return a>0?a:-a;}

using namespace std;
#define MAXN 45
#define MAXM 1605

LL xr[MAXN][MAXN],gs[MAXN];
int m,n,num[MAXN],mp[MAXN][MAXN];

void Gauss() {
    int p,now=1;
    for(int i=1;i<=m;i++) {
        p=now;
        while(p<=m&&!(gs[p]>>i&1)) p++;
        if(p>m) continue;
        swap(gs[p],gs[i]);swap(num[p],num[i]);
        for(int j=1;j<=m;j++) if(i!=j&&gs[j]>>i&1) 
            gs[j]^=gs[i];
        now++;
    }
    p=m;
    for(int i=m;i;i--) {
        if(!(gs[i]>>i&1)) {mp[1][i]=1;continue;}
        for(int j=i+1;j<=m;j++) if(gs[i]>>j&1) mp[1][i]^=mp[1][j];
    }
    for(int i=2;i<=n;i++) 
        for(int j=1;j<=m;j++) 
            mp[i][j]=mp[i-1][j]^mp[i-1][j-1]^mp[i-1][j+1]^mp[i-2][j];
}

int main() {
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++) xr[1][i]=1ll<<i;
    for(int i=2;i<=n+1;i++) 
        for(int j=1;j<=m;j++) 
            xr[i][j]^=xr[i-1][j]^xr[i-1][j-1]^xr[i-1][j+1]^xr[i-2][j];
    for(int i=1;i<=m;i++) gs[i]=xr[n+1][i],num[i]=i;
    Gauss();
    for(int i=1;i<=n;i++) {
        printf("%d",mp[i][1]);
        for(int j=2;j<=m;j++) printf(" %d",mp[i][j]);
        printf("\n");
    }
    return 0;
}