Hdu 3255 Farming(长方形体积并)

传送门:Hdu 3255 Farming

---

题意：有N块农田，每块农田中种一种作物，每种作物都有一个价格，当在同一区域内种植了两种不同的作物时，作物价格大的生存下来，作物价格小的死亡。求最后的所有作物的能买的总钱数。

---

思路:将价格转化为长方形的高,利用长方形体积并

---

#include <map>
#include <set>
#include <stack>
#include <queue>
#include <cmath>
#include <ctime>
#include <vector>
#include <cstdio>
#include <cctype>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
using namespace std;
#define INF 0x3f3f3f3f
#define inf -0x3f3f3f3f
#define lson l,m,rt<<1
#define rson m+1,r,rt<<1|1
#define mem0(a) memset(a,0,sizeof(a))
#define mem1(a) memset(a,-1,sizeof(a))
#define mem(a, b) memset(a, b, sizeof(a))
typedef long long ll;
const int maxn=50010;
struct node{
    int x,y,h,d,price;
}a[2*maxn],b[2*maxn];
int sum[maxn*8],mark[maxn*8],t[maxn*2],price[maxn];

int read(){
   char ch=getchar(),last=' ';
   while(ch<'0' || ch>'9')last=ch,ch=getchar();
   int ans=0;
   while(ch>='0' && ch<='9')ans=ans*10+ch-'0',ch=getchar();
   if(last=='-')ans=-ans;
   return ans;
}

bool cmp1(node u,node v){
    if(u.h==v.h)
        return u.d>v.d;
    return u.h<v.h;
}

int L,d,R;

void pushup(int rt,int l,int r){
    if(mark[rt])
        sum[rt]=t[r+1]-t[l];
    else if(l==r)
        sum[rt]=0;
    else
        sum[rt]=sum[rt<<1]+sum[rt<<1|1];
}

void update(int l,int r,int rt){
    if(L<=l&&R>=r){
        mark[rt]+=d;
        pushup(rt,l,r);
        return ;
    }
    int m=(l+r)>>1;
    if(L<=m)
        update(lson);
    if(R>m)
        update(rson);
    pushup(rt,l,r);
}

bool cmp(int x,int y){
    return x>y;
}

int main(){
    int _=read(),n,m,id,x1,x2,y1,y2;
    for(int case1=1;case1<=_;case1++){
        n=read(),m=read();
        for(int i=1;i<=m;i++)
            price[i]=read();
        for(int i=1;i<=n;i++){
            x1=read(),y1=read(),x2=read(),y2=read(),id=read();
            a[i*2-1]=(node){x1,x2,y1,1,price[id]};
            a[i*2]=(node){x1,x2,y2,-1,price[id]};
            t[i*2-1]=x1,t[i*2]=x2;
        }
        sort(price+1,price+m+1,cmp);
        sort(t+1,t+2*n+1);
        int M=unique(t+1,t+2*n+1)-t-1;
        long long ans=0,pre=0;
        for(int i=1;i<=m;i++){  //使用几个种子
            int tot=0,Count=0;
            for(int j=1;j<=2*n;j++)
                if(a[j].price>=price[i])
                    b[++tot]=a[j];
            sort(b+1,b+tot+1,cmp1);
            memset(sum,0,sizeof(sum));
            memset(mark,0,sizeof(mark));
            long long tmp=0;
            for(int j=1;j<tot;j++){
                L=lower_bound(t+1,t+M+1,b[j].x)-t,R=lower_bound(t+1,t+M+1,b[j].y)-t-1,d=b[j].d;
                update(1,M,1);
                tmp+=1LL*sum[1]*(b[j+1].h-b[j].h);
            }
            ans+=(tmp-pre)*price[i];
            pre=tmp;
        }
        printf("Case %d: %lld\n",case1,ans);
    }
    return 0;
}