xdoj-1032-找规律

最新推荐文章于 2024-12-04 16:31:55 发布

Ares2333

最新推荐文章于 2024-12-04 16:31:55 发布

阅读量457

点赞数

CC 4.0 BY-SA版权

分类专栏： Lucas 快速幂大组合数

本文链接：https://blog.youkuaiyun.com/Ares2333/article/details/53102809

Lucas 同时被 3 个专栏收录

1 篇文章

订阅专栏

快速幂

1 篇文章

订阅专栏

大组合数

1 篇文章

订阅专栏

本文介绍了一种高效计算快速幂的方法，并在此基础上实现了组合数的计算，还提供了一个求解 Lucas 定理的实现。这些算法对于解决涉及大量整数运算的问题非常有用。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

#include<cstdio>
typedef long long ll;
const ll mod=10007;

ll feima(ll a,ll b)
{
    ll c=1;
    while(b)
    {
        if(b&1)
            c=c*a%mod;
        a=a*a%mod;
        b=b>>1;
    }
    return c;
}
ll C(ll a,ll b)
{
    if(a<b) return 0;
    if(a==b)    return 1;
    if(b>a-b)   b=a-b;
    ll ca=1,cb=1;
    for(ll i=0;i<b;i++)
    {
        ca=ca*(a-i)%mod;
        cb=cb*(b-i)%mod;
    }
    return ca*feima(cb,mod-2)%mod;
}
/*LL Lucas(LL n,LL m)  
{  
    if(m==0) return 1;  
    return C(n%p,m%p)*Lucas(n/p,m/p)%p;  
}*/ 
ll Lucas(ll n,ll m)
{
    ll ans=1;
    while(n&&m&&ans)
    {
        ans=ans%mod*C(n%mod,m%mod)%mod;
        n/=mod;
        m/=mod;
    }
    return ans;
}

int main()
{
    ll n,m;
    while(~scanf("%lld%lld",&n,&m))
        printf("%lld\n",Lucas(m,n));
    return 0;
}