联考20200730 T3 小B的农场


分析:
没挖掘出性质可以\(O(n^2logn)\)得到60分。。。
发现\(1*n\)的农田是一定划得出来的,所以答案最小为\(2(max(W,H)+1)\)
之后我们发现如果要比这个答案大,农田一定过横向或者纵向的中线
假设过纵向中线(横向的同样处理)
每一次下边界向下拓展时,新加入的会对左右限制的点用线段维护一下就好了
复杂度\(O(nlogn)\)

#include<cstdio>
#include<cmath>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<vector>
#include<string>

#define maxn 300005
#define INF 0x3f3f3f3f
#define MOD 998244353

using namespace std;

inline long long getint()
{
	long long num=0,flag=1;char c;
	while((c=getchar())<'0'||c>'9')if(c=='-')flag=-1;
	while(c>='0'&&c<='9')num=num*10+c-48,c=getchar();
	return num*flag;
}

int n,W,H,ans;
struct node{
	int x,y;
}p[maxn];
int mx[maxn<<2],lz[maxn<<2];
int s1[maxn],t1,s2[maxn],t2;

inline void build(int i,int l,int r)
{
	lz[i]=0,mx[i]=-p[l].y;
	if(l==r)return;
	int mid=(l+r)>>1;
	build(i<<1,l,mid),build(i<<1|1,mid+1,r);
}

inline void update(int i,int l,int r,int ql,int qr,int x)
{
	if(qr<l||r<ql)return;
	if(ql<=l&&r<=qr){lz[i]+=x,mx[i]+=x;return;}
	int mid=(l+r)>>1;
	update(i<<1,l,mid,ql,qr,x),update(i<<1|1,mid+1,r,ql,qr,x);
	mx[i]=max(mx[i<<1],mx[i<<1|1])+lz[i];
}
inline bool cmp(node x,node y)
{return x.y<y.y||(x.y==y.y&&x.x<y.x);}

inline void solve()
{
	sort(p+1,p+n+1,cmp);
	build(1,1,n),t1=t2=0;
	for(int i=1;i<=n;i++)
	{
		if(i>1)update(1,1,n,i-1,i-1,W);
		ans=max(ans,(p[i].y+mx[1]));
		if(p[i].x<=W/2)
		{
			update(1,1,n,s1[t1],i-1,-p[i].x);
			while(t1&&p[s1[t1]].x<p[i].x)update(1,1,n,s1[t1-1],s1[t1]-1,p[s1[t1]].x-p[i].x),--t1;
			s1[++t1]=i;
		}
		else
		{
			update(1,1,n,s2[t2],i-1,p[i].x-W);
			while(t2&&p[s2[t2]].x>p[i].x)update(1,1,n,s2[t2-1],s2[t2]-1,p[i].x-p[s2[t2]].x),--t2;
			s2[++t2]=i;
		}
	}
}

int main()
{
	W=getint(),H=getint(),n=getint();
	for(int i=1;i<=n;i++)p[i].x=getint(),p[i].y=getint();
	p[++n]=(node){0,0},p[++n]=(node){W,H};
	solve();
	for(int i=1;i<=n;++i)swap(p[i].x,p[i].y);
	swap(W,H);
	solve();
	printf("%d\n",ans*2); 
}

原文地址:https://www.cnblogs.com/IzayoiDoyo/p/13405994.html