[CQOI2014]和谐矩阵
嘟嘟嘟
遇到这种(看似)构造的题,我好像一般都做不出来……
然而这题正解是高斯消元解异或方程组……
首先我们容易列出式子a[i][j] ^ a[i - 1][j] ^ a[i + 1][j] ^ a[i][j - 1] ^ a[i][j + 1] = 0。于是我们列出所有像这样的\(n * m\)个式子,然后\(O((nm) ^ 3)\)高斯消元加bitset优化就过了。
讲真我还不会高斯消元解异或方程组,就现学了一下。其实就是把运算改成了异或,然后bitset可以把一个一个消改成一行和一行消。
#include<cstdio>
#include<iostream>
#include<cmath>
#include<algorithm>
#include<cstring>
#include<cstdlib>
#include<cctype>
#include<vector>
#include<stack>
#include<queue>
#include<ctime>
#include<bitset>
using namespace std;
#define enter puts("")
#define space putchar(' ')
#define Mem(a, x) memset(a, x, sizeof(a))
#define In inline
typedef long long ll;
typedef double db;
const int INF = 0x3f3f3f3f;
const db eps = 1e-8;
const int maxn = 1605;
inline ll read()
{ll ans = 0;char ch = getchar(), last = ' ';while(!isdigit(ch)) last = ch, ch = getchar();while(isdigit(ch)) ans = (ans << 1) + (ans << 3) + ch - '0', ch = getchar();if(last == '-') ans = -ans;return ans;
}
inline void write(ll x)
{if(x < 0) x = -x, putchar('-');if(x >= 10) write(x / 10);putchar(x % 10 + '0');
}int n, m;
const int dx[] = {-1, 0, 1, 0, 0}, dy[] = {0, 1, 0, -1, 0};
bitset<maxn> f[maxn];In int num(int x, int y)
{return (x - 1) * m + y;
}int ans[maxn];
In void Gauss(int n)
{for(int i = 1; i <= n; ++i){int pos = i;while(pos <= n && !f[pos][i]) ++pos;if(pos == n + 1) continue;swap(f[i], f[pos]);for(int j = i + 1; j <= n; ++j) if(f[j][i]) f[j] ^= f[i];}for(int i = n; i; --i)if(!f[i][i]) ans[i] = 1;else for(int j = i + 1; j <= n; ++j) if(f[i][j]) ans[i] ^= ans[j];
}int main()
{n = read(), m = read();for(int i = 1; i <= n; ++i)for(int j = 1; j <= m; ++j)for(int k = 0; k <= 4; ++k){int x = i + dx[k], y = j + dy[k];if(x < 1 || x > n || y < 1 || y > m) continue;f[num(i, j)][num(x, y)] = 1;}// for(int i = 1; i <= n; ++i) cout << f[i].to_string() << endl;Gauss(n * m);for(int i = 1; i <= n; ++i){for(int j = 1; j <= m; ++j) write(ans[num(i, j)]), space;enter;}return 0;
}
转载于:https://www.cnblogs.com/mrclr/p/10461668.html
[CQOI2014]和谐矩阵相关推荐
- bzoj3503: [Cqoi2014]和谐矩阵
高斯消元解异或方程组.学了bitset.对比如下 #include<cstdio> #include<cstring> #include<iostream> #in ...
- bzoj 3503: [Cqoi2014]和谐矩阵(高斯消元)
3503: [Cqoi2014]和谐矩阵 Time Limit: 10 Sec Memory Limit: 128 MBSec Special Judge Submit: 1101 Solved ...
- bzoj3503【CQOI2014】和谐矩阵
3503: [Cqoi2014]和谐矩阵 Time Limit: 10 Sec Memory Limit: 128 MBSec Special Judge Submit: 841 Solv ...
- JZOJ 3596. 【CQOI2014】和谐矩阵
Description 我们称一个由0和1组成的矩阵是和谐的,当且仅当每个元素都有偶数个相邻的1.一个元素相邻的元素包括它本身,及他上下左右的4个元素(如果存在). 给定矩阵的行数和列数,请计算并输出 ...
- 【BZOJ3503】【Cqoi2014】和谐矩阵 高斯消元,解异或方程组
#include <stdio.h> int main() {puts("转载请注明出处");puts("地址:blog.csdn.net/vmurder/a ...
- [暑假的bzoj刷水记录]
(这篇我就不信有网站来扣) 这个暑假打算刷刷题啥的 但是写博客好累啊 堆一起算了 隔一段更新一下. 7月27号之前刷的的就不写了 , 写的累 代码不贴了,可以找我要啊.. 2017.8.27upd ...
- bzoj 乱刷计划 50/50
前言 话说第一个板刷计划由于种种原因而告一段落了..其实那一版还有很多题想做,那就只能放一放了 附上效果图一张(几乎每一题都在我博客有题解): 打算 可以复习,重做自己做过的题,不局限于没做过的 乱刷 ...
- 一句话题解(20170801~20170125)
8.1 bzoj 4720 noip2016 换教室 floyd预处理+期望(薛定谔的猫) bzoj 4318 OSU! 三次函数期望值 从一次.二次推得 8.2 bzoj 1076 状压+期望DP ...
- linux系统数据库类型,linux系统下调度数据库类型资源库中的kettle job
已经存在kettle的一个资源库enfo,在目录/works/wxj下面有一个job (testmailsuccess.kjb)如何实现手工在kettle外部执行此job和让系统每天定时的调用此job ...
最新文章
- print_r php encode,详细介绍PHP在调试时echo print() print_r() var_dump()的区别分享
- bash文件外传入参数
- Winform中通过代码设置DevExpress的TextEdit的类型为Numbernic
- java实现资源监视器_实现Java监视的12个步骤程序存在缺陷
- Linux内核源代码分析——插入之后会发生什么?
- ios之JavaScript
- matlab序列码,Matlab7 R14(sn 序列号)
- sqlserver之创建视图
- php圆的面积计算公式,圆的面积怎么计算如何求解
- mPass 微服务开发平台
- 美丽诗文背诵-未完待续
- 隐马尔科夫模型(HMM)等文章记录
- SMT贴片元件拆卸的妙招分享
- 对 Access 数据库中日期字段操作时 OleDbParameter 的问题及解决
- 神经网络中的偏置值什么意思
- threejs+vue房子模型,实现BSP模型剪切,导入obj+mtl文件
- 用html5实现一个播放器,6 个基于 HTML5 实现的多媒体播放器
- php选择星座,PHP日期得到星座
- 【渝粤题库】陕西师范大学151214 纳税筹划 作业(专升本)
- preparedStatement.setDate()丢失时间问题