hdu1542(矩形面积并)

链接：点击打开链接

题意 ：给出n个矩形的左下角和右上角的坐标，求矩形面积的并

代码：

#include <stdio.h>#include <string.h>#include <iostream>#include <algorithm>using namespace std;const int SIZE=505;int add[SIZE<<2];                               //add为区间标记，与懒惰标记类似double x[SIZE<<2],sum[SIZE<<2];struct node{    int ss;                                     //ss=1为下边，ss=-1为上边    double l,r,h;                               //分别记录线段的左端，右端和高度    node(){}    node(double a,double b,double c,int d):l(a),r(b),h(c),ss(d){}    friend bool operator<(node p,node q){       //因为是从下到上扫描因此，需要按照高度排序        return p.h<q.h;    }}s[SIZE];void pushup(int rt,int l,int r){    if(add[rt])                                 //这就是如何解决的线段重叠的办法,如果当前区间    sum[rt]=x[r+1]-x[l];                        //有标记用x数组更新，而不是左右儿子更新    else if(l==r)    sum[rt]=0;    else    sum[rt]=sum[rt<<1]+sum[rt<<1|1];            //与懒惰标记的区别在于标记不会消失(除非遇到上边)}                                               //这也就是线段不会重复相加减的关键void update(int L,int R,int c,int l,int r,int rt){    int m;    if(L<=l&&r<=R){        add[rt]+=c;                             //update与普通线段树基本相同        pushup(rt,l,r);        return;    }    m=(l+r)>>1;    if(L<=m)    update(L,R,c,l,m,rt<<1);    if(R>m)    update(L,R,c,m+1,r,rt<<1|1);    pushup(rt,l,r);}int main(){                                     //用的是从下到上扫描，也就是扫描线，又因为    int n,i,k,l,m,r,cas;                        //x坐标可能过多，因此需要对x轴进行离散化,从    double a,b,c,d,ans;                         //下到上扫描也就是遇到下边进行成段加1，遇到    cas=1;                                      //上边就是成段减1，所以问题就变成了线段树成    while(scanf("%d",&n)!=EOF&&n){              //段更新，但我们还需要解决线段重叠的问题，只        k=1,ans=m=0;                            //要修改pushup就能巧妙的解决这个问题        for(i=0;i<n;i++){            scanf("%lf%lf%lf%lf",&a,&b,&c,&d);            x[m]=a;            s[m++]=node(a,c,b,1);            x[m]=c;            s[m++]=node(a,c,d,-1);        }        sort(x,x+m);        sort(s,s+m);        for(i=1;i<m;i++)                        //对x坐标进行离散化，也可以用set        if(x[i]!=x[i-1])        x[k++]=x[i];        memset(add,0,sizeof(add));        memset(sum,0,sizeof(sum));        for(i=0;i<m-1;i++){            l=lower_bound(x,x+k,s[i].l)-x;            r=lower_bound(x,x+k,s[i].r)-x-1;    //找出更新的区间，这是左必右开区间，因此            update(l,r,s[i].ss,0,k-1,1);        //r需要减1，这也是与普通线段树的区别,普通            ans+=(sum[1]*(s[i+1].h-s[i].h));    //的线段树其实是更新的节点中保存的点,假如        }                                       //对[1,2]和[2,3]区间都加1，普通线段树中2        printf("Test case #%d\nTotal explored area: %.2lf\n\n",cas++,ans);    }                                           //那个点会变成2，而这个是表示的线段,因此    return 0;                                   //2应该为1，因此用左闭右开区间}