51nod 1486 大大走格子[容斥原理]

题面传送门

题解

好早以前就考过,当时根本不会

$O(nm)$的dp很简单,但是不够用啊。考虑通过障碍物转移,从$(x1,y1)$到$(x2,y2)$的方案数就是$C(x2-x1,x2-x1+y2-y1)$,就是想想总共走多少步中选几步往下走,把起点终点加到障碍物里一起转移。

先按坐标排个序,设$f[i]$为前i个障碍物只经过第i个的方案数,$f[i]=C(x_i-1,x_i+y_i-2)\Sigma{f[i]*C(x_i-x_j,x_i-x_j+y_i-y_j)}$$(j<i)$。

逆元预处理一下就好

#include<cstdio>
#include<algorithm>
#define P 1000000007
#define LL long long
using namespace std;

int h,w,n,fac[200005],inv[200005];
LL f[2005];
struct node{
    int x,y;
    bool operator < (const node &b) const {
        return x<b.x||(x==b.x&&y<b.y);
    }
}a[2005];

inline void Pre(int n){
    int i;
    for(fac[0]=1,i=1;i<=n;++i) fac[i]=1ll*fac[i-1]*i%P;
    for(inv[0]=inv[1]=1,i=2;i<=n;++i) inv[i]=1ll*(P-P/i)*inv[P%i]%P;
    for(i=1;i<=n;++i) inv[i]=1ll*inv[i]*inv[i-1]%P;
}

inline int C(int x,int y){return 1ll*fac[y]*inv[x]%P*inv[y-x]%P;}

inline int calc(int i,int j){
    return C(a[i].x-a[j].x,a[i].x-a[j].x+a[i].y-a[j].y);
}

int main(){
    scanf("%d%d%d",&h,&w,&n);
    Pre(h+w);
    for(int i=1;i<=n;++i) scanf("%d%d",&a[i].x,&a[i].y);
    a[++n]=(node){h,w};a[0]=(node){1,1};
    sort(a+1,a+n+1);
    for(int i=1;i<=n;++i){
        f[i]=calc(i,0);
        for(int j=1;j<i;++j){
                if(a[j].y<=a[i].y)
                f[i]=(1ll*f[i]-f[j]*calc(i,j)%P+P)%P;
        }
    }
    printf("%lld\n",f[n]);
    return 0;
}