bzoj2982: combination

最新推荐文章于 2020-07-16 21:43:51 发布

weixin_30699235

最新推荐文章于 2020-07-16 21:43:51 发布

阅读量98

点赞数

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/Bloodline/p/6230947.html

本文介绍了一种利用Lucas定理求组合数C(n,m)模10007的算法实现。通过递归分解问题规模，结合快速幂取模运算，解决了大整数组合数的高效计算问题。

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

Lucas定理模版题

T行，每行一个数，为C(n, m) mod 10007的答案。

 1 #include<bits/stdc++.h>
 2 #define rep(i,l,r) for(int i=l;i<=r;++i)
 3 using namespace std;
 4 typedef long long ll;
 5 const ll p=10007;
 6 ll n,m,T;
 7 ll Pow(int x,int y){
 8     if(!y)return 1;
 9     ll t=Pow(x,y/2)%p; (t*=t)%=p;
10     if(y&1)return t*x;else return t;
11 }
12 ll c(ll n,ll m){
13     if(n<m)return 0;
14     ll ans=1;
15     rep(i,1,m){
16         ans=ans*(((n-m+i)%p)*Pow(i,p-2)%p)%p;  
17     }
18     return ans;
19 }
20 ll lucas(ll n,ll m){
21     if(!m)return 1;
22     return lucas(n/p,m/p)*c(n%p,m%p)%p;
23 }
24 int main(){
25     cin>>T;
26     while(T--){
27         cin>>n>>m;
28         cout<<lucas(n,m)<<endl;
29     }
30 }