[BZOJ]1923: [Sdoi2010]外星千足虫 高斯消元解xor方程组

Description

题解：

高斯消元解xor方程组的入门题，其实挺好学的，也是和高斯消元解普通方程一样的思路，将最后的方程变成一个倒三角形的形式，然后用bitset可以加速。

代码：

#include<cstdio>
#include<cstring>
#include<bitset>
#include<iostream>
#include<algorithm>
using namespace std;
const int Maxn=1010;
bitset<Maxn>a[Maxn*2];
int ans[Maxn],n,m,Ans=0;
bool no=false;
void work()
{
    for(int i=1;i<=n;i++)
    {
        int j;
        for(j=i;j<=m&&!a[j][i];j++);//这里不能写j=i+1! 
        if(j==m+1){no=true;return;}
        //j为第一个a[j][i]=1
        swap(a[i],a[j]);
        Ans=max(Ans,j);//因为至少到j才能确定i是什么 
        for(int k=i+1;k<=m;k++)
        if(a[k][i])a[k]^=a[i];//因为要变成倒三角，所以必须让某些位置变为0
    }
    for(int i=n;i;i--)
    {
        for(int j=i;j<=n;j++)a[i][n+1]=a[i][n+1]^(ans[j]*a[i][j]);
        ans[i]=a[i][n+1];
    }
}
int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=m;i++)
    {
        char str[Maxn];
        scanf("%s",str+1);
        for(int j=1;j<=n;j++)a[i][j]=str[j]-'0';
        scanf("%s",str+1);
        a[i][n+1]=str[1]-'0';
    }
    work();
    if(no){puts("Cannot Determine");return 0;}
    printf("%d\n",Ans);
    for(int i=1;i<=n;i++)puts((ans[i]==1)?"?y7M#":"Earth");
}