NOI 2012 魔幻棋盘 | 二维差分 + 二维线段树

题目：luogu 2086

二维线段树，按套路差分原矩阵，gcd( x₁, x₂, ……, x_n ) = gcd( x_i, x₂ - x₁, ……, x_n - x_n-1 )，必须要有一个原数 x_i，恰好每次询问都包含一个固定点 ( X , Y )，差分以它为中心就可以保证它是原值。以 e 为中心的二维差分如图。

对于一维序列，修改区间 [ l , r ] 只需修改差分后的 l 和 r + 1 两点。那么对于二维，差分后的修改如下所示：

中间的灰色格子为守卫者所在地（为方便表示多个区域把它拆成了 9 格），以它为差分中心，若修改矩阵在它的右下方，即黄色区域，设它的左上角坐标为 ( x1, y1)，右下角坐标 ( x2, y2 )，则需在 ( x1, y1 ) 处 + c，( x1, y2+1 ) 处 - c，( x2+1, y1 ) 处 - c，( x2 + 1, y2 + 1) 处 +c。

其他区域同理。再加上守卫者正上方、正下方、正左方、正右方的区域，一共 8 种情况。如果一个修改矩阵包含多个区域，就将其划分为若干个子矩阵处理。

这样就将巨慢的区间修改转换成了单点修改，每修改完一个点，pushup( ) 维护 gcd 即可。

另外读入的时候，500000 * 500000 的二维数组是开不了的，需要动态分配内存；或者开一个 500000 的一维数组读入，判断每一行的末端在哪里。

讨论的时候有些麻烦，下面列一下大致的情况：（黄色表示修改的区域，绿色加，红色减）

仔细观察后，发现可分为三种情况：① 包含守卫者 ( X , Y )；② 包含 X 行或 Y 列；③ 只包含左上或左下、右上或右下区域。

细节诸多：

差分的时候不能用原数组减，因为前面的值会变，需要另开一个数组保存；

返回 gcd 的时候要取绝对值；二维线段树内层空间要开大。

 #include <cstdio>

 #include <string>

 const int N = ;

 typedef long long ll;

 int n, m; ll a[N], b[N];

 ll read() {

     ll x = , f = ;

     char c = getchar();

     while (!isdigit(c)) {

         if (c == '-') f = -;

         c = getchar();

     }

     while (isdigit(c)) {

         x = (x << ) + (x << ) + (c ^ );

         c = getchar();

     }

     return x * f;

 }

 ll abs(ll x) {

     if (x >= ) return x; return -x;

 }

 ll gcd(ll x, ll y) {

     if (y) return gcd(y, x % y); return abs(x);

 }

 struct inNode {

     ll val;

     inNode *ls, *rs;

     inNode() : val(), ls(NULL), rs(NULL) {}

 } inPool[N << ];

 inNode *newInNode () {

     static int cnt = ;

     return &inPool[++cnt];

 }

 void buildY(inNode *&cur, int l, int r, int x) {

     if (!cur) cur = newInNode();

     if (l == r) cur->val = a[(x - ) * m + l];

     else {

         int mid = l + ((r - l) >> );

         buildY(cur->ls, l, mid, x);

         buildY(cur->rs, mid + , r, x);

         cur->val = gcd(cur->ls->val, cur->rs->val);

     }

 }

 void updateY(inNode *&cur, int l, int r, int y, ll val) {

     if (l == r) cur->val += val;

     else {

         int mid = l + ((r - l) >> );

         if (y <= mid) updateY(cur->ls, l, mid, y, val);

         else updateY(cur->rs, mid + , r, y, val);

         cur->val = gcd(cur->ls ? cur->ls->val : , cur->rs ? cur->rs->val : );

     }

 }

 ll queryY(inNode *&cur, int l, int r, int y1, int y2) {

     if (y1 <= l && r <= y2) return cur->val;

     int mid = l + ((r - l) >> ); ll res = ;

     if (y1 <= mid) res = gcd(res, queryY(cur->ls, l, mid, y1, y2));

     if (mid < y2) res = gcd(res, queryY(cur->rs, mid + , r, y1, y2));

     return res;

 }

 struct outNode {

     inNode *root;

     outNode *ls, *rs;

     outNode() : root(NULL), ls(NULL), rs(NULL) {}

 } outPool[N << ], *root;

 outNode *newOutNode() {

     static int cnt = ;

     return &outPool[++cnt];

 }

 void maintain(inNode *&cur, inNode *&lc, inNode *&rc, int l, int r) {

     if (!cur) cur = newInNode();

     if (l == r) cur->val = gcd(lc->val, rc->val);

     else {

         int mid = l + ((r - l) >> );

         maintain(cur->ls, lc->ls, rc->ls, l, mid);

         maintain(cur->rs, lc->rs, rc->rs, mid + , r);

         cur->val = gcd(cur->ls->val ,cur->rs->val);

     }

 }

 void buildX(outNode *&cur, int l, int r) {

     if (!cur) cur = newOutNode();

     if (l == r) buildY(cur->root, , m, l);

     else {

         int mid = l + ((r - l) >> );

         buildX(cur->ls, l, mid);

         buildX(cur->rs, mid + , r);

         maintain(cur->root, cur->ls->root, cur->rs->root, , m);

     }

 }

 void change(inNode *&cur, int l, int r, int y, ll val) {

     if (l == r) cur->val = val;

     else {

         int mid = l + ((r - l) >> );

         if (y <= mid) change(cur->ls, l, mid, y, val);

         else change(cur->rs, mid + , r, y, val);

         cur->val = gcd(cur->ls ? cur->ls->val : , cur->rs ? cur->rs->val : );

     }

 }

 void updateX(outNode *&cur, int l, int r, int x, int y, ll val) {

     if (x <=  || y <=  || x > n || y > m) return;

     if (l == r) updateY(cur->root, , m, y, val);

     else {

         int mid = l + ((r - l) >> ); ll lv = , rv = ;

         if (x <= mid) updateX(cur->ls, l, mid, x, y, val);

         else updateX(cur->rs, mid + , r, x, y, val);

         if (cur->ls) lv = queryY(cur->ls->root, , m, y, y);

         if (cur->rs) rv = queryY(cur->rs->root, , m, y, y);

         change(cur->root, , m, y, gcd(lv, rv));

     }

 }

 ll queryX(outNode *&cur, int l, int r, int x1, int x2, int y1, int y2) {

     if (x1 <= l && r <= x2) return queryY(cur->root, , m, y1, y2);

     int mid = l + ((r - l) >> ); ll res = ;

     if (x1 <= mid) res = gcd(res, queryX(cur->ls, l, mid, x1, x2, y1, y2));

     if (mid < x2) res = gcd(res, queryX(cur->rs, mid + , r, x1, x2, y1, y2));

     return res;

 }

 int main() {

     n = read(), m = read();    int X = read(), Y = read(), T = read();

     for (int i = ; i <= n * m; ++ i) a[i] = read();

     for (int i = ; i <= n * m; ++ i) {

         if ((i - ) % m +  < Y) b[i] = a[i] - a[i + ];

         else if ((i - ) % m +  > Y) b[i] = a[i] - a[i - ];

         else b[i] = a[i];

     }

     for (int i = ; i <= n * m; ++ i) {

         if ((i - ) / m +  < X) a[i] = b[i] - b[i + m];

         else if ((i - ) / m +  > X) a[i] = b[i] - b[i - m];

         else a[i] = b[i];

     }

     buildX(root, , n);

     while (T--) {

         int opt = read(), x1 = read(), y1 = read(), x2 = read(), y2 = read();

         if (opt == ) {

             x1 = X - x1, x2 = X + x2, y1 = Y - y1, y2 = Y + y2;

             printf("%lld\n", queryX(root, , n, x1, x2, y1, y2));

         } else {

             ll val = read();

             if (x1 <= X && x2 >= X && y1 <= Y && y2 >= Y) {    //包含(X,Y)

                 updateX(root, , n, X, Y, val);

                 updateX(root, , n, x1 - , y1 - , val);

                 updateX(root, , n, x1 - , y2 + , val);

                 updateX(root, , n, x2 + , y1 - , val);

                 updateX(root, , n, x2 + , y2 + , val);

                 updateX(root, , n, x1 - , Y, -val);

                 updateX(root, , n, x2 + , Y, -val);

                 updateX(root, , n, X, y1 - , -val);

                 updateX(root, , n, X, y2 + , -val);

             } else if (y1 <= Y && y2 >= Y) {                //包含Y列

                 if (x1 < X) {                                //在(X,Y)上方

                     updateX(root, , n, x2, Y, val);

                     updateX(root, , n, x1 - , y1 - , val);

                     updateX(root, , n, x1 - , y2 + , val);

                     updateX(root, , n, x1 - , Y, -val);

                     updateX(root, , n, x2, y1 - , -val);

                     updateX(root, , n, x2, y2 + , -val);

                 } else {                                    //在(X,Y)下方

                     updateX(root, , n, x1, Y, val);

                     updateX(root, , n, x2 + , y1 - , val);

                     updateX(root, , n, x2 + , y2 + , val);

                     updateX(root, , n, x2 + , Y, -val);

                     updateX(root, , n, x1, y1 - , -val);

                     updateX(root, , n, x1, y2 + , -val);

                 }

             } else if (x1 <= X && x2 >= X) {                //包含X行

                 if (y1 < Y) {                                //在(X,Y)左边

                     updateX(root, , n, X, y2, val);

                     updateX(root, , n, x1 - , y1 - , val);

                     updateX(root, , n, x2 + , y1 - , val);

                     updateX(root, , n, X, y1 - , -val);

                     updateX(root, , n, x1 - , y2, -val);

                     updateX(root, , n, x2 + , y2, -val);

                 } else {                                    //在(X,Y)右边

                     updateX(root, , n, X, y1, val);

                     updateX(root, , n, x1 - , y2 + , val);

                     updateX(root, , n, x2 + , y2 + , val);

                     updateX(root, , n, X, y2 + , -val);

                     updateX(root, , n, x1 - , y1, -val);

                     updateX(root, , n, x2 + , y1, -val);

                 }

             } else if (x1 < X && y1 < Y) {                    //左上区域

                 updateX(root, , n, x2, y2, val);

                 updateX(root, , n, x1 - , y1 - , val);

                 updateX(root, , n, x1 - , y2, -val);

                 updateX(root, , n, x2, y1 - , -val);

             } else if (x1 < X && y1 > Y) {                    //右上区域

                 updateX(root, , n, x2, y1, val);

                 updateX(root, , n, x1 - , y2 + , val);

                 updateX(root, , n, x1 - , y1, -val);

                 updateX(root, , n, x2, y2 + , -val);

             } else if (x1 > X && y1 < Y) {                    //左下区域

                 updateX(root, , n, x1, y2, val);

                 updateX(root, , n, x2 + , y1 - , val);

                 updateX(root, , n, x1, y1 - , -val);

                 updateX(root, , n, x2 + , y2, -val);

             } else {                                        //右下区域

                 updateX(root, , n, x1, y1, val);

                 updateX(root, , n, x2 + , y2 + , val);

                 updateX(root, , n, x1, y2 + , -val);

                 updateX(root, , n, x2 + , y1, -val);

             }

         }

     }

     return ;

 }

NOI 2012 魔幻棋盘 | 二维差分 + 二维线段树的更多相关文章

【坐标变换】【二维偏序】【线段树】Gym - 100820G - Racing Gems
题意:第一象限有n个点,你从x正半轴任选一个位置出发,vy恒定,vx可以任意变化,不过只能在-vy/r到vy/r之间变化,问你最多能经过多少个点. 暴力dp是n^2,不可取. 注意到,一个点,所能到达 ...
[BZOJ3196] 二逼平衡树 [权值线段树套位置平衡树]
题面洛咕题面思路没错我就是要不走寻常路! 看看那些外层位置数据结构,必须二分的,$O(n\log^3 n)$的做法吧! 看看那些cdq分治/树状数组套线段树的,空间$O(n\log^2 n)$挤 ...
二叉搜索树TREE（线段树，区间DP）
前言线段树+区间DP题,线段树却不是优化DP的,是不是很意外? 题面二叉搜索树是一种二叉树,每个节点都有一个权值,并且一个点的权值比其左子树里的点权值都大,比起右子树里的点权值都小. 一种朴素的向 ...
[BZOJ3307]:雨天的尾巴（LCA+树上差分+权值线段树）
题目传送门题目描述: N个点,形成一个树状结构.有M次发放,每次选择两个点x,y对于x到y的路径上(含x,y)每个点发一袋Z类型的物品.完成所有发放后,每个点存放最多的是哪种物品. 输入格式: 第一 ...
洛谷1083（差分+二分 or 线段树）
第一种方法:可以二分最大天数订单的答案然后通过差分求一下是否可行. ; int n, m, a[maxn], ans; struct section { int cnt, l, r; }b[maxn] ...
2019杭电多校6 hdu6638 Snowy Smile（二维最大矩阵和线段树）
http://acm.hdu.edu.cn/showproblem.php?pid=6638 题意:给你一些点的权值,让找一个矩形圈住一部分点,问圈住点的最大权值和分析:由于是稀疏图,明显要先把x, ...
刷题总结——二逼平衡树（bzoj3224线段树套splay）
题目: Description 您需要写一种数据结构(可参考题目标题),来维护一个有序数列,其中需要提供以下操作:1.查询k在区间内的排名2.查询区间内排名为k的值3.修改某一位值上的数值4.查询k在 ...
2018HDU多校二 -F 题 Naive Operations(线段树)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6315 In a galaxy far, far away, there are two integer ...
NOIp 2012 #2 借教室 Label：区间修改线段树
题目描述在大学期间,经常需要租借教室.大到院系举办活动,小到学习小组自习讨论,都需要向学校申请借教室.教室的大小功能不同,借教室人的身份不同,借教室的手续也不一样. 面对海量租借教室的信息,我们自然 ...

随机推荐

python 读取共享内存
测试环境 centos7 python3.6.5 首先使用c创建内存,这里的方法是:作为参数读一个二进制数据文件进去,把文件的内容作为共享内存的内容定义块 #include <stdio.h& ...
Testing Is the Engineering Rigor of Software Development
Testing Is the Engineering Rigor of Software Development Neal Ford DEVELOPERS LOVE TO USE TORTURED M ...
sql的一些知识_where
简单相同查询,查询age=22的三列信息,并按照降序排列 ORDER BY weight DESC 此外,where还支持一下判断操作符值得注意的是,如果比较的值为数值,则不需要加引号,如果是字符, ...
树莓派 Zero作为飞控图传
前言原创文章,转载引用务必注明链接,水平有限,如有疏漏,欢迎指正. 本文使用Markdown写成,为获得更好的阅读体验和正常的链接.图片显示,请访问我的博客原文: http://www.cnblog ...
hql 多对多查询
这种查询,hibernate 建议用 From Dealer s inner join fetch s.carSerieses cs 实现,注意这种实现只支持b.c,不支持b.cs. 如果要用b.c ...
Android 音频的播放之二MediaPlayer
MediaPlayer类可用于控制音频/视频文件或流的播放.关于怎样使用这个类的方法还能够阅读VideoView类的文档. 1．状态图对播放音频/视频文件和流的控制是通过一个状态机来管理的. 下图显 ...
shell脚本编写-自动部署及监控
1.编写脚本自动部署反向代理.web.nfs: I.部署nginx反向代理两个web服务,调度算法使用加权轮询 II.所有web服务使用共享存储nfs,保证所有web都对其有读写权限,保证数据一致性: ...
gulp（基础篇）
今天在写项目的时候用到了gulp构建工具,虽然一年前就有用过,但是一直只存在于我的“有道云笔记”里,今天又一次用到,固然是巩固一下,这里来记录一下吧:这里我主要想要记录的就是初学者在第一次使用gulp ...
用DD-WRT自建计费WiFi热点
架设无线网络(Wlan)向周围的用户有偿共享网络. 传统的方法,能够使用专业的无线与宽带计费网关设备和软件.比方MikroTik等. 只是,对于个人架设WiFi热点来说.这些方案太过昂贵,安装设置也很 ...
C语言--函数篇
1-1.函数简单调用 1 #include<stdio.h> 2 #include<string.h> 3 #include<windows.h> 4 int ...

NOI 2012 魔幻棋盘 | 二维差分 + 二维线段树

NOI 2012 魔幻棋盘 | 二维差分 + 二维线段树的更多相关文章

随机推荐

热门专题