二维平面鱼群模拟-优快云博客

本文链接：https://blog.youkuaiyun.com/shiyongyang/article/details/78342183

题目大意

池塘里有一些鱼和一个渔网，池塘可以看成一个二维的平面，而渔网可以看成一个与坐标轴平行的矩形。每条鱼都被给予了一个标号，分别从1到n标号，n表示池塘里鱼的总数。鱼的游动可以概括为两个动作：
1 l r d : 表示标号在[l，r]这个区间内的鱼向x轴正方向游动了d个单位长度。
2 l r d：表示标号在[l，r]这个区间内的鱼向y轴正方向游动了d个单位长度。
还有询问操作：3 l r 表示询问现在标号在[l，r]这个区间内的鱼有多少在渔网内。

输入格式

第一行包含一个整数T，表示测试数据组数。
对于每组测试数据：
第一行包含一个整数n表示鱼的总数。
第二行包含四个整数x1,y1,x2,y2，表示渔网的左下角坐标和右上角坐标。
接下来n行每行两个整数xi,yi，表示标号为i的鱼初始时刻的坐标。
接下来一行包含一个整数m，表示后面的事件数目。
接下来的m行每行为以下三种类型的一种：
1 l r d : 表示标号在[l，r]这个区间内的鱼向x轴正方向游动了d个单位长度。
2 l r d：表示标号在[l，r]这个区间内的鱼向y轴正方向游动了d个单位长度。
3 l r : 表示询问现在标号在[l，r]这个区间内的鱼有多少在渔网内。

输出格式

对于每组数据的每个询问，输出一个整数表示对应的答案。

数据范围

对于30%的数据1≤n,m≤1000
对于100%的数据1≤T≤10,1≤n,m≤30000,1≤l≤r≤n. 1≤d≤10^9 ,x1≤x2,y1≤y2。保证任意时刻所有涉及的坐标值在[−10^9,10^9]范围内。

/*
hdu 5283
用鱼的编号为下标建立线段树，每一个节点维护当前区间内小于x1（y1）的最大x（y）值和小于x2（y2）点最大x（y）的值。
因为鱼的运动方向具有单调性，所以一个鱼一但出去了就再也回不来了 ，那么久删除这些节点，不断更行1号节点的值。注意x，y
各开一棵线段树。
*/
#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)//x,x<<1,x<<1|1
{
    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 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;//a代表小于x2的最大坐标 
            b[0][e].a1=b[1][e].a1=l;//a1代表坐标表示的鱼的编号 
            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;//i表示当前区间内坐标小于x1点最大坐标，i1表示编号 
            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 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 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+10,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;
}