[SDOI2010]外星千足虫(高斯消元)

洛谷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;
}
  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值