MAZE(2019年牛客多校第二场E题+线段树+矩阵乘法)
题目链接
题意
在一张\(n\times m\)的矩阵里面,你每次可以往左右和下三个方向移动(不能回到上一次所在的格子),\(1\)表示这个位置是墙,\(0\)为空地。
现在有\(q\)次操作,操作一是将\((x,y)\)这个位置的状态取反,操作二问你从\((1,x)\)走到\((n,y)\)的方案数。
思路
首先我们考虑不带修改操作时求方案数:
我们发现从第\(i-1\)行到第\(i\)行的\(j\)这个位置只能通过\((i-1,j)\)到达,因此可以从第\(i-1\)行到\((i,j)\)的位置只能是与\((i-1,j)\)的路径上不能有墙的点,从而我们可以得知\(dp[x][y]=\sum\limits_{i=L}^{R}dp[x-1][i]\),其中\(dp[x][y]\)表示到达\((x,y)\)这个位置的方案数,\(L,R\)表示与\((x-1,y)\)联通的左右端点。
我们发现这是一个递推式,因此我们可以用矩阵乘法来维护这个东西,我们用一个例子来理解:
假设要从第\(i-1\)行到达第\(i\)行,且第\(i-1\)行的状态为\("10010"\),那么将递推式表示成矩阵乘法为:
\begin{matrix}
dp[i][1] & dp[i][2] & dp[i][3] & dp[i][4] & dp[i][5]
\end{matrix}
\right)=
\left(
\begin{matrix}
dp[i-1][1] & dp[i-1][2] & dp[i-1][3] & dp[i-1][4] & dp[i-1][5]
\end{matrix}
\right)
\times
\left(
\begin{matrix}
0 & 0 & 0 & 0 & 1\\
0 & 1 & 1 & 0 & 1\\
0 & 1 & 1 & 0 & 1\\
0 & 0 & 0 & 0 & 1\\
0 & 0 & 0 & 0 & 1
\end{matrix}
\right)
\]
得到了相邻两项的递推式那么从第\(1\)行到第\(n\)行的答案那么答案就是\(dp[n+1][y]\),为什么是\(n+1\)而不是\(n\)呢?因为如果是\(n\),那么得到的只有从\(n-1\)行到达这个位置的方案数,缺少了从第\(n\)行的其他位置到达这个位置的方案数。
那么待修改操作的我们该怎么处理呢?
我们发现修改一个点的位置只会影响当前行与下一行的系数矩阵,并不会影响其他行直接的系数矩阵,如果我们暴力修改然后暴力求解递推式那么对于每次操作都要从第\(1\)行递推到第\(n\)行,那么每次修改操作复杂度为\(O(1)\),求解复杂度为\(O(nm^3)\),很明显不能满足题目给的时限。
我们发现如果我们用线段树来维护这个东西,那么每次修改的复杂度为\(O(m^3long(n))\),求解复杂度为\(O(1)\),那么总体复杂度就比上面上了一个\(n\)。
而维护方式也很简单,定义线段树的每个结点都是一个系数矩阵,表示从结点代表的左端点\((l)\)到右端点\((r)\)\(+1\)的递推式中的矩阵相乘,而且系数矩阵中的\(sum[i][j]\)表示的是从\(l\)的第\(i\)列到\(r+1\)的第\(j\)列的方案数,最后答案为根结点的\(sum[x][y]\)。
代码
#include <set>
#include <map>
#include <deque>
#include <queue>
#include <stack>
#include <cmath>
#include <ctime>
#include <bitset>
#include <cstdio>
#include <string>
#include <vector>
#include <cassert>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;
typedef pair<LL, LL> pLL;
typedef pair<LL, int> pLi;
typedef pair<int, LL> pil;;
typedef pair<int, int> pii;
typedef unsigned long long uLL;
#define lson (rt<<1)
#define rson (rt<<1|1)
#define lowbit(x) x&(-x)
#define name2str(name) (#name)
#define bug printf("*********\n")
#define debug(x) cout<<#x"=["<<x<<"]" <<endl
#define FIN freopen("/home/dillonh/CLionProjects/Dillonh/in.txt","r",stdin)
#define IO ios::sync_with_stdio(false),cin.tie(0)
const double eps = 1e-8;
const int mod = 1000000007;
const int maxn = 50000 + 7;
const double pi = acos(-1);
const int inf = 0x3f3f3f3f;
const LL INF = 0x3f3f3f3f3f3f3f3fLL;
int n, m, q, op, x, y;
int mp[maxn][11];
struct node {
int l, r, sum[11][11];
node operator * (const node& a) const {
node c;
for(int i = 1; i <= m; ++i) {
for(int j = 1; j <= m; ++j) {
c.sum[i][j] = 0;
for(int k = 1; k <= m; ++k) {
c.sum[i][j] = (c.sum[i][j] + 1LL * sum[i][k] * a.sum[k][j] % mod) % mod;
}
}
}
return c;
}
}segtree[maxn<<2];
void push_up(int rt) {
int l = segtree[rt].l, r = segtree[rt].r;
segtree[rt] = segtree[lson] * segtree[rson];
segtree[rt].l = l, segtree[rt].r = r;
}
void build(int rt, int l, int r) {
segtree[rt].l = l, segtree[rt].r = r;
if(l == r) {
for(int i = 1; i <= m; ++i) {
for(int j = 1; j <= m; ++j) segtree[rt].sum[i][j] = 0;
for(int j = i; j <= m && !mp[l][j]; ++j) segtree[rt].sum[i][j] = 1;
for(int j = i; j >= 1 && !mp[l][j]; --j) segtree[rt].sum[i][j] = 1;
}
return;
}
int mid = (l + r) >> 1;
build(lson, l, mid);
build(rson, mid + 1, r);
push_up(rt);
}
void update(int rt, int pos) {
if(segtree[rt].l == segtree[rt].r) {
for(int i = 1; i <= m; ++i) {
for(int j = 1; j <= m; ++j) segtree[rt].sum[i][j] = 0;
for(int j = i; j <= m && !mp[pos][j]; ++j) segtree[rt].sum[i][j] = 1;
for(int j = i; j >= 1 && !mp[pos][j]; --j) segtree[rt].sum[i][j] = 1;
}
return;
}
int mid = (segtree[rt].l + segtree[rt].r) >> 1;
if(pos <= mid) update(lson, pos);
else update(rson, pos);
push_up(rt);
}
int main() {
#ifndef ONLINE_JUDGE
FIN;
#endif
scanf("%d%d%d", &n, &m, &q);
for(int i = 1; i <= n; ++i) {
for(int j = 1; j <= m; ++j) {
scanf("%1d", &mp[i][j]);
}
}
build(1, 1, n);
while(q--) {
scanf("%d%d%d", &op, &x, &y);
if(op == 1) {
mp[x][y] ^= 1;
update(1, x);
} else {
printf("%d\n", segtree[1].sum[x][y]);
}
}
return 0;
}
MAZE(2019年牛客多校第二场E题+线段树+矩阵乘法)的更多相关文章
- 2019年牛客多校第二场 H题Second Large Rectangle
题目链接 传送门 题意 求在\(n\times m\)的\(01\)子矩阵中找出面积第二大的内部全是\(1\)的子矩阵的面积大小. 思路 处理出每个位置往左连续有多少个\(1\),然后对每一列跑单调栈 ...
- Kth Minimum Clique(2019年牛客多校第二场D题+k小团+bitset)
目录 题目链接 题意 思路 代码 题目链接 传送门 题意 找第\(k\)小团. 思路 用\(bitset\)来标记每个结点与哪些结点直接有边,然后进行\(bfs\),在判断新加入的点与现在有的点是否都 ...
- 2019年牛客多校第二场 F题Partition problem 爆搜
题目链接 传送门 题意 总共有\(2n\)个人,任意两个人之间会有一个竞争值\(w_{ij}\),现在要你将其平分成两堆,使得\(\sum\limits_{i=1,i\in\mathbb{A}}^{n ...
- 2019年牛客多校第一场B题Integration 数学
2019年牛客多校第一场B题 Integration 题意 给出一个公式,求值 思路 明显的化简公式题,公式是分母连乘形式,这个时候要想到拆分,那如何拆分母呢,自然是裂项,此时有很多项裂项,我们不妨从 ...
- 2019牛客多校第二场H题(悬线法)
把以前的题补补,用悬线求面积第二大的子矩形.我们先求出最大子矩阵的面积,并记录其行三个方向上的悬线长度.然后排除这个矩形,记得还得特判少一行或者少一列的情况 #include <bits/std ...
- Cutting Bamboos(2019年牛客多校第九场H题+二分+主席树)
题目链接 传送门 题意 有\(n\)棵竹子,然后有\(q\)次操作,每次操作给你\(l,r,x,y\),表示对\([l,r]\)区间的竹子砍\(y\)次,每次砍伐的长度和相等(自己定砍伐的高度\(le ...
- 2019年牛客多校第一场 I题Points Division 线段树+DP
题目链接 传送门 题意 给你\(n\)个点,每个点的坐标为\((x_i,y_i)\),有两个权值\(a_i,b_i\). 现在要你将它分成\(\mathbb{A},\mathbb{B}\)两部分,使得 ...
- 2019年牛客多校第一场 H题XOR 线性基
题目链接 传送门 题意 求\(n\)个数中子集内所有数异或为\(0\)的子集大小之和. 思路 对于子集大小我们不好维护,因此我们可以转换思路变成求每个数的贡献. 首先我们将所有数的线性基的基底\(b\ ...
- 2018牛客多校第二场a题
一个人可以走一步或者跳x步,但不能连着跳,问到这个区间里有几种走法 考虑两种状态 对于这一点,我可以走过来,前面是怎么样的我不用管,也可以跳过来但是,跳过来必须保证前一步是走的 dp[i][0]表示 ...
随机推荐
- NOI2016优秀的拆分
一种想法是枚举分割位置, 然后考虑前面部分有多少种可行的AA拆分方式, 后面部分有多少种可行的BB拆分方式, 然后乘法原理即可 那么问题是如何快速求出合法方案 解法是首先枚举长度len, 然后将序列分 ...
- [LeetCode] 727. Minimum Window Subsequence 最小窗口序列
Given strings S and T, find the minimum (contiguous) substring W of S, so that T is a subsequence of ...
- Ubuntu安装微信、钉钉等各种windows软件
详见这个博客,用sudo dpkg -i 安装软件时,如果出现错误,是因为缺少安装依赖关系,用下面的命令解决: sudo apt-get install -f
- Supervisor-守护进程工具
前言 手头上的服务器程序1.0版本终于要进入线上测试阶段了
- Python连载31-threading的使用
一. 例子:我们对传参是有要求的必须传入一个元组,否则报错 二. import _thread as thread import time def loop1(in1): print("St ...
- linux vsftp查看ftp账号信息的方法
linux vsftp查看ftp账号信息的方法 查看注册的FTP账号 在/etc/vsftpd/chroot_list 密码看不到 只能重置密码 passwd username
- 深入学习OpenCV中图像灰度化原理,图像相似度的算法
最近一段时间学习并做的都是对图像进行处理,其实自己也是新手,各种尝试,所以我这个门外汉想总结一下自己学习的东西,图像处理的流程.但是动起笔来想总结,一下却不知道自己要写什么,那就把自己做过的相似图片搜 ...
- window下使用curl操作elasticsearch
1.下载curlzip,https://curl.haxx.se/download.html; 2.解压,在bin文件夹中找到curl.exe,右键“以管理员身份运行”,cmd e: 换盘符:出现E: ...
- Linq分批次,每组1000条
/// <summary> /// 分组插入每次插入1000 /// </summary> /// <param name="data">< ...
- webapi 集成NLog
参考项目代码:SwaggerDemoApi 安装 打开nuget管理器--->搜索nlog,安装箭头所指的两个文件到你的项目中,config安装到你的API项目即可,nlog文件安装到你用得到n ...