hdu 3340 线段树思路活用

开始看到这道题目的时候，神题啊。。  想了一个下午，没思路。。。

后面到网上一看这位大哥的解说，豁然开朗http://blog.csdn.net/gotoac/article/details/7588254

每个区间分割多边形是可能是三角形或者梯形或者四边形，都是上底+下底乘高除2，这个小学就会的东西，看来现在还是不能够活用的啊。。。。

#include<iostream>
#include<stdlib.h>
#include<cstdio>
#include<algorithm>
using std::swap;
using std::sort;
using std::unique;
const int N = 25005;
struct order
{
       int x[6],y[6],s;
}Order[N];
int X[N*5];
double sum[N<<4],lh[N<<4],rh[N<<4],len[N<<4];
int Fin(int k,int len)
{
    int l=0,r=len;
    while(l<r)
    {
        int m=(r+l)>>1;
        if(X[m]==k)return m;
        if(X[m]>k)r=m-1;
        else l=m+1;
    }
    return l;
}
double Cal(double ly,double ry,int d,int l,int r)
{
       return ly+(1.0*(r-l)*(ry-ly)/d);
}
void build(int t,int l,int r)
{
     len[t]=X[r+1]-X[l];
     sum[t]=lh[t]=rh[t]=0;
     if(l==r)return ;
     int m=(r+l)>>1;
     build(t<<1,l,m);
     build(t<<1|1,m+1,r);
}
void PushDown(int t,int l,int r)
{
     if(lh[t]||rh[t])
     {
        int t1=t<<1;
        int t2=t1|1;
        int m=(r+l)>>1;
        double mh=Cal(lh[t],rh[t],len[t],X[l],X[m+1]);
        lh[t1]+=lh[t];
        rh[t1]+=mh;
        lh[t2]+=mh;
        rh[t2]+=rh[t];
        sum[t1]+=len[t1]*(lh[t]+mh)/2;
        sum[t2]+=len[t2]*(rh[t]+mh)/2;
        lh[t]=rh[t]=0;
     }
}
void PushUp(int t)
{
    sum[t]=sum[t<<1]+sum[t<<1|1];
}
void update(int t,int l,int r,int L,int R,double ly,double ry)
{
     if(L<=l&&r<=R)
     {
        double ld=Cal(ly,ry,X[R+1]-X[L],X[L],X[l]);
        double rd=Cal(ly,ry,X[R+1]-X[L],X[L],X[r+1]);
        lh[t]+=ld;
        rh[t]+=rd;
        sum[t]+=len[t]*(ld+rd)/2;
        return ;
     }
     PushDown(t,l,r);
     int m=(r+l)>>1;
     if(L<=m)update(t<<1,l,m,L,R,ly,ry);
     if(R>m)update(t<<1|1,m+1,r,L,R,ly,ry);
     PushUp(t);
}
double query(int t,int l,int r,int L,int R)
{
     if(L<=l&&r<=R)return sum[t];
     PushDown(t,l,r);
     int m=(r+l)>>1;
     double ans=0;
     if(L<=m)ans+=query(t<<1,l,m,L,R);
     if(R>m)ans+=query(t<<1|1,m+1,r,L,R);
     return ans;
}
int main()
{
    int t,n;
    char ord[5];
    scanf("%d",&t);
    while(t--)
    {
          scanf("%d",&n);
          int topx=0;
          for(int i=0;i<n;i++)
          {
              scanf("%s",ord);
              if(ord[0]=='Q')
              {
                 Order[i].s=0;
                 scanf("%d%d",&Order[i].x[0],&Order[i].y[0]);
                 X[topx++]=Order[i].x[0];
                 X[topx++]=Order[i].y[0];
              }
              else
              {
                 scanf("%d",&Order[i].s);
                 for(int j=0;j<Order[i].s;j++)
                 {
                     scanf("%d%d",&Order[i].x[j],&Order[i].y[j]);
                     X[topx++]=Order[i].x[j];
                 }
                 Order[i].x[Order[i].s]=Order[i].x[0];
                 Order[i].y[Order[i].s]=Order[i].y[0];
                 Order[i].s++;
              }
          }
          sort(X,X+topx);
          topx=unique(X,X+topx)-X;
          build(1,0,topx-1);

          for(int i=0;i<n;i++)
          {
              if(Order[i].s==0)
              {
                 if(Order[i].x[0]>Order[i].y[0])swap(Order[i].x[0],Order[i].y[0]);
                 int L=Fin(Order[i].x[0],topx-1);
                 int R=Fin(Order[i].y[0],topx-1)-1;
                 printf("%.3lf\n",L>R?0:query(1,0,topx-1,L,R));
              }
              else
              {
                 for(int j=1;j<Order[i].s;j++)
                 {
                     int L=Fin(Order[i].x[j-1],topx-1);
                     int R=Fin(Order[i].x[j],topx-1);
                     if(R>L)update(1,0,topx-1,L,R-1,-1.0*Order[i].y[j-1],-1.0*Order[i].y[j]);
                     if(L>R)update(1,0,topx-1,R,L-1,1.0*Order[i].y[j],1.0*Order[i].y[j-1]);
                 }
              }
          }
    }
    return 0;
}