【组合数学 && DP && 向左向上走 (1, 1) 到 (n, m) 的路径数】CodeForces - 559C Gerald and Giant Chess

Step1 Problem:

给你 n*m 的矩阵，矩阵中有 k 个黑色的块，你只能向上或者向右走，求从 (1,1) 走到 (n,m) 不能经过黑色的块的方案数。
数据范围：
1<=h, w<=1e5, 1<=n<=2000.

Step2 Ideas:

学习v5zsq
h,w 很大， n 才 2000，我们可以从黑块入手。
前置技能：(x1, y1) 到 (x2, y2) 的所有路径方案数：x = x2-x1, y = y2-y1, C(x+y, x) 就是方案数。
从 (x1, y1) 到 (x2, y2) 一共走 x+y 步，选择其中任意的 x 步为向右走。
状态 dp[i]：(1, 1) 到第 i 个黑块的合法路径数量，我们设起点和终点也是黑块，那么 dp[n+1] 就是结果。
get_num(i, j)：代表 i 到 j 的所有路径方案数。
dp[i] = get_num(0, i) - sum{dp[j]*num[i][j]}, j = 1,2,…,i-1.

Step3 Code:

#include<bits/stdc++.h>
using namespace std;
#define ll long long
const int N = 2e5+5;
const int M = 2005;
const int MOD = 1e9+7;
struct node
{
    int x, y;
    bool operator < (const node &b) const {
        if(x == b.x) return y < b.y;
        else return x < b.x;
    }
};
node a[N];
ll A[N], inv[N];
ll dp[M];
ll Pow(ll x, int m)
{
    ll ans = 1;
    while(m)
    {
        if(m&1) ans *= x, ans %= MOD;
        x = x*x, x %= MOD;
        m >>= 1;
    }
    return ans;
}
void get_inv()
{
    A[0] = 1;
    inv[0] = 1;
    for(int i = 1; i <= N-5; i++)
    {
        A[i] = A[i-1] * i, A[i] %= MOD;
        inv[i] = Pow(A[i], MOD-2);
    }
}
ll C(int n, int m)
{
    return A[n]*inv[m]%MOD * inv[n-m]%MOD;
}
ll get_num(int x1, int y1, int x2, int y2)
{
    int x = x2-x1;
    int y = y2-y1;
    return C(x+y, x);
}
int main()
{
    int n, m, k;
    get_inv();
    scanf("%d %d %d", &n, &m, &k);
    for(int i = 1; i <= k; i++)
        scanf("%d %d", &a[i].x, &a[i].y);
    a[k+1] = (node){n, m};
    sort(a, a+k+1);
    for(int i = 1; i <= k+1; i++)
    {
        dp[i] = get_num(1, 1, a[i].x, a[i].y);
    }
    for(int i = 1; i <= k+1; i++)
    {
        for(int j = i+1; j <= k+1; j++)
        {
            if(a[i].y > a[j].y) continue;// j 在 i 的右上角才行（包括正上和正右）
            dp[j] = ((dp[j] - (dp[i] * get_num(a[i].x, a[i].y, a[j].x, a[j].y))%MOD)%MOD + MOD)%MOD;
        }
    }
    printf("%lld\n", dp[k+1]);
    return 0;
}