hdu 3255 Farming 线段树空间体积

刚刚敲完空间体积相交，但是这道题竟然一开始没有思路，果然还是该打，总结线段树的做题经验就是要仔细看数据的范围。

解题思路：

     因为每一种植物的价值最大不过到100，那么我们需要的价值就是面积*价值，那么我们何不把价值作为第三个坐标轴z轴栏运算呢~

     接下来不说了，，看代码

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<string.h>
using std::unique;
using std::sort;
const int N = 30005;
struct Re
{
       int x1,x2,y1,y2,z1,z2;
}Re[N];
struct line
{
       int y,x1,x2,s;
       bool operator < (const line & tmp)const
       {
           return y<tmp.y;
       }
}Line[N<<2];
int Z[N<<2],X[N<<3],cover[N<<3],Price[N<<2];
__int64 sum[N<<3];
int Fin(int k,int len)
{
    int l=0,r=len;
    while(l<r)
    {
        int m=(l+r)>>1;
        if(X[m]==k)return m;
        if(X[m]>k)r=m-1;
        else l=m+1;
    }
    return l;
}
void PushUp(int t,int l,int r)
{
     if(cover[t]>0)sum[t]=(__int64)(X[r+1]-X[l]);
     else
       if(l==r)sum[t]=0;
       else  sum[t]=sum[t<<1]+sum[t<<1|1];
}
void update(int t,int l,int r,int L,int R,int val)
{
     if(L<=l&&r<=R)
     {
        cover[t]+=val;
        PushUp(t,l,r);
        return ;
     }
     int m=(r+l)>>1;
     if(L<=m)update(t<<1,l,m,L,R,val);
     if(R>m)update(t<<1|1,m+1,r,L,R,val);
     PushUp(t,l,r);
}
int main()
{
    int t,n,m,k,T=0;
    scanf("%d",&t);
    while(t--)
    {
          scanf("%d%d",&n,&m);
          int topz=0;
          Z[topz++]=0;
          for(int i=1;i<=m;i++){scanf("%d",&Price[i]);Z[topz++]=Price[i];}
          for(int i=0;i<n;i++)
          {
              scanf("%d%d%d%d%d",&Re[i].x1,&Re[i].y1,&Re[i].x2,&Re[i].y2,&k);
              Re[i].z1=0,Re[i].z2=Price[k];
          }
          sort(Z,Z+topz);
          topz=unique(Z,Z+topz)-Z;
          __int64 ans=0;
          for(int i=0;i<topz-1;i++)
          {
              int topx=0,topl=0;
              for(int j=0;j<n;j++)
              {
                  if(Re[j].z1<=Z[i]&&Re[j].z2>Z[i])
                  {
                      Line[topl].y=Re[j].y1;Line[topl].x1=Re[j].x1;Line[topl].x2=Re[j].x2;Line[topl].s=1;topl++;
                      Line[topl].y=Re[j].y2;Line[topl].x1=Re[j].x1;Line[topl].x2=Re[j].x2;Line[topl].s=-1;topl++;
                      X[topx++]=Re[j].x1;X[topx++]=Re[j].x2;
                  }
              }
              sort(X,X+topx);
              sort(Line,Line+topl);
              topx=unique(X,X+topx)-X;
              memset(sum,0,sizeof(sum));
              memset(cover,0,sizeof(cover));
              __int64 ret=0;
              for(int j=0;j<topl-1;j++)
              {
                  int l=Fin(Line[j].x1,topx-1);
                  int r=Fin(Line[j].x2,topx-1)-1;
                  if(l<=r)update(1,0,topx-1,l,r,Line[j].s);
                  ret+=(sum[1]*(__int64)(Line[j+1].y-Line[j].y));
              }
              ans+=ret*(__int64)(Z[i+1]-Z[i]);
          }
          printf("Case %d: %I64d\n",++T,ans);
    }
    return 0;
}