最近遇到的各种DP 2.0

题目链接

wls的每日一题，容斥板子。

http://oj.daimayuan.top/problem/467

题意

思路

$n=10^6$ , 显然我们不能依照$n$来做，考虑根据障碍物$m$进行转移。

我们记录 $dp[i]$ 为 从 位置$（1，1）$到第$ i $个障碍物的合法走法。

记第$i$个障碍物的位置是$(x_i,y_i)$ ,

1.当第$i$个障碍物前没有障碍物时，显然有：

$dp[i]$ = $C_{x_i+y_i-2}^{x_i-1}$

2.假设第$i$个障碍物前有一个障碍物$j$，那么我们就减去障碍物$j$对答案的贡献：

$dp[i]$ = $dp[i]$ - $dp[j]$ * (从$j$号障碍物到$i$号障碍物的走法)

即 $$dp[i]=dp[i]-\sum _{j=1}^{i-1}dp[j]\ast C_{x_i+y_i-x_j-y_j}^{x_i-y_i}$$
$$=C_{x_i+y_i-2}^{x_i-1}-\sum _{j=1}^{i-1}dp[j]\ast C_{x_i+y_i-x_j-y_j}^{x_i-y_i}(其中x_j<=x_i,y_j<=y_i)$$

如图 ， $dp[i]$ = 总的走法 - 蓝色区域走法*绿色区域走法

我想这个题的时候曾经有一个疑问：

“从$j$号障碍物到$i$号障碍物的走法” 这里可能会经过其他障碍物，不会造成重复吗？

想想其实是不会的，我们假设 $j$ 和 $i$ 中间插着一个障碍物 $z$ ，
那么一定不会存在一条 $j-->z-->i$ 的路线被统计两次，
因为$z$会先于$i$被访问，而 $dp[z]$ 中一定不存在经过j的情况。

3.关于答案统计：

这个简单，创造第$n+1$个障碍物在$(n,n)$处 , 答案即 $dp[n+1]$

代码：

#include<bits/stdc++.h>
using namespace std;
#define int long long
const int N = 2e6+100;
const int M = 3030;
const int mod = 1e9+7;
int fac[N],inv[N];
int dp[M];
struct node{
    int x;
    int y;
    bool operator< (node B)const{
        if(x==B.x) return y<B.y;
        return x<B.x;
    }
};
vector<node> v;
int ksm(int a,int b,int p){
    int ans = 1;
    while(b){
        if(b&1) ans = ans*a%p;
        a = a*a%p;
        b >>= 1;
    }
    return ans%p;
}
void pre(){
    const int mx = 2e6;
    fac[0] = 1;
    for(int i=1;i<=mx;i++){
        fac[i] = fac[i-1]*i%mod;
    }
    inv[mx] = ksm(fac[mx],mod-2,mod)%mod;  // 逆元
    for(int i=mx-1;i>=0;i--){
        inv[i] = inv[i+1]*(i+1)%mod;
    }
}
int C(int n,int m){
    return fac[n]*inv[n-m]%mod*inv[m]%mod;   // 组合数
}
signed main(){
    ios::sync_with_stdio(false);
    pre();    // 预处理逆元和累乘
    int n,m;
    cin>>n>>m;
    for(int i=1;i<=m;i++){
        int x,y;
        cin>>x>>y;
        v.push_back({x,y});
    }
    v.push_back({n,n});    // 创造第n+1个障碍物在(n,n)处 
    sort(v.begin(),v.end());
    for(int i=0;i<v.size();i++){
        node now = v[i];
        dp[i] = C(now.x+now.y-2,now.x-1)%mod;
    }
    for(int i=0;i<v.size();i++){
        node now = v[i];
        for(int j=0;j<i;j++){
            if(v[j].x<=v[i].x&&v[j].y<=v[i].y){   // 保证j在i前面
                node tmp = v[j];
                int cnt = C(now.x-tmp.x+now.y-tmp.y,now.y-tmp.y)%mod;
                dp[i] = (dp[i] - cnt*dp[j]%mod+mod)%mod;
            }
        }
    }
    cout<<dp[v.size()-1]<<endl;
}

上课了，待更

…