POJ 1151 Atlantis

需求：区间更新（覆盖次数，总长度），整体查询（总长度）。

因为是总体查询，只需要看len[1]即可，而且正因为如此不需要精细到每个叶子节点的len，只要更新到当前需要的区间，回溯到根即可。若是下一次有更底层的更新在这层可以滤过，不需要更新到底所以从整体上省去了lazy操作。

解法：离散化之后查询每次更新前整体总长度与上一次的纵坐标差相乘，累加即可。

G++交的话有小trick，printf("%.2f")而不是printf(".2lf")

#include <algorithm>
#include <iostream>
#include <fstream>
#include <sstream>
#include <iomanip>

#include <map>
#include <set>
#include <list>
#include <stack>
#include <queue>
#include <deque>
#include <vector>
#include <string>
#include <bitset>
#include <memory>
#include <complex>
#include <numeric>

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <time.h>
#include <ctype.h>
#include <locale.h>

using namespace std;

#pragma pack(4)

const double  eps = 1e-8;
const double   pi = acos(-1.0);
const int     inf = 0x7f7f7f7f;
const __int64 INF = 0x7f7f7f7f7f7f7f7fll;

#define loop(a,n)                            \
    for(int i=0;n>i;i++)                     \
        cout<<a[i]<<(i!=n-1?' ':'\n')
#define loop2(a,n,m)                         \
    for(int i=0;n>i;i++)                     \
        for(int j=0;m>j;j++)                 \
            cout<<a[i][j]<<(j!=m-1?' ':'\n')

#define   at(a,i) ((a)&(1<<(i)))
#define   nt(a,i) ((a)^(1<<(i)))
#define set1(a,i) ((a)|(1<<(i)))
#define set0(a,i) ((a)&(~(1<<(i))))

#define gret(a,b) (((a)-(b))>eps)
#define less(a,b) (((b)-(a))>eps)
#define greq(a,b) (((a)-(b))>-eps)
#define leeq(a,b) (((b)-(a))>-eps)
#define equl(a,b) (fabs((a)-(b))<eps)

#define lmax(a,b) ((a)>(b)?(a):(b))
#define lmin(a,b) ((a)<(b)?(a):(b))
#define fmax(a,b) (gret(a,b)?(a):(b))
#define fmin(a,b) (less(a,b)?(a):(b))

const int MAXV = 1200;

struct info
{
    int s;
    double x1,x2,y;
    info(double _x1=0,double _x2=0,double _y=0,int _s=0):
        x1(_x1),x2(_x2),y(_y),s(_s){}
    bool operator < (const info argu) const
    {
        return less(y,argu.y);
    }
}a[MAXV<<1];

#define lc(i) ((i)<<1)
#define rc(i) ((i)<<1|1)

int n,m,times=1,cov[MAXV<<2];
double u1,v1,u2,v2,len[MAXV<<2],xx[MAXV<<1];

void up(int i,int l,int r)
{
    if(cov[i]) len[i]=xx[r-1]-xx[l-2];
    else if(l==r) len[i]=0;
    else len[i]=len[lc(i)]+len[rc(i)];
}
void change(int x,int y,int w,int l=1,int r=MAXV-1,int i=1)
{
    if(x<=l&&r<=y) cov[i]+=w;
    else
    {
        if(x<=(l+r)/2) change(x,y,w,l,(l+r)/2,lc(i));
        if((l+r)/2<y) change(x,y,w,(l+r)/2+1,r,rc(i));
    }
    up(i,l,r);
}
int bin(double v)
{
    return lower_bound(xx,xx+m,v)-xx+1;
}

int main()
{
    #ifndef ONLINE_JUDGE
    freopen("Atlantis[2].txt","r",stdin);
    #endif

    while(scanf("%d",&n),n)
    {
        memset(len,0,sizeof(len));
        memset(cov,0,sizeof(cov));
        for(int i=0;n>i;i++)
        {
            scanf("%lf %lf %lf %lf",&u1,&v1,&u2,&v2);
            xx[lc(i)]=u1;
            xx[rc(i)]=u2;
            a[lc(i)]=info(u1,u2,v1,+1);
            a[rc(i)]=info(u1,u2,v2,-1);
        }
        sort(a,a+2*n);
        sort(xx,xx+2*n);
        for(int i=m=1;2*n>i;i++)
        {
            if(!equl(xx[i],xx[i-1])) xx[m++]=xx[i];
        }
        double cnt=0;
        change(bin(a[0].x1)+1,bin(a[0].x2),a[0].s);
        for(int i=1;2*n>i;i++)
        {
            cnt+=len[1]*(a[i].y-a[i-1].y);
            change(bin(a[i].x1)+1,bin(a[i].x2),a[i].s);
        }
        printf("Test case #%d\n",times++);
        printf("Total explored area: %.2f\n\n",cnt);
    }
    return 0;
}