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;
}