Description

一个黑白网格,点一次会改变这个以及与其连通的其他方格的颜色,求最少点击次数使得所有全部变成黑色.

Sol

高斯消元解异或方程组.

先建立一个方程组.

\(x_i\) 表示这个点是否被用过.

因为第二次使用同一个点,这个点的贡献就被消除了,所以每个点只会被用 0/1 次.

\(a_{ij}\) 表示 \(j\) 点对 \(i\) 是否有影响,有影响为 1 否则为 0.

最后的一位表示最后的状态^最初的状态.

这样就列出来了 \(n*m\) 个方程组,一共 \(n*m\) 个未知数.

然后求解,就跟高斯消元一样,用异或的方式消掉就可以了.

Code

#include <cstdio>
#include <cstring>
#include <iostream>
using namespace std; #define H(i,j) ((i-1)*m+j)
const int N = 35; int n,m;
int a[N][N]; void print(){
cout<<"----------------------------"<<endl;
for(int i=1;i<=30;i++){
for(int j=1;j<=31;j++) cout<<a[i][j];
cout<<endl;
}
} void init(){
n=5,m=6;
for(int i=1;i<=n;i++){
char s[N];memset(s,0,sizeof(s));
scanf("%s",s+1);
for(int j=1;j<=m;j++) a[H(i,j)][n*m+1]=1^(s[j]-'0');
// cout<<a[H(i,j)][n*m+1];
// cout<<endl<<s+1<<endl<<endl;
} for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) {
a[H(i,j)][H(i,j)]=1;
if(i+1<=n) a[H(i,j)][H(i+1,j)]=1;
if(i-1>=1) a[H(i,j)][H(i-1,j)]=1;
if(j+1<=m) a[H(i,j)][H(i,j+1)]=1;
if(j-1>=1) a[H(i,j)][H(i,j-1)]=1;
}
// print();
} void gauss(int n){
for(int i=1,r,j;i<n;i++){
for(r=j=i;j<n;j++) if(a[j][i]) { r=j;break; }
// cout<<i<<" "<<r<<endl;
if(r!=i) for(int k=i;k<=n;k++) swap(a[i][k],a[r][k]);
for(int k=i+1;k<n;k++) if(a[k][i]) {
for(j=i;j<=n;j++) a[k][j]^=a[i][j];
}
// print();
}
// print();
for(int i=n-1;i;i--){
for(int j=i+1;j<n;j++) a[i][n]^=a[j][n]*a[i][j];
} } void out(){
int ans=0;
for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(a[H(i,j)][n*m+1]) ans++;
cout<<ans<<endl;
for(int i=1;i<=n;i++) for(int j=1;j<=m;j++) if(a[H(i,j)][n*m+1]) cout<<i<<" "<<j<<endl;
} int main(){
// freopen("in.in","r",stdin);
// freopen("log.out","w",stdout);
init();
gauss(31);
out(); return 0;
}

  

hihoCoder 1196 高斯消元·二的更多相关文章

  1. hihocoder 1196 高斯消元.二

    传送门 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 在上一回中,小Hi和小Ho趁着便利店打折,买了一大堆零食.当他们结账后,看到便利店门口还有其他的活动. 店主:买了 ...

  2. hihoCoder#1196 : 高斯消元&#183;二(开关灯问题)

    传送门 高斯消元解异或方程组 小Ho在游戏板上忙碌了30分钟,任然没有办法完成,于是他只好求助于小Hi. 小Ho:小Hi,这次又该怎么办呢? 小Hi:让我们来分析一下吧. 首先对于每一个格子的状态,可 ...

  3. hiho #1196 : 高斯消元&#183;二

    #1196 : 高斯消元·二 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 在上一回中,小Hi和小Ho趁着便利店打折,买了一大堆零食.当他们结账后,看到便利店门口还有其 ...

  4. hihocoder 第五十二周 高斯消元&#183;二【高斯消元解异或方程 难点【模板】】

    题目地址:http://hihocoder.com/contest/hiho57/problem/1 输入 第1..5行:1个长度为6的字符串,表示该行的格子状态,1表示该格子是亮着的,0表示该格子是 ...

  5. HihoCoder 1195 高斯消元&#183;一(高斯消元)

    题意 https://hihocoder.com/problemset/problem/1195 思路 高斯消元是解决高元方程的一种算法,复杂度 \(O(n^3)\) . 过程大致是: 构造一个未知数 ...

  6. hihoCoder 1195 高斯消元.一

    传送门 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Ho:喂不得了啦,那边便利店的薯片半价了! 小Hi:啥?! 小Ho:那边的便利店在打折促销啊. 小Hi:走走走, ...

  7. hihoCoder #1195 高斯消元&#183;一

    题意:便利店老板为了促销,推出了组合包的形式,将不同数量的各类商品打包成一个组合.比如2袋薯片,1听可乐的组合只要5元,而1袋薯片,2听可乐的组合只要4元.通过询问老板知道:一共有N种不同的商品和M种 ...

  8. [HIHO1196]高斯消元&#183;二(高斯消元、枚举自由变元)

    题目链接:http://hihocoder.com/problemset/problem/1196 #include <bits/stdc++.h> using namespace std ...

  9. [hihoCoder] 高斯消元&#183;一 [TPLY]

    高斯消元一 题目链接 : http://hihocoder.com/problemset/problem/1195?sid=1269842 很"好aoaoaoaoaoaoa"的高斯 ...

随机推荐

  1. 嵌入式Linux驱动学习之路(二十三)NAND FLASH驱动程序

    NAND FLASH是一个存储芯片. 在芯片上的DATA0-DATA7上既能传输数据也能传输地址. 当ALE为高电平时传输的是地址. 当CLE为高电平时传输的是命令. 当ALE和CLE都为低电平时传输 ...

  2. C#如何获取物理网卡,虚拟网卡,以及无线网卡

    就不废话了,直接上代码 /// <summary></summary> /// 显示本机各网卡的详细信息 /// <summary></summary> ...

  3. HTML5 &lt;a&gt;标签的ping属性用法

    随着移动互联网用户井喷式的增长,web前端开发中的HTML5在近几年备受瞩目,越来越多的人从事html5开发相关工作.今天小编也来凑个热闹,和大家一起来谈谈HTML5中<a>标签的ping ...

  4. c语言之内存的申请malloc() 和释放free()

    c语言之内存的申请malloc() 和释放free() 1.如何使用 malloc 函数 malloc是一个函数,专门用来从堆上分配内存.使用malloc函数需要几个要求: 内存分配给谁?分配多大内存 ...

  5. OpenNebula Restfull 接口请求示例

    Fri Jun 20 07:28:20 2014 [I]: 10.0.2.2 - - [20/Jun/2014 07:28:20] "POST /vmtemplate HTTP/1.1&qu ...

  6. webservice主流框架Axis、Axis2、XFire、CXF的比较

    http://ws.apache.org/axis/ http://axis.apache.org/axis2/java/core/ http://xfire.codehaus.org/ http:/ ...

  7. 如何定义自己的ViewGroup

    在发展中,有时会遇到一些要求.布局和控制系统不仅提供使用,以满足我们的发展,所以这一次就行,通常是你自己的自定义布局(ViewGroup)并控制(View)该.我在这里,我们将用一个简单的例子,当他们 ...

  8. maven GroupId 和ArtifactId的含义

    GroupID是项目组织唯一的标识符,实际对应Java的包的结构,是main目录里java的目录结构. ArtifactID就是项目的唯一的标识符,实际对应项目的名称,就是项目根目录的名称.

  9. Minimum Inversion Number~hdu 1394

    The inversion number of a given number sequence a1, a2, ..., an is the number of pairs (ai, aj) that ...

  10. Gitlab_ansible_jenkins三剑客②Gitlab的后台管理

    系统信息和日志 健康状态 使用gitlab的用户管理和审批功能 创建用户 创建一个lead普通账号 进入test-repo仓库 这样就把dev添加到了test-repo这个项目中,并且有了develo ...