线段树扫描线模板

#include <cstdio>//hdu1542 求平面矩形面积交
#include <algorithm>
#include <vector>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
const int N=1222;
using namespace std;
int n,cnt;
struct node
{
    double l,r,h;
    int f;
    node(double _l=0,double _r=0,double _h=0,int _f=0):l(_l),r(_r),h(_h),f(_f){};
}a[N*2];
bool cmp(node a,node b){return a.h<b.h;}
double sum[N<<2];
int flag[N<<2];
vector<double> vt;
int Getid(double x) {return lower_bound(vt.begin(),vt.end(),x)-vt.begin()+1;}
void Pushup(int l,int r,int rt)
{
    if(flag[rt]) sum[rt]=vt[r]-vt[l-1];
    else if(l==r) sum[rt]=0;
    else sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}
void Update(int L,int R,int x,int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        flag[rt]+=x;
        Pushup(l,r,rt);
        return;
    }
    int m=l+r>>1;
    if(L<=m) Update(L,R,x,lson);
    if(m<R) Update(L,R,x,rson);
    Pushup(l,r,rt);
}
int main()
{
    int i,cas=0;
    while(~scanf("%d",&n)&&n)
    {
        vt.clear();
        cnt=0;
        for(i=1;i<=n;i++)
        {
            double x1,x2,y1,y2;
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            vt.push_back(x1);vt.push_back(x2);
            a[++cnt]=node(x1,x2,y1,1);
            a[++cnt]=node(x1,x2,y2,-1);
        }
        sort(vt.begin(),vt.end());
        vt.erase(unique(vt.begin(),vt.end()),vt.end());
        sort(a+1,a+1+cnt,cmp);
        double ans=0;
        for(i=1;i<=cnt;i++)
        {
            int l=Getid(a[i].l);
            int r=Getid(a[i].r)-1;
            Update(l,r,a[i].f,1,(int)vt.size(),1);
            ans+=sum[1]*(a[i+1].h-a[i].h);
        }
        printf("Test case #%d\n",++cas);
        printf("Total explored area: %.2f\n\n",ans);
    }
}
/*********************************************/
#include <cstdio>//hdu1828求矩形周长交
#include <algorithm>
#include <vector>
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
const int N=100010;
using namespace std;
int n,cnt,sum[N<<2],flag[N<<2],num[N<<2],nl[N<<2],nr[N<<2];
struct node
{
    int l,r,h,f;
    node(int _l=0,int _r=0,int _h=0,int _f=0):l(_l),r(_r),h(_h),f(_f){}
}a[N*2];
bool Cmp(node a,node b)
{
    if(a.h!=b.h) return a.h<b.h;
    return a.f>b.f;
}
vector<int> vt;
int Getid(int x) {return lower_bound(vt.begin(),vt.end(),x)-vt.begin()+1;}
void pushup(int l,int r,int rt)
{
    if(flag[rt])
    {
        sum[rt]=vt[r]-vt[l-1];
        num[rt]=nl[rt]=nr[rt]=1;
    }
    else if(l==r) sum[rt]=num[rt]=nl[rt]=nr[rt]=0;
    else
    {
        sum[rt]=sum[rt<<1]+sum[rt<<1|1];
        nl[rt]=nl[rt<<1];
        nr[rt]=nr[rt<<1|1];
        num[rt]=num[rt<<1]+num[rt<<1|1]-nl[rt<<1|1]*nr[rt<<1];
    }
}
void Update(int L,int R,int x,int l,int r,int rt)
{
    if(L<=l&&r<=R)
    {
        flag[rt]+=x;
        Pushup(l,r,rt);
        return;
    }
    int m=l+r>>1;
    if(L<=m) Update(L,R,x,lson);
    if(m<R) Update(L,R,x,rson);
    Pushup(l,r,rt);
}
int main()
{
    int i;
    while(~scanf("%d",&n))
    {
        vt.clear();
        cnt=0;
        int x1,x2,y1,y2;
        for(i=1;i<=n;i++)
        {
            scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
            vt.push_back(x1),vt.push_back(x2);
            a[++cnt]=node(x1,x2,y1,1);
            a[++cnt]=node(x1,x2,y2,-1);
        }
        sort(vt.begin(),vt.end());
        vt.erase(unique(vt.begin(),vt.end()),vt.end());
        sort(a+1,a+1+cnt,Cmp);
        int ans=0,last=0;
        for(i=1;i<=cnt;i++)
        {
            int l=Getid(a[i].l);
            int r=Getid(a[i].r)-1;
            Update(l,r,a[i].f,1,(int)vt.size(),1);
            ans+=abs(sum[1]-last)+num[1]*(a[i+1].h-a[i].h)*2;
            last=sum[1];
        }
        printf("%d\n",ans);
    }
}