解析:
计数类DP。
首先明显不可过。。。
考虑用总方案数减去至少经过一个黑格子的路线数量即为答案。
重点是如何做到不重不漏。把所有格子按照横纵坐标排序。假设左上角是第 0 个黑格子,右下角是第 n + 1 个格子。设 f [ i ] 表示从左上角走到排序后的第 i 个点,并且中途不经过其他黑格子的路线数。则有:
枚举的 j 作为从左上角到第 i 个点经过的第一个黑格子,总路线数为 f [ i ] ,然后到第 i 个点就可以任意走,再根据加法原理以及乘法原理可得所有非法方案,并且可以发现第一个经过得黑格子不同,则不同得 j 对应得路线一定不同,而 j 又取遍了所有 i 之前得黑格子,所有没有遗漏。
代码:
#include <bits/stdc++.h>
#define int long long
#define x first
#define y second
using namespace std;
const int mod=1e9+7;
const int Max=2010;
int n,m,k;
int mul[200005],c[200005],f[Max];
pair<int,int>num[Max];
inline int get_int()
{
int x=0,f=1;
char c;
for(c=getchar();(!isdigit(c))&&(c!='-');c=getchar());
if(c=='-') f=-1,c=getchar();
for(;isdigit(c);c=getchar()) x=(x<<3)+(x<<1)+c-'0';
return x*f;
}
inline int ksm(int a,int b)
{
int ans=1;
while(b)
{
if(b&1) ans=(ans*a)%mod;
b>>=1;
a=(a*a)%mod;
}
return ans;
}
inline int C(int n,int m)
{
if(n>=0&&m>=0&&n>=m) return (mul[n]*c[n-m]%mod*c[m]%mod)%mod;
return 0;
}
signed main()
{
n=get_int(),m=get_int(),k=get_int();
mul[0]=f[0]=c[0]=1;
for(int i=1;i<=200000;i++) mul[i]=(mul[i-1]*i)%mod,c[i]=ksm(mul[i],mod-2);
for(int i=1;i<=k;i++)
{
int x=get_int(),y=get_int();
num[i]=make_pair(x,y);
}
sort(num+1,num+k+1);
num[k+1]=make_pair(n,m);
for(int i=1;i<=k+1;i++)
{
f[i]=C(num[i].x+num[i].y-2,num[i].x-1);
for(int j=1;j<i;j++)
{
if(num[j].x>num[i].x||num[j].y>num[i].y) continue;
f[i]=(f[i]-f[j]*C(num[i].x+num[i].y-num[j].x-num[j].y,num[i].x-num[j].x))%mod;
}
}
cout<<(f[k+1]%mod+mod)%mod;
return 0;
}