zoukankan      html  css  js  c++  java
  • Luogu2447/BZOJ1923 [SDOI2010]外星千足虫

    题目传送门

    算法分析

    不难看出这相当于是解\(n\)元一次异或方程组,类比\(n\)元一次方程组的解法,不难想到使用高斯消元。

    发现数据范围比较大,\(n^3\)是肯定过不去的,因此想到01矩阵bitset存储,用bitset Xor运算实现对其他方程的转化,时间复杂度\(\frac{n^3}{32}\)可以通过。

    求最早确定结果的时间,只需要再记录一个\(id[i]\)表示当前第\(i\)个方程是输入的第\(id[i]\)个方程,对每一次使用的方程的\(id\)取最大值即可。

    代码实现

    #include<bits/stdc++.h>
    using namespace std;
    #define maxn 1005
    #define maxm 2005
    template<typename Tp>
    void read(Tp &x){
    	x=0;int f=1;char c=getchar();
    	while(c>'9'||c<'0'){if(c=='-'){f=-1;}c=getchar();}
    	while(c>='0'&&c<='9'){x=(x<<1)+(x<<3)+(c^48);c=getchar();}x*=f;
    }
    int n,m;
    bitset<maxn>a[maxm];
    int id[maxm],ans;
    int main(){
    	read(n);read(m);int x;
    	for(int i=1;i<=m;i++){
    		for(int j=1;j<=n+1;j++){
    			scanf("%1d",&x);
    			a[i][j]=x;
    		}
    		id[i]=i;
    	}
    	for(int i=1;i<=n;i++){
    		int p=i;
    		while(!a[p][i]&&p<=m)p++;
    		ans=max(ans,id[p]);
    		if(p==m+1){
    			puts("Cannot Determine");
    			return 0;
    		}
    		if(i!=p){
    			swap(a[i],a[p]);swap(id[i],id[p]);
    		}
    		for(int j=1;j<=m;j++){
    			if(i==j)continue;
    			if(a[j][i])a[j]=a[j]^a[i];
    		}
    	}
    	printf("%d\n",ans);
    	for(int i=1;i<=n;i++){
    		puts(a[i][n+1]?"?y7M#":"Earth");
    	}
    	return 0;
    }
    
    
  • 相关阅读:
    实验0 了解和熟悉操作系统
    学习进度条
    0302软件构建与教学
    评论任务
    学习进度条
    sprint3个人总结
    软件工程学期总结
    6.3 学术诚信与职业道德
    阅读《构建之法》第8、9、10章
    nodejs学习心得
  • 原文地址:https://www.cnblogs.com/ZigZagKmp/p/11491105.html
Copyright © 2011-2022 走看看