洛谷2447
模2相当于异或运算,将高斯消元的板子改一改就行了。
注意到n的范围到1000,于是可以考虑用bitset优化,复杂度为数组的1/32。
高斯消元查找交换的两行时,若发现某一列上(即对于某一只虫子)每行均为0,可直接判断无解,否则记录下在当前行之后最小的为1的一行,再对每一列记录的数据求最大值,即为所需方程个数。
#include<bitset>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int n,m;
bitset<1010> f[2005];
inline void swap(bitset<1010> &a,bitset<1010> &b){
bitset<1010> t=a;
a=b,b=t;
}
inline int max(int a,int b){
return a>b?a:b;
}
int main()
{
scanf("%d %d",&n,&m);
for(int i=1;i<=m;++i)//读入
{
int x;
for(int j=1;j<=n;++j)
{
scanf("%1d",&x);
f[i][j]=x;
}
scanf("%1d",&x);
f[i][n+1]=x;
}
int w,ans=0;
for(int i=1;i<=n;++i)//Gauss-Jordan模板
{
w=i;
while(w<=m&&!f[w][i]) ++w;
if(w==m+1){
printf("Cannot Determine"); return 0;
}
ans=max(ans,w);//记录所需方程最大值
swap(f[w],f[i]);
for(int j=1;j<=m;++j)
{
if(i==j||!f[j][i]) continue;
f[j]^=f[i];
}
}
printf("%d\n",ans);
for(int i=1;i<=n;++i)
if(f[i][n+1]) printf("?y7M#\n");
else printf("Earth\n");
return 0;
}