codeforces 560E dp+组合数学

最新推荐文章于 2024-04-19 20:40:32 发布

three trees

最新推荐文章于 2024-04-19 20:40:32 发布

阅读量178

点赞数

CC 4.0 BY-SA版权

分类专栏：组合数学动态规划

本文链接：https://blog.csdn.net/qq_43101466/article/details/102786828

组合数学同时被 2 个专栏收录

13 篇文章

订阅专栏

动态规划

3 篇文章

订阅专栏

本文介绍了一种在二维坐标系中，从起点(1,1)到终点(h,w)，避开坏点，只允许向上、向右或斜向上的路径计数算法。通过动态规划和组合数学方法，计算出所有可能路径的数量。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

题意

在一个二维坐标系中，起点在 $(1, 1)$ ，终点在 $(h, w)$ ，期间只能向上，向右，斜向上
并且不能经过坏点的方案数。坏点有n个

分析

我们分析两个坏点之间的方案数： $C_{x_i-x_j+y_i-y_j}^{x_i-x_j}$
我们定义 $d p [i]$ 为从 $(1, 1)$ 点到第 $i$ 个坏点之间不经过坏点的方案数
那我们知道 $dp[0]=C_{x_0-1+y_0-1}^{x_0-1}$

$dp[i]=Cxi−1+yi−1xi−1+∑j=1i−1dp[j]∗Cxi−xj+yi−yjxi−xjdp[i]=C_{x_i-1+y_i-1}^{x_i-1}+\sum_{j=1}^{i-1}dp[j]*C_{x_i-x_j+y_i-y_j}^{x_i-x_j}$

我们对所有的坏点进行排序（点 $(h, w)$ 也算一个坏点），然后跑到最后一个坏点的dp值即为答案

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
template <typename T>
void out(T x) { cout << x << endl; }
ll fast_pow(ll a, ll b, ll p) {ll c = 1; while(b) { if(b & 1) c = c * a % p; a = a * a % p; b >>= 1;} return c;}
ll exgcd(ll a, ll b, ll &x, ll &y) { if(!b) {x = 1; y = 0; return a; } ll gcd = exgcd(b, a % b, y, x); y-= a / b * x; return gcd; }
const int N = 2e3 + 5;
const int NN = 2e5 + 10;
const int mod = 1e9 + 7;
struct node 
{
    int x, y;
    bool operator < (const node b) const 
    {
        if(x == b.x)
            return y < b.y;
        return x < b.x;
    }
}pp[N];
ll fac[NN], inv[NN];
void init()
{
    inv[0] = fac[0] = 1;
    for(int i = 1; i < NN; i ++)
        fac[i] = fac[i - 1] * i % mod;
    inv[NN - 1] = fast_pow(fac[NN - 1], mod - 2, mod);
    for(int i = NN - 2; i >= 1; i --)
        inv[i] = inv[i + 1] * (i + 1) % mod;
}
ll C(int n, int m)
{ 
    return fac[n] * inv[m] % mod * inv[n - m] % mod;
}
ll dp[N];
int main()
{
    ios::sync_with_stdio(false);
    init();
    int h, w, n;
    cin >> h >> w >> n;
    for(int i = 0; i < n; i ++)
        cin >> pp[i].x >> pp[i].y;
    pp[n].x = h;
    pp[n].y = w;
    n ++;
    sort(pp, pp + n);
    dp[0] = C(pp[0].x + pp[0].y - 2, pp[0].x - 1);
    for(int i = 1; i <= n; i ++)
    {
        dp[i] = C(pp[i].x + pp[i].y -2, pp[i].x - 1);
        for(int j = 0; j < i; j ++)
        {
            if(pp[i].x >= pp[j].x && pp[i].y >= pp[j].y)//只能向上，向下或者斜上
                dp[i] = (dp[i] - dp[j] * C(pp[i].x - pp[j].x + pp[i].y - pp[j].y, pp[i].x - pp[j].x) % mod + mod) % mod;
        }
    }
    cout << dp[n - 1] << endl;
}