ARC063F すぬけ君の塗り絵 2 / Snuke's Coloring 2

题面

一句话题面：给你一些点，求这些点之中夹的最大的矩形周长。（考虑边界）

Solution

首先是一个结论,答案矩形一定经过(x=frac{w}{2})或经过(y=frac{h}{2}),不然答案一定不优.

怎么说?因为答案一定(ge 2*max(h,w)+1),这个可以通过左右|上下显然得出.

接下来我们考虑扫描线,对于从左往右的(p_i.x),令(p_i.x)为右边界,单调栈维护上下边界然后左边界直接每一次-就行了.

唯一的问题在于弹栈时的一些小操作,代码中都有注释.

Code

/*
  mail: mleautomaton@foxmail.com
  author: MLEAutoMaton
  This Code is made by MLEAutoMaton
*/
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<math.h>
#include<algorithm>
#include<queue>
#include<set>
#include<map>
#include<iostream>
using namespace std;
#define ll long long
#define REP(a,b,c) for(int a=b;a<=c;a++)
#define re register
#define file(a) freopen(a".in","r",stdin);freopen(a".out","w",stdout)
#define int ll
inline int gi(){
	int f=1,sum=0;char ch=getchar();
	while(ch>'9' || ch<'0'){if(ch=='-')f=-1;ch=getchar();}
	while(ch>='0' && ch<='9'){sum=(sum<<3)+(sum<<1)+ch-'0';ch=getchar();}
	return f*sum;
}
const int N=2000010;
int w,h,n;
struct node{
	int x,y;
	bool operator<(const node &b)const{return x<b.x || (x==b.x && y<b.y);}
}p[N];
typedef pair<int,int> pii;
#define mp make_pair
pii a[N],b[N];//stack
int mx[N],tag[N],ans;
void pushup(int o){mx[o]=max(mx[o<<1],mx[o<<1|1]);}
void pushdown(int o){
	if(!tag[o])return;
	tag[o<<1]+=tag[o];tag[o<<1|1]+=tag[o];
	mx[o<<1]+=tag[o];mx[o<<1|1]+=tag[o];
	tag[o]=0;
}
void modify(int o,int l,int r,int posl,int posr,int v){
	if(posl<=l && r<=posr){tag[o]+=v;mx[o]+=v;return;}
	pushdown(o);int mid=(l+r)>>1;
	if(posl<=mid)modify(o<<1,l,mid,posl,posr,v);
	if(mid<posr)modify(o<<1|1,mid+1,r,posl,posr,v);
	pushup(o);
}
void work(){
	memset(mx,0,sizeof(mx));memset(tag,0,sizeof(tag));
	sort(p+1,p+n+1);int l=0,r=0;
	for(int i=1;i<n;i++){
		if(p[i].y<=h/2){
			//单调栈维护
			int lst=i-1;
			while(l && p[i].y>a[l].second){
				modify(1,1,n,a[l].first,lst,a[l].second-p[i].y);
				lst=a[l--].first-1;
			}
			if(lst!=i-1)
				a[++l]=mp(lst+1,p[i].y);
		}
		else{
			//单调栈维护
			int lst=i-1;
			while(l && p[i].y<b[r].second){
				modify(1,1,n,b[r].first,lst,p[i].y-b[r].second);
				lst=b[r--].first-1;
			}
			if(lst!=i-1)b[++r]=mp(lst+1,p[i].y);
		}
		modify(1,1,n,i,i,h-p[i].x);//为了便于计算答案.
		a[++l]=mp(i,0);b[++r]=mp(i,h);//为了减去上面或者下面.
		ans=max(ans,mx[1]+p[i+1].x);//计算答案.
	}
}
signed main(){
	w=gi();h=gi();n=gi();
	for(int i=1;i<=n;i++)p[i].x=gi(),p[i].y=gi();
	p[++n]=(node){0,0};p[++n]=(node){w,h};
	work();
	for(int i=1;i<=n;i++)swap(p[i].x,p[i].y);swap(h,w);
	work();
	printf("%lld
",ans*2);
	return 0;
}