51nod1486 大大走格子

题解:

这是一道DP题。
如果没有那个障碍的条件，就是简单的组合数C(n+m-2,n-1)+逆元了。
但是有这个条件，我们不能直接减掉经过那些障碍的方案数，因为会有重复。
仔细思考一下，会发现答案其实就是f[0]-f[1]+f[2]-f[3]+…+f[n]这样的（f[i]表示经过i个障碍物的方案数）但是这样我们就要写一个n^3的DP，复杂度会爆炸。
我们可以换一种想法，f[i]表示从出发点到当前这个障碍，中间不经过其他障碍的方案数。那么我们只要最后把总方案数-(f[1]*c[1]++f[n]*c[n])就好了。c[i]表示从i这个障碍物到终点的方案数。这样写就不会有重复，转移也很容易。
代码如下：

#include<cstdio>
#include<string>
#include<algorithm>
using namespace std;
const int maxn=2005,tt=1000000007;
int h,w,n;
long long ans,inv[1000005],f[maxn];
struct dyt{
    int x,y;
    bool operator <(const dyt &b) const{return x<b.x||(x==b.x&&y<b.y);}
}a[maxn];
long long qsm(long long x,int y){
    long long ret=1ll,p=x;
    while (y){
        if (y%2==1) ret=(ret*p)%tt;
        p=(p*p)%tt; y/=2;
    }
    return ret;
}
long long c(int x,int y) {return inv[x]*qsm(inv[y]*inv[x-y]%tt,tt-2)%tt;}
int main(){
    inv[0]=1ll; for (int i=1;i<=1000000;i++) inv[i]=(inv[i-1]*i)%tt;
    scanf("%d %d %d",&h,&w,&n); ans=c(h+w-2,w-1);
    for (int i=1;i<=n;i++) scanf("%d %d",&a[i].x,&a[i].y);
    sort(a+1,a+1+n);
    for (int i=1;i<=n;i++) {
        f[i]=c(a[i].x+a[i].y-2,a[i].x-1);
        for (int j=1;j<i;j++)
        if (a[i].x>=a[j].x&&a[i].y>=a[j].y)
            f[i]=(f[i]+tt-f[j]*c(a[i].x-a[j].x+a[i].y-a[j].y,a[i].x-a[j].x))%tt;
    }
    for (int i=1;i<=n;i++) ans=(ans+tt-f[i]*c(h-a[i].x+w-a[i].y,h-a[i].x)%tt)%tt;
    printf("%lld\n",ans);
    return 0;
}