hdu3642(三维转二维扫描线)

 比较有意思的一题，刚开始做这题的时候想到的low方法是扫描面，但这种肯定不对。。。想了一段时间，突然发现。。。枚举每个z点，就可以转化为二维上的问题。

这题还有一点需要注意，不能简单的枚举每个z点，算每个面上的重叠面积，这点是错的，具体操作请看代码

#include<iostream>
#include<stdio.h>
#include<cstring>
#include<algorithm>
#define MAX 1200
using namespace std;
typedef struct
{
    int x,y1,y2,z1,z2,f;
    void set(int xx,int yy1,int yy2,int zz1,int zz2,int ff)
    {
        x=xx;
        y1=yy1;
        y2=yy2;
        z1=zz1;
        z2=zz2;
        f=ff;
    }
}Line;
Line line[2*MAX];
Line temp[MAX*2];
bool cmp(const Line &a,const Line &b)
{
    return a.x<b.x;
}
int y[2*MAX],z[2*MAX];
typedef struct
{
    int left,right,realleft,realright;
    int one,two,three,cover;
    void set(int l,int r)
    {
        left=l;
        right=r;
        realleft=y[l];
        realright=y[r];
        three=two=one=cover=0;
    }
}P;
P p[MAX<<2];
void build(int k,int l,int r)
{
    p[k].set(l,r);
    if(r-l==1)
        return;
    int mid=(r+l)>>1;
    build(k<<1,l,mid);
    build(k<<1|1,mid,r);
}
void push_up(int k)
{
    if( p[k].cover>2)
    {
        p[k].three= p[k].realright- p[k].realleft;
        p[k].two= p[k].one=0;
    }
    if( p[k].cover==2)
    {
        p[k].one=0;
        if( p[k].right- p[k].left==1)
        {
            p[k].three=0;
            p[k].two= p[k].realright- p[k].realleft;
        }
        else
        {
            int lt=k<<1,rt=lt+1;
            p[k].three = p[lt].one+ p[lt].two+ p[lt].three+ p[rt].one+ p[rt].two+ p[rt].three;
            p[k].two= p[k].realright- p[k].realleft -p[k].three;
        }
    }
    if( p[k].cover==1)
    {
        if( p[k].right- p[k].left==1)
        {
            p[k].three= p[k].two=0;
            p[k].one= p[k].realright- p[k].realleft;
        }
        else
        {
            int lt=k<<1,rt=lt+1;
            p[k].two= p[lt].one+ p[rt].one;
            p[k].three= p[lt].two+ p[rt].two+ p[lt].three+ p[rt].three;
            p[k].one= p[k].realright- p[k].realleft -p[k].two- p[k].three;
        }
    }
    if( p[k].cover==0)
    {
        if( p[k].right- p[k].left==1)
        {
            p[k].one= p[k].two= p[k].three=0;
        }
        else
        {
            int lt=k<<1,rt=lt+1;
            p[k].one= p[lt].one+ p[rt].one;
            p[k].two= p[lt].two+ p[rt].two;
            p[k].three= p[lt].three+ p[rt].three;
        }
    }
}
void update(int k,Line a)
{
    if(a.y1==p[k].realleft&& a.y2==p[k].realright)
    {
        p[k].cover+=a.f;
        push_up(k);
        return;
    }
    if( p[k<<1].realright>= a.y2)
        update(k<<1,a);
    else
        if( p[k<<1|1].realleft<=a.y1)
        update(k<<1|1,a);
    else
    {
        Line m=a;
        m.y2= p[k<<1].realright;
        update(k<<1,m);
        m=a;
        m.y1= p[k<<1|1].realleft;
        update(k<<1|1,m);
    }
    push_up(k);
}
int main()
{
    long long sum,avg;
    int t,x1,x2,y1,y2,z1,z2,num=0;
    scanf("%d",&t);
    while(t--)
    {
        num++;
        avg=0;
        int n,t=0,t1=0;
        scanf("%d",&n);
        for(int i=0;i<n;i++)
        {
            scanf("%d%d%d%d%d%d",&x1,&y1,&z1,&x2,&y2,&z2);
            line[t++].set(x1,y1,y2,z1,z2,1);
            line[t++].set(x2,y1,y2,z1,z2,-1);
            y[t1]=y1;
            z[t1++]=z1;
            y[t1]=y2;
            z[t1++]=z2;
        }
        sort(line,line+t,cmp);
        sort(y,y+t1);
        sort(z,z+t1);
        int len1=unique(y,y+t1)-y;
        int len2=unique(z,z+t1)-z;
        build(1,0,len1-1);
        for(int i=0;i<len2-1;i++)
        {
            sum=0;
            int m=0;
            for(int j=0;j<t;j++)
                if( line[j].z1<=z[i]&& line[j].z2>z[i])
                    temp[m++]=line[j];
            if(m==0)
                continue;
            update(1,temp[0]);
            for(int j=1;j<m;j++)
            {//cout<<p[1].three<<endl;
                sum+=(long long) p[1].three*(temp[j].x- temp[j-1].x);
                update(1,temp[j]);
            }
            avg+=sum*(z[i+1]-z[i]);//cout<<avg<<endl;
        }
        cout<<"Case "<<num<<": "<<avg<<endl;
    }
}