HDU 3265 Posters

矩形面积并，一个拆成四个

#include<cstdio>
#include<cstring>
#include<cmath>
#include<map>
#include<algorithm>
using namespace std;

const long long maxn=+;
struct Seg
{
    long long x;
    long long Y1,Y2;//离散化之后的坐标
    long long flag;
}s[*maxn];
long long X1[*maxn],X2[*maxn],Y1[*maxn],Y2[*maxn];
map<long long ,long long>m;
long long M[*maxn];long long k;
long long n,tot,N;
long long sum,ans;

struct SegTree
{
    long long len;
    long long cover;
} segTree[maxn*];

bool cmp(const Seg&a,const Seg&b)
{
    return a.x<b.x;
}

void lsh()
{
    k=;
    m.clear();
    for(long long i=; i<N; i++)
    {
        if(m[Y1[i]]==) M[k++]=Y1[i],m[Y1[i]]=;
        if(m[Y2[i]]==) M[k++]=Y2[i],m[Y2[i]]=;
    }
    sort(M,M+k);
    m.clear();
    for(long long i=; i<k; i++) m[M[i]]=i;
}

void build(long long l,long long r,long long rt)
{
    segTree[rt].cover=;
    segTree[rt].len=;
    if(l==r) return;
    long long m=(l+r)/;
    build(l,m,*rt);
    build(m+,r,*rt+);
}

void pushUp(long long rt,long long l,long long r)
{
    if(segTree[rt].cover) segTree[rt].len=M[r]-M[l-];
    else segTree[rt].len=segTree[*rt].len+segTree[*rt+].len;
}

void update(long long info,long long L,long long R,long long l,long long r,long long rt)
{
    if(L<=l&&r<=R)
    {
        segTree[rt].cover=segTree[rt].cover+info;
        pushUp(rt,l,r);
        return;
    }

    long long m=(l+r)/;
    if(L<=m) update(info,L,R,l,m,*rt);
    if(R>m)  update(info,L,R,m+,r,*rt+);
    pushUp(rt,l,r);
}

int main()
{

    while(~scanf("%lld",&n))
    {
        if(n==) break;
        N=;
        for(long long i=; i<=n; i++)
        {
            long long a1,a2,a3,a4,b1,b2,b3,b4;
            scanf("%lld%lld%lld%lld%lld%lld%lld%lld",&a1,&b1,&a2,&b2,&a3,&b3,&a4,&b4);
            if(b3>b1&&a2>a1){
            X1[N]=a1,Y1[N]=b1;
            X2[N]=a2,Y2[N]=b3;
            N++;}

            if(b2>b4&&a2>a1){
            X1[N]=a1,Y1[N]=b4;
            X2[N]=a2,Y2[N]=b2;
            N++;}

            if(a3>a1&&b4>b3){
            X1[N]=a1,Y1[N]=b3;
            X2[N]=a3,Y2[N]=b4;
            N++;}

            if(a2>a4&&b4>b3){
            X1[N]=a4,Y1[N]=b3;
            X2[N]=a2,Y2[N]=b4;
            N++;}
        }
        if(N==)
        {
            printf("0\n");
            continue;
        }
        lsh();

        tot=;
        for(long long i=; i<N; i++)
        {
            s[tot].x=X1[i],s[tot].Y1=m[Y1[i]],s[tot].Y2=m[Y2[i]],s[tot].flag=,tot++;
            s[tot].x=X2[i],s[tot].Y1=m[Y1[i]],s[tot].Y2=m[Y2[i]],s[tot].flag=-,tot++;
        }
        sort(s,s+tot,cmp);

        ans=;

        build(,k,);
        update(s[].flag,s[].Y1+,s[].Y2,,k,);
        for(long long i=; i<tot; i++)
        {
            sum=segTree[].len;
            ans=ans+sum*(s[i].x-s[i-].x);
            update(s[i].flag,s[i].Y1+,s[i].Y2,,k,);
        }

        printf("%lld\n",ans);
    }
    return ;
}