#241. Lefkaritika

把点按横坐标排序,然后枚举正方形的底边所在行,可以对于横坐标相同的点,纵坐标大的没有卵用。
对于每个正方形,假设在它上方的、纵坐标最小的前提下横坐标最小的点可以管辖它。
那么很显然有些正方形是不会被管辖的。这些正方形夹在两个相邻的点中间,可以被O(1)算出来。
然后剩下的正方形都有且只有一个点可以管辖。那么可以通过直接枚举管辖它的点来搞事情。
枚举右边纵坐标小于它的点,左边横坐标小于等于它的点,然后得到一个矩形,其中的正方形个数可以O(1)算出来。
∑起来就是答案了。
到simpleoj上一看会发现这题没人A,其实这题的满分就只有86分233。
#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<vector>
#include<queue>
#include<cmath>
#include<map>
using namespace std;
#define rep(i,j,k) for(i=j;i<=k;++i)
#define per(i,j,k) for(i=j;i>=k;--i)
#define sqr(x) ((x)*(x))
#define G getchar()
#define LL long long
#define pll pair<LL,LL>
#define mkp make_pair
#define X first
#define Y second
#define N 200005
#define NN 1005
#define inf 1061109568
LL n,m,k,ans,key[NN],cnt,q[NN],Ft,Rr,le[NN],ri[NN];pll a[NN];
int read(){
	int x=0;char ch=G;bool flg=0;
	while((ch<48||ch>57)&&ch!=45)ch=G;
	if(ch==45)flg=1,ch=G;
	for(;ch>47&&ch<58;ch=G)x=x*10+ch-48;
	return flg?-x:x;
}

void cal(LL L,LL R,LL U,LL D){
	R=R-L,D=D-U;
	if(R<=0||D<=0)return;
	if(R<=D)ans+=R*(R+1)>>1;
	else ans+=D*(R+R-D+1)>>1;
}
void cal2(LL L,LL R,LL U,LL D){
	R=R-L,D=D-U;
	if(R<=0||D<=0)return;
	if(R<=D)ans-=R*(R+1)>>1;
	else ans-=D*(R+R-D+1)>>1;
}
int main(){
	LL x,y,i,j;
	n=read();m=read();k=read();
	rep(i,1,k){
		x=read();y=read();a[i]=mkp(x,y);
	}
	sort(a+1,a+k+1);
	rep(i,1,m){
		cnt=0;
		rep(j,1,k)if(a[j].Y>=i&&a[j].X!=a[key[cnt]].X)key[++cnt]=j;
		if(!cnt){
			cal(1,n,i,m);continue;
		}
		if(a[key[1]].X>1)cal(1,a[key[1]].X-1,i,m);
		rep(j,2,cnt)cal(a[key[j-1]].X+1,a[key[j]].X-1,i,m);
		if(a[key[cnt]].X<n)cal(a[key[cnt]].X+1,n,i,m);
		Rr=1;a[key[q[Ft=0]=0]=0]=mkp(0,i);
		rep(j,1,cnt){
			while(Ft<Rr&&a[key[q[Rr-1]]].Y>a[key[j]].Y)Rr--;
			le[j]=q[Rr-1];q[Rr++]=j;
		}
		Rr=1;a[key[q[Ft=0]=cnt+1]=k+1]=mkp(n+1,i-1);
		per(j,cnt,1){
			while(Ft<Rr&&a[key[q[Rr-1]]].Y>=a[key[j]].Y)Rr--;
			ri[j]=q[Rr-1];q[Rr++]=j;
		}
		rep(j,1,cnt){
			cal(a[key[le[j]]].X+1,a[key[ri[j]]].X-1,i,a[key[j]].Y-1);
			cal2(a[key[le[j]]].X+1,a[key[j]].X-1,i,a[key[j]].Y-1);
			cal2(a[key[j]].X+1,a[key[ri[j]]].X-1,i,a[key[j]].Y-1);
		}
	}
	printf("%lld\n",ans);
	return 0;
}

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值