Lucas

最新推荐文章于 2020-09-29 17:00:31 发布

原创最新推荐文章于 2020-09-29 17:00:31 发布 · 494 阅读

0 ·

CC 4.0 BY-SA版权

模板专栏收录该内容

46 篇文章

订阅专栏

#include"iostream"
using namespace std;
const int MOD=1e9+7;
typedef long long LL;
const int maxn=25;
int N,K;
LL mul(LL a,LL b,LL mod)
{
    double x=a;
    LL res=a*b-(LL)(x*b/mod)*mod;
    return (res+mod)%mod;
}
LL ksm(LL a,LL b,LL mod)
{
    LL res=1,base=a;
    while(b)
    {
        if(b&1)res=mul(res,base,MOD);
        b>>=1;
        base=mul(base,base,MOD);
    }
    return res;
}
LL C(LL n,LL m)
{
    if(m<0||m>n)return 0;
    if(m==0)return 1;
    LL res=1,b=1;
    for(LL i=1;i<=m;i++)
    {
        res=mul(res,n-m+i,MOD);
        b=mul(b,i,MOD);
    }
    res=mul(res,ksm(b,MOD-2,MOD),MOD);
    return res;
}
LL Lucas(LL n,LL m)
{
    if(m==0)return 1;
    if(m<0||m>n)return 0;
    return mul(C(n%MOD,m%MOD),Lucas(n/MOD,m/MOD),MOD);
}

int main()
{
    int N,M;
    while(cin>>N>>M)
    {
        cout<<"C="<<C(N,M)<<"\n";
        cout<<"Lucas="<<Lucas(N,M)<<"\n";
    }
}