JZOJ 4694 火神的鱼

来源：互联网发布：linux开机启动管理编辑：程序博客网时间：2024/05/01 02:06

题目大意

在池塘里用一张网（左下角x1,y1,右上角x2,y2）捕鱼，网和池塘都是一个矩形，池塘里有n条鱼，这些鱼会四处游动，火神会在鱼游动的间隙问你现在网里有多少条鱼。
鱼的游动可以概括为两个动作：
1 l r d : 表示标号在[l，r]这个区间内的鱼向x轴正方向游动了d个单位长度。
2 l r d：表示标号在[l，r]这个区间内的鱼向y轴正方向游动了d个单位长度。

多组数据。
T<=10,n,m<=30000,0<=d<=1e9,所有涉及的坐标在[−10^9,10^9]范围内。
时间限制 1s
空间限制 256M

解题思路

建两棵线段树，由于d是正数，所以从网里面钻出来的鱼不可能再进去，所以我们可以用线段树上区间[l,r]储存的信息是：
（以x轴上的线段树为例）
1、标号为[l,r]的鱼中有多少条；
2、标号为[l,r]的鱼中坐标小于x1的最大坐标；
3、标号为[l,r]的鱼中坐标在[x1,x2]范围内的最大坐标；

于是鱼的游动就可以用区间加减解决，加减之后就维护上述的第二、三项，及时把坐标脱离范围的鱼加入或移出网外就可以了。

#include<cstdio>#include<cstring>#include<algorithm>#define maxn 100006#define fr(i,a,b) for(i=a;i<=b;i++)using namespace std;typedef long long ll;const ll ding=2147483647;struct nod{    int x,y,c,s;} tr1[maxn*3],tr2[maxn*3];int T,i,n,m,x1,x2,y1,y2,z,l,r,d,s,a[maxn],b[maxn];void merge(int v){    tr1[v].s=tr1[v+v].s+tr1[v+v+1].s;    tr1[v].x=max(tr1[v+v].x,tr1[v+v+1].x);    tr1[v].y=max(tr1[v+v].y,tr1[v+v+1].y);    tr2[v].s=tr2[v+v].s+tr2[v+v+1].s;    tr2[v].x=max(tr2[v+v].x,tr2[v+v+1].x);    tr2[v].y=max(tr2[v+v].y,tr2[v+v+1].y);    tr1[v].c=tr2[v].c=0;    return;}void maket(int v,int st,int en){    if (st==en)    {        tr1[v].c=tr2[v].c=0;        tr1[v].s=tr2[v].s=0;        tr1[v].x=tr1[v].y=tr2[v].x=tr2[v].y=-ding;        if (a[st]<x1) tr1[v].x=a[st];        else if (a[st]<=x2) tr1[v].y=a[st];        if (b[st]<y1) tr2[v].x=b[st];        else if (b[st]<=y2) tr2[v].y=b[st];        if (a[st]>=x1 && a[st]<=x2 && b[st]>=y1 && b[st]<=y2)             tr1[v].s=1;        return;    }    int m=(st+en) >> 1;    maket(v+v,st,m);    maket(v+v+1,m+1,en);    merge(v);    return;}void update1(int v){    if (!tr1[v].c) return;    if (tr1[v+v].x>-ding) tr1[v+v].x+=tr1[v].c;    if (tr1[v+v].y>-ding) tr1[v+v].y+=tr1[v].c;    if (tr1[v+v+1].x>-ding) tr1[v+v+1].x+=tr1[v].c;    if (tr1[v+v+1].y>-ding) tr1[v+v+1].y+=tr1[v].c;    tr1[v+v].c+=tr1[v].c;    tr1[v+v+1].c+=tr1[v].c;    tr1[v].c=0;    return;}void update2(int v){    if (!tr2[v].c) return;    if (tr2[v+v].x>-ding) tr2[v+v].x+=tr2[v].c;    if (tr2[v+v].y>-ding) tr2[v+v].y+=tr2[v].c;    if (tr2[v+v+1].x>-ding) tr2[v+v+1].x+=tr2[v].c;    if (tr2[v+v+1].y>-ding) tr2[v+v+1].y+=tr2[v].c;    tr2[v+v].c+=tr2[v].c;    tr2[v+v+1].c+=tr2[v].c;    tr2[v].c=0;    return;}void modify1(int v,int st,int en,int l,int r,int x){    if (st==l && en==r)    {        if (tr1[v].x>-ding) tr1[v].x+=x;        if (tr1[v].y>-ding) tr1[v].y+=x;        tr1[v].c+=x;        return;    }    update1(v);    update2(v);    int m=(st+en) >> 1;    if (r<=m) modify1(v+v,st,m,l,r,x);    else if (l>m) modify1(v+v+1,m+1,en,l,r,x);    else    {        modify1(v+v,st,m,l,m,x);        modify1(v+v+1,m+1,en,m+1,r,x);    }    merge(v);    return;}void modify2(int v,int st,int en,int l,int r,int x){    if (st==l && en==r)    {        if (tr2[v].x>-ding) tr2[v].x+=x;        if (tr2[v].y>-ding) tr2[v].y+=x;        tr2[v].c+=x;        return;    }    update1(v);    update2(v);    int m=(st+en) >> 1;    if (r<=m) modify2(v+v,st,m,l,r,x);    else if (l>m) modify2(v+v+1,m+1,en,l,r,x);    else    {        modify2(v+v,st,m,l,m,x);        modify2(v+v+1,m+1,en,m+1,r,x);    }    merge(v);    return;}void add(int v,int st,int en,int x){    if (st==en)    {        if (x==1)        {            if (tr2[v].y>=y1 && tr2[v].y<=y2) tr1[v].s=1;            else tr1[v].s=0;            tr1[v].y=tr1[v].x,tr1[v].x=-ding;        } else        {            if (tr1[v].y>=x1 && tr1[v].y<=x2) tr1[v].s=1;            else tr1[v].s=0;            tr2[v].y=tr2[v].x,tr2[v].x=-ding;        }        return;    }    update1(v);    update2(v);    int m=(st+en) >> 1;    if (x==1)    {        if (tr1[v+v].x>=x1) add(v+v,st,m,x);        else add(v+v+1,m+1,en,x);    } else    {        if (tr2[v+v].x>=y1) add(v+v,st,m,x);        else add(v+v+1,m+1,en,x);    }    merge(v);    return;}void del(int v,int st,int en,int x){    if (st==en)    {        tr1[v].x=tr1[v].y=tr2[v].x=tr2[v].y=-ding;        tr1[v].s=0;        return;    }    update1(v);    update2(v);    int m=(st+en) >> 1;    if (x==1)    {        if (tr1[v+v].y>x2) del(v+v,st,m,x);        else del(v+v+1,m+1,en,x);    } else    {        if (tr2[v+v].y>y2) del(v+v,st,m,x);        else del(v+v+1,m+1,en,x);    }    merge(v);    return;}void findd(int v,int st,int en,int l,int r){    if (st==l && en==r)    {        s+=tr1[v].s;        return;    }    update1(v);    update2(v);    int m=(st+en) >> 1;    if (r<=m) findd(v+v,st,m,l,r);    else if (l>m) findd(v+v+1,m+1,en,l,r);    else    {        findd(v+v,st,m,l,m);        findd(v+v+1,m+1,en,m+1,r);    }    merge(v);    return;}int main(){    scanf("%d",&T);    while (T--)    {        scanf("%d",&n);        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);        fr(i,1,n) scanf("%d%d",&a[i],&b[i]);        maket(1,1,n);        scanf("%d",&m);        fr(i,1,m)        {            scanf("%d%d%d",&z,&l,&r);            if (z!=3)            {                scanf("%d",&d);                if (z==1)                {                    modify1(1,1,n,l,r,d);                    while (tr1[1].x>=x1) add(1,1,n,1);                    while (tr1[1].y>x2) del(1,1,n,1);                } else if (z==2)                {                    modify2(1,1,n,l,r,d);                    while (tr2[1].x>=y1) add(1,1,n,2);                    while (tr2[1].y>y2) del(1,1,n,2);                }            } else            {                s=0;                findd(1,1,n,l,r);                printf("%d\n",s);            }        }    }    return 0;}

0 0