17/10/12训练日记

死乞白赖的a了一个题线段树这个地方是真的不好理解，尤其是牵扯到离散化

给定平面上若干矩形,求出被这些矩形覆盖过至少两次的区域的面积.

 

Input

输入数据的第一行是一个正整数T(1<=T<=100),代表测试数据的数量.每个测试数据的第一行是一个正整数N(1<=N<=1000),代表矩形的数量,然后是N行数据,每一行包含四个浮点数,代表平面上的一个矩形的左上角坐标和右下角坐标,矩形的上下边和X轴平行,左右边和Y轴平行.坐标的范围从0到100000. 注意:本题的输入数据较多,推荐使用scanf读入数据.

 

Output

对于每组测试数据,请计算出被这些矩形覆盖过至少两次的区域的面积.结果保留两位小数.

 

Sample Input

251 1 4 21 3 3 72 1.5 5 4.53.5 1.25 7.5 46 3 10 730 0 1 11 0 2 12 0 3 1

 

Sample Output

7.630.00
典型的线段树求矩阵面积交的题，但还是觉得题目有点坑，先离散化然后枚举x坐标。此题用lazy不合适，必须跟新至单元，
因为每个节点既受父节点影响，也受子节点影响。
#include<iostream>
#include<stdio.h>#include<algorithm>#include<string.h>#include<iostream>using namespace std;const int MAXN=2020;struct Node{    int l,r;    int c;    double lf,rf;    double cnt;    double more;}segTree[MAXN*3];struct Line{    double x,y1,y2;    double f;}line[MAXN];double y[MAXN];bool cmp(Line a,Line b){    return a.x<b.x;}void Build(int i,int l,int r){    segTree[i].l=l;    segTree[i].r=r;    segTree[i].cnt=0;    segTree[i].more=0;    segTree[i].lf=y[l];    segTree[i].rf=y[r];    if(l+1==r)return;    int mid=(l+r)/2;    Build(i*2,l,mid);    Build(i*2+1,mid,r);}void calen(int i){    if(segTree[i].c>=2)    {        segTree[i].more=segTree[i].cnt=segTree[i].rf-segTree[i].lf;        return;    }    else if(segTree[i].c==1)    {        segTree[i].cnt=segTree[i].rf-segTree[i].lf;        if(segTree[i].l+1==segTree[i].r)segTree[i].more=0;        else segTree[i].more=segTree[i<<1].cnt+segTree[(i<<1)|1].cnt;    }    else    {        if(segTree[i].l+1==segTree[i].r)        {            segTree[i].cnt=segTree[i].more=0;        }        else        {            segTree[i].cnt=segTree[i<<1].cnt+segTree[(i<<1)|1].cnt;            segTree[i].more=segTree[i<<1].more+segTree[(i<<1)|1].more;        }    }}void update(int i,Line e){    if(e.y1==segTree[i].lf&&segTree[i].rf==e.y2)    {        segTree[i].c+=e.f;        calen(i);        return;    }    if(e.y2<=segTree[i<<1].rf) update(i<<1,e);    else if(e.y1>=segTree[(i<<1)|1].lf) update((i<<1)|1,e);    else    {        Line temp=e;        temp.y2=segTree[i<<1].rf;        update(i<<1,temp);        temp=e;        temp.y1=segTree[(i<<1)|1].lf;        update((i<<1)|1,temp);    }    calen(i);}int main(){    int T;    int n;    double x1,y1,x2,y2;    scanf("%d",&T);    while(T--)    {        scanf("%d",&n);        int t=1;        for(int i=1;i<=n;i++)        {            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);            line[t].x=x1;            line[t].y1=y1;            line[t].y2=y2;            line[t].x=x1;            line[t].f=1;            y[t]=y1;            t++;            line[t].y1=y1;            line[t].y2=y2;            line[t].x=x2;            line[t].f=-1;            y[t]=y2;            t++;        }        sort(line+1,line+t,cmp);        sort(y+1,y+t);        Build(1,1,t-1);        update(1,line[1]);        double ans=0;        for(int i=2;i<t;i++)        {            ans+=segTree[1].more*(line[i].x-line[i-1].x);            update(1,line[i]);        }        printf("%.2lf\n",ans);    }    return 0;}