把点按横坐标排序,然后枚举正方形的底边所在行,可以对于横坐标相同的点,纵坐标大的没有卵用。
对于每个正方形,假设在它上方的、纵坐标最小的前提下横坐标最小的点可以管辖它。
那么很显然有些正方形是不会被管辖的。这些正方形夹在两个相邻的点中间,可以被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;
}