HDU - 1828(矩形周长并 扫描线算法 )

主题算法思想参考notonlysuccess算法；

自己写了一遍，有很多地方比较巧妙；

线段树维护的值为len（扫到当前线段时的所有覆盖的线段总长），sumseg（扫到当前线段时的所有覆盖的线段的端点总数目） 需另外两个辅助参数lbd，rbd来计算；

本代码风格并未采用下推标记线段树写法，更简洁；因为只使用len【1】和sumseg【1】；所以只向上维护即可；

    #define lson l,m,rt<<1
    #define rson m+1,r,rt<<1|1

    const int maxn = 31111;
    int sumseg[maxn<<2],col[maxn<<2],len[maxn<<2],lbd[maxn<<2],rbd[maxn<<2];

    void build(int l,int r,int rt){
    sumseg[rt]=len[rt]=col[rt]=lbd[rt]=rbd[rt]=0;
    if(l==r) return ;
    int m=(l+r)>>1;
    build(lson);
    build(rson);
    }

    void pushup(int l,int r,int rt){
    if(col[rt]){
         sumseg[rt]=2;
         len[rt]=r-l+1;
         lbd[rt]=l; rbd[rt]=r;
    }
    else if(l==r){
            len[rt] = sumseg[rt]=0;
    }
    else{
         len[rt] = len[rt<<1] + len[rt<<1|1];
         sumseg[rt] = sumseg[rt<<1] +sumseg[rt<<1|1];
         if(sumseg[rt<<1] && sumseg[rt<<1|1] && rbd[rt<<1]+1==lbd[rt<<1|1])
            sumseg[rt]-=2;
         lbd[rt]=lbd[rt<<1]; rbd[rt]=rbd[rt<<1|1];
    }
    }

    int ql,qr,val;
    void update(int l,int r,int rt){
    if(ql<=l && r<=qr){
         col[rt]+=val;
         pushup(l,r,rt);
         return ;
    }
    int m=(l + r)>>1;
    if(ql<=m) update(lson);
    if(qr> m) update(rson);
    pushup(l,r,rt);
    }

    #include <cstdio>
    #include <vector>
    #include <cstdlib>
    #include <cstring>
    #include <iostream>
    #include <algorithm>
    using namespace std;

    const int add = 10002;
    int L=1,R=25000,n;

    struct segment {
    int l,r,x,flag;
    segment(int a,int b,int c,int d):l(a),r(b),x(c),flag(d){}
    bool operator < (const segment& rhs) const {
    return x < rhs.x || x == rhs.x && flag > rhs.flag;
    }
    };

    vector<segment> a;
    int main()
    {
       while(scanf("%d",&n)==1){
           a.clear();
           for(int i=1;i<=n;i++){
               int aa,b,c,d;
               scanf("%d %d %d %d",&aa,&b,&c,&d);
               aa+=add; b+=add; c+=add; d+=add;
               a.push_back(segment(b,d,aa,1));
               a.push_back(segment(b,d,c,-1));
           }
           sort(a.begin(),a.end());

           build(L,R,1);
           int ans=0,last=0;
           for(int i=0;i<a.size();i++){
               segment& seg = a[i];
               ql = seg.l; qr = seg.r-1; val = seg.flag;
               update(L,R,1);
               ans += abs(len[1] - last);
               last = len[1];
               ans += sumseg[1]*(i==a.size()-1 ? 0: a[i+1].x-a[i].x);
           }
           printf("%d\n",ans);
       }
       return 0;
    }