扫描线 n块菜地最大利润种菜（重叠区域种最贵） hdu3255 Farming

题目链接   hdu3255 Farming

题目大意：有n块蔬菜地，每块蔬菜地中种一种作物，每种作物都有一个价格，当在同一区域内种植了两种不同的作物时，作物价格大的生存下来，作物价格小的死亡。问你最后能得到的最大利润是多少。

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

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

#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;
}