hdu1255 矩阵的交 线段树+扫描线

/*
    不是叶子节点 ，且cnt=1.注意这里，cnt=1确切的意义是什么，
    应该是，可以确定，这个区间被完全覆盖了1次，
    而有没有被完全覆盖两次或以上则不知道无法确定，那么怎么怎么办了，
    只要加上t[lch].s + t[rch].s  即，看看左右孩子区间被覆盖了一次或以上的长度，
    那么叠加在双亲上就是双亲被覆盖两次或以上的长度
*/

#include<stdio.h>
#include<string.h>
#include<algorithm>
using namespace std;
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define maxn 1005
struct segment
{
    double l,r,h;
    int f;
}s[maxn*];
struct node
{
    int cnt;
    double len1;//记录一次及以上的
    double len2;//记录二次及以上的
}tree[maxn*];
double mark[maxn<<];
bool cmp(segment a,segment b)
{
    return a.h<b.h;
}
void build(int l,int r,int rt)
{
    tree[rt].cnt=;
    tree[rt].len1=tree[rt].len2=;
    if(l==r)
    {
        return ;
    }
    int m=(l+r)/;
    build(lson);
    build(rson);
}
/*
    重点。
    若该节点的cnt >= 2,说明被至少两条线段覆盖，那么len1=len2=区间长度。
    若该节点的cnt == 1，说明该区间被一条线段覆盖，len1=区间长度，只要左右节点的len1有值，
    那么那些长度一定是至少被覆盖两次的，因此len2为左右节点的len1之和。
    若该节点的cnt = 0，说明没被完全覆盖，直接用其左右节点更新。
    还要注意特判叶子节点。
*/
void getlen(int l,int r,int rt)
{
    if(tree[rt].cnt>=)//>=2时很好理解，就是覆盖2次或以上，直接减一减
    {
        tree[rt].len1=mark[r+]-mark[l];
        tree[rt].len2=mark[r+]-mark[l];
    }
    else if(tree[rt].cnt==)//此处需要注意
                            //由于cnt大于0,所以覆盖一次或以上的部分直接得到，而由于完全覆盖一次
                            //可能其他部分有覆盖，此时就可以根据左右孩子来，由于len1记录覆盖一次或
                            //以上的，因为已经完全覆盖了一次，只要加上左右孩子的一次或以上覆盖的值
                            //那么就是2次或以上覆盖的值。
    {
        tree[rt].len1=mark[r+]-mark[l];//因为覆盖一次，len1还要继续加
        if(l == r)//由于覆盖一次，又只有一个l==r，所以len2为0
        {
            tree[rt].len2=;
        }
        else
        {
            tree[rt].len2=tree[rt<<].len1+tree[rt<<|].len1;
        }
    }
    else if(tree[rt].cnt==)//为完全覆盖，根据孩子来。此处与上面等于1出相似，可以发现，如果为完全覆盖
                            //就可以根据左右孩子得到覆盖的，那cnt等于1时，要求覆盖2次的时候就也可以
                            //根据孩子来。
    {
        if(l == r)
            tree[rt].len1=tree[rt].len2=;
        else
        {
            tree[rt].len1=tree[rt<<].len1+tree[rt<<|].len1;
            tree[rt].len2=tree[rt<<].len2+tree[rt<<|].len2;
        }
    }
}
void updata(int L,int R,int c,int l,int r,int rt)
{
    if(l>=L&&R>=r)
    {
        tree[rt].cnt+=c;
        getlen(l,r,rt);
        return ;
    }
    int m=(l+r)/;
    if(m>=L)
        updata(L,R,c,lson);
    if(R>m)
        updata(L,R,c,rson);
    getlen(l,r,rt);
}
int find(double val,int x,int y)
{
    int l=x,r=y,m;
    while(l<=r)
    {
        m=(l+r)/;
        if(mark[m]==val)
            return m;
        else if(mark[m]>val)
            r=m-;
        else l=m+;
    }
    return -;
}
int main()
{
    int t,n,i;
    double x1,x2,y1,y2;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%d",&n);
        int m=;
        for(i=;i<n;i++)
        {
            scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
            s[m].l=x1;s[m].r=x2;s[m].h=y1;s[m].f=;mark[m++]=x1;
            s[m].l=x1;s[m].r=x2;s[m].h=y2;s[m].f=-;mark[m++]=x2;
        }
        sort(s,s+m,cmp);
        sort(mark,mark+m);
        int k=;
        for(i=;i<m;i++)
        {
            if(mark[i]!=mark[i-])
                mark[k++]=mark[i];
        }
        /*
        for(i=0;i<k;i++)
            printf("%.2lf ",mark[i]);
        printf("\n");
        */
        build(,k-,);
        double ans=;
        for(i=;i<m;i++)
        {
            int ll=find(s[i].l,,k-);
            int rr=find(s[i].r,,k-)-;
            updata(ll,rr,s[i].f,,k-,);
            ans+=(s[i+].h-s[i].h)*tree[].len2;
        }
        printf("%.2lf\n",ans);
    }
}