BZOJ3476 : [Usaco2014 Mar]The Lazy Cow

旋转坐标系后转化为正方形，$x'=x+y$，$y'=x-y+1000001$，$k'=2k-1$

两根扫描线从左往右扫

f[i]表示y坐标下边界为i时的价值和

每次加入/删除一个点等价于一段区间加减

然后查询全局最大值

线段树维护扫描线之间的点

 

 

+ View Code?

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26

#include<cstdio> #include<algorithm> const int E=1000001,N=E*2-1,M=N*4+10; struct P{int x,y,g;}a[100010]; inline bool cmp(P a,P b){return a.x<b.x;} int n,m,k,i,j,x,y,ans,v[M],tag[M]; inline void add1(int x,int p){v[x]+=p,tag[x]+=p;} void add(int x,int a,int b,int c,int d,int p){   if(c<=a&&b<=d){add1(x,p);return;}   if(tag[x])add1(x<<1,tag[x]),add1(x<<1|1,tag[x]),tag[x]=0;   int mid=(a+b)>>1;   if(c<=mid)add(x<<1,a,mid,c,d,p);   if(d>mid)add(x<<1|1,mid+1,b,c,d,p);   v[x]=v[x<<1]>v[x<<1|1]?v[x<<1]:v[x<<1|1]; } int main(){   scanf("%d%d",&n,&k);k=k*2+1;   for(i=1;i<=n;i++)scanf("%d%d%d",&a[i].g,&x,&y),a[i].x=x+y,a[i].y=x-y+E;   std::sort(a+1,a+n+1,cmp);   for(i=j=1;i<=n;i++){     while(a[i].x-a[j].x+1>k)add(1,1,N,a[j].y-k+1,a[j].y,-a[j].g),j++;     add(1,1,N,a[i].y-k+1,a[i].y,a[i].g);     if(ans<v[1])ans=v[1];   }   return printf("%d",ans),0; }