#241. Lefkaritika

把点按横坐标排序,然后枚举正方形的底边所在行,可以对于横坐标相同的点,纵坐标大的没有卵用。
对于每个正方形,假设在它上方的、纵坐标最小的前提下横坐标最小的点可以管辖它。
那么很显然有些正方形是不会被管辖的。这些正方形夹在两个相邻的点中间,可以被O(1)算出来。
然后剩下的正方形都有且只有一个点可以管辖。那么可以通过直接枚举管辖它的点来搞事情。
枚举右边纵坐标小于它的点,左边横坐标小于等于它的点,然后得到一个矩形,其中的正方形个数可以O(1)算出来。
∑起来就是答案了。
到simpleoj上一看会发现这题没人A,其实这题的满分就只有86分233。
#include
#include
#include
#include
#include
#include
#include
#include
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
#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]].Xa[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=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;
}

你可能感兴趣的:(simpleOJ)