三角阵中的数(lucas)

## 三角阵中的数

问题描述

输入
有多组测试数据， 每组数据占一行， 其上有两个整数m 、n
（1<=m<=n<=10^18）。

输入样例
2 3
3 5
输出样例
3
10

思路：
明显规律题，规律是第m行的n个数的值是组合数Cn m ,就是组合数取模即可。
知道组合数取模后就可以快乐上lucas的板子了。

AC code：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll mod = 10007;

ll fast_pow(ll a,ll b,ll p)
{
    ll ans = 1;
    a %= p;
    while(b > 0)
    {
        if(b & 1) ans = (ans * a) % p;
        a = (a * a) % p;
        b >>= 1;
    }
    return ans;
}

ll comb(ll n,ll m,ll p)
{
    if(n < m){
        return 0;
    }
    if(n == m) return 1;
    if(m > n - m) m = n - m;
    ll up = 1,down = 1;
    for(ll i = 1;i <= m; ++i)
    {
        up = (up * (n - i + 1)) % p;
        down = (down * i) % p;
    }
    return (up * fast_pow(down, p - 2,p)) % p;
}

ll lucas(ll n, ll m, ll p)
{
    ll ans = 1;
    while(m > 0){
        ans = ans * comb(n % p, m % p, p) % p;
        n /= p, m /= p;
    }
    return ans;
}

int main()
{
    ll n,m;
    while(cin >> m >> n)
    {
        cout << lucas(n, m, mod) << endl;
    }
    return 0;
}