BZOJ1923:[Sdoi2010]外星千足虫 (高斯消元+二进制压位)

题目传送门:http://www.lydsy.com/JudgeOnline/problem.php?id=1923


题目分析:一开始还很好奇什么是异或高斯消元,后来发现就是用高消解异或方程组。

普通的n3高消理论上来说是过不了这题的,虽然实际上能过。观察一下,发现高斯消元的主要时间在于选定某一行作为主元后,将后面几行的对应位置变成0。而本题由于是异或操作,可以用bitset或者位运算加速(由于我不会用bitset,所以写了后者QAQ)。高消的时候记录一下主元原本来自第几行,取个max即可。时间复杂度O(n2m60)

高斯消元的时候,如果信息不足,会出现无解(常数不为0),或者无限解的情况。如果信息多余(即存在无用的线性方程组),也可能因为不符合当前解而导致无解。虽然本题保证了数据不会有矛盾,但还是要考虑清楚这些问题。

最后说一点:纯高斯消元的题,过了样例基本就能A了。


CODE:

#include<iostream>
#include<string>
#include<cstring>
#include<cmath>
#include<cstdio>
#include<cstdlib>
#include<stdio.h>
#include<algorithm>
using namespace std;

const int maxn=2010;
const int maxm=40;
typedef long long LL;

LL a[maxn][maxm];
int ans[maxn];
int n,m;

int Gauss()
{
    int k=0;
    int R,r=1,c=0;
    while ( r<=m && c<n )
    {
        int x=c/60;
        int y=c%60;
        R=m+1;
        for (int i=r; i<=m; i++)
            if (a[i][x]&(1LL<<y))
            {
                R=i;
                break;
            }

        if (R==m+1) r--;
        else
        {
            k=max(k,R);
            for (int i=0; i<=n/60; i++) swap(a[R][i],a[r][i]);
            for (int i=r+1; i<=m; i++)
                if (a[i][x]&(1LL<<y))
                    for (int j=0; j<=n/60; j++) a[i][j]^=a[r][j];
        }

        r++;
        c++;
    }

    if (r<=n) return 0;
    int x=n/60;
    int y=n%60;
    for (int i=n; i>=1; i--)
    {
        ans[i]=((bool)( a[i][x]&(1LL<<y) ));
        for (int j=i; j<n; j++)
            ans[i]^=( ((bool)( a[i][j/60]&(1LL<<(j%60)) ))*ans[j+1] );
    }
    return k;
}

int main()
{
    freopen("1923.in","r",stdin);
    freopen("1923.out","w",stdout);

    scanf("%d%d",&n,&m);
    for (int i=1; i<=m; i++)
        for (int j=0; j<=n; j++)
        {
            char c=getchar();
            while ( c!='0' && c!='1' ) c=getchar();
            int x=c-'0';
            int p=j/60;
            int q=j%60;
            if (x) a[i][p]|=(1LL<<q);
        }

    int temp=Gauss();
    if (!temp) printf("Cannot Determine\n");
    else
    {
        printf("%d\n",temp);
        for (int i=1; i<=n; i++)
            if (ans[i]) printf("?y7M#\n");
            else printf("Earth\n");
    }

    return 0;
}
发布了160 篇原创文章 · 获赞 76 · 访问量 10万+
展开阅读全文

没有更多推荐了,返回首页

©️2019 CSDN 皮肤主题: 技术黑板 设计师: CSDN官方博客

分享到微信朋友圈

×

扫一扫,手机浏览