[BOI2007]Mokia 摩基亚

Description:

摩尔瓦多的移动电话公司摩基亚（Mokia）设计出了一种新的用户定位系统。和其他的定位系统一样，它能够迅速回答任何形如“用户C的位置在哪？”的问题，精确到毫米。但其真正高科技之处在于，它能够回答形如“给定区域内有多少名用户？”的问题。

在定位系统中，世界被认为是一个W×W的正方形区域，由1×1的方格组成。每个方格都有一个坐标(x,y)，1<=x,y<=W。坐标的编号从1开始。对于一个4×4的正方形，就有1<=x<=4,1<=y<=4（如图）：

请帮助Mokia公司编写一个程序来计算在某个矩形区域内有多少名用户。

Hint:

\(n \le 2*10^5\)

Solution:

考虑把答案拆成四个矩阵,前缀和运算一下,就是cdq裸题了

#include <map>
#include <set>
#include <stack>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#define ls p<<1 
#define rs p<<1|1
using namespace std;
typedef long long ll;
const int mxn=2e6+5;
int n,m,cnt,tot,t[mxn],hd[mxn];

inline int read() {
	char c=getchar(); int x=0,f=1;
	while(c>'9'||c<'0') {if(c=='-') f=-1;c=getchar();}
	while(c<='9'&&c>='0') {x=(x<<3)+(x<<1)+(c&15);c=getchar();}
	return x*f;
}
inline void chkmax(int &x,int y) {if(x<y) x=y;}
inline void chkmin(int &x,int y) {if(x>y) x=y;}

struct Q {
	int id,x,y,val,opt;
}q[mxn<<2];

int cmp(Q x,Q y) {
	return x.id<y.id;
}

int cmp1(Q a,Q b) {
	return a.x==b.x?a.y<b.y:a.x<b.x;
}
 
int lb(int x) {return x&-x;}

void mod(int x,int y) {
	while(x<=n) t[x]+=y,x+=lb(x);
}

int query(int x) {
	int res=0;
	while(x) res+=t[x],x-=lb(x);
	return res;
}

void cdq(int l,int r) {
	if(l==r) return ;
	int mid=(l+r)>>1; 
	cdq(l,mid); cdq(mid+1,r);
	sort(q+l,q+mid+1,cmp1); 
	sort(q+mid+1,q+r+1,cmp1);
	int lp=l;
	for(int i=mid+1;i<=r;++i) {
		while(q[lp].x<=q[i].x&&lp<=mid) { //这里没取等调了半年
			if(!q[lp].opt) mod(q[lp].y,q[lp].val);
			++lp;
		}
		if(q[i].opt) q[i].val+=query(q[i].y);
	}
	for(int i=l;i<lp;++i) 
		if(!q[i].opt) mod(q[i].y,-q[i].val);
}


int main()
{
	int opt,x,y,w,z;
	n=read(); n=read()+1; opt=read();
	while(opt!=3) {
		if(opt==1) 
			++tot, q[tot].id=tot,q[tot].x=read()+1,q[tot].y=read()+1,q[tot].val=read();
		else {
			x=read(); y=read(); w=read()+1,z=read()+1;
			++tot, q[tot]=(Q){tot,x,y,0,1};
			++tot, q[tot]=(Q){tot,w,z,0,1};
			++tot, q[tot]=(Q){tot,w,y,0,1};
			++tot, q[tot]=(Q){tot,x,z,0,1};
		}
		opt=read();
	}
	cdq(1,tot); sort(q+1,q+tot+1,cmp);
	for(int i=1;i<=tot;++i) 
		if(q[i].opt) 
			printf("%d\n",q[i].val+q[i+1].val-q[i+2].val-q[i+3].val),i+=3;
    return 0;
}