[BZOJ]5001: 搞事情 高斯消元解xor方程组

Description

给定一个NM的01矩阵（原矩阵）。每次可以选择原矩阵中任意一个点“搞事情”：具体来说会把原矩阵中选中点与其上下左右共五个点取反。您要通过最少操作次数使得原矩阵数值全部变成0！然而，您输出的不是最少操作次数，而是一个NM的答案矩阵。答案矩阵中每个点代表对于原矩阵该点被“搞事情”操作的次数。也就是答案矩阵所有点数值加起来和为总操作最少次数。由于在相同操作次数下的答案矩阵可能有多个，出题人又懒的不想写SPJ！所以您需要输出字典序最小的答案矩阵。

题解：

高斯消元解xor方程组的裸题啊，判断无解的方法是一个方程组所有的系数都为0但结果为1，字典序最小判断一下就可以了，理论极限复杂度很高，但还是跑的很快。拿了暂时的rank1，不知道能持续多久……（现在已经不是了）

代码:

#include<bits/stdc++.h>
using namespace std;
#define LL long long
const int Maxn=410;
int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return x*f;
}
int n,m,ans[Maxn],Ans=500;
bitset<Maxn>a[Maxn];
bool no=false;
void gauss()
{
    for(int i=1;i<=n*m;i++)
    {
        int j;
        for(j=i;j<=n*m&&!a[j][i];j++);
        if(j>n*m)continue;
        swap(a[i],a[j]);
        for(int k=i+1;k<=n*m;k++)
        if(a[k][i])a[k]=a[k]^a[i];
    }
    for(int i=1;i<=n*m;i++)
    {
        if(!a[i][n*m+1])continue;
        bool all_zero=true;
        for(int j=1;j<=n*m;j++)
        if(a[i][j]){all_zero=false;break;}
        if(all_zero){no=true;return;}
    }
}
int mp[Maxn];
void dfs(int x,int cnt)
{
    if(cnt>Ans)return;
    if(!x)
    {
        if(cnt<Ans){Ans=cnt;for(int i=1;i<=n*m;i++)mp[i]=ans[i];return;}
        bool tf=false;
        for(int i=1;i<=n*m;i++)
        if(ans[i]<mp[i]){tf=true;break;}
        else if(ans[i]>mp[i])break;
        if(tf)
        {
            for(int i=1;i<=n*m;i++)
            mp[i]=ans[i];
        }return;
    }
    if(a[x][x])
    {
        int t=a[x][n*m+1];
        for(int i=x+1;i<=n*m;i++)t=t^(ans[i]*a[x][i]);
        ans[x]=t;
        dfs(x-1,cnt+t);
    }
    else
    {
        ans[x]=0;dfs(x-1,cnt);
        ans[x]=1;dfs(x-1,cnt+1);
    }
}
int num[25][25],cnt=0;
int tx[]={1,-1,0,0};
int ty[]={0,0,1,-1};
bool in(int x,int y){return(x<=n&&y<=m&&x&&y);}
int main()
{
    n=read();m=read();
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    num[i][j]=++cnt;
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    a[num[i][j]][n*m+1]=read(),a[num[i][j]][num[i][j]]=1;
    for(int i=1;i<=n;i++)
    for(int j=1;j<=m;j++)
    for(int k=0;k<4;k++)
    {
        int nx=i+tx[k],ny=j+ty[k];
        if(in(nx,ny))a[num[i][j]][num[nx][ny]]=1;
    }
    gauss();
    if(no){puts("IMPOSSlBLE");return 0;}
    dfs(n*m,0);
    for(int i=1;i<=n*m;i++)printf("%d%c",mp[i],((i%m==0)?'\n':' '));
}