HDU 3944 DP? （ lucas+预处理

题意：

求从$(0, 0)$开始计算的杨辉三角到给定的点$(a,b)$路径上最小的权值$mod(p)$ 要求只能斜着或者竖着走

题意：

观察一下 我们发现 我们需要尽可能多的走边上 也就是$１$这样　就会只有两种情况可以考虑了
case1:
从斜边上下去　然后竖着到达终点　
$ans = 0$ $ans += m+1$
考虑$ans += C_{m+1}^{m}+C_{m+2}^{m}＋.....+C_{n}^{m}$　
考虑加入$C_{m+1}^{m+1}$　这样就可以将公式合并为$C_{n+1}^{m+1}$
则最终的结果为$ans = m+C_{n+1}^{m+1}$
同理Case2 :
竖直边下去　然后斜着到达终点
$ans = 0, ans += n-m$
考虑$C_{m}^{0}+C_{m+1}^{1}+C_{m+2}^{2}+....+C_{n}^{m}$
$C_{m+1}^{0} == C_{m}^{0}$　所以合并二项式为
$C_{n+1}^{m}$
最终结果为$n-m+C_{n+1}^{m}$
数据量非常大，　我们需要预处理模数的阶乘表，然后使Lucas才能得到结果

#include <bits/stdc++.h>

using namespace std;

#define ll long long 
#define pb push_back
#define ls o<<1
#define rs o<<1|1
#define fi first
#define se second
#define CLR(a, b) memset(a, (b), sizeof(a))
const int INF = 0x3f3f3f3f;
const int mod = 1e9+7;
const int MAXN = 1e4+10;

void F() {
    #ifndef ONLINE_JUDGE
        freopen("in.txt", "r", stdin);
    #endif
}

ll fac[MAXN][MAXN], invv[MAXN][MAXN];
ll prime[MAXN], vis[MAXN];
ll k = 0;
ll pow_mod(ll a, ll b, ll p) {
    ll res = 1;
    while(b) {
        if(b&1) res = res*a%p;
        a = a*a%p;
        b >>= 1;
    }
    return res;
}
ll inv(ll x, ll p) {
    return pow_mod(x, p-2, p);
}
ll C(ll a, ll b, ll p) {
    if(a < b) return 0;
    if(a == b) return 1;
    ll x = vis[p];
    return fac[x][a]*(invv[x][a-b]*invv[x][b]%p)%p;
}
ll lucas(ll a, ll b, ll p) {
    if(a<p && b<p) return C(a, b, p);
    return C(a%p, b%p, p)*lucas(a/p, b/p, p)%p;
}

void pri() {
    CLR(prime, 1);
    for(int i = 2; i <= MAXN; ++i) {
        if(prime[i]) {
            prime[++k] = i;
            vis[i] = k;
            for(int j = i*2; j < MAXN; j+=i) {
                prime[j] = false;
            }
        }
    }
}
void init() {
    pri();
    for(int i = 1; i <= k; ++i) {
        fac[i][0] = 1;
        invv[i][0] = 1;
        for(int j = 1; j < prime[i]; ++j) {
            fac[i][j] = (fac[i][j-1]*j)%prime[i];
            invv[i][j] = inv(fac[i][j], prime[i]);
        }
    }
}

int main() {
    F();
    init();
    ll a, b, p; int kk = 0;
    while(cin >> a >> b >> p) {
        cout << "Case #" << ++kk << ": ";
        if(2*b < a) 
            cout << (lucas(a+1, b, p)+a-b)%p << endl;
        else
            cout << (lucas(a+1, b+1, p)+b)%p << endl;
    }
    return 0;
}