Codeforces Round #313(div.2) E. Gerald and Giant Chess

题目链接

题目大意：有一个$h*w$的棋盘，一个兵要从$S(1,1)$走到$T(h,w)$，且只能向下或向右走，还有$n$个格子是不能走的，问从$S$走到$T$的路径数。

分析：如果每个格子都能走，那么从$(1,1)$走到$(p,q)$的答案为$C_{p+q-2}^{p-1}$（要向下走$p-1$步，向左走$q-1$步，我们从$p+q-2$个位置中选出$p-1$个向下走，剩下的就是向右走的啦）。

简单题解：
我们设$A_0,A_1,A_2,\cdots ,A_{n-1}$表示排序后的不能走的格子。
设$G(A, B)$表示从$A(x1,y1)$到$B(x2,y2)$的方案数，由上面的分析很容易得到

$$G=C_{x2-x1+y2-y1}^{x2-x1}$$ 设 $F_i$表示走到 $A_i$不经过 $A_0, A_1,\cdots, A_{i-1}$的答案，那么有转移方程
$$F_i =G(S,A_i) - \sum_{j=1}^{i-1} G(A_j,A_i)*F_j$$ 把 $(h,w)$也看成不能走的格子，记为 $A_n$，那么 $F_n$就是答案。

算组合数的话预处理阶乘及逆元即可。

时间复杂度：$O(n^2)$

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
typedef long long LL;

inline int read()
{
    int x = 0, f = 1, t = getchar();
    while(t < 48 || t > 57) t == 45 ? f = -1 : 0, t = getchar();
    while(t > 47 && t < 58) x = (x << 1) + (x << 3) + t - 48, t = getchar();
    return x * f;
}

const int maxn = 2005;
const int maxl = 200000;
const int MOD = 1e9 + 7;

struct Black
{
    int x, y;
    bool operator < (const Black &r) const {
        return x < r.x || (x == r.x && y < r.y);
    }
};
int h, w, n;
Black data[maxn];
LL fac[maxl + 5], inv[maxl + 5];
LL f[maxn];

LL power_mod(LL x, int y)
{
    LL res = 1;
    for(; y; y >>= 1, x = x * x % MOD)
        if(y & 1) res = res * x % MOD;
    return res;
}
void init_fact()
{
    fac[0] = 1; for(int i = 1; i <= maxl; ++i) fac[i] = fac[i - 1] * i % MOD;
    inv[maxl] = power_mod(fac[maxl], MOD - 2); for(int i = maxl; i >= 1; --i) inv[i - 1] = inv[i] * i % MOD;
}
LL C(int x, int y)
{
    if(x < y) return 0;
    return fac[x] * inv[y] % MOD * inv[x - y] % MOD;
}
LL calc(int x, int y)
{
    if(x < 1 || y < 1) return 0;
    return C(x + y - 2, x - 1);
}

int main()
{
#ifndef ONLINE_JUDGE
    freopen("input.txt", "r", stdin);
    freopen("output.txt", "w", stdout);
#endif

    init_fact();

    h = read(), w = read(), n = read();
    for(int i = 0; i < n; ++i) data[i] = (Black){read(), read()};
    data[n] = (Black){h, w};
    sort(data, data + n + 1);

    for(int i = 0; i <= n; ++i)
    {
        LL res = 0;
        for(int j = 0; j < i; ++j)
            res = (res + f[j] * calc(data[i].x - data[j].x + 1, data[i].y - data[j].y + 1)) % MOD;
        f[i] = (calc(data[i].x, data[i].y) - res + MOD) % MOD;
    }
    int ans = f[n];
    printf("%d\n", ans);

#ifndef ONLINE_JUDGE
    fclose(stdin), fclose(stdout);
#endif
    return 0;
}