DP——51nod1486 大大走格子

题面：51nod1486
至于CF原题是什么，我们不去管它
这题的DP思路很有趣。首先如果没有不能走的格子的话，$n*m$的棋盘的走法数就是$C_{n+m-2}^{m-1}$，因为通过转移方程可以很直观地发现是一个杨辉三角
现在考虑到有不能走的格子的问题，我们发现只有这个格子左上的不能走的格子会对走法数产生影响，所以我们只要考虑这些格子就可以了。
然而我们不能直接枚举行和列，所以一个很显然的方案是枚举不能走的格子
状态：$f[i]$表示不经过任何不能走的格子走到这个“不能走“的格子的走法数
我们把右下角的点也当做是一个“不能走”的格子，这样就可以DP了
对不能走的格子进行双目排序之后，我们可以得到一个转移方程

$$f[i]=C_{r[i]+c[i]-2}^{c[i]-1}-{\sum_{j=1}^{i-1}f[j]*C_{r[i]-r[j]+c[i]-c[j]}^{c[i]-c[j]}(r[j]<r[i]且c[j]<c[i])}$$
什么意思呢？前一个很好理解吧，后面是把有不能走的格子的减掉，而且这样是不重复不遗漏的
因为到 $j$这个不能走的格子前没有经过任何的不能走的点，后面的话随你怎么走都可以。统计的时候呢这个格子只会被同样在它之前的格子计数，而不会被后面的计数，这样就可以保证不重复不遗漏了（他们说这个是容斥，我觉得就是普通的去重而已）
最后答案就是$f[最后一个点]$了
因为这题要取模，所以要用到逆元。我是在求组合数的时候暴力求逆元，时间复杂度 $O(n^2logn)$，需要卡常。如果直接把逆元筛出来的话就是 $O(n^2)$了

#include <stdio.h>
#include <algorithm>
#define int long long
using namespace std;
const int MOD=1e9+7;
inline int read(){
    int k=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){k=k*10+ch-'0';ch=getchar();}
    return k*f;
}
struct ppap{int x,y;}a[100010];
inline bool cmp(ppap a,ppap b){return a.x==b.x?a.y<b.y:a.x<b.x;}
int n,m,p,jc[200010],f[200010];
inline int mi(int a,int b){
    int x=1,y=a;
    while(b){
        if(b&1)x=x*y%MOD;y=y*y%MOD;b>>=1;
    }
    return x;
}
inline int C(int x,int y){return jc[x+y]*mi(jc[y],MOD-2)%MOD*mi(jc[x],MOD-2)%MOD;}//这里实际上在求C(x+y,y)
signed main()
{
    n=read();m=read();p=read();
    for(int i=1;i<=p;i++)a[i].x=read(),a[i].y=read();
    a[++p].x=n;a[p].y=m;
    sort(a+1,a+p+1,cmp);
    jc[0]=1;for(int i=1;i<=200010;i++)jc[i]=jc[i-1]*i%MOD;
    for(int i=1;i<=p;i++){
        f[i]=C(a[i].x-1,a[i].y-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]-f[j]*C(a[i].x-a[j].x,a[i].y-a[j].y)%MOD+MOD)%MOD;
    }
    printf("%lld",f[p]);
    return 0;
}