【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#");
}

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

相关文章推荐

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

执笔生宣,浅墨入画,绘一幅丹青风华

【高斯消元解XOR方程组】各种题

MT神牛的ppt很酷 QW又无节操的蒯了道经典高斯消元的题(Orz) 于是准备给自己普及普及这方面知识 只要以后有这方面题就会写了贴上来 我认为比较重要或是难理解的东西有以下几个 1...

BZOJ 1923 [Sdoi2010] 外星千足虫

高斯消元解异或方程组

[题解]bzoj1923(SDOI2010)外星千足虫

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

【高斯消元解异或方程组】[CQOI2014][HYSBZ/BZOJ3503]和谐矩阵

题目 分析:一个元素只可能被它上下左右即自己影响,设mat[i][j]为当且元素的值。那么mat[i][j]^mat[i+1][j]^mat[i][j+1]^mat[i-1][j]^a[i][j-1...

BZOJ 1770 USACO 2009 Nov Lights 灯 高斯消元异或方程组

%%% http://wenku.baidu.com/link?url=Y55O7adUuZ1biB4m_9wdrJhts9KqqLwnMlI25xgxcpUok9Gv5CrRK9kV9Rxci2Eq...

【高斯消元】【图论】[Wc2011][HYSBZ/BZOJ2155]Xor

题目大意找1到n的路径中异或和最大的路径。分析首先,我们考虑这样一个问题:先看另一道题N个点M条边的边带权的无向图,求一个回路使XOR和最大(回路中的路径可以走多次)。另一道题的分析有这样一个结论(从...

[高斯消元 线性基 生成树 随机化权值Xor] BZOJ 3569 DZY Loves Chinese II

很好的建图姿势: 我们找到这个图的任意一棵生成树 然后对于每条非树边将其的权值赋为一个随机数 对于每条树边 我们将这条树边的权值设为所有覆盖这条树边的边权的异或和 那么图不连通当且仅当删除一条...

bzoj2466高斯消元求解XOR方程

http://www.lydsy.com/JudgeOnline/problem.php?id=2466 不会做,暴力- - 所以T掉不解释(n 正解是高斯消元。 预备知识:矩阵乘法,行列式的...

高斯消元解方程组

  • 2017年01月11日 16:15
  • 441B
  • 下载
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:【BZOJ1923】外星千足虫,高斯消元解xor方程组
举报原因:
原因补充:

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