首页 > 代码库 > 组合数学lucas定理 BZOJ2982 combination

组合数学lucas定理 BZOJ2982 combination

2982: combination

Time Limit: 1 Sec  Memory Limit: 128 MB
Submit: 597  Solved: 357
[Submit][Status][Discuss]

Description

LMZn个不同的基友,他每天晚上要选m个进行[河蟹],而且要求每天晚上的选择都不一样。那么LMZ能够持续多少个这样的夜晚呢?当然,LMZ的一年有10007天,所以他想知道答案mod 10007的值。(1<=m<=n<=200,000,000)

Input

  第一行一个整数t,表示有t组数据。(t<=200)
  接下来t行每行两个整数n, m,如题意。

Output

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

Sample Input

4
5 1
5 2
7 3
4 2

Sample Output

5
10
35
6

HINT

Source

此题是lucas定理的裸题

lucas定理的证明改日贴上来

 

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 using namespace std;
 6 const int mod=10007;
 7 int t,n,m;
 8 int qpow(int a,int b){
 9     int ans;
10     for(ans=1;b;b>>=1,a=a*a%mod)
11         if(b&1) ans=ans*a%mod;
12     return ans;
13 }
14 int zuhe(int nn,int mm){
15     if(nn<mm) return 0;
16     if(mm>nn-mm) mm=nn-mm;
17     long long s1=1,s2=1;
18     for(int i=0;i<mm;i++){
19         s1=s1*(nn-i)%mod;
20         s2=s2*(i+1)%mod;
21     }
22     return s1*qpow(s2,mod-2)%mod;//这里在int范围内 
23 }
24 int lucas(int nn,int mm){
25     if(!mm) return 1;
26     return zuhe(nn%mod,mm%mod)*lucas(nn/mod,mm/mod)%mod;
27 }
28 int main(){
29     scanf("%d",&t);
30     while(t--){
31         scanf("%d%d",&n,&m);
32         printf("%d\n",lucas(n,m));
33     } 
34     return 0;
35 }

 

组合数学lucas定理 BZOJ2982 combination