题解 洛谷 P2086 【[NOI2012]魔幻棋盘】
先考虑只有一维的情况,要求支持区间加和求区间 \(\gcd\),根据 \(\gcd\) 的性质,发现:
\]
其中 \(a_i\) 为原序列 \(a\) 中的任意一个元素,其与序列 \(a\) 的差分序列的 \(\gcd\) 即为原序列的 \(\gcd\)。根据该性质,对于一维的情况,就可以通过线段树单点修改维护差分序列,区间查询 \(\gcd\) 来实现了。对于 \(a_i\),根据修改维护其当前值即可。
然后考虑二维的情况,比如对于这样的一个 \(4 \times 4\) 的矩形:
a_{1,1}\quad a_{1,2}\quad a_{1,3}\quad a_{1,4} \\
a_{2,1}\quad a_{2,2}\quad a_{2,3}\quad a_{2,4} \\
a_{3,1}\quad a_{3,2}\quad a_{3,3}\quad a_{3,4} \\
a_{4,1}\quad a_{4,2}\quad a_{4,3}\quad a_{4,4}
\end{aligned}
\]
设棋盘守护者的位置为 \((x,y)\),先对每一行进行差分,得:
a_{1,y}\quad a_{1,2}-a_{1,1}\quad a_{1,3}-a_{1,2}\quad a_{1,4}-a_{1,3} \\
a_{2,y}\quad a_{2,2}-a_{2,1}\quad a_{2,3}-a_{2,2}\quad a_{2,4}-a_{2,3} \\
a_{3,y}\quad a_{3,2}-a_{3,1}\quad a_{3,3}-a_{3,2}\quad a_{3,4}-a_{3,3} \\
a_{4,y}\quad a_{4,2}-a_{4,1}\quad a_{4,3}-a_{4,2}\quad a_{4,4}-a_{4,3}
\end{aligned}
\]
对每一行来维护差分序列,根据每一行的 \(\gcd\),就能求得整个矩形的 \(\gcd\)。发现这样实现复杂度依然无法接受,于是对每一列也进行差分,也就是进行二维差分,得:
&a_{x,y}\qquad\qquad\qquad\quad a_{x,2}-a_{x,1}\qquad\qquad\qquad\quad a_{x,3}-a_{x,2}\qquad\qquad\qquad\quad a_{x,4}-a_{x,3} \\\\
&a_{2,y}-a_{1,y}\qquad a_{2,2}-a_{2,1}-a_{1,2}+a_{1,1}\qquad a_{2,3}-a_{2,2}-a_{1,3}+a_{1,2}\qquad a_{2,4}-a_{2,3}-a_{1,4}+a_{1,3} \\\\
&a_{3,y}-a_{2,y}\qquad a_{3,2}-a_{3,1}-a_{2,2}+a_{2,1}\qquad a_{3,3}-a_{3,2}-a_{2,3}+a_{2,2}\qquad a_{3,4}-a_{3,3}-a_{2,4}+a_{2,3} \\\\
&a_{4,y}-a_{3,y}\qquad a_{4,2}-a_{4,1}-a_{3,2}+a_{3,1}\qquad a_{4,3}-a_{4,2}-a_{3,3}+a_{3,2}\qquad a_{4,4}-a_{4,3}-a_{3,4}+a_{3,3}
\end{aligned}
\]
发现二维差分后的这个矩形的 \(\gcd\) 即为原矩形的 \(\gcd\)。然后考虑维护,对于第一行和第一列,去掉位置 \((1,1)\) 的 \(a_{x,y}\) 后剩下的序列,其为一维上的问题,用两个一维线段树维护即可。对于去掉第一行和第一列的剩下的矩形,用二维线段树维护即可。两者都是只用支持单点修改和区间查询。对于 \(a_{x,y}\),根据修改维护其当前值即可。
二维线段树可以用树套树或者四分树来实现,因为树套树的空间不好处理,并且题目中也说明了查询是随机的,所以我这里用四分树来实现二维线段树了。
\(code:\)
#include<bits/stdc++.h>
#define maxn 2000010
#define maxm 32000010
#define a(i,j) a[(i-1)*m+j]
#define b(i,j) b[(i-1)*m+j]
#define ls (cur<<1)
#define rs (cur<<1|1)
#define mid ((l+r)>>1)
#define midx ((u+d)>>1)
#define midy ((l+r)>>1)
using namespace std;
typedef long long ll;
template<typename T> inline void read(T &x)
{
x=0;char c=getchar();bool flag=false;
while(!isdigit(c)){if(c=='-')flag=true;c=getchar();}
while(isdigit(c)){x=(x<<1)+(x<<3)+(c^48);c=getchar();}
if(flag)x=-x;
}
int n,m,sx,sy,q,root,rt=1,tree_cnt;
int uls[maxm],urs[maxm],dls[maxm],drs[maxm];
ll w;
ll a[maxn],b[maxn],val[maxm];
ll gcd(ll a,ll b)
{
return b?gcd(b,a%b):abs(a);
}
void build(int u,int d,int l,int r,int &cur)
{
if(u>d||l>r) return;
if(!cur) cur=++tree_cnt;
if(u==d&&l==r)
{
val[cur]=b(u,l);
return;
}
build(u,midx,l,midy,uls[cur]),build(u,midx,midy+1,r,urs[cur]);
build(midx+1,d,l,midy,dls[cur]),build(midx+1,d,midy+1,r,drs[cur]);
val[cur]=gcd(gcd(val[uls[cur]],val[urs[cur]]),gcd(val[dls[cur]],val[drs[cur]]));
}
void modify(int u,int d,int l,int r,int x,int y,ll v,int cur)
{
if(l==r&&u==d)
{
val[cur]+=v;
return;
}
if(x<=midx)
{
if(y<=midy) modify(u,midx,l,midy,x,y,v,uls[cur]);
else modify(u,midx,mid+1,r,x,y,v,urs[cur]);
}
else
{
if(y<=midy) modify(midx+1,d,l,midy,x,y,v,dls[cur]);
else modify(midx+1,d,midy+1,r,x,y,v,drs[cur]);
}
val[cur]=gcd(gcd(val[uls[cur]],val[urs[cur]]),gcd(val[dls[cur]],val[drs[cur]]));
}
ll query(int U,int D,int L,int R,int u,int d,int l,int r,int cur)
{
if(U>D||L>R) return 0;
if(U<=u&&D>=d&&L<=l&&R>=r) return val[cur];
ll v=0;
if(U<=midx)
{
if(L<=midy) v=gcd(v,query(U,D,L,R,u,midx,l,midy,uls[cur]));
if(R>midy) v=gcd(v,query(U,D,L,R,u,midx,midy+1,r,urs[cur]));
}
if(D>midx)
{
if(L<=midy) v=gcd(v,query(U,D,L,R,midx+1,d,l,midy,dls[cur]));
if(R>midy) v=gcd(v,query(U,D,L,R,midx+1,d,midy+1,r,drs[cur]));
}
return v;
}
struct Segment_Tree
{
ll a[maxn],val[maxn];
void build(int l,int r,int cur)
{
if(l==r)
{
val[cur]=a[l];
return;
}
build(l,mid,ls),build(mid+1,r,rs);
val[cur]=gcd(val[ls],val[rs]);
}
void modify(int l,int r,int pos,ll v,int cur)
{
if(l==r)
{
val[cur]+=v;
return;
}
if(pos<=mid) modify(l,mid,pos,v,ls);
else modify(mid+1,r,pos,v,rs);
val[cur]=gcd(val[ls],val[rs]);
}
ll query(int L,int R,int l,int r,int cur)
{
if(L>R) return 0;
if(L<=l&&R>=r) return val[cur];
ll v=0;
if(L<=mid) v=gcd(v,query(L,R,l,mid,ls));
if(R>mid) v=gcd(v,query(L,R,mid+1,r,rs));
return v;
}
}T1,T2;
int main()
{
read(n),read(m),read(sx),read(sy),read(q);
for(int i=1;i<=n;++i)
for(int j=1;j<=m;++j)
read(a(i,j));
for(int i=1;i<n;++i)
for(int j=1;j<m;++j)
b(i,j)=a(i+1,j+1)-a(i,j+1)-a(i+1,j)+a(i,j);
for(int i=1;i<n;++i) T1.a[i]=a(i+1,sy)-a(i,sy);
for(int i=1;i<m;++i) T2.a[i]=a(sx,i+1)-a(sx,i);
w=a(sx,sy),build(0,n,0,m,root),T1.build(0,n,root),T2.build(0,m,root);
while(q--)
{
int u,d,l,r,opt;
ll v;
read(opt),read(u),read(l),read(d),read(r);
if(!opt)
{
u=sx-u,l=sy-l,d=sx+d,r=sy+r;
v=gcd(w,query(u,d-1,l,r-1,0,n,0,m,root));
v=gcd(v,gcd(T1.query(u,d-1,0,n,rt),T2.query(l,r-1,0,m,rt)));
printf("%lld\n",v);
}
else
{
read(v);
if(sx>=u&&sx<=d&&sy>=l&&sy<=r) w+=v;
if(sy>=l&&sy<=r) T1.modify(0,n,u-1,v,rt),T1.modify(0,n,d,-v,rt);
if(sx>=u&&sx<=d) T2.modify(0,m,l-1,v,rt),T2.modify(0,m,r,-v,rt);
modify(0,n,0,m,u-1,l-1,v,root),modify(0,n,0,m,d,r,v,root);
modify(0,n,0,m,u-1,r,-v,root),modify(0,n,0,m,d,l-1,-v,root);
}
}
return 0;
}
题解 洛谷 P2086 【[NOI2012]魔幻棋盘】的更多相关文章
- BZOJ2877:[NOI2012]魔幻棋盘
浅谈树状数组与主席树:https://www.cnblogs.com/AKMer/p/9946944.html 题目传送门:https://lydsy.com/JudgeOnline/problem. ...
- 题解 洛谷P5018【对称二叉树】(noip2018T4)
\(noip2018\) \(T4\)题解 其实呢,我是觉得这题比\(T3\)水到不知道哪里去了 毕竟我比较菜,不大会\(dp\) 好了开始讲正事 这题其实考察的其实就是选手对D(大)F(法)S(师) ...
- 题解 洛谷 P3396 【哈希冲突】(根号分治)
根号分治 前言 本题是一道讲解根号分治思想的论文题(然鹅我并没有找到论文),正 如论文中所说,根号算法--不仅是分块,根号分治利用的思想和分块像 似却又不同,某一篇洛谷日报中说过,分块算法实质上是一种 ...
- 题解-洛谷P5410 【模板】扩展 KMP(Z 函数)
题面 洛谷P5410 [模板]扩展 KMP(Z 函数) 给定两个字符串 \(a,b\),要求出两个数组:\(b\) 的 \(z\) 函数数组 \(z\).\(b\) 与 \(a\) 的每一个后缀的 L ...
- 题解-洛谷P4229 某位歌姬的故事
题面 洛谷P4229 某位歌姬的故事 \(T\) 组测试数据.有 \(n\) 个音节,每个音节 \(h_i\in[1,A]\),还有 \(m\) 个限制 \((l_i,r_i,g_i)\) 表示 \( ...
- 题解-洛谷P4724 【模板】三维凸包
洛谷P4724 [模板]三维凸包 给出空间中 \(n\) 个点 \(p_i\),求凸包表面积. 数据范围:\(1\le n\le 2000\). 这篇题解因为是世界上最逊的人写的,所以也会有求凸包体积 ...
- 题解-洛谷P4859 已经没有什么好害怕的了
洛谷P4859 已经没有什么好害怕的了 给定 \(n\) 和 \(k\),\(n\) 个糖果能量 \(a_i\) 和 \(n\) 个药片能量 \(b_i\),每个 \(a_i\) 和 \(b_i\) ...
- 题解-洛谷P5217 贫穷
洛谷P5217 贫穷 给定长度为 \(n\) 的初始文本 \(s\),有 \(m\) 个如下操作: \(\texttt{I x c}\),在第 \(x\) 个字母后面插入一个 \(c\). \(\te ...
- 洛谷 P1169||bzoj1057 [ZJOI2007]棋盘制作
洛谷P1169 bzoj1057 这个题目跟最大全0子矩阵是类似的.正方形的话,只要把任意极大子正方形(”极大“定义见后面的”论文“)当成把某个极大子矩形去掉一块变成正方形即可,容易解决. 解法1:看 ...
随机推荐
- Python3-算法-冒泡排序
冒泡排序 它重复地走访过要排序的数列,一次比较两个元素,如果他们的顺序错误就把他们交换过来,走访数列的工作是重复地进行直到没有再需要交换,也就是说该数列已经排序完成,这个算法的名字由来是因为越大的元素 ...
- BOM问题-对于php的影响
甲.BOM说明 BOM(Byte Order Mark),是UTF编码方案里用于标识编码的标准标记.这个标记是可选的,UTF-8不需要BOM来表明字节顺序,但可以用BOM来表明当前编码方式.但如果文件 ...
- Myeclipse启动WebLogic 总是报账号密码无效<Authentication denied: Boot identity not valid
在MyEclipse下配置了Weblogic 11后,每次启动从报错: Critical> 看了下描述,是用户名及密码什么的问题,我想起来,配置Weblogic 的域的时候将密码改成了12345 ...
- Redis系列(八):发布与订阅
Redis的发布与订阅,有点类似于消息队列,发送者往频道发送消息,频道的订阅者接收消息. 1. 发布与订阅示例 首先,在本机开启第1个Redis客户端,执行如下命令订阅blog.redis频道: SU ...
- springMVC中的HttpSession与Model
目录 1.1 spring的@MODELATTRIBUTE 2.1 session的概念 3.1 示例 4.1 为什么springmvc框架要使用model这个对象呢? 突然发问:相信很多人在做WEB ...
- CSS中link和@import的使用区别
我们都知道在html中引入外部的CSS 有2种方式,link标签和@import,他们又什么区别呢? 1.从属关系区别@import是 CSS 提供的语法规则,只有导入样式表的作用:link是HTML ...
- elasticsearch 单节点搭建与爬坑记录
elasticsearch 单节点搭建与爬坑记录 prepare 虚拟机或者云服务器(这里用的是阿里云ECS) linux---centos7 安装完毕的jdk 相应的安装包(在https:/ ...
- 云小课 | “VPC连接”知多少
摘要:华为云提供了丰富的网络服务,可满足多种网络互连场景. 同Region的两个VPC怎么连通?” “跨Region的两个VPC又怎么连通?” “VPC内的ECS搭建了一个应用,需要访问Interne ...
- 关于小白一天摆弄Axure产品原型的心得
因为项目所需,被迫快速设计产品原型,制作app使用视频,这里简单分享一些一个小白第一次触碰产品原型设计的一些心得: 原型设计工具的选用 我这里针对的是贫穷学生党,分享的是自己尝试的几款.有钱大佬自行绕 ...
- sql 增加修改表格属性字段
1.修改字段默认值 alter table 表名 drop constraint 约束名字 ------注解:删除表的字段的原有约束 alter table 表名 add constraint 约束名 ...