数据结构(二维线段树,差分): NOI2012 魔幻棋盘
貌似想复杂了……
#include <iostream>
#include <cstring>
#include <cstdio>
#define mid ((l+r)>>1)
using namespace std;
const int maxn=;
const int maxm=;
int N,M,Q;
struct Array{
long long num[maxm];
long long *operator [](int x){
return &num[(x-)*M];
}
}a,b; long long ABS(long long x){return x>?x:-x;}
long long GCD(long long x,long long y){
return y?GCD(y,x%y):ABS(x);
} int cnt,ch[maxn][];
long long w[maxn]; void Merge(int &x,int p1,int p2,int l,int r){
if(!x)x=++cnt;
w[x]=GCD(w[p1],w[p2]);
if(l!=r){
Merge(ch[x][],ch[p1][],ch[p2][],l,mid);
Merge(ch[x][],ch[p1][],ch[p2][],mid+,r);
}
} void Build(int &x,int l,int r,long long t[]){
x=++cnt;
if(l==r){
w[x]=t[l];
return;
}
Build(ch[x][],l,mid,t);
Build(ch[x][],mid+,r,t);
w[x]=GCD(w[ch[x][]],w[ch[x][]]);
} int rt[maxm<<];
void Build(int x,int l,int r){
if(l==r){
Build(rt[x],,M,a[l]);
return;
}
Build(x<<,l,mid);
Build(x<<|,mid+,r);
Merge(rt[x],rt[x<<],rt[x<<|],,M);
} int x,y,tx,ty;
int tp,x1,y1,x2,y2;
long long d; void Change(int x,int l,int r,long long d){
if(l==r){w[x]+=d;return;}
if(mid>=ty)Change(ch[x][],l,mid,d);
else Change(ch[x][],mid+,r,d);
w[x]=GCD(w[ch[x][]],w[ch[x][]]);
} void Update(int x,int p1,int p2,int l,int r){
w[x]=GCD(w[p1],w[p2]);
if(l==r)return;
if(mid>=ty)Update(ch[x][],ch[p1][],ch[p2][],l,mid);
else Update(ch[x][],ch[p1][],ch[p2][],mid+,r);
} void Modify(int x,int l,int r,long long d){
if(l==r){
Change(rt[x],,M,d);
return;
}
if(mid>=tx)Modify(x<<,l,mid,d);
else Modify(x<<|,mid+,r,d);
Update(rt[x],rt[x<<],rt[x<<|],,M);
} long long Que(int x,int l,int r){
if(l>=y1&&r<=y2)
return w[x];
long long ret=;
if(mid>=y1)ret=Que(ch[x][],l,mid);
if(mid<y2)ret=GCD(ret,Que(ch[x][],mid+,r));
return ret;
} long long Query(int x,int l,int r){
if(l>=x1&&r<=x2)
return Que(rt[x],,M);
long long ret=;
if(mid>=x1)ret=Query(x<<,l,mid);
if(mid<x2)ret=GCD(ret,Query(x<<|,mid+,r));
return ret;
} int main(){
freopen("chessa.in","r",stdin);
freopen("chessa.out","w",stdout);
scanf("%d%d",&N,&M);
scanf("%d%d",&x,&y);
scanf("%d",&Q);
for(int i=;i<=N;i++)
for(int j=;j<=M;j++)
scanf("%lld",&b[i][j]);
for(int i=;i<=N;i++)
for(int j=;j<=M;j++){
a[i][j]=b[i][j];
if(i<x)a[i][j]-=b[i+][j];
if(i>x)a[i][j]-=b[i-][j];
if(j<y)a[i][j]-=b[i][j+];
if(j>y)a[i][j]-=b[i][j-];
if(i!=x&&j!=y)a[i][j]+=b[i+(i<x?:-)][j+(j<y?:-)];
} Build(,,N);
while(Q--){
scanf("%d",&tp);
scanf("%d%d",&x1,&y1);
scanf("%d%d",&x2,&y2);
if(tp==){
x1=x-x1;x2=x+x2;
y1=y-y1;y2=y+y2;
printf("%lld\n",Query(,,N));
}
if(tp==){
scanf("%lld",&d); {
if(x1<=x&&y1<=y){
tx=x1-;ty=y1-;
if(tx&&ty)Modify(,,N,d);
} if(x1<=x&&y1>y){
tx=x1-;ty=y1;
if(tx)Modify(,,N,-d);
} if(x1>x&&y1<=y){
tx=x1;ty=y1-;
if(ty)Modify(,,N,-d);
} if(x1>x&&y1>y){
tx=x1;ty=y1;
Modify(,,N,d);
}
} {
if(x1<=x&&y2<y){
tx=x1-;ty=y2;
if(tx)Modify(,,N,-d);
} if(x1<=x&&y2>=y){
tx=x1-;ty=y2+;
if(tx&&ty<=M)Modify(,,N,d);
} if(x1>x&&y2<y){
tx=x1;ty=y2;
Modify(,,N,d);
} if(x1>x&&y2>=y){
tx=x1;ty=y2+;
if(ty<=M)Modify(,,N,-d);
}
} {
if(x2<x&&y1<=y){
tx=x2;ty=y1-;
if(ty)Modify(,,N,-d);
} if(x2<x&&y1>y){
tx=x2;ty=y1;
Modify(,,N,d);
} if(x2>=x&&y1<=y){
tx=x2+;ty=y1-;
if(tx<=N&&ty)Modify(,,N,d);
} if(x2>=x&&y1>y){
tx=x2+;ty=y1;
if(tx<=N)Modify(,,N,-d);
}
} {
if(x2<x&&y2<y){
tx=x2;ty=y2;
Modify(,,N,d);
} if(x2<x&&y2>=y){
tx=x2;ty=y2+;
if(ty<=M)Modify(,,N,-d);
} if(x2>=x&&y2<y){
tx=x2+;ty=y2;
if(tx<=N)Modify(,,N,-d);
} if(x2>=x&&y2>=y){
tx=x2+;ty=y2+;
if(tx<=N&&ty<=M)Modify(,,N,d);
}
} if(x1<=x&&x2>=x){
if(y1<=y){
tx=x;ty=y1-;
if(ty)Modify(,,N,-d);
}
if(y1>y){
tx=x;ty=y1;
Modify(,,N,d);
} if(y2<y){
tx=x;ty=y2;
Modify(,,N,d);
}
if(y2>=y){
tx=x;ty=y2+;
if(ty<=M)Modify(,,N,-d);
}
} if(y1<=y&&y2>=y){
if(x1<=x){
tx=x1-;ty=y;
if(tx)Modify(,,N,-d);
}
if(x1>x){
tx=x1;ty=y;
Modify(,,N,d);
} if(x2<x){
tx=x2;ty=y;
Modify(,,N,d);
}
if(x2>=x){
tx=x2+;ty=y;
if(tx<=N)Modify(,,N,-d);
}
} if(x1<=x&&x2>=x&&y1<=y&&y2>=y){
tx=x;ty=y;
Modify(,,N,d);
}
}
}
return ;
}
数据结构(二维线段树,差分): NOI2012 魔幻棋盘的更多相关文章
- POJ 2155 Matrix (二维线段树)
Matrix Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 17226 Accepted: 6461 Descripti ...
- HDU 4819 Mosaic (二维线段树)
Mosaic Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 102400/102400 K (Java/Others)Total S ...
- BZOJ2877 NOI2012魔幻棋盘(二维线段树)
显然一个序列的gcd=gcd(其差分序列的gcd,序列中第一个数).于是一维情况直接线段树维护差分序列即可. 容易想到将该做法拓展到二维.于是考虑维护二维差分,查询时对差分矩阵求矩形的gcd,再对矩形 ...
- NOI 2012 魔幻棋盘 | 二维差分 + 二维线段树
题目:luogu 2086 二维线段树,按套路差分原矩阵,gcd( x1, x2, ……, xn ) = gcd( xi , x2 - x1 , ……, xn - xn-1 ),必须要有一个原数 xi ...
- [BZOJ2877][NOI2012]魔幻棋盘(二维线段树)
https://blog.sengxian.com/solutions/bzoj-2877 注意二维线段树的upd()也是一个O(log n)的函数(pushdown()应该也是但没写过). #inc ...
- tyvj P1716 - 上帝造题的七分钟 二维树状数组区间查询及修改 二维线段树
P1716 - 上帝造题的七分钟 From Riatre Normal (OI)总时限:50s 内存限制:128MB 代码长度限制:64KB 背景 Background 裸体就意味着 ...
- [JZOJ3615]【NOI2014模拟】数列(平面几何+二维线段树)
Description 给定一个长度为n的正整数数列a[i]. 定义2个位置的f值为两者位置差与数值差的和,即f(x,y)=|x-y|+|a[x]-a[y]|. 你需要写一个程序支持2种操作(k都是正 ...
- poj 2155 matrix 二维线段树 线段树套线段树
题意 一个$n*n$矩阵,初始全为0,每次翻转一个子矩阵,然后单点查找 题解 任意一种能维护二维平面的数据结构都可以 我这里写的是二维线段树,因为四分树的写法复杂度可能会退化,因此考虑用树套树实现二维 ...
- BZOJ 3489 A simple rmq problem 可持久化KDtree/二维线段树
题目链接:https://www.lydsy.com/JudgeOnline/problem.php?id=3489 题意概述: 给出一个序列,每次询问一个序列区间中仅出现了一次的数字最大是多少,如果 ...
- bzoj 1513 POI2006 Tet-Tetris 3D 二维线段树+标记永久化
1511: [POI2006]OKR-Periods of Words Time Limit: 5 Sec Memory Limit: 64 MBSubmit: 351 Solved: 220[S ...
随机推荐
- CentOS 6.7 编译安装Nginx 1.8.0
1.配置编译环境 yum update && yum upgrade yum groupinstall "Development Tools" 或者 yum ins ...
- rabbitmq 消息持久化
rabbitmq 消息持久化 2016-02-18 11:19 224人阅读 评论(0) 收藏 举报 分类: 综合(15) 版权声明:本文为博主原创文章,未经博主允许不得转载. 二: 任务分发 & ...
- 我是一块cpu 《转载》
我是一块cpu,原装intel,在一台普通的台式计算机里供职.我有个小弟是内存,我要靠他时时刻刻陪伴我工作,其实有时候我并不是没有某某地址的资料,而是懒得翻--麻烦. 还有一个老大哥叫bios,每次那 ...
- Activity间的跳转,startActivity与startActivityForResult
JreduCh04 2016-07-30跳转 (由一个画面跳转到另一个画面)两种方法:Intent中 startActivity.startActivityForResult.后者可设置request ...
- Activity Threa创建Window和View分析
http://blog.csdn.net/ljsbuct/article/details/7094580 1. 入口. 以前一直都说Activity的人口是onCreate方法.其实android上一 ...
- Android开发app如何设定应用图标下的应用名称为汉字以及自定义图标
一.应用名称为汉字 二.自定义图标
- 【转】 KVC/KVO原理详解及编程指南
原文地址:http://blog.csdn.net/wzzvictory/article/details/9674431 前言: 1.本文基本不讲KVC/KVO的用法,只结合网上的资料说说对这种技术的 ...
- javaScript 删除数组中指定元素
Array.prototype.indexOf = function(val) { for (var i = 0; i < this.length; i++) { if (this[i] == ...
- PHP算法 《树形结构》 之 伸展树(1) - 基本概念
伸展树的介绍 1.出处:http://dongxicheng.org/structure/splay-tree/ A. 概述 二叉查找树(Binary Search Tree,也叫二叉排序树,即Bin ...
- Sqoop import加载HBase案例详解
简单写一下如何将订单表sqoop到hbase表中的步骤. 下表: 1.通过hbase shell 打开hbase. 2.创建一个hbase表 create 'so','o' 3.将so表的数据导入到h ...