【HDU 3255 Farming 】 线段树之扫描线（体积并）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=3255

题目大意：给你多块蔬菜地，蔬菜地有重叠部分，每部分种的蔬菜可能相同也可能不同，重叠部分种植蔬菜价值最贵的，最后让你求最大利润。

解题思路：

     这是我做的最坑爹的一题，几乎让我快哭了，Wrong answer了三十多次， 重拍了五六次代码。

     每块蔬菜地种植蔬菜收获的利润为 val=x*y*price。  面积乘以价格，题目的重点转换在于如何确定重叠区域怎么让它种植最贵的蔬菜。

     观察利润计算公式 :   x*y*price <==> x*y*h    可以转换为求体积并。

     说实话，这题一看真的很简单，一交就是一直错，代码每个部分都有调试N遍。 昨天一个就一直在调试，五个小时无果。和别人一对照，思路一点都没错，今天下午再重拍了一遍代码，又开始了Wrong answer之路，然后带着那颗受伤的心去改数组范围。点最大3W， 我线段树开12W多，存X的点开6W多。 点开4W错， 6WTLE， 索性所有数组范围开150005，错。 150101 AC。 我想哭了.....

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;

#define lson l,mid,ID<<1
#define rson mid+1,r,ID<<1|1
const int maxn=150101;
typedef long long lld;
int flag[maxn];
lld sum[maxn];
int X[maxn], Z[maxn];

struct Node   ///这题快把我整哭了
{
    int lx, rx, y, z, s;
    Node(){}
    Node(int lx_, int rx_ , int y_, int z_, int s_)
    {
        lx=lx_, rx=rx_, y=y_, z=z_, s=s_;
    }
    bool operator<(const Node &S) const
    {
        if(y==S.y) return s>S.s;
        else return y<S.y;
    }
}line[maxn], tmp[maxn];


int find(int x, int M)
{
    int l,r,m;
    l=1;
    r=M;
    while(l<=r)
    {
        m=(l+r)>>1;
        if(X[m]==x)
            return m;
        if(X[m]<x)
            l=m+1;
        else
            r=m-1;
    }
}
void Push_up(int ID,int l,int r)
{
    if(flag[ID])sum[ID]=X[r+1]-X[l];
    else if(l==r)sum[ID]=0;
    else sum[ID]=sum[ID<<1]+sum[ID<<1|1];
}
void Update(int x,int y,int z,int l,int r,int ID)
{
    int mid;
    if(x<=l&&r<=y)
    {
        flag[ID]+=z;
        Push_up(ID,l,r);
        return ;
    }
    mid=(l+r)>>1;
    if(x<=mid)
        Update(x,y,z,lson);
    if(y>mid)
        Update(x,y,z,rson);
    Push_up(ID,l,r);
}

int main()
{
    int n, m, T, tcase=0;
    cin >> T;
    while(T--)
    {
        cin >> n >> m;
        Z[0]=0;
        for(int i=1; i<=m; i++)
            cin >> Z[i];
        int num=0;
        for(int i=0; i<n; i++)
        {
            int x1, x2, y1, y2, id;
            scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&id);
            line[++num]=Node(x1,x2,y1,Z[id],1);
            X[num]=x1;
            line[++num]=Node(x1,x2,y2,Z[id],-1);
            X[num]=x2;
        }
        sort(Z,Z+m+1);
        sort(X+1,X+num+1);
        sort(line+1,line+num+1);
        int ep=1;
        for(int i=2; i<=num; i++)
            if(X[i]!=X[ep]) X[++ep]=X[i];
        lld ans=0;
        for(int i=0; i<m; i++)
        {
            memset(sum,0,sizeof(sum));
            memset(flag,0,sizeof(flag));
            lld tp=0, cnt=0;
            for(int j=1; j<=num; j++)
               if(line[j].z>Z[i]) tmp[++cnt]=line[j];
            for(int j=1; j<cnt; j++)
            {
                int l=find(tmp[j].lx,ep);
                int r=find(tmp[j].rx,ep)-1;
                Update(l,r,tmp[j].s,1,ep-1,1);
                tp+=sum[1]*(lld)(tmp[j+1].y-tmp[j].y);
            }
            ans+=tp*(lld)(Z[i+1]-Z[i]);
        }
        printf("Case %d: %I64d\n",++tcase,ans);
    }
    return 0;
}