2877: [Noi2012]魔幻棋盘

最新推荐文章于 2019-08-30 21:20:00 发布
最新推荐文章于 2019-08-30 21:20:00 发布
阅读量642 收藏
点赞数
分类专栏： 线段树 数学题
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/CRZbulabula/article/details/51013473
版权

2877: [Noi2012]魔幻棋盘

Time Limit: 50 Sec   Memory Limit: 512 MB
Submit: 808   Solved: 317
[ Submit][ Status][ Discuss]

Description

Input

第一行为两个正整数N,M，表示棋盘的大小。 第二行为两个正整数X,Y，表示棋盘守护者的位置。 第三行仅有一个正整数T，表示棋盘守护者将进行次操作。 接下来N行，每行有M个正整数，用来描述初始时棋盘上每个位置的数。 接下来T行，按操作的时间顺序给出T次操作。每行描述一次操作，以一个数字0或1开头： 若以数字0开头，表示此操作为询问，随后会有四个非负整数x1,y1,x2,y2，表示询问的区域是以棋盘守护者的位置为基础向上扩展x1行，向下扩展y1行，向左扩展x2列，向右扩展y2列得到的矩形区域（详见样例）。 若以数字1开头，表示此操作为修改，随后会有四个正整数x1,y1,x2,y2和一个整数c，表示修改区域的上、下边界分别为第x1,x2行，左、右边界分别为第y1,y2列（详见样例），在此矩形区域内的所有数统一加上c（注意c可能为负数）。

Output


 对于每次询问操作，每行输出一个数，表示该区域内所有数的最大公约数。

Sample Input

2 2
1 1
4
6 12
18 24
0 0 0 1 0
1 1 1 1 2 6
1 2 1 2 2 6
0 0 0 1 1

Sample Output

6 6

HINT

Source

[ Submit][ Status][ Discuss]

还是去看noi2012原题吧  wuliwala传送门

大概就是给你一个二维棋盘然后在里面修改区间查询区间gcd，查询操作有一个固定中心
gcd(a,b) = gcd(a,a-b)--更相减损术
对于一维的操作，由上式得gcd(a1,a2,a3,...,an) = gcd(a1,a2-a1,a3-a2,...,an-an-1)
所以。。原序列gcd=原序列的差分序列的gcd
线段树维护区间gcd,修改操作只改2个点

对于二维，二维差分定义为a[i][j] = a[i][j]-a[i-1][j]-a[i][j-1]+a[i-1][j-1]
当然朝四个方向都可以
利用题目给的中心，二维线段树维护四个方向的差分，然后x行与y列特别维护指向中心的一维差分
每次修改最多8个点bulabulabulabula(写的时候N多地方出错)

查询返回值应先定义为0而不是1
涉及gcd操作应使用绝对值
需要修改的单点被多次修改！！！(额外附加判断数组)

万恶的编写错误--Modify函数调用que函数时节点o应为当前节点的左右子节点，Modify函数赋值出错，Judge函数竟然写了两个dx
而且而且。。。getLL居然没考虑负数？！！！！！ 



#include<iostream>      
#include<cstdio>      
#include<cstring>      
#include<cmath>    
#include<vector>      
using namespace std;      
         
const int maxn = 5E5 + 10;     
const int dx[4][4] = {{0,0,1,1},{0,0,1,1},{0,0,-1,-1},{0,0,-1,-1}};    
const int dy[4][4] = {{0,1,0,1},{0,-1,0,-1},{0,1,0,1},{0,-1,0,-1}};     
typedef long long LL;      
         
LL c[maxn*25];      
int n,m,t,x,y,X1,X2,Y1,Y2,cnt = 0,tt,lc[maxn*20],rc[maxn*20];   
int jx[10],jy[10],cur;     
         
vector <LL> v[maxn];    
vector <int> bo[maxn];
         
LL gcd(LL a,LL b) {return b?gcd(b,a%b):a;}      
     
LL ABS(LL x) {return x>0?x:x*-1;}  
     
void PRE(int l,int r,int o)      
{      
    cnt = max(cnt,o);      
    if (l == r) return;      
    int mid = (l+r) >> 1;      
    PRE(l,mid,2*o); PRE(mid+1,r,2*o+1);      
}  
         
LL que(int o,int l,int r,int ql,int qr)    
{    
    if (ql <= l && r <= qr) return ABS(c[o]);    
    LL ret = 0;    
    int mid = (l+r) >> 1;    
    if (ql <= mid) ret = gcd(ret,que(lc[o],l,mid,ql,qr));    
    if (qr > mid) ret = gcd(ret,que(rc[o],mid+1,r,ql,qr));    
    return ret;    
}  
       
LL Modi(int l,int r,int o,int pos,LL va,int typ)      
{      
    if (l == r) {if (typ) c[o] = va; else c[o] += va; return ABS(c[o]);}      
    int mid = (l+r) >> 1; LL ret;  
    if (pos <= mid) {      
        if (!lc[o]) lc[o] = ++cnt;      
        ret = Modi(l,mid,lc[o],pos,va,typ);      
    }      
    else {      
        if (!rc[o]) rc[o] = ++cnt;      
        ret = Modi(mid+1,r,rc[o],pos,va,typ);      
    }      
    c[o] = gcd(ABS(c[lc[o]]),ABS(c[rc[o]]));  
    return ret;  
}      
         
LL Modify(int l,int r,int o,int pn,int pm,LL va)      
{      
    if (!pn || !pm || pn > n || pm > m) return -1;    
    if (l == r) return Modi(1,m,o,pm,va,0);      
    int mid = (l+r) >> 1;      
    LL ret1,ret2;  
    if (pn <= mid) {  
        ret1 = Modify(l,mid,2*o,pn,pm,va);  
        ret2 = que(2*o+1,1,m,pm,pm);   
    }     
    else {  
        ret1 = Modify(mid+1,r,2*o+1,pn,pm,va);  
        ret2 = que(2*o,1,m,pm,pm);      
    }  
    return Modi(1,m,o,pm,gcd(ret1,ret2),1);  
}      
         
LL getLL()    
{    
    LL ret = 0;  LL xx = 1;  
    char ch = getchar();    
    while (ch < '0' || ch > '9') {  
        if (ch == '-') xx = -1;  
        ch = getchar();   
    }  
    while ('0' <= ch && ch <= '9') ret = ret*10 + 1LL*(ch-'0'),ch = getchar();    
    return ret*xx;    
}    
         
void Judge(int px,int py,LL c)    
{    
    if (!px || !py || px > n || py > m || px == x || py == y) return;   
	if (bo[px][py] == tt) return; 
    int po;    
    if (px < x && py < y) po = 0;    
    else if (px < x && py > y) po = 1;    
    else if (px > x && py < y) po = 2;    
    else po = 3;    
    LL tot = 0;    
    for (int l = 0; l < 4; l++) {    
        int xx = px + dx[po][l];    
        int yy = py + dy[po][l];    
        LL t;    
        if ((!dx[po][l] || !dy[po][l]) && dx[po][l] != dy[po][l]) t = -1;  else t = 1;    
        if (X1 <= xx && xx <= X2 && Y1 <= yy && yy <= Y2) tot += c*t;    
    }    
    if (tot != 0) Modify(1,n,1,px,py,tot),bo[px][py] = tt,jx[cur] = px,jy[cur++] = py;    
}     
         
LL query(int o,int l,int r,int nl,int nr,int ml,int mr)    
{    
    if (nl <= l && r <= nr) return que(o,1,m,ml,mr);    
    LL ret = 0;    
    int mid = (l+r) >> 1;    
    if (nl <= mid) ret = gcd(ret,query(2*o,l,mid,nl,nr,ml,mr));    
    if (nr > mid) ret = gcd(ret,query(2*o+1,mid+1,r,nl,nr,ml,mr));    
    return ret;    
}      
       
int main()      
{      
    #ifdef YZY      
        //freopen("chess1.in","r",stdin);      
        //freopen("test.txt","w",stdout);  
    #endif      
             
    cin >> n >> m >> x >> y >> t; PRE(1,n,1);      
    for (int i = 1; i <= n; i++) {    
        v[i].push_back(0);    bo[i].push_back(0);
        for (int j = 1; j <= m; j++) {    
            LL x; x = getLL();    
            v[i].push_back(x);    bo[i].push_back(0);
        }    
    }    
           
    for (int i = 1; i <= n; i++)    
        for (int j = 1; j <= m; j++) {    
            if (i == x && j == y) Modify(1,n,1,i,j,v[i][j]);    
            else if (i == x && j < y) Modify(1,n,1,i,j,v[i][j]-v[i][j+1]);    
            else if (i == x && j > y) Modify(1,n,1,i,j,v[i][j]-v[i][j-1]);    
            else if (j == y && i < x) Modify(1,n,1,i,j,v[i][j]-v[i+1][j]);    
            else if (j == y && i > x) Modify(1,n,1,i,j,v[i][j]-v[i-1][j]);    
            else if (i < x && j < y) Modify(1,n,1,i,j,v[i][j]-v[i+1][j]-v[i][j+1]+v[i+1][j+1]);    
            else if (i < x && j > y) Modify(1,n,1,i,j,v[i][j]-v[i+1][j]-v[i][j-1]+v[i+1][j-1]);    
            else if (i > x && j < y) Modify(1,n,1,i,j,v[i][j]-v[i-1][j]-v[i][j+1]+v[i-1][j+1]);    
            else Modify(1,n,1,i,j,v[i][j]-v[i][j-1]-v[i-1][j]+v[i-1][j-1]);    
        }  
            
    while(t--) {      
        int typ = getLL();    
        if (typ == 0) {    
            X1 = getLL(); Y1 = getLL();     
            X2 = getLL(); Y2 = getLL();    
            X1 = x - X1; X2 = x + X2;    
            Y1 = y - Y1; Y2 = y + Y2;    
            printf("%lld\n",query(1,1,n,X1,X2,Y1,Y2));    
        }    
        else {    
        	++tt;
            X1 = getLL(); Y1 = getLL();     
            X2 = getLL(); Y2 = getLL();    
            LL c; c = getLL();  cur = 0;  
            Judge(X1-1,Y1-1,c); Judge(X1,Y1-1,c); Judge(X1,Y1,c); Judge(X1-1,Y1,c);    
            Judge(X1-1,Y2,c); Judge(X1-1,Y2+1,c); Judge(X1,Y2+1,c); Judge(X1,Y2,c);    
            Judge(X2,Y1,c);  Judge(X2+1,Y1,c); Judge(X2+1,Y1-1,c); Judge(X2,Y1-1,c);   
            Judge(X2,Y2,c); Judge(X2,Y2+1,c); Judge(X2+1,Y2,c); Judge(X2+1,Y2+1,c);    
            if (X1 <= x && x <= X2 && Y1 <= y && y <= Y2) Modify(1,n,1,x,y,c);    
            if (X1 <= x && x <= X2) {    
                if (Y1 <= y && y <= Y2) Modify(1,n,1,x,Y2+1,-c),Modify(1,n,1,x,Y1-1,-c);    
                else if (y > Y2) Modify(1,n,1,x,Y2,c),Modify(1,n,1,x,Y1-1,-c);    
                else Modify(1,n,1,x,Y1,c),Modify(1,n,1,x,Y2+1,-c);    
            }    
            if (Y1 <= y && y <= Y2) {    
                if (X1 <= x && x <= X2) Modify(1,n,1,X1-1,y,-c),Modify(1,n,1,X2+1,y,-c);    
                else if (x > X2) Modify(1,n,1,X2,y,c),Modify(1,n,1,X1-1,y,-c);    
                else Modify(1,n,1,X1,y,c),Modify(1,n,1,X2+1,y,-c);    
            }    
        }    
    }      
    return 0;      
}



CRZbulabula
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
【NOI2012】魔幻棋盘
weixin_30244681的博客
08-25 210
Description 将要读二年级的小 Q 买了一款新型益智玩具——魔幻棋盘，它是一个N行M列的网格棋盘，每个格子中均有一个正整数。棋盘守护者在棋盘的第X行Y列（行与列均从1开始编号） 并且始终不会移动。棋盘守护者会进行两种操作： （a）询问：他会以自己所在位置为基础，向四周随机扩展出一块大小不定的矩形区域，向你询问这一区域内所有数的最大公约数是多少。 （b）修改：他会随意挑选棋盘上...
BZOJ 2877 NOI2012 魔幻棋盘 二维线段树
世界
01-07 4639
题目大意：给定一个矩阵，支持两种操作： 1.将某个子矩阵中的每个值增加一个数 2.询问某个子矩阵中的所有数的GCD 已知所有询问恒过定点(x,y) 算了BZOJ没有原题我还是把原题发上来吧- - 《论代码长度与注释长度以及题目简单程度的比例失调关系以及日本饮用水资源的解决方案》 《10K+代码是怎样炼成的》 《GCD与修改标记的正确用法》 《出题人我*你吗系列》 《
[NOI2012]魔幻棋盘
aonsfobv0165995的博客
08-30 156
这一道题的话就我而言可能更像数学里的分类讨论？ 就是讨论一下修改矩形的四个角在守护者的哪边就好了，只是我的代码可能有点麻烦，如果有神犇看到的话能给点建议吗，拜谢！ 然后很多东西我会放在代码里面，大佬们就看代码吧 #include <iostream> #include <cstdio> #include <cstdlib> #include <a...
bzoj 2877: [Noi2012]魔幻棋盘
fyc的博客
10-12 337
口头AC
bzoj2877: [Noi2012]魔幻棋盘 树套数+差分
IED98的专栏
07-13 1364
神题，只能膜拜题解。然后发现了差分这个东西，所谓差分就是sigma{a[i]}=b[i]; a[i]=b[i]-b[i-1];(是不是和辗转相减很像。。。。。) 我们再来看一下二维的怎么做。。然后我们换一种维护的方法。因为题目要求一定会查询到守护者a[x, y]，所以我们以(x, y)为原点，建立直角座标系。二维的一样是得差分，对于第四象限的点：b[i, j] = a[i, j] - a[i-1
ACM:NOI:CSP比赛经验分享与代码程序资源
最新发布
05-14
ACM:NOI:CSP比赛经验分享与代码程序资源
NOI2012 测试数据
07-10
NOI 2012 测试数据 day1: random, bicycling, chess day2: park, delicacy, tritown
oi:NOIP，NOI，IOI
05-01
NOI收录了1001~2009年的试题、解析与测试数据。NOIP收录了1995~2018年的初赛试题与参考答案。收录了1995~2018年的复赛试题、测试数据与部分参考程序。书《全国青少年信息学（计算机）奥林匹克分区联赛试题解析（中学...
noi2012浙江省选初测成绩试题数据
04-14
浙江noi2012省队选拔（ZJOI2012）一试初测成绩、试题及测试数据
noi2012试题
12-01
noi比赛试题如何取得题解，大家可以共同学习的PASCAL利器
[BZOJ2877][Noi2012]魔幻棋盘（差分+二维线段树）
女装的你，如此好看！
09-22 353
Address 洛谷P2086 BZOJ2877 LOJ#2672 Solution 调了一整个下午 发现每次询问的矩形一定包含 (X,Y)(X,Y)(X,Y) ，所以我们自然地想到把棋盘拆成 444 个部分，分别为 (X,Y)(X,Y)(X,Y) 的左上部分（横坐标小于等于 XXX 且纵坐标小于等于 YYY ），右上部分（横坐标小于等于 XXX 且纵坐标大于等于 YYY ），左下部分，右下部分。...
bzoj2877 [Noi2012]魔幻棋盘 二维线段树
Todobe
07-02 611
题目大意： 维护一个矩阵支持两个操作： 0、查询子矩形gcd。 1、子矩形整体加一个值。 所有的询问都会一定包含同一个点。题目分析： PoPoQQQ讲解的非常好啊这个…… 二维线段树可以资磁子矩形查询，但是它不资磁子矩形修改啊orz 由于gcd(a,b)=gcd(a-b,b)，所以我们可以考虑差分。 但是上面这个式子使用是有前提的，就是把一些数差分之后，其中必须有一个数时原数才行，恰
[BZOJ2877][Noi2012]魔幻棋盘 && 二维线段树+差分
shiyukun1998的专栏
06-03 1535
二维线段树原来就是线段树的节点上加个线段树- -  差分以(x, y)为中心向四周进行差分 这样就不用分别维护一阶差分和二阶差分什么的（如果以(1, 1) 为中心进行差分还要分开维护 而且涉及区间修改根本无法实现嘛） 还有hlq神犇提出了8二维树状数组的方法（太可怕了） 膜膜膜 这代码就6K啦？ 有点不敢相信 #include #include #include #include #inc
【二维线段树（二维区间GCD）】[NOI2012]魔幻棋盘
weixin_30446197的博客
07-09 173
题目描述 分析 这是经典的区间gcd（最大公约数）问题。 差分GCD gcd(a,b)gcd(a,b,c)=gcd(a,b−ka)=gcd(gcd(a,b),c)=gcd(gcd(a,b−a),c)=gcd(gcd(a,b−a),c−k(gcd(a,b−a)))=gcd(gcd(a,b−a),c−b) 所以区间[l,r]的最大公约数等于其差分和第一个数的最...
[BZOJ2877][NOI2012]魔幻棋盘(二维线段树)
weixin_30745553的博客
03-20 90
https://blog.sengxian.com/solutions/bzoj-2877 注意二维线段树的upd()也是一个O(log n)的函数（pushdown()应该也是但没写过）。 1 #include<cstdio> 2 #include<algorithm> 3 #define Ls (x<<1) 4 #define Rs (...
BZOJ2877 NOI2012魔幻棋盘（二维线段树）
weixin_30236595的博客
02-01 124
　　显然一个序列的gcd=gcd(其差分序列的gcd，序列中第一个数)。于是一维情况直接线段树维护差分序列即可。 　　容易想到将该做法拓展到二维。于是考虑维护二维差分，查询时对差分矩阵求矩形的gcd，再对矩形的两个边界求一下原本的gcd即可。 　　但这样大概需要三个二维线段树，空间可能不太够。由于查询区域是由一个给定点拓展的，可以改为以该点为中心建差分矩阵，这样剩下部分是一个十字形，可以直接一...
【LOJ】#2672. 「NOI2012」魔幻棋盘
weixin_30746117的博客
06-13 79
题解 代码7.1KB，人傻代码长 恶心死我了这代码真的把我写恶心了= = 想一想就知道这个东西……维护到中心的差分，然后用二维线段树维护一下矩形的gcd 嗯，我说完了，你写吧。 首先这个二维线段树的单点修改啊，要把第一维遍历过的节点也修改了，还不能直接修改，需要把第一维的左右子树里套的线段树拿出来合并 然后这个差分啊，我为了好搞一点，拆成了9个矩阵，也就是 中心，中心上方，中心下方，中心...
bzoj2877 [Noi2012]魔幻棋盘 [二维线段树]
Shion的专栏
06-02 1227
http://www.lydsy.com/JudgeOnline/problem.php?id=2877 题解 人生苦短。 #include #define p(i, j) (((i) - 1) * m + (j)) #define g(i, j) (((i) - 1) * ((m << 2) + 100) + (j)) #define MAXN 555555 #define lch
noi excel公式
12-30
NOI是资产净收入的缩写，是指净营业收入减去直接经营成本和间接经营成本后的净额。在财务报表中，NOI扮演着衡量企业盈利能力的关键角色，特别是在房地产领域。 计算NOI的公式如下： NOI = 净营业收入 - 直接经营成本 - 间接经营成本 净营业收入是指企业在一定时期内从营业活动中所获得的净收入，包括租金、销售收入等。直接经营成本是与直接生产或提供服务相关的成本，如物料和人工成本。间接经营成本是与间接经营活动相关的成本，如行政费用、市场营销费用等。 通过计算NOI，企业可以了解其核心业务的盈利能力。这对于投资者和管理者来说都非常重要。投资者可以通过比较不同企业的NOI来评估其收益情况，从而进行投资决策。对于企业管理者来说，NOI是评估经营绩效的重要指标，可以用于监控和改进业务运营。 总之，NOI是计算企业净收入的重要指标，在财务分析和投资决策中具有重要意义。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值