【BZOJ1923】外星千足虫,高斯消元解xor方程组

原创 2016年08月29日 20:27:50

Time:2016.08.29
Author:xiaoyimi
转载注明出处谢谢


传送门
思路:
原本以为是高斯消元解取模方程,后来发现这题意不就和异或方程一样吗

【异或(XOR)运算由于与“奇偶性”密切相关,经常出现在有关奇偶性和二进制的题目中】
很多异或问题都要涉及到解异或方程组,因此首先要搞懂异或方程组的解法。

(1)异或方程组的格式(设有n个未知数,m个方程): a11*x1 XOR a12*x2 XOR … XOR a1n*xn = b1
a21*x1 XOR a22*x2 XOR … XOR a2n*xn = b2 …………………………………………………… am1*x1
XOR am2*x2 XOR … XOR amn*xn = bm
其中的所有a、b、x的值均为0或1。解异或方程组就是给出了所有的系数a和b之后,求出解(x1, x2 …
xn)。一般来说,题目的要求无非就是如下几种:<1>求任意一组解;<2>求解的总数;<3>求出最优解(比如字典序最小的解或者加权以后权值最大/小的解等)。

n个未知数,m个方程,上高斯消元就可以了
但是我们发现直接做的复杂度是O(n2m)
所以又要引入一种神奇的STL东东——bitset
用它做位运算好像能降低一个指数
这样就可以满足要求了
关于无解的情况:只要有一个未知数在所有的方程中系数均为0,就是无解了
代码:

#include<cstdio>
#include<iostream>
#include<cstring>
#include<bitset>
#include<string>
using namespace std;
int n,m,mx,ans[1003];
bool flag;
bitset<1003> f[2003];
char get()
{
    char ch=getchar();
    while (ch!='0'&&ch!='1') ch=getchar();
    return ch;
}
void gauss()
{
    int t;
    for (int i=1;i<=n;++i)
    {
        flag=0;
        for (int j=i;j<=m;++j)
            if (flag=f[j][i])
            {t=j;mx=max(mx,j);break;}
        if (!flag) return;
        if (t!=i) swap(f[t],f[i]);
        for (int j=i+1;j<=m;++j)
            if (f[j][i])
                f[j]^=f[i];
    }
    for (int i=n;i;--i)
    {
        t=f[i][n+1];
        for (int j=i+1;j<=n;++j) t^=ans[j]*f[i][j];
        ans[i]=t;
    }
}
main()
{
    scanf("%d%d",&n,&m);
    for (int i=1;i<=m;++i)
        for (int j=1;j<=n+1;++j)
            if (get()=='1') f[i].set(j);
    gauss();
    if (!flag) {printf("Cannot Determine");return 0;}
    printf("%d\n",mx);
    for (int i=1;i<=n;++i)
        if (!ans[i])
            puts("Earth");
        else
            puts("?y7M#");
}

版权声明:不转不是中国人!(滑稽)

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

题目描述传送门题解高斯消元求解异或方程组 和普通消元的思路差不多,同样是找到一个可用的交换,然后把原来的加减运算全部改为异或 如果a(i,i)无法交换成0的话,说明这里有可能出现无解或多解,那么回...
  • Clove_unique
  • Clove_unique
  • 2017年01月12日 08:56
  • 370

Bzoj 1923: [Sdoi2010]外星千足虫(bitset优化高斯消元)

传送门: 题意:给你一些虫子,每个有奇数条足或者偶数条足,有n个虫子,进行m次试验, 每次试验有一个字符串,后面跟一个字符,字符串中位置为0表示这个虫子没有 被放入,为1则相反,最后这几...
  • acm_fighting
  • acm_fighting
  • 2016年09月13日 20:15
  • 291

bzoj 1923: [Sdoi2010]外星千足虫 (高斯消元)

1923: [Sdoi2010]外星千足虫 Time Limit: 10 Sec  Memory Limit: 64 MB Submit: 753  Solved: 474 [Submit][S...
  • clover_hxy
  • clover_hxy
  • 2016年04月23日 10:58
  • 638

BZOJ 1923: [Sdoi2010]外星千足虫

高斯消元。。。记录一下最大用到哪个就好 有自由元就多组解咯。
  • cgh_Andy
  • cgh_Andy
  • 2017年04月25日 21:14
  • 162

【bzoj1923】[Sdoi2010]外星千足虫 高斯消元解xor方程组

DescriptionInput第一行是两个正整数 N, M。 接下来 M行,按顺序给出 Charles 这M次使用“点足机”的统计结果。每行 包含一个“01”串和一个数字,用一个空格隔开。“01”串...
  • LOI_DQS
  • LOI_DQS
  • 2016年03月17日 08:53
  • 558

bzoj1923: [Sdoi2010]外星千足虫

传送门:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1923 裸裸的高斯消元,bitset大法好... #include #in...
  • thy_asdf
  • thy_asdf
  • 2015年07月29日 16:34
  • 477

[Sdoi2010]外星千足虫

解方程的东西 无疑几个算法:变成图进行网络流/模拟高斯消元/并查集(逻辑方程) 这题用的是朴素高斯消元 实际上没有网络上说什么位压缩运算那么复杂 直接把加法改成xor就行了 附代码 //dsy1...
  • Laishao_yuan
  • Laishao_yuan
  • 2012年07月27日 09:40
  • 761

【bzoj 1923】[Sdoi2010]外星千足虫(高斯消元)

执笔生宣,浅墨入画,绘一幅丹青风华
  • reverie_mjp
  • reverie_mjp
  • 2016年04月24日 07:17
  • 264

洛谷 P2447 [SDOI2010]外星千足虫(高斯消元bitset优化)

传送门题解上午才学了高斯,下午准备做几个题实践一下,然后就看见了这个题。。 我看完题,就觉得,这题好裸啊,然后就去看数据范围,发现n...
  • stone41123
  • stone41123
  • 2017年09月09日 15:14
  • 128

bzoj1923 [Sdoi2010]外星千足虫 高斯消元

Description 对于 20%的数据,满足 N=M≤20; 对于 40%的数据,满足 N=M≤500; 对于 70%的数据,满足 N≤500,M≤1,000; 对于 100%的数...
  • jpwang8
  • jpwang8
  • 2018年01月20日 09:52
  • 30
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【BZOJ1923】外星千足虫,高斯消元解xor方程组
举报原因:
原因补充:

(最多只允许输入30个字)