汕头市队赛 SRM1X T2 ——扫描线

绵津见-终 SRM 13

背景

“西瓜也是可以种在海上的！”——绵津见

然而种在海上的西瓜最需要防范的，是时不时会涌向瓜田的阵阵海浪。

幸好，身为海神的绵津见可以释放魔法“水平如镜”来阻止海浪拍打西瓜。

然而，当西瓜一个接一个成熟之时，它们就要离开瓜田，飘向遥远的彼岸。绵津见的魔法无法保护离开瓜田的西瓜们，但至少，也得知道西瓜们遭遇了多大的风浪啊。

描述

我们用一个坐标系来描述大海，绵津见的瓜田位于x轴下方，每当有一个西瓜成熟时，它会从x轴上一点出发，沿一条平行y轴的直线往y轴正方向前进。

某个时刻，海上会有一个海浪生成。每个海浪都是一条平行于x轴的线段，并且会往y轴负方向前进。当它达到x轴下方时，它会受到“水平如镜”的影响而消失。

所有西瓜、海浪每个时刻会前进一个单位距离，西瓜与西瓜之间，海浪与海浪之间互不影响，相互重叠的情况也是允许发生的。

当一个西瓜与海浪重叠时，认定为该西瓜遭遇该海浪的拍打，西瓜与海浪的运动不受此次拍打的影响，原路前进。

每个时刻的流程如下：

①在该时刻成熟的西瓜被摆放在x轴上，该时刻生成的海浪出现。

②进行一次判定操作：如果有西瓜和某海浪上的一点重合，判定为发生了一次拍打。

③所有西瓜前进一个单位距离，即y坐标加1。

④进行一次判定操作。

⑤所有海浪前进一个单位距离，即y坐标减1。y坐标为-1的海浪消失。

你需要回答每个西瓜、每个海浪共涉及多少次拍打。

输入格式

第一行两个整数n,m表示西瓜的数量和海浪的数量。

接下来n行，每行两个整数ti,xi描述一个西瓜，分别表示西瓜成熟的时间ti,成熟时出现在(xi,0)位置。

接下来m行，每行四个整数Ti,li,ri,yi描述一个海浪，分别表示海浪生成的时间ti,生成时两个端点为(li,yi)，(ri,yi)。

输出格式

输出第一行n个数，第i个数表示第i个西瓜被拍打的次数。

输出第二行m个数，第i个数表示第i个海浪拍打到的西瓜数。

样例输入	样例输出
5 3 2 3 2 5 2 6 4 6 4 7 1 1 6 0 2 1 6 0 2 6 7 3	1 1 2 1 1 0 3 3

数据范围与约定

对于初：  ,,

对于续：  ,,

对于终：  ,,

样例解释

​不解释

————————————————————————————————————

这道题分析之后 我们发现西瓜和海浪相碰需要满足

设西瓜出发时间T1 位置x 海浪T2 l r 高度 y

l<=x<=r  T2-y<=T1<=T2+y

很明显这是符合扫描线的性质 所以我们可以将T排序 将 l x r 离散化一波

将询问拆为两个 相减得到答案

按T扫一波然后进行区间求和 因为点和矩阵相互影响是对称的

所以询问和求和是相反的 这样就可以直接维护一波答案了

代码简短但是略复杂 还是要花时间慢慢体会的

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int M=1e5+;
int read(){
    int ans=,f=,c=getchar();
    while(c<''||c>''){if(c=='-') f=-; c=getchar();}
    while(c>=''&&c<=''){ans=ans*+(c-''); c=getchar();}
    return ans*f;
}
int n,m;
int xs[M*],xp;
int qp,ep;
int s1[*M],s2[*M],ans1[M],ans2[M];
int lowbit(int x){return x&-x;}
void add(int x,int v,int s[]){
    while(x<=xp){
        s[x]+=v;
        x+=lowbit(x);
    }
}
int query(int x,int s[]){
    int ans=;
    while(x){ans+=s[x]; x-=lowbit(x);}
    return ans;
}
struct Q{
    int l,r,y,id,s;
    bool operator <(const Q& h)const{return y<h.y;}
    void cal(){
        ans1[id]+=(query(r,s1)-query(l-,s1))*s;
        add(l,-s,s2);
        add(r+,s,s2);
    }
}e[*M];
struct pos{
    int x,T,id;
    bool operator <(const pos& h)const{return T<h.T;}
    void cal(){
        add(x,,s1);
        ans2[id]=query(x,s2);
    }
}q[*M];
void $(int&x){x=lower_bound(xs,xs+xp,x)-xs+;}
int main()
{
    int l,r,k,x;
    n=read(); m=read();
    for(int i=;i<n;i++){
        k=read(); x=read();
        q[qp++]=(pos){xs[xp++]=x,k,i};
    }
    for(int i=;i<m;i++){
        k=read(); l=read(); r=read(); x=read();
        e[ep++]=(Q){l,r,k-x-,i,-};
        e[ep++]=(Q){l,r,k+x,i,};
        xs[xp++]=l;
        xs[xp++]=r;
    }
    sort(xs,xs+xp);
    for(int i=;i<qp;++i)$(q[i].x);
    for(int i=;i<ep;++i)$(e[i].l),$(e[i].r);
    sort(q,q+qp);
    sort(e,e+ep);
    for(int i=,j=;i<ep;i++){
        while(j<qp&&q[j].T<=e[i].y) q[j++].cal();
        e[i].cal();
    }
    for(int i=;i<n;i++) printf("%d ",ans2[i]); printf("\n");
    for(int i=;i<m;i++) printf("%d ",ans1[i]); printf("\n");
    return ;
}