#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 1061109568LL 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;}