BZOJ2877 NOI2012魔幻棋盘（二维线段树）
显然一个序列的gcd=gcd(其差分序列的gcd，序列中第一个数)。于是一维情况直接线段树维护差分序列即可。
容易想到将该做法拓展到二维。于是考虑维护二维差分，查询时对差分矩阵求矩形的gcd，再对矩形的两个边界求一下原本的gcd即可。
但这样大概需要三个二维线段树，空间可能不太够。由于查询区域是由一个给定点拓展的，可以改为以该点为中心建差分矩阵，这样剩下部分是一个十字形，可以直接一维线段树维护，就只需要一个二维线段树了。
注意题面有锅，详见discuss，被坑了一年。
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<vector>
#include<map>
using namespace std;
#define ll long long
#define N 500010
#define pii pair<int,int>
#define PII pair< pii , pii >
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
ll gcd(ll n,ll m){return m==0?n:gcd(m,n%m);}
ll read()
{
ll x=0,f=1;char c=getchar();
while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}
while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();
return x*f;
}
int n,m,T,X,Y,root[2],cnt[2],CNT,ROOT;
ll BIT[2][N];
vector<ll> a[N];
map<PII,int> id;
struct data{int l,r,L,R;ll gcd;
}tree[2][N<<1],TREE[N*20];
void BIT_add(int op,int n,int x,ll y){while (x<=n) BIT[op][x]+=y,x+=x&-x;}
ll BIT_query(int op,int x){ll s=0;while (x) s+=BIT[op][x],x-=x&-x;return s;}
ll calc(int i,int j)
{
if (i<X&&j<Y) return a[i][j]-a[i+1][j]-a[i][j+1]+a[i+1][j+1];
if (i<X&&j>Y) return a[i][j]-a[i+1][j]-a[i][j-1]+a[i+1][j-1];
if (i>X&&j<Y) return a[i][j]-a[i-1][j]-a[i][j+1]+a[i-1][j+1];
if (i>X&&j>Y) return a[i][j]-a[i-1][j]-a[i][j-1]+a[i-1][j-1];
return 0;
}
void add(int &k,int op,int l,int r,int x,ll p)
{
if (!k) k=++cnt[op];
if (l==r) {tree[op][k].gcd+=p;return;}
int mid=l+r>>1;
if (x<=mid) add(tree[op][k].l,op,l,mid,x,p);
else add(tree[op][k].r,op,mid+1,r,x,p);
tree[op][k].gcd=gcd(tree[op][tree[op][k].l].gcd,tree[op][tree[op][k].r].gcd);
}
ll query(int k,int op,int l,int r,int x,int y)
{
if (x>y||!k) return 0;
if (l==x&&r==y) return tree[op][k].gcd;
int mid=l+r>>1;
if (y<=mid) return query(tree[op][k].l,op,l,mid,x,y);
else if (x>mid) return query(tree[op][k].r,op,mid+1,r,x,y);
else return gcd(query(tree[op][k].l,op,l,mid,x,mid),query(tree[op][k].r,op,mid+1,r,mid+1,y));
}
void BUILD(int &k,int l,int r,int x)
{
id[make_pair(make_pair(x,x),make_pair(l,r))]=k=++CNT;
if (l==r) {TREE[k].gcd=calc(x,l);return;}
int mid=l+r>>1;
BUILD(TREE[k].l,l,mid,x);
BUILD(TREE[k].r,mid+1,r,x);
TREE[k].gcd=gcd(TREE[TREE[k].l].gcd,TREE[TREE[k].r].gcd);
}
void BUILD2(int &k,int l,int r,int u,int x,int y)
{
if (!k) id[make_pair(make_pair(l,r),make_pair(x,y))]=k=++CNT;
TREE[k].L=id[make_pair(make_pair(l,u),make_pair(x,y))];
TREE[k].R=id[make_pair(make_pair(u+1,r),make_pair(x,y))];
TREE[k].gcd=gcd(TREE[TREE[k].L].gcd,TREE[TREE[k].R].gcd);
if (x==y) return;
int mid=x+y>>1;
BUILD2(TREE[k].l,l,r,u,x,mid);
BUILD2(TREE[k].r,l,r,u,mid+1,y);
}
void update(int k,int l,int r,int x)
{
TREE[k].gcd=gcd(TREE[TREE[k].L].gcd,TREE[TREE[k].R].gcd);
if (l==r) return;
int mid=l+r>>1;
if (x<=mid) update(TREE[k].l,l,mid,x);
else update(TREE[k].r,mid+1,r,x);
}
void build(int &k,int l,int r)
{
if (l==r) {BUILD(k,0,m+1,l);return;}
id[make_pair(make_pair(l,r),make_pair(0,m+1))]=k=++CNT;
int mid=l+r>>1;
build(TREE[k].L,l,mid);
build(TREE[k].R,mid+1,r);
BUILD2(k,l,r,(l+r>>1),0,m+1);
}
void ADD(int &k,int l,int r,int x,ll p)
{
if (l==r) {TREE[k].gcd+=p;return;}
int mid=l+r>>1;
if (x<=mid) ADD(TREE[k].l,l,mid,x,p);
else ADD(TREE[k].r,mid+1,r,x,p);
TREE[k].gcd=gcd(TREE[TREE[k].l].gcd,TREE[TREE[k].r].gcd);
}
void Add(int k,int l,int r,int x,int y,ll p)
{
if (l==r) {ADD(k,0,m+1,y,p);return;}
int mid=l+r>>1;
if (x<=mid) Add(TREE[k].L,l,mid,x,y,p);
else Add(TREE[k].R,mid+1,r,x,y,p);
update(k,0,m+1,y);
}
ll QUERY(int k,int l,int r,int x,int y)
{
if (l==x&&r==y) return TREE[k].gcd;
int mid=l+r>>1;
if (y<=mid) return QUERY(TREE[k].l,l,mid,x,y);
else if (x>mid) return QUERY(TREE[k].r,mid+1,r,x,y);
else return gcd(QUERY(TREE[k].l,l,mid,x,mid),QUERY(TREE[k].r,mid+1,r,mid+1,y));
}
ll Query(int k,int l,int r,int xl,int xr,int yl,int yr)
{
if (xl>xr||yl>yr) return 0;
if (l==xl&&r==xr) return QUERY(k,0,m+1,yl,yr);
int mid=l+r>>1;
if (xr<=mid) return Query(TREE[k].L,l,mid,xl,xr,yl,yr);
else if (xl>mid) return Query(TREE[k].R,mid+1,r,xl,xr,yl,yr);
else return gcd(Query(TREE[k].L,l,mid,xl,mid,yl,yr),Query(TREE[k].R,mid+1,r,mid+1,xr,yl,yr));
}
signed main()
{
#ifndef ONLINE_JUDGE
freopen("bzoj2877.in","r",stdin);
freopen("bzoj2877.out","w",stdout);
const char LL[]="%I64d\n";
#else
const char LL[]="%lld\n";
#endif
n=read(),m=read(),X=read(),Y=read(),T=read();
for (int j=0;j<=m+1;j++) a[0].push_back(0);
for (int i=1;i<=n;i++)
{
a[i].push_back(0);
for (int j=1;j<=m;j++)
a[i].push_back(read());
a[i].push_back(0);
}
for (int j=0;j<=m+1;j++) a[n+1].push_back(0);
build(ROOT,0,n+1);
for (int i=1;i<=n;i++)
BIT_add(0,n,i,a[i][Y]-a[i-1][Y]),add(root[0],0,1,n,i,a[i][Y]-a[i-1][Y]);
for (int j=1;j<=m;j++)
BIT_add(1,m,j,a[X][j]-a[X][j-1]),add(root[1],1,1,m,j,a[X][j]-a[X][j-1]);
while (T--)
{
int op=read();
if (op==0)
{
int up=read(),left=read(),down=read(),right=read();
ll ans=gcd(BIT_query(0,X-up),query(root[0],0,1,n,X-up+1,X+down));
ans=gcd(ans,gcd(BIT_query(1,Y-left),query(root[1],1,1,m,Y-left+1,Y+right)));
ans=gcd(ans,Query(1,0,n+1,X-up,X-1,Y-left,Y-1));
ans=gcd(ans,Query(1,0,n+1,X-up,X-1,Y+1,Y+right));
ans=gcd(ans,Query(1,0,n+1,X+1,X+down,Y-left,Y-1));
ans=gcd(ans,Query(1,0,n+1,X+1,X+down,Y+1,Y+right));
printf(LL,abs(ans));
}
else
{
int xl=read(),yl=read(),xr=read(),yr=read();ll c=read();
if (yl<=Y&&Y<=yr)
{
BIT_add(0,n,xl,c),add(root[0],0,1,n,xl,c);
if (xr<n) BIT_add(0,n,xr+1,-c),add(root[0],0,1,n,xr+1,-c);
}
if (xl<=X&&X<=xr)
{
BIT_add(1,m,yl,c),add(root[1],1,1,m,yl,c);
if (yr<m) BIT_add(1,m,yr+1,-c),add(root[1],1,1,m,yr+1,-c);
}
if (xl<=X&&yl<=Y)
{
if (xr<X&&yr<Y) Add(1,0,n+1,xr,yr,c);
if (xr<X) Add(1,0,n+1,xr,yl-1,-c);
if (yr<Y) Add(1,0,n+1,xl-1,yr,-c);
Add(1,0,n+1,xl-1,yl-1,c);
}
if (xr>=X&&yl<=Y)
{
if (xl>X&&yr<Y) Add(1,0,n+1,xl,yr,c);
if (xl>X) Add(1,0,n+1,xl,yl-1,-c);
if (yr<Y) Add(1,0,n+1,xr+1,yr,-c);
Add(1,0,n+1,xr+1,yl-1,c);
}
if (xl<=X&&yr>=Y)
{
if (xr<X&&yl>Y) Add(1,0,n+1,xr,yl,c);
if (xr<X) Add(1,0,n+1,xr,yr+1,-c);
if (yl>Y) Add(1,0,n+1,xl-1,yl,-c);
Add(1,0,n+1,xl-1,yr+1,c);
}
if (xr>=X&&yr>=Y)
{
if (xl>X&&yl>Y) Add(1,0,n+1,xl,yl,c);
if (xl>X) Add(1,0,n+1,xl,yr+1,-c);
if (yl>Y) Add(1,0,n+1,xr+1,yl,-c);
Add(1,0,n+1,xr+1,yr+1,c);
}
}
}
return 0;
}
BZOJ2877 NOI2012魔幻棋盘（二维线段树）的更多相关文章
- [BZOJ2877][NOI2012]魔幻棋盘(二维线段树)
https://blog.sengxian.com/solutions/bzoj-2877 注意二维线段树的upd()也是一个O(log n)的函数(pushdown()应该也是但没写过). #inc ...
- NOI 2012 魔幻棋盘 | 二维差分 + 二维线段树
题目:luogu 2086 二维线段树,按套路差分原矩阵,gcd( x1, x2, ……, xn ) = gcd( xi , x2 - x1 , ……, xn - xn-1 ),必须要有一个原数 xi ...
- BZOJ2877 [Noi2012]魔幻棋盘
本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000 作者博客:http://www.cnblogs.com/ljh2000-jump/ ...
- UVA 11297 线段树套线段树(二维线段树)
题目大意: 就是在二维的空间内进行单个的修改,或者进行整块矩形区域的最大最小值查询 二维线段树树,要注意的是第一维上不是叶子形成的第二维线段树和叶子形成的第二维线段树要 不同的处理方式,非叶子形成的 ...
- POJ2155 Matrix二维线段树经典题
题目链接 二维树状数组 #include<iostream> #include<math.h> #include<algorithm> #include<st ...
- HDU 1823 Luck and Love（二维线段树）
之前只知道这个东西的大概概念,没具体去写,最近呵呵,今补上. 二维线段树 -- 点更段查 #include <cstdio> #include <cstring> #inclu ...
- poj 2155:Matrix（二维线段树，矩阵取反，好题）
Matrix Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 17880 Accepted: 6709 Descripti ...
- poj 1195:Mobile phones（二维线段树，矩阵求和）
Mobile phones Time Limit: 5000MS Memory Limit: 65536K Total Submissions: 14391 Accepted: 6685 De ...
- POJ 2155 Matrix （二维线段树）
Matrix Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 17226 Accepted: 6461 Descripti ...
随机推荐
- 二十：让行内元素在div中垂直居中
(1)使用display:table-cell配合vertical-align:center(淘宝也是这样用的) <div class="method4"> <s ...
- H5 68-伪元素选择器
<!DOCTYPE html> <html lang="en"> <head> <meta charset="UTF-8&quo ...
- @RequestBody ajax 415 400
使用springmvc和Ajax进行数据交互时使用标签@RequestBody时我报了这两个错,刚开始对springmvc的使用和注解有点迷,然后踩坑上了. 先说下怎么才会踩上去.首先@Request ...
- Python—包介绍
包(Package) 当你的模块文件越来越多,就需要对模块文件进行划分,比如把负责跟数据库交互的都放一个文件夹,把与页面交互相关的放一个文件夹, . └── my_proj ├── crm #代码目录 ...
- 【RSYSLOG】The Property Replacer【转】
最近在调整日志平台的日志格式,一下是RSYSLOG的 Property Replacer 说明.鉴于RSYSLOG官网略坑,转发一下,原地址忘记了- - ||| The property replac ...
- 1168: mxh对lfx的询问（前缀和+素数表）
题目描述: AS WE ALL KNOW, lfx是咱们组的神仙,但是mxh想考一考lfx一个简单的问题,以此看一下lfx到底是不是神仙.但是lfx要准备补考,于是请你来帮忙回答问题: 给定一个整数N ...
- 1171: lfx捧杯稳啦！
escription Lfx在复习离散的时候突然想到了一个算法题,毕竟是lfx, 算法题如下: 他想知道这样的问题,先定义1~n中即是3的倍数,又是11的倍数的那些数的和sum, 他想知道sum有多少 ...
- 正则表达式验证input文本框
方便以后的查找,直接copy代码在这里了. eg: //公司邮箱验证 if ($("#Email").val() != "") { var myreg = /^ ...
- hadoop和java 配置环境变量的的tar
第一步:打开工具上传tar包 如下图 第二步:在文件路径下查看是否上传成功 第三步:解压tar包 tar -zxvf hadoop.2.6.5.tar.gz 第四步:配置环 ...
- Git之项目使用
现在最为盛行的版本控制器,非git莫属了, 那就看看在项目中我们是如何使用它的吧 一. 在已经存在秘钥对的情况下,我们需要在本地进行相关配置 git config --global user.name ...