51nod 1486

http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1486

1486 大大走格子

题目来源： CodeForces

基准时间限制：1 秒 空间限制：131072 KB 分值: 160 难度：6级算法题

收藏

关注

有一个h行w列的棋盘，里面有一些格子是不能走的，现在要求从左上角走到右下角的方案数。

Input

单组测试数据。
第一行有三个整数h, w, n(1 ≤ h, w ≤ 10^5, 1 ≤ n ≤ 2000)，表示棋盘的行和列，还有不能走的格子的数目。
接下来n行描述格子，第i行有两个整数ri, ci (1 ≤ ri ≤ h, 1 ≤ ci ≤ w)，表示格子所在的行和列。
输入保证起点和终点不会有不能走的格子。

Output

输出答案对1000000007取余的结果。

Input示例

3 4 2
2 2
2 3

Output示例

2

 跟前面做过的很类似，但是增加了一些障碍格子，很闹心啊。。。容斥也不太懂，看了dp的作法还挺好懂的把，写的太丑被卡时一早上。。。

我们不妨计算出所有的方案然后减去有障碍格子的方案，令f[i]表示从(1,1)到第i个障碍点的合法路径，也就是说除了i点，路径上得点都合法。

这样的话不难得出方程  f[i]=C(xi+yi-2,xi-1)-Σ_j=1^i-1f[j]*C(xi+yi-xj-yj,xi-xj) ;

这个方程相当于枚举了所有不合法的路径中第一个障碍格子，显然之后的所有路径都至少包含一个障碍格，所以都是非法的，

这个方程前面就是所有的路径方案，后面就是所有的包含障碍的路径方案，减去就是要求的状态了，

显然对于(a,b)-(c,d)格子间所有的路径数我们根据组合数能轻易的得到，打表就好了，注意并不是i之前的所有点都一定能到达i，记得判断一下。

如果把(h,w)看做是一个最后的障碍点的话，ans=f[n+1].

 #include<bits/stdc++.h>
 using namespace std;
 #define LL long long
 LL mod=1e9+;
 LL dp[];
 LL f[],inv[];
 struct node{int x,y;}P[];
 int read(){
     int x=;char c=getchar();
     while(!isdigit(c)) c=getchar();
     while(isdigit(c)) x=x*+c-'',c=getchar();
     return x;
 }
 bool cmp(node A,node B)
 {
     if(A.x==B.x) return A.y<B.y;
     return A.x<B.x;
 }
 LL qpow(LL a,LL b)
 {
     LL r=;
     while(b){
         if(b&) r=r*a%mod;
         a=a*a;
         b>>=;
     }
     return r;
 }
 inline LL C(LL N,LL M) {return f[N]*inv[M]%mod*inv[N-M]%mod;}
 int main()
 {
     LL h,w,n,i,j,k;
     scanf("%lld%lld%lld",&h,&w,&n);
     for(i=;i<=n;++i)  P[i].x=read(),P[i].y=read();
     sort(P+,P++n,cmp);
     P[n+].x=h; P[n+].y=w;
     inv[]=f[]=;
     for(i=;i<=h+w;++i)
     {
         f[i]=i*f[i-]%mod;
         inv[i]=qpow(f[i],mod-);
     }
     dp[]=C(P[].x+P[].y-,P[].x-);
     for(i=;i<=n+;++i)
     {
         LL s=;
         for(j=;j<i;++j)
         {
           if(P[j].x<=P[i].x&&P[j].y<=P[i].y)
           s=(s+dp[j]*C(P[i].x+P[i].y-P[j].x-P[j].y,P[i].x-P[j].x))%mod;
         }
         dp[i]=(C(P[i].x+P[i].y-,P[i].x-)-s+mod)%mod;
     }
     printf("%lld\n",dp[n+]);
     return ;
 }