这容斥真是写的我心态爆炸…
考虑用至少0个坏点的-至少1个坏点的+至少两个坏点的-至少三个坏点的+至少四个坏点的。
我们发现对于斜着的正方形,可以直接在框住它的大正方形处计数,边长为i的大正方形内就有i个正方形。
且我们发现每个点出现且仅出现在一个正方形上。因此对于至少一个坏点的正方形,我们只需要统计过一个坏点的正方形框架个数即可。然后就是毒瘤的分类讨论!
至少两个三个四个的可以通过枚举两个坏点,算出另外两点的可能坐标,Hash来判断。
还要注意去重。
复杂度
O(n+k2)
O
(
n
+
k
2
)
#include <bits/stdc++.h>
using namespace std;
#define ll long long
#define inf 0x3f3f3f3f
#define N 1000003
#define mod 100000007
inline char gc(){
static char buf[1<<16],*S,*T;
if(T==S){T=(S=buf)+fread(buf,1,1<<16,stdin);if(T==S) return EOF;}
return *S++;
}
inline ll read(){
ll x=0,f=1;char ch=gc();
while(ch<'0'||ch>'9'){if(ch=='-') f=-1;ch=gc();}
while(ch>='0'&&ch<='9') x=x*10+ch-'0',ch=gc();
return x*f;
}
int n,m,K,ans0,ans1,ans2,ans3,ans4,h[N],num=0;
inline void inc(int &x,int y){x+=y;if(x>=mod) x-=mod;}
struct P{
int x,y;
}a[2010];
struct Hash_table{
ll key;int next;
}data[2010];
inline void ins(ll key){
int x=key%N;
for(int i=h[x];i;i=data[i].next)
if(data[i].key==key) return;
data[++num].key=key;data[num].next=h[x];h[x]=num;
}
inline bool hs(ll key){
int x=key%N;
for(int i=h[x];i;i=data[i].next)
if(data[i].key==key) return 1;return 0;
}
inline int cal(int l1,int l2,int h){
h--;if(l1>l2) swap(l1,l2);
if(h<=l1) return (ll)h*(h+1)/2%mod;
if(h<=l2) return ((ll)l1*(l1+1)/2+(ll)l1*(h-l1))%mod;
if(h>=l1+l2) return (ll)l1*l2%mod;
return ((ll)l1*(l1+1)/2+(ll)l1*(l2-l1)+(ll)(l1-1+l1-h+l2)*(h-l2)/2)%mod;
}
inline void check(int x1,int y1,int x2,int y2){
if(x1<0||x1>n||x2<0||x2>n||y1<0||y1>m||y2<0||y2>m) return;
int res=hs((ll)x1*(m+1)+y1)+hs((ll)x2*(m+1)+y2);
ans2++;ans3+=res;ans4+=res==2;
}
int main(){
// freopen("a.in","r",stdin);
n=read();m=read();K=read();
for(int i=1;i<=min(n,m);++i) inc(ans0,(ll)i*(n-i+1)*(m-i+1)%mod);
for(int i=1;i<=K;++i){
a[i].x=read(),a[i].y=read();ins((ll)a[i].x*(m+1)+a[i].y);
inc(ans1,min(a[i].x,a[i].y)+min(a[i].x,m-a[i].y)+min(n-a[i].x,a[i].y)+min(n-a[i].x,m-a[i].y));
inc(ans1,cal(a[i].x,n-a[i].x,a[i].y));
inc(ans1,cal(a[i].x,n-a[i].x,m-a[i].y));
inc(ans1,cal(a[i].y,m-a[i].y,a[i].x));
inc(ans1,cal(a[i].y,m-a[i].y,n-a[i].x));
}for(int i=1;i<=K;++i)
for(int j=i+1;j<=K;++j){
int x1=a[i].x,y1=a[i].y,x2=a[j].x,y2=a[j].y;
check(x1+y2-y1,y1+x1-x2,x2+y2-y1,y2+x1-x2);
check(x1+y1-y2,y1+x2-x1,x2+y1-y2,y2+x2-x1);
if((x1+x2+y1+y2)&1) continue;
check(x1+x2+y2-y1>>1,y1+y2+x1-x2>>1,x1+x2+y1-y2>>1,y1+y2+x2-x1>>1);
}printf("%d\n",((ans0-ans1+ans2-ans3/3+ans4/6)%mod+mod)%mod);
return 0;
}