HDU 3944 DP？（Lucas求组合排列）

题目大意：给定一个坐标n，k，求这个坐标到0，0，路过的数字的最小和。

思路：因为只能向下走和向右下走，所以只有两种情况。

要么一开始向下走，然后向右下走。

要么一开始往右下走，然后垂直往下走。

其实两种走法都一样，因为杨辉三角是对称的 ，所以我们只考虑第一种情况。第二种情况（2*k>n)的时候，k = n-k就可以归为第一种情况。

然后我们得到表达式n−k+C(n,k)+C(n−1,k−1)+……+C(n−k+1,1)+C(n−k,0)

又因为化简公式
$$C_n^r=C_{n-1}^r+C_{n-1}^{r-1}$$
得到最后公式为n-k+C（n+1，k)。

又由于数据量很大，如果所有阶乘都存下来的话肯定爆内存，所以先把素数筛出来，单独求阶乘就好了。

因为n,k很大，mod很小，在1e5以内，所以使用Lucas算法。

$$C_n^r\%mod=C_{n/mod}^{r/mod}\ast C_{n\%mod}^{r\%mod}\%mod$$

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

typedef long long ll;
const ll maxn = 2000;
ll Jc[maxn][10000+10];
int prime[maxn];
int len = 0;
void calJc()    //求maxn以内的数的阶乘
{
    len = 0;
    for(int i=2;i<10000;i++)
    {
        int flag = 1;
        for(int j = 2;j*j<=i;j++)
        {
            if(i %j == 0)
            {
                flag = 0;
                break;
            }
        }
        if(flag)prime[len++] = i;
    }

    for(int i = 0;i<len;i++)
            Jc[i][0] = Jc[i][1] = 1;
    for(int i = 0;i<len;i++)
    {
        for(int j = 2;j<prime[i]+7;j++)
        {
            Jc[i][j] = Jc[i][j-1] *j %prime[i];
        }
    }
}
ll pow(ll a, ll n, ll mod)    //快速幂 a^n % p
{
    ll ans = 1;
    while(n)
    {
        if(n & 1) ans = ans * a % mod;
        a = a * a % mod;
        n >>= 1;
    }
    return ans;
}
ll niYuan(ll a,ll mod)   //费马小定理求逆元
{
    return pow(a, mod - 2, mod);
}
ll C(ll n, ll r,ll mod)    //计算C(n, r)
{
    if(r>n)return 0;
    int local = lower_bound(prime,prime+len,mod)-prime;
    return Jc[local][n] * niYuan(Jc[local][r], mod) % mod
        * niYuan(Jc[local][n - r], mod) ;
}

ll Lucas(ll n,ll r,ll mod)
{
    if(n<mod && r<mod)return C(n,r, mod);
    else return C(n%mod, r%mod, mod)*Lucas(n/mod, r/mod, mod)%mod;
}
int main()
{
    ll n,k;
int mod;
    calJc();
    int kase = 0;
    while(scanf("%lld%lld%d",&n, &k, &mod)!=EOF)
    {
        if(2*k >n)
            k = n-k;

        ll ans = n-k + Lucas(n+1,k,mod) ;
        //cout<<Lucas(2,1,2)<<"??"<<endl;
        printf("Case #%d: %d\n",++kase,ans%mod);
    }
    return 0;
}