考虑容斥：计数没有任意一个点的矩形个数。
首先长宽分别为 \(x,y\) 的矩形的所有子矩形个数为：
\[\frac {x (x+1) y(y+1)} 4
\]
按横坐标排序，如果限制左边的取值在相邻两个点之间，并且限制上下界，其右边的范围可以确定。
如上图，限制左边在第零个点与第一个点之间。发现蓝绿色矩形间会有重（矩形 \((1,1),(7,7)\) 和矩形 \((1,8),(7,10)\)），直接容斥减掉。
当两个点的时候，我们发现左限制有重（矩形 \((5,1),(7,4)\) 和矩形 \((5,6),(7,8)\) 等），容斥减掉。
#include<bits/stdc++.h>
#define int long long
using namespace std;
namespace Fast_OI{char buf[1000000],*p1=buf,*p2=buf,obuf[1000000],*p3=obuf;#define getchar() (p1==p2&&(p2=(p1=buf)+fread(buf,1,1000000,stdin),p1==p2)?EOF:*p1++)#define putchar(x) (p3-obuf<1000000?*p3++=x:(fwrite(obuf,1,p3-obuf,stdout),p3=obuf,*p3++=x))int read(){int x=0;bool f=1;char c=getchar();while(!isdigit(c)){if(c=='-')f=0;c=getchar();}while(isdigit(c))x=(x<<3)+(x<<1)+(c^48),c=getchar();return f?x:-x;}void puts(const char* str){while(*str!='\0')putchar(*str ++); putchar('\n'); }void write(int x){if(x<0)putchar('-'),x=-x;if(x>9)write(x/10);putchar(x%10+48);} void flush() { fwrite(obuf,1,p3-obuf,stdout); }
}using namespace Fast_OI;
const int N=1e3+5;
const int mod=1e9+7;
int n,m,ans;
int k;
pair<int,int> a[N];
set<pair<int,pair<int,int> > > s;
inline int calc(int x,int y){if(x<=0||y<=0) return 0;return (x*(x+1)/2%mod)*(y*(y+1)/2%mod)%mod;
}
signed main(){freopen("count.in","r",stdin);freopen("count.out","w",stdout);n=read(),m=read(),k=read();for(int i=1;i<=k;i++)a[i]=make_pair(read(),read());sort(a+1,a+1+k);a[k+1]=make_pair(n,0);for(int i=0;i<=k;i++){s.clear();s.insert(make_pair(1,make_pair(m,a[i].first)));for(int j=i+1;j<=k;j++){auto ir=s.upper_bound(make_pair(a[j].second,make_pair(m+1,m+1)));auto il=ir;if(il!=s.begin()) il--;else continue;int l=il->first,r=il->second.first,v=il->second.second;if(l>a[j].second||r<a[j].second)continue;s.erase(il);ans-=calc(v-a[i].first,r-l+1);ans+=calc(v-a[i+1].first,r-l+1);ans+=calc(a[j].first-1-a[i].first,r-l+1);ans-=calc(a[j].first-1-a[i+1].first,r-l+1);ans%=mod;if(ans<0) ans+=mod;if(l<a[j].second)s.insert(make_pair(l,make_pair(a[j].second-1,a[j].first-1)));if(r>a[j].second)s.insert(make_pair(a[j].second+1,make_pair(r,a[j].first-1)));}for(auto it=s.begin();it!=s.end();it++){int l=it->first,r=it->second.first,v=it->second.second;ans-=calc(v-a[i].first,r-l+1);ans+=calc(v-a[i+1].first,r-l+1);ans+=calc(n-a[i].first,r-l+1);ans-=calc(n-a[i+1].first,r-l+1);ans%=mod;if(ans<0) ans+=mod;}}ans=(calc(n,m)-ans+mod)%mod;printf("%lld",ans);return 0;
}