【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。我想哭了.....

View Code

  1 #include <iostream>
  2 #include <cstdio>
  3 #include <cstring>
  4 #include <algorithm>
  5 using namespace std;
  6 
  7 #define lson l,mid,ID<<1
  8 #define rson mid+1,r,ID<<1|1
  9 const int maxn=150101;
 10 typedef long long lld;
 11 int flag[maxn];
 12 lld sum[maxn];
 13 int X[maxn], Z[maxn];
 14 
 15 struct Node   ///这题快把我整哭了
 16 {
 17     int lx, rx, y, z, s;
 18     Node(){}
 19     Node(int lx_, int rx_ , int y_, int z_, int s_)
 20     {
 21         lx=lx_, rx=rx_, y=y_, z=z_, s=s_;
 22     }
 23     bool operator<(const Node &S) const
 24     {
 25         if(y==S.y) return s>S.s;
 26         else return y<S.y;
 27     }
 28 }line[maxn], tmp[maxn];
 29 
 30 
 31 int find(int x, int M)
 32 {
 33     int l,r,m;
 34     l=1;
 35     r=M;
 36     while(l<=r)
 37     {
 38         m=(l+r)>>1;
 39         if(X[m]==x)
 40             return m;
 41         if(X[m]<x)
 42             l=m+1;
 43         else
 44             r=m-1;
 45     }
 46 }
 47 void Push_up(int ID,int l,int r)
 48 {
 49     if(flag[ID])sum[ID]=X[r+1]-X[l];
 50     else if(l==r)sum[ID]=0;
 51     else sum[ID]=sum[ID<<1]+sum[ID<<1|1];
 52 }
 53 void Update(int x,int y,int z,int l,int r,int ID)
 54 {
 55     int mid;
 56     if(x<=l&&r<=y)
 57     {
 58         flag[ID]+=z;
 59         Push_up(ID,l,r);
 60         return ;
 61     }
 62     mid=(l+r)>>1;
 63     if(x<=mid)
 64         Update(x,y,z,lson);
 65     if(y>mid)
 66         Update(x,y,z,rson);
 67     Push_up(ID,l,r);
 68 }
 69 
 70 int main()
 71 {
 72     int n, m, T, tcase=0;
 73     cin >> T;
 74     while(T--)
 75     {
 76         cin >> n >> m;
 77         Z[0]=0;
 78         for(int i=1; i<=m; i++)
 79             cin >> Z[i];
 80         int num=0;
 81         for(int i=0; i<n; i++)
 82         {
 83             int x1, x2, y1, y2, id;
 84             scanf("%d%d%d%d%d",&x1,&y1,&x2,&y2,&id);
 85             line[++num]=Node(x1,x2,y1,Z[id],1);
 86             X[num]=x1;
 87             line[++num]=Node(x1,x2,y2,Z[id],-1);
 88             X[num]=x2;
 89         }
 90         sort(Z,Z+m+1);
 91         sort(X+1,X+num+1);
 92         sort(line+1,line+num+1);
 93         int ep=1;
 94         for(int i=2; i<=num; i++)
 95             if(X[i]!=X[ep]) X[++ep]=X[i];
 96         lld ans=0;
 97         for(int i=0; i<m; i++)
 98         {
 99             memset(sum,0,sizeof(sum));
100             memset(flag,0,sizeof(flag));
101             lld tp=0, cnt=0;
102             for(int j=1; j<=num; j++)
103                if(line[j].z>Z[i]) tmp[++cnt]=line[j];
104             for(int j=1; j<cnt; j++)
105             {
106                 int l=find(tmp[j].lx,ep);
107                 int r=find(tmp[j].rx,ep)-1;
108                 Update(l,r,tmp[j].s,1,ep-1,1);
109                 tp+=sum[1]*(lld)(tmp[j+1].y-tmp[j].y);
110             }
111             ans+=tp*(lld)(Z[i+1]-Z[i]);
112         }
113         printf("Case %d: %I64d\n",++tcase,ans);
114     }
115     return 0;
116 }

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

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