BZOJ3503 [Cqoi2014]和谐矩阵（洛谷P3164）

高斯消元

BZOJ题目传送门
洛谷题目传送门

一道高斯消元解异或方程组的题目。

我们把每一个位置看成一个未知数，一共$nm$个。每个未知数都有一个条件，一共$nm$个。用这$nm$个建一个方程组求解，时间复杂度$O(n^3m^3)$。可以通过。

可以发现只要知道到第一行，根据条件可以递推出所有行。而当且仅当第$m+1$行全为0时才“和谐”（不为0的话去掉这行就不“和谐”了）。

这种方法一共$n$个未知数，每一列对应的$m+1$行都有一个方程。高斯消元求解即可。

（最后我因为不会异或方程组去翻了题解。。。）

对于这道题构造异或方程组的方法大体就是给第$j$列一个$1<<j-1$的编号，递推完后第$i$行第$j$列的系数为第$m+1$行$i$列编号第$j$列。

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define N 45
using namespace std;
typedef long long LL;
const int t1[4]={-1,-1,-1,-2};
const int t2[4]={-1,0,1,0};
int n,m;
LL a[N][N],b[N][N],ans[N][N];
inline void gauss(){
    for (int i=1,x=1;i<=m;x=++i){
        while (x<=m&&!a[x][i]) x++;
        if (x>m) continue; swap(a[i],a[x]);
        for (int j=i+1;j<=m;j++)
            if (a[j][i])
                for (int k=i;k<=m;k++)
                    a[j][k]^=a[i][k];
    }
    for (int i=m;i;i--){
        ans[1][i]=a[i][i]?a[i][m+1]:1;
        if (ans[1][i])
            for (int j=1;j<=i-1;j++)
                if (a[j][i]) a[j][m+1]^=1;
    }
}
int main(){
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;i++)
        b[1][i]=(1ll<<i-1);
    for (int i=2;i<=n+1;i++)
        for (int j=1;j<=m;j++)
            for (int k=0;k<4;k++)
                b[i][j]^=b[i+t1[k]][j+t2[k]];
    for (int i=1;i<=m;i++)
        for (int j=1;j<=m;j++)
            a[i][j]=(b[n+1][i]>>j-1)&1;
    gauss();
    for (int i=2;i<=n;i++)
        for (int j=1;j<=m;j++){
            ans[i][j]=ans[i-1][j-1];
            for (int k=1;k<4;k++)
                ans[i][j]^=ans[i+t1[k]][j+t2[k]];
        }
    for (int i=1;i<=n;i++){
        for (int j=1;j<m;j++)
            printf("%lld ",ans[i][j]);
        printf("%lld\n",ans[i][m]);
    }
    return 0;
}