[SDOI2010]外星千足虫（高斯消元）

高斯消元裸题。。。

方法一：暴力，O（2^n）20分

方法二：直接Gauss,加点玄学技巧搞得好的话70分

方法三：使用bitset优化，复杂度：$O(\frac{n^3}{ω})$

 

不会的同学看一下代码吧，跟裸Gauss差不多...只不过消元地方用异或消

真的是代码一看就懂了呢...

 

#include<cstdio>
#include<iostream>
#include<cstring>
#include<bitset>
#define writeln(x)  write(x),puts("")
#define writep(x)   write(x),putchar(' ')
using namespace std;
inline int read(){
    int ans=0,f=1;char chr=getchar();
    while(!isdigit(chr)){if(chr=='-') f=-1;chr=getchar();}
    while(isdigit(chr)){ans=(ans<<3)+(ans<<1)+chr-48;chr=getchar();}
    return ans*f;
}void write(int x){
    if(x<0) putchar('-'),x=-x;
    if(x>9) write(x/10);
    putchar(x%10+'0');
}const int M=2005;
bitset<M> a[M];
int n,m,ans;
inline void Gauss(){
    for(int i=1,now;i<=n;i++){now=i;
        while(!a[now][i]&&now<=m) ++now;
        if(now==m+1){return (void)(puts("Cannot Determine"));}
        ans=max(ans,now);
        if(now!=i) swap(a[now],a[i]);
        for(int j=1;j<=m;j++)
            if(i==j||!a[j][i]) continue;
            else a[j]^=a[i];
    }writeln(ans);
    for(int i=1;i<=n;i++){
        if(!a[i][n+1]) puts("Earth");
        else puts("?y7M#");
    }
}
int main(){
    n=read(),m=read();
    for(int i=1;i<=m;i++)
        for(int j=1;j<=n+1;j++){
            char chr;
            while((chr=getchar())<'0'||chr>'1');
            a[i][j]=chr-48;
        }Gauss();
    return 0;
}