【线段树专题】poj1177

这回变成了矩形周长并，

思想和矩形面积并是一样的。

但是有一点区别的是周长并需要判断当前子节点是不是已经被覆盖了。如果被覆盖了那么ans不需要改动，

如果当前节点没有 （被覆盖的上级节点）那么这个节点的sum直如果变化的话就要加到答案里去了。

还是老问题，交c++WA交g++AC。。。。

对了。。上下方向和左右方向都要做一遍

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<algorithm>
#include<vector>
using namespace std;
int ans=0,n,tail,tmp_tail,s_tot,q[20100],tmp_q[20100],tot;
struct rec
{
    int l,r,lch,rch,sum,da,biao;
}tree[20100];

struct rec2
{
    int h,l,r,biao;
    int h2,l2,r2,biao2;
}segment[20100];

int cmp(rec2 x,rec2 y) {return x.h<y.h;}

void clear()
{
    memset(tree,0,sizeof(tree));
    memset(q,0,sizeof(q));
    memset(tmp_q,0,sizeof(tmp_q));
    tail=0;
    tmp_tail=0;
    tot=0;
}

void maketree(int l,int r)
{
    int now=++tot;
    tree[now].l=l;
    tree[now].r=r;
    if (l==r) 
    {
        tree[now].da=q[l+1]-q[l];
        return;
    }
    int mid=(l+r)>>1;
    tree[now].lch=tot+1;maketree(l,mid);
    tree[now].rch=tot+1;maketree(mid+1,r);
    tree[now].da=tree[tree[now].lch].da+tree[tree[now].rch].da;
}
void Input()
{
    int x1,y1,x2,y2;
    for (int i=1;i<=n;i++)
    {
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        s_tot++;
        segment[s_tot].h=y1,segment[s_tot].l=x1,segment[s_tot].r=x2,segment[s_tot].biao=1;
        segment[s_tot].h2=x1,segment[s_tot].l2=y1,segment[s_tot].r2=y2,segment[s_tot].biao2=1;
        s_tot++;
        segment[s_tot].h=y2,segment[s_tot].l=x1,segment[s_tot].r=x2,segment[s_tot].biao=-1;
        segment[s_tot].h2=x2,segment[s_tot].l2=y1,segment[s_tot].r2=y2,segment[s_tot].biao2=-1;
    }
}
void pre_work()
{
    for (int i=1;i<=s_tot;i++)
    {
        tmp_q[++tmp_tail]=segment[i].l;
        tmp_q[++tmp_tail]=segment[i].r;
    }
    sort(tmp_q+1,tmp_q+tmp_tail+1);
    for (int i=1;i<=tmp_tail;i++)
        if (tmp_q[i-1]!=tmp_q[i])
            q[++tail]=tmp_q[i];
    maketree(1,tail-1);
    sort(segment+1,segment+s_tot+1,cmp);
}

void change(int x,int l,int r,int w,int f)
{
    if (tree[x].l>=l&&tree[x].r<=r)
    {
        int tmp=tree[x].sum;
        tree[x].biao+=w;
        if (tree[x].biao!=0) tree[x].sum=tree[x].da;
        else 
            if (tree[x].l!=tree[x].r) 
                tree[x].sum=tree[tree[x].lch].sum+tree[tree[x].rch].sum;
            else
                tree[x].sum=0;
        if (f==0) ans+=abs(tree[x].sum-tmp);
        return ;
    }
    int nowf=0;
    if (tree[x].biao>0) nowf=1;
    int mid=(tree[x].l+tree[x].r)>>1;
    if (l<=mid) change(tree[x].lch,l,r,w,f+tree[x].biao);
    if (r>mid) change(tree[x].rch,l,r,w,f+tree[x].biao);
    int tmp=tree[x].sum;
    if (tree[x].biao!=0) tree[x].sum=tree[x].da;
        else 
            if (tree[x].l!=tree[x].r) 
                tree[x].sum=tree[tree[x].lch].sum+tree[tree[x].rch].sum;
            else
                tree[x].sum=0;
}
void work()
{
    int l,r;
    for (int i=1;i<=s_tot;i++)
    {
        l=lower_bound(q+1,q+tail+1,segment[i].l)-q;
        r=lower_bound(q+1,q+tail+1,segment[i].r)-q;
        change(1,l,r-1,segment[i].biao,0);
    }    
}

void change_work()
{
    for (int i=1;i<=s_tot;i++)
    {
        segment[i].h=segment[i].h2;
        segment[i].l=segment[i].l2;
        segment[i].r=segment[i].r2;
        segment[i].biao=segment[i].biao2;
    }
}
int main()
{
    freopen("poj1177.in","r",stdin);
    freopen("poj1177.out","w",stdout);
    memset(segment,0,sizeof(segment));
    s_tot=0;
    scanf("%d",&n);
    clear();
    Input();
    pre_work();
    work();
    clear();
    change_work();
    pre_work();
    work();
    printf("%d",ans);
}