线段树 面积并 模板题

hdu 1542

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=1542

面积并不需要pushdown操作，因为不需要查询，但要pushup向上更新

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
#define maxn 100030
#define lson step<<1
#define rson step<<1|1
#define getMID (l+r)>>1
using namespace std;

struct Tree
{
    int l,r,date;
    double tsum;
}tree[maxn*3];

struct Node
{
    double l,r,h;
    int date;//记录出边，和入边，分别为1和-1 
}s[maxn];
int lazy[maxn*3],cnt;//用lazy数组来记录某段区间是否要纳入长度进行计算 
double pos[maxn];//离散数组 
bool cmp(Node p,Node q)
{
    return p.h<q.h;
}
void getlen(int root)
{
    if (lazy[root]>0)
     tree[root].tsum=pos[tree[root].r+1]-pos[tree[root].l];//如果lazy数组为正数，更新长度 
    else if (tree[root].l==tree[root].r)//更新到叶子节点，置零 
     tree[root].tsum=0;
    else tree[root].tsum=tree[root<<1].tsum+tree[root<<1|1].tsum; 
}
void build(int root,int l,int r)
{
    tree[root].l=l;
    tree[root].r=r;
    tree[root].tsum=0;
    if (l==r) return;
    int mid=(l+r)>>1;
    build(root<<1,l,mid);
    build(root<<1|1,mid+1,r);
}
void update(int root,int left,int right,int l,int r,int val)
{
    if (left==l && right==r)
    {
        lazy[root]+=val;
        getlen(root);
        return;
    }
    int mid=(left+right)>>1;
    if (r<=mid) update(root<<1,left,mid,l,r,val);
    else if (l>mid) update(root<<1|1,mid+1,right,l,r,val);
    else
    {
       update(root<<1,left,mid,l,mid,val);
       update(root<<1|1,mid+1,right,mid+1,r,val);
    }
    getlen(root);
}
void getline(double x1,double x2,double y1,int date) 
{
    s[cnt].l=x1;
    s[cnt].r=x2;
    s[cnt].h=y1;
    s[cnt].date=date;
    cnt++;
}
int main()
{
    int n,num,cas=0;
    while (scanf("%d",&n)!=EOF && n)
    {
        cnt=0;
        num=0;
        memset(lazy,0,sizeof(lazy));
        for (int i=0;i<n;i++)
        {
            double x1,x2,y1,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            getline(x1,x2,y1,1);
            getline(x1,x2,y2,-1);
            pos[num++]=x1;//离散化 
            pos[num++]=x2;
        }

        sort(s,s+cnt,cmp);//按高度y排好序 
        sort(pos,pos+num);
        int tem=1;
        for (int i=1;i<num;i++)//去重 
        {
            if (pos[i]!=pos[i-1])
             pos[tem++]=pos[i];
        }
        /*int tem;去重 
		int len=unique(pos,pos+num)-pos;
		tem=len; */
        build(1,0,tem-1);
        double res=0;
        for (int i=0;i<cnt-1;i++)//注意是cnt-1 
        {
            int l=lower_bound(pos,pos+tem,s[i].l)-pos;
            int r=lower_bound(pos,pos+tem,s[i].r)-pos-1;//由于是以节点来维护线段长度，可让r减去1 
            update(1,0,tem-1,l,r,s[i].date);
            res+=(s[i+1].h-s[i].h)*tree[1].tsum;//逐次计算每次扫描时的面积，累加 
            //注意，每次更新后就应计算面积 
        }
        printf("Test case #%d\n",++cas);
        printf("Total explored area: %.2f\n\n",res);
    }
}