3503: [Cqoi2014]和谐矩阵

3503: [Cqoi2014]和谐矩阵

链接

分析：

　　对于每个点，可以列出一个方程a[i][j]=a[i][j-1]^a[i][j+1]^a[i-1][j]^a[i+1][j]，于是可以列出n*m个方程，高斯消元，复杂度$O(n^3m^3)$。可以再bitset优化一下。

　　还有一种复杂度更优的做法：如果知道了第一行，那么整个矩阵都可以推出来了，即每个点可以有第一行的几个位置异或得到。

　　所以可以推出每一行每个点，与第一行的那些点有关系，推得时候转化以下上面的式子，使得每一个点只与上面的行有关系。

　　由于第n+1行只能全部是0，所以可以推出第n+1行每个点与第一行每个点的关系，然后列出m个方程，高斯消元，复杂度$O(m^3)$。

代码：

#include<cstdio>
#include<algorithm>
#include<cstring>
#include<iostream>
#include<cmath>
#include<cctype>
#include<set>
#include<queue>
#include<vector>
#include<map>
using namespace std;
typedef long long LL;

inline int read() {
    int x=,f=;char ch=getchar();for(;!isdigit(ch);ch=getchar())if(ch=='-')f=-;
    for(;isdigit(ch);ch=getchar())x=x*+ch-'';return x*f;
}

const int N = ;
LL a[N][N], b[N][N], c[N][N];
int n, m;

void Gauss() {
    for (int k = ; k <= m; ++k) {
        int r = k;
        while (r <= m && !a[r][k]) r ++;
        if (r > m) continue;
        if (r != k) for (int j = ; j <= m; ++j) swap(a[k][j], a[r][j]);
        for (int i = k + ; i <= m; ++i)
            if (a[i][k]) for (int j = ; j <= m; ++j) a[i][j] ^= a[k][j];
    }
    for (int i = m; i; --i) {
        if (!a[i][i]) { c[][i] = ; continue; }
        c[][i] = a[i][m + ];
        for (int j = i + ; j <= m; ++j) if (a[i][j]) c[][i] ^= c[][j];
    }
}
int main() {
    n = read(), m = read();
    for (int i = ; i <= m; ++i) b[][i] = (1ll << (i - ));
    for (int i = ; i <= n + ; ++i)
        for (int j = ; j <= m; ++j)
            b[i][j] = b[i - ][j - ] ^ b[i - ][j] ^ b[i - ][j + ] ^ b[i - ][j];
    for (int i = ; i <= m; ++i)
        for (int j = ; j <= m; ++j) a[i][j] = (b[n + ][i] >> (j - )) & ;
    Gauss();
    for (int i = ; i <= n; ++i)
        for (int j = ; j <= m; ++j)
            c[i][j] = c[i - ][j - ] ^ c[i - ][j] ^ c[i - ][j + ] ^ c[i - ][j];
    for (int i = ; i <= n; ++i, puts(""))
        for (int j = ; j <= m; ++j) printf("%lld ", c[i][j]);
    return ;
}