bzoj2683&&bzoj4066
题解:
前一题不是强制在线,后一题是强制在线
树套树空间会炸
说一下cdq分治+树状数组
首先我们利用cdq分治使得查询和操作保证先后关系
然后矩阵查询变成4个矩阵的差
那么我们就可以运用扫描线的方法来维护了
时间nlogn^2,空间O(n)
后一题是kd-tree
查询的方法和线段树基本一样
如果矩阵被包含就返回答案,如果不被包含就直接退出
否则递归下去
然后修改的话和替罪羊树一样
达到一定时候就重构
注意一下x相同insert和build要保证y的顺序区分左右
另外学习了一下map中用struct作为key的方法 虽然这题完全不用用
代码:
#include <bits/stdc++.h>
using namespace std;
#define rint register int
#define IL inline
#define rep(i,h,t) for (rint i=h;i<=t;i++)
#define dep(i,t,h) for (rint i=t;i>=h;i--)
#define mid ((h+t)/2)
#define me(x) memset(x,0,sizeof(x))
const int INF=2e9;
const int N=6e6+1e4;
char ss[<<],*A=ss,*B=ss;
IL char gc()
{
return A==B&&(B=(A=ss)+fread(ss,,<<,stdin),A==B)?EOF:*A++;
}
template<class T>void read(T &x)
{
rint f=,c; while (c=gc(),c<||c>) if (c=='-') f=-; x=(c^);
while (c=gc(),c>&&c<) x=(x<<)+(x<<)+(c^); x*=f;
}
IL void umax(int &x,int y)
{
if (x<y) x=y;
}
IL void umin(int &x,int y)
{
if (x>y) x=y;
}
struct re{
int d[],v;
}p[N];
int cmp_d,ans,rt,num;
bool cmp(re x,re y)
{
return x.d[cmp_d]<y.d[cmp_d]||(x.d[cmp_d]==y.d[cmp_d]&&x.d[cmp_d^]<y.d[cmp_d^]);
}
struct kd
{
int Mx[N],My[N],Nx[N],Ny[N],count2[N],ls[N],rs[N];
IL void clear()
{
me(count2);
me(ls); me(rs);
}
IL void updata(int x)
{
count2[x]=count2[ls[x]]+count2[rs[x]]+p[x].v;
if (ls[x])
{
umax(Mx[x],Mx[ls[x]]);
umax(My[x],My[ls[x]]);
umin(Nx[x],Nx[ls[x]]);
umin(Ny[x],Ny[ls[x]]);
}
if (rs[x])
{
umax(Mx[x],Mx[rs[x]]);
umax(My[x],My[rs[x]]);
umin(Nx[x],Nx[rs[x]]);
umin(Ny[x],Ny[rs[x]]);
}
}
int build(int h,int t,int o)
{
cmp_d=o; nth_element(p+h,p+mid,p+t+,cmp);
int x=mid;
Mx[x]=Nx[x]=p[x].d[];
My[x]=Ny[x]=p[x].d[];
count2[x]=p[x].v;
if (h!=x) ls[x]=build(h,mid-,o^); else ls[x]=;
if (x!=t) rs[x]=build(mid+,t,o^); else rs[x]=;
updata(x);
return x;
}
void insert(int &k,int x,int y,int z,int o)
{
if (!k)
{
k=++num;
Mx[num]=Nx[num]=p[num].d[];
My[num]=Ny[num]=p[num].d[];
count2[num]=p[num].v=z;
return;
}
if (p[k].d[]==x&&p[k].d[]==y)
{
count2[k]+=z; p[k].v+=z; return;
}
if (!o)
if (x<p[k].d[]||(x==p[k].d[]&&y<p[k].d[])) insert(ls[k],x,y,z,o^);
else insert(rs[k],x,y,z,o^);
else if (y<p[k].d[]||(y==p[k].d[]&&x<p[k].d[])) insert(ls[k],x,y,z,o^);
else insert(rs[k],x,y,z,o^);
updata(k);
}
void query(int k,int x1,int x2,int y1,int y2)
{
if (!k||Mx[k]<x1||Nx[k]>x2||My[k]<y1||Ny[k]>y2) return;
if (x2>=Mx[k]&&Nx[k]>=x1&&y2>=My[k]&&Ny[k]>=y1)
{
ans+=count2[k]; return;
}
if (x1<=p[k].d[]&&p[k].d[]<=x2&&y1<=p[k].d[]&&p[k].d[]<=y2)
ans+=p[k].v;
query(ls[k],x1,x2,y1,y2); query(rs[k],x1,x2,y1,y2);
}
}kd;
struct re1{
int a,b;
};
struct cmp2
{
bool operator() (const re1 x,const re1 y)
{
if (x.a<y.a||(x.a==y.a&&x.b<y.b)) return(); else return();
}
};
map<re1,int,cmp2> M;
int main()
{
int m;
freopen("1.in","r",stdin);
freopen("1.out","w",stdout);
read(m);
int cnt=,rt=,i=;
while ()
{
i++;
ans=;
int k;
read(k);
if (k==) break;
if (k==)
{
int x,y,z;
read(x); read(y); read(z);
x^=ans; y^=ans; z^=ans;
re1 k3=(re1){x,y};
int kk=M[k3];
if (!kk) p[num+].d[]=x,p[num+].d[]=y,kk=num+;
M[k3]=kk;
kd.insert(rt,x,y,z,);
cnt++;
if (cnt%==)
{
kd.clear();
rt=kd.build(,num,);
}
} else
{
int x1,x2,y1,y2;
read(x1),read(y1),read(x2),read(y2);
x1^=ans,x2^=ans,y1^=ans,y2^=ans;
ans=;
kd.query(rt,x1,x2,y1,y2);
cout<<ans<<endl;
}
} return ;
}
bzoj2683&&bzoj4066的更多相关文章
- [BZOJ2683][BZOJ4066]简单题
[BZOJ2683][BZOJ4066]简单题 试题描述 你有一个N*N的棋盘,每个格子内有一个整数,初始时的时候全部为0,现在需要维护两种操作: 命令 参数限制 内容 1 x y A 1<=x ...
- [基本操作] kd 树
概念就不说了吧,网上教程满天飞 学了半天才知道,kd 树实质上只干了两件事情: 1.快速定位一个点 / 矩形 2.有理有据地优化暴力 第一点大概是可以来做二维平面上给点/矩形打标记的问题 第二点大概是 ...
- Bzoj4066 简单题
Time Limit: 50 Sec Memory Limit: 20 MBSubmit: 2185 Solved: 581 Description 你有一个N*N的棋盘,每个格子内有一个整数,初 ...
- BZOJ4066 简单题(KD-Tree)
板子题. #include<iostream> #include<cstdio> #include<cmath> #include<cstdlib> # ...
- 【bzoj4066】 简单题
http://www.lydsy.com/JudgeOnline/problem.php?id=4066 (题目链接) 题意 维护一个矩阵,两个操作,给某一个元素加上A,求其中一个子矩阵的元素之和.强 ...
- [BZOJ2683]简单题/[BZOJ1176][BalkanOI2007]Mokia
[BZOJ2683]简单题 题目大意: 一个\(n\times n(n\le5\times10^5)\)的矩阵,初始时每个格子里的数全为\(0\).\(m(m\le2\times10^5)\)次操作, ...
- 【BZOJ1176】[Balkan2007]Mokia/【BZOJ2683】简单题 cdq分治
[BZOJ1176][Balkan2007]Mokia Description 维护一个W*W的矩阵,初始值均为S.每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=1600 ...
- 【BZOJ4066】简单题(KD-Tree)
[BZOJ4066]简单题(KD-Tree) 题面 BZOJ 题解 如果这题不卡空间,并且不强制在线的话 显然可以用\(CDQ\)分治做 但是它又卡空间又强制在线,于是我们欢快的来用\(KD-Tree ...
- 【bzoj1176】[Balkan2007]Mokia/【bzoj2683】简单题 CDQ分治+树状数组
bzoj1176 题目描述 维护一个W*W的矩阵,初始值均为S(题目描述有误,这里的S没有任何作用!).每次操作可以增加某格子的权值,或询问某子矩阵的总权值.修改操作数M<=160000,询问数 ...
随机推荐
- 使用chttpfile的一个错误
先贴一部分代码 CString strHttpName="http://localhost/TestReg/RegForm.aspx"; // 需要提交数据的页面 CString ...
- kerberos简单介绍
重要术语 1. KDC 全称:key distributed center 作用:整个安全认证过程的票据生成管理服务,其中包含两个服务,AS和TGS 2. AS 全称:authentication s ...
- Bootstrap -- 模态框实现拖拽移动
### 这里实现这个效果 需要引入 jquery-ui.min.js类库 jquery-ui.min.css样式 使用它提供的draggable()方法实现 ### 菜鸟教程 http://www.r ...
- [MySQL]InnoDB引擎的行锁和表锁
1.行锁和表锁 在mysql 的 InnoDB引擎支持行锁,与Oracle不同,mysql的行锁是通过索引加载的,即是行锁是加在索引响应的行上的,要是对应的SQL语句没有走索引,则会全表扫描, 行锁则 ...
- PyJWT 使用
最近要用 Falsk 开发一个大点的后端,为了安全考虑,弃用传统的Cookie验证.转用JWT. 其实 Falsk 有一个 Falsk-JWT 但是我觉得封装的太高,还是喜欢通用的 PyJWT . J ...
- 虚拟机时间同步14 Aug 04:09:18 ntpdate[2941]: no server suitable for synchronization found
因为虚拟机经常挂起,所以需要时间同步 [root@slave1 /root]$ cp /usr/share/zoneinfo/Asia/Shanghai /etc/localtime cp: over ...
- Nginx(./configure --help)
# ./configure --help --help print this message --prefix=PATH set installation prefix --sbin-path=PAT ...
- Oracle12C_____处理数据库01033连接错误问题.sql
Microsoft Windows [版本 6.1.7601] 版权所有 (c) 2009 Microsoft Corporation.保留所有权利. C:\Users\Administrator&g ...
- python基础之迭代器与生成器
一.什么是迭代器: 迭代是Python最强大的功能之一,是访问集合元素的一种方式. 迭代器对象从集合的第一个元素开始访问,直到所有的元素被访问完结束. 迭代器是一个可以记住遍历的位置的对象. 迭代器的 ...
- poj3107树的重心
/*树的重心求法:两次dfs,第一次dfs处理出每个结点的size,以此求每个结点大儿子的size,第二次dfs将每个结点大儿子的size和余下结点数进行比较,所有结点里两个值之间差值最小的那个点就是 ...