HDU-3255-Farming

来源：互联网发布：淘宝床上四件套特价编辑：程序博客网时间：2024/05/22 23:34
求体积的一个比较好的题，由于m很小，直接枚举就行。
代码：
#include<cstdio>#include<cstring>#include<iostream>#include<algorithm>using namespace std;const int inf=1<<29;const int maxn=3e4+1000;const int maxm=maxn*8;struct Line{    int x1;    int x2;    int y;    int flag;    int pri;    Line(){}    Line(int sx1,int sx2,int sy,int spri,int sflag)    {        x1=sx1;        x2=sx2;        y=sy;        pri=spri;        flag=sflag;    }    bool operator < (const Line &l)const    {        if(y==l.y)            return x1<l.x1;        return y<l.y;    }}L[maxn*2];struct Node{    int l;    int r;    int len;    int cover;}t[maxm];int n,m,cnt,X[maxn*2],index[maxn*2],cntx,pri[4];void Build(int l,int r,int index){    t[index].l=l;    t[index].r=r;    t[index].len=t[index].cover=0;    if(l==r)        return ;    int mid=(l+r)>>1;    Build(l,mid,index<<1);    Build(mid+1,r,index<<1|1);}void PushUp(int index){    if(t[index].cover>0)        t[index].len=X[t[index].r+1]-X[t[index].l];    else if(t[index].l==t[index].r)        t[index].len=0;    else        t[index].len=t[index<<1].len+t[index<<1|1].len;}void Update(int l,int r,int index,int val){    if(l>r)        return;    if(t[index].l==l&&t[index].r==r)    {        t[index].cover+=val;        PushUp(index);        return;    }    int mid=(t[index].l+t[index].r)>>1;    if(r<=mid)        Update(l,r,index<<1,val);    else if(l>mid)        Update(l,r,index<<1|1,val);    else    {        Update(l,mid,index<<1,val);        Update(mid+1,r,index<<1|1,val);    }    PushUp(index);}int main(){    int T,cas=1;    scanf("%d",&T);    while(T--)    {        scanf("%d%d",&n,&m);        for(int i=1;i<=m;i++)            scanf("%d",&pri[i]);        cnt=cntx=0;        for(int i=0;i<n;i++)        {            int x1,x2,y1,y2,s;            scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&s);            L[cnt++]=Line(x1,x2,y1,pri[s],1);            L[cnt++]=Line(x1,x2,y2,pri[s],-1);            X[cntx++]=x1;            X[cntx++]=x2;        }        sort(X,X+cntx);        cntx=unique(X,X+cntx)-X;        sort(pri+1,pri+m+1);        sort(L,L+cnt);        long long ans=0;        for(int i=1;i<=m;i++)        {            int scnt=0;            Build(0,cntx,1);            for(int j=0;j<cnt;j++)                if(L[j].pri>=pri[i])                    index[scnt++]=j;            for(int j=0;j<scnt-1;j++)            {                int sl=lower_bound(X,X+cntx,L[index[j]].x1)-X;                int sr=lower_bound(X,X+cntx,L[index[j]].x2)-X;                sr--;                Update(sl,sr,1,L[index[j]].flag);                ans+=(long long)t[1].len*(L[index[j+1]].y-L[index[j]].y)*(pri[i]-pri[i-1]);            }        }        printf("Case %d: %I64d\n",cas++,ans);    }    return 0;}
0 0