题目

P2447 [SDOI2010]外星千足虫

解析

sol写到自闭,用文字描述描述了半个小时没描述出来,果然还是要好好学语文

用高斯消元求解异或方程组。

因为

  • \(奇数\bigoplus奇数=偶数\)
  • \(偶数\bigoplus偶数=偶数\)
  • \(奇数\bigoplus偶数=奇数\)

\(0\)为偶数,\(1\)为奇数,

  • \((奇数+奇数)\mod 2=0\)
  • \((偶数+偶数)\mod 2=0\)
  • \((奇数+偶数)\mod 2=1\)

若把第一个里面的奇偶数分别换成\(1\)和\(0\),则对于\((x_1+x_2)\bmod 2\)的操作,可以看做异或操作(\(x_1\bigoplus x_2\))。

易证,\((x_1+x_2+x_3+···+x_n)\mod 2 = x_1\bigoplus x_2\bigoplus x_3 \bigoplus ···\bigoplus x_n\)。

对于主元所在的列,我们只让主元行上的数为\(1\),其余的为\(0\),于是我们让每一行与当前主元行比较,若某一行的这个数为\(1\),就让这一行异或主元行。

因为我们之前处理当前主元行以上的内容时,把除了当时主元行上的所有当时主元所在列上的数都异或成了\(0\),所以我们当前主元行主元之前的数都为0,根据异或的性质,发现当前主元行前面的\(0\)对之前处理的行没有影响,这样更新到最后,我们会得到一个单位矩阵,

其实手动一模拟就出来了。。

如$$\begin{bmatrix}

0&1&1&\mid&1\

1&0&1&\mid&0\

0&0&1&\mid&1

\end{bmatrix}

\to\begin{bmatrix}

1&0&1&\mid&0\

0&1&1&\mid&1\

0&0&1&\mid&1

\end{bmatrix}

\to\begin{bmatrix}

1&0&0&\mid&1\

0&1&0&\mid&0\

0&0&1&\mid&1

\end{bmatrix}\$$

于是就得到了答案。

代码

#include <bits/stdc++.h>
using namespace std;
const int N = 2010;
bitset<N> a[N];
char s[N];
int n, m, ans; template<class T>inline void read(T &x) {
x = 0;int f = 0;char ch = getchar();
while(!isdigit(ch)) f |= (ch == '-'),ch = getchar();
while(isdigit(ch)) x = x * 10 + ch -'0',ch = getchar();
x = f ? -x : x;
return ;
} void Gauss() {
for (int i = 1; i <= n; ++i) {
int k = i;
while (!a[k][i] && k <= m) k++;
if (k == m + 1) {
ans = -1;
return;
}
ans = max(ans, k);
if (k != i) swap(a[k], a[i]);
for (int j = 1; j <= m; ++j)
if (j != i && a[j][i]) a[j] ^= a[i];
}
return;
} int main() {
read(n), read(m);
for (int i = 1, x; i <= m; ++i) {
scanf("%s", s);
for (int j = 0; j < n; ++j) a[i][j + 1] = s[j] - '0';
read(x);
a[i][n + 1] = x;
}
Gauss();
if (ans == -1) printf("Cannot Determine\n");
else {
printf("%d\n", ans);
for (int i = 1; i <= n; ++i)
printf(a[i][n + 1] == 1 ? "?y7M#\n" : "Earth\n");
}
return 0;
}

P2447 [SDOI2010]外星千足虫 (高斯消元)的更多相关文章

  1. Luogu P2447 [SDOI2010]外星千足虫 高斯消元

    链接 给出的条件是异或类型的方程,可以直接用bitset优化高斯消元. 至于求K,在高斯消元时记录用到的最大的方程的编号即可. 代码: // luogu-judger-enable-o2 #inclu ...

  2. BZOJ 1923: [Sdoi2010]外星千足虫 [高斯消元XOR]

    1923: [Sdoi2010]外星千足虫 对于 100%的数据,满足 N≤1,000,M≤2,000. 裸高斯消元解异或方程组 给定方程顺序要求用从上到下最少的方程,那么找主元时记录一下最远找到哪个 ...

  3. 【BZOJ1923】[Sdoi2010]外星千足虫 高斯消元

    [BZOJ1923][Sdoi2010]外星千足虫 Description Input 第一行是两个正整数 N, M. 接下来 M行,按顺序给出 Charles 这M次使用“点足机”的统计结果.每行 ...

  4. BZOJ1923:[SDOI2010]外星千足虫(高斯消元)

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

  5. BZOJ.1923.[SDOI2010]外星千足虫(高斯消元 异或方程组 bitset)

    题目链接 m个方程,n个未知量,求解异或方程组. 复杂度比较高,需要借助bitset压位. 感觉自己以前写的(异或)高斯消元是假的..而且黄学长的写法都不需要回代. //1100kb 324ms #i ...

  6. BZOJ 1923: [Sdoi2010]外星千足虫 高斯消元+bitset

    高斯消元求解异或方程组,可以多学一下 $bitset$ 在位运算中的各种神奇操作. #include <cstdio> #include <bitset> #define N ...

  7. 【BZOJ-1923】外星千足虫 高斯消元 + xor方程组

    1923: [Sdoi2010]外星千足虫 Time Limit: 10 Sec  Memory Limit: 64 MBSubmit: 766  Solved: 485[Submit][Status ...

  8. LG2447/BZOJ1923 「SDOI2010」外星千足虫 高斯消元

    问题描述 LG2447 BZOJ1923 题解 显然是一个高斯消元,但是求的东西比较奇怪 发现这个方程组只关心奇偶性,于是可以用一个\(\mathrm{bitset}\)进行优化,用xor来进行消元操 ...

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

    高斯消元解异或方程组 #include <iostream> #include <algorithm> #include <cstdio> #include < ...

随机推荐

  1. [.NET] 使用 Senparc.Weixin 接入微信公众号开发:简单实现自动回复

    使用 Senparc.Weixin 接入微信公众号开发:简单实现自动回复 目录 一.前提 二.基本配置信息简析 三.配置服务器地址(URL) 四.请求处理 一.前提 先申请微信公众号的授权,找到或配置 ...

  2. asp.net mvc前台显示带htm标签的解决办法(Razor —@Html.Raw())

    数据是从后台富文本编辑后丢在数据库后取出的,不加Html.Raw(),前台就会把Html标签一同显示

  3. 通用mapper认识和用法

    目录 0. 认识 1. 导包 2. mybatis的config文件:mybatis-mapper-config.xml 3. spring与mybatis整合配置文件:mybatis.xml 4. ...

  4. Spring 中 IoC 容器简介

    IoC 是一种通过描述来生成或者获取对象的技术,可以说 Spring 是一种基于 IoC 容器编程的框架 在一个系统中可以生成各种对象,并且这些对象都需要进行管理.为了描述这些对象关系,我们需要一个容 ...

  5. Spring注解IOC/DI(4)

    2019-03-08/11:10:17 演示:使用注解的方式完成注入对象中的效果 注解参考链接:https://www.cnblogs.com/szlbm/p/5512931.html Spring中 ...

  6. 日志分析工具Log Parser介绍

    摘要: 微软动态CRM专家罗勇 ,回复321或者20190322可方便获取本文,同时可以在第一间得到我发布的最新博文信息,follow me! 分析Dynamics 365 Customer Enga ...

  7. Android View的重绘过程之Layout

    博客首页:http://www.cnblogs.com/kezhuang/p/ View绘制的三部曲,测量,布局,绘画现在我们分析布局部分测量部分在上篇文章中已经分析过了.不了解的可以去我的博客里找一 ...

  8. kotlin 第一个Android项目

    一.创建过程 二.TextView点击事件 class MainActivity : AppCompatActivity() { lateinit var tv:TextView; //初始化Text ...

  9. Canvas 绘图学习笔记2

    1 绘制文本 fillText(string,x,y,maxWidth)  //填充试绘制文本 strokeText(string,x,y,maxWidth) 画线试绘制文本 设置字体样式: cont ...

  10. 使用免费 mongodb数据库 + 免费node.js服务器搭建小程序接口

    由于微信的小程序只支持不带端口的域名接口,不支持IP地址和接口,所以我们需要映射到80端口并绑定备案过的域名才能被微信小程序访问到.简单点就是接口需要 https 协议才行,找了许久的免费的数据库与n ...