#include<bits/stdc++.h>
#define int   long long
#define ll long long
#define down cout<<'\n';
#define debug cout<<" cucuucucuuu",down
#define NHP     ios_base::sync_with_stdio(0);cout.tie(0);cin.tie(0);
#define modwwe  int t;cin>>t; while(t--)
#define bit(i,j) (i>>j&1)
#define sobit(a) __builtin_popcountll(a)
#define task "test"
#define fin(x) freopen(x".inp","r",stdin)
#define fou(x) freopen(x".out","w",stdout)
#define pb push_back
#define mask(k) (1<<k)
#define mp make_pair
#define checktime   cerr << (double)clock() / CLOCKS_PER_SEC * 1000  << " ms";
using namespace std;
#define getchar_unlocked getchar

inline int scan()
{
    char c = getchar_unlocked();
    int x = 0;
    while (c < '0' || c > '9')
    {
        c = getchar_unlocked();
    }
    while (c >= '0' && c <= '9')
    {
        x = (x << 1) + (x << 3) + c - '0';
        c = getchar_unlocked();
    }
    return x;
}
void phongbeo();
const int inf = 1e18;
const ll mod2 = 1e9+7;
const int  mod1 = 998244353;
const ll base=67;
int add(int x,int y)
{
    if(x+y>=mod2) x-=mod2;
    if(x+y<0)x+=mod2;
    return x+y;
}
struct icd
{
    long double a;
    int b;
};
struct ib
{
    int a;
    int b;
};
struct ic
{
    int a, b, c;
};
struct id
{
    int a, b, c, d;
};
struct ie
{
    int a, b, c, d, e;

};
int n, m, s1, s2, s4, s3, sf, k, s5, s6, mx, s7, s8, s9, mx2, res, dem2 = 0, dem = 0, s33, dem3, dem4, mid, l2, r2, center;
int  i, s10, s12,k1,k2,k3,s11,lim,w,l,r ;
int kk;
int el = 19;
main()
{
    if(fopen(task".inp","r"))
    {
        fin(task);
        fou(task);
    }
    NHP
    /// cin>>s1;
    // modwwe
    phongbeo();
    // checktime
}
int dp[3001][3001];
int fac[6001];
int invf[6001];
int ipow(int x,int y)
{
    if(y==0) return 1;
    int res=ipow(x,y/2);
    res=res*res%mod2;
    if(y%2==1) res=res*x%mod2;
    return res;
}
int rmin,rmax=0,cmin,cmax=0;
void phongbeo()
{
    cin>>n>>m;
    cin>>k;
    invf[0]=1;
    fac[0]=1;
    for(int i=1; i<=6000; i++)
        fac[i]=fac[i-1]*i%mod2,
               invf[i]=ipow(fac[i],mod2-2);
    rmin=n+1;
    cmin=m+1;
    for(int i=1; i<=k; i++)
    {
        cin>>l>>r;
        rmin=min(rmin,l);
        rmax=max(rmax,l);
        cmin=min(cmin,r);
        cmax=max(cmax,r);
    }
    s6=(rmax-rmin+1)*(cmax-cmin+1);
    s6-=k;
    s5=1;
    for(int i=1; i<=s6; i++)
        s5=s5*i%mod2;
    s6=s5;
    s2=n-rmax+rmin-1;
    s4=fac[s2]*invf[n-rmax]%mod2*invf[rmin-1]%mod2;
    s3=m-cmax+cmin-1;
    s4=s4*(fac[s3]*invf[m-cmax]%mod2*invf[cmin-1]%mod2)%mod2;
    dp[0][0]=1;
    s4=s4*s6%mod2;
    for(int i=0; i<=s2; i++)
    {
        for(int j=1; j<=s3; j++)
            dp[i][j]=add(dp[i][j],dp[i][j-1]*fac[rmax-rmin+1+i]%mod2);
        if(i!=s2)
        {
            for(int j=0; j<=s3; j++)
                dp[i+1][j]=add(dp[i+1][j],dp[i][j]*fac[cmax-cmin+j+1]%mod2);
        }
    }
    cout<<dp[s2][s3]*s4%mod2;
}