【JZOJ 4694】【HDU 5283】火神的鱼

Description

火神最爱的就是吃鱼了，所以某一天他来到了一个池塘边捕鱼。池塘可以看成一个二维的平面，而他的渔网可以看成一个与坐标轴平行的矩形。
池塘里的鱼不停地在水中游动，可以看成一些点。有的时候会有鱼游进渔网，有的时候也会有鱼游出渔网。所以火神不知道什么时候收网才可以抓住最多的鱼，现在他寻求你的帮助。
他对池塘里的每条鱼都给予了一个标号，分别从1到n标号，n表示池塘里鱼的总数。鱼的游动可以概括为两个动作：
1 l r d : 表示标号在[l，r]这个区间内的鱼向x轴正方向游动了d个单位长度。
2 l r d：表示标号在[l，r]这个区间内的鱼向y轴正方向游动了d个单位长度。
在某些时刻火神会询问你现在有多少条他关心的鱼在渔网内(边界上的也算)，请你来帮助他吧。

Solution

先把题目的两个维度拆开来分析，
每个数在变化的过程中，只会不停的增加，也就是每个点最多只会进一次区间，出一次区间，
所有每个点进入区间和出区间的操作可以一个一个来做，
用线段树，线段树内存mi,mi1,分别表示x<X1的最大的x值以及位置，ma,ma1表示x<X2的最大的x以及位置，（注意，一个点可能即是ma也是mi）
y也开一颗差不多的，
同时还要记录一下当前区间x,y都合法的个数（答案），
mi,mi1的初始值为-maxlongint，
预处理时，如果一个点的x,y其中一个越过了合法的区间，很显然这个点就没有用了，它的mi,ma就没有值，
每次区间修改x，ma就有可能>X2，这样就不合法所以，因此我们要找到是哪个不合法，然后把它同时在两颗树中删掉，也有可能mi>X1，把它从mi中删掉，同时看一下它的x,y是否合法，合法就加进s，如果有很多都要做操作，则一个一个来，
修改y一样，

复杂度：O(nlog(n))

Code

#include<iostream>#include<cstdio>#include<cstdlib>#include<cmath>#include<cstring>#include<algorithm>#define fo(i,a,b) for(int i=a;i<=b;i++)#define fod(i,a,b) for(int i=a;i>=b;i--)#define OK(q,w) (X1<=(q)&&(q)<=X2&&Y1<=(w)&&(w)<=Y2)#define B b[E]using namespace std;const int N=30500,INF=2147483640;int read(int &n){    char ch=' ';int q=0,w=1;    for(;(ch!='-')&&((ch<'0')||(ch>'9'));ch=getchar());    if(ch=='-')w=-1,ch=getchar();    for(;ch>='0' && ch<='9';ch=getchar())q=q*10+ch-48;n=q*w;return n;}int m,n,ans,X1,Y1,X2,Y2;struct qwqw{int x,y;}a[N];struct qqww{int i,i1,a,a1,s,la;}b[2][3*N];void merge(int q,int w,int e){    fo(E,0,1)    {        B[q]=B[w];        B[q].s=B[w].s+B[e].s;        if(B[q].i<B[e].i)B[q].i=B[e].i,B[q].i1=B[e].i1;        if(B[q].a<B[e].a)B[q].a=B[e].a,B[q].a1=B[e].a1;    }}void doit(int l,int r,int e){    fo(E,0,1)if(B[e].la)    {        if(B[e].i!=-INF)B[e].i+=B[e].la;        B[e].a+=B[e].la;        if(l!=r)B[e*2].la+=B[e].la,B[e*2+1].la+=B[e].la;        B[e].la=0;    }}void build(int l,int r,int e){    if(l==r)    {        if(a[l].x<=X2&&a[l].y<=Y2)        {            b[0][e].a=a[l].x,b[1][e].a=a[l].y;            b[0][e].a1=b[1][e].a1=l;            if(OK(b[0][e].a,b[1][e].a))b[0][e].s=1;            if(a[l].x<X1)b[0][e].i=a[l].x,b[0][e].i1=l;            if(a[l].y<Y1)b[1][e].i=a[l].y,b[1][e].i1=l;        }        return;    }    int mid=(l+r)/2;    build(l,mid,e*2),build(mid+1,r,e*2+1);    merge(e,e*2,e*2+1);}void change(int l,int r,int e,int E,int l1,int r1,int l2){    if(l==l1&&r==r1){B[e].la+=l2;doit(l,r,e);    return;}    int mid=(l+r)/2;    doit(l,mid,e*2),doit(mid+1,r,e*2+1);    if(r1<=mid)change(l,mid,e*2,E,l1,r1,l2);        else if(mid<l1)change(mid+1,r,e*2+1,E,l1,r1,l2);            else change(l,mid,e*2,E,l1,mid,l2),change(mid+1,r,e*2+1,E,mid+1,r1,l2);    merge(e,e*2,e*2+1);}void delt(int l,int r,int e,int l2){    if(l==r){b[0][e]=b[1][e]=b[0][0];return;}    int mid=(l+r)/2;    doit(l,mid,e*2),doit(mid+1,r,e*2+1);    if(l2<=mid)delt(l,mid,e*2,l2);        else delt(mid+1,r,e*2+1,l2);    merge(e,e*2,e*2+1);}void join(int l,int r,int e,int l2){    if(l==r)    {        if(OK(b[0][e].a,b[1][e].a))b[0][e].s=1;        if(b[0][1].i>=X1)b[0][e].i=-INF,b[0][e].i1=0;        if(b[1][1].i>=Y1)b[1][e].i=-INF,b[1][e].i1=0;        return;    }    int mid=(l+r)/2;    doit(l,mid,e*2),doit(mid+1,r,e*2+1);    if(l2<=mid)join(l,mid,e*2,l2);        else join(mid+1,r,e*2+1,l2);    merge(e,e*2,e*2+1);}int find(int l,int r,int e,int l1,int r1){    if(l==l1&&r==r1)return b[0][e].s;    int mid=(l+r)/2;    doit(l,mid,e*2),doit(mid+1,r,e*2+1);    if(r1<=mid)return find(l,mid,e*2,l1,r1);        else if(mid<l1)return find(mid+1,r,e*2+1,l1,r1);            else return find(l,mid,e*2,l1,mid)+find(mid+1,r,e*2+1,mid+1,r1);} int main(){    int q,w,_,l,r;    b[0][0].i=-INF,b[0][0].a=-INF;      fo(i,1,3*N-10)fo(j,0,1)b[j][i]=b[0][0];    read(_);    while(_--)    {        read(n);        read(X1),read(Y1),read(X2),read(Y2);        fo(i,1,n)read(a[i].x),read(a[i].y);        build(1,n,1);        read(m);        fo(i,1,m)        {            read(q),read(l),read(r);            if(q==1)            {                read(w);                change(1,n,1,0,l,r,w);                while(b[0][1].a>X2)delt(1,n,1,b[0][1].a1);                while(b[0][1].i>=X1)                    join(1,n,1,b[0][1].i1);            }            if(q==2)            {                read(w);                change(1,n,1,1,l,r,w);                while(b[1][1].a>Y2)delt(1,n,1,b[1][1].a1);                while(b[1][1].i>=Y1)join(1,n,1,b[1][1].i1);            }            if(q==3)printf("%d\n",find(1,n,1,l,r));        }        fo(i,1,3*n)fo(j,0,1)b[j][i]=b[0][0];    }    return 0;}