Bzoj 1923: [Sdoi2010]外星千足虫(bitset优化高斯消元)

传送门:

---

题意:给你一些虫子,每个有奇数条足或者偶数条足,有n个虫子,进行m次试验,
每次试验有一个字符串,后面跟一个字符,字符串中位置为0表示这个虫子没有
被放入,为1则相反,最后这几个虫子的足之和对2取模为后面跟着的字符,
求最后每个虫子的奇偶或者不确定,最早在第几次试验的时候可以确定

---

思路:高斯消元

---

#include<bitset>
#include<algorithm>
#include<cstdio>
using namespace std;
bitset<1010>a[2010];
int x[2010],ans=0;
//如果方程有m行,n列,调用Gauss(m,n)

int Gauss(int equ,int var){
    int i,j,k,col,temp;
    int max_r;//当前绝对值最大的行
    for(i=0;i<=var;i++)
        x[i]=0;
    for(k=0,col=0;k<equ&&col<var;k++,col++){
        max_r=k;
        for(i=k;i<equ;i++)
            if(a[i][col]){
                max_r=i;
                break;
            }
        ans=max(max_r,ans);
        // 与第k行交换.
        if(max_r!=k)
            swap(a[k],a[max_r]);
        // 说明该col列第k行以下全是0了，则处理当前行的下一列.
        if(a[k][col]==0){
            k--;
            continue;
        }
        for(i=k+1;i<equ;i++)
            if(a[i][col])   a[i]^=a[k];
    }
    // 1. 无解的情况: 化简的增广阵中存在(0, 0, ..., a)这样的行(a != 0).
    for(i=k;i<equ;i++){
        if(a[i][col]!=0)
            return -1;
    }
    //2.无穷解的情况: 在var * (var + 1)的增广阵中出现(0, 0, ..., 0)这样的行，即说明没有形成严格的上三角阵.
    // 且出现的行数即为自由变元的个数.
    if(k<var)
        return var-k;
    for(i=var-1;i>=0;i--){
        temp=a[i][var];
        for(j=i+1;j<var;j++)
            if(a[i][j]!=0)
                temp-=a[i][j]*x[j];
        temp=(temp%2+2)%2;
        x[i]=temp/a[i][i];
    }
    return 0;
}

char s[2010];

int main(){
    int n,m;
    scanf("%d%d",&n,&m);
    for(int i=0;i<m;i++){
        scanf("%s",s);
        for(int j=0;j<n;j++)    a[i][j]=s[j]-'0';
        scanf("%s",s);
        a[i][n]=s[0]-'0';
    }
    if(Gauss(m,n)>0)    printf("Cannot Determine");
    else{
        printf("%d\n",ans+1);
        for(int i=0;i<n;i++)
            if(x[i])puts("?y7M#");
            else  puts("Earth");
    }
    return 0;
}