【bzoj2982】【combination】【Lucas定理】

Description

LMZ有n个不同的基友，他每天晚上要选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

题解：直接上lucas好了。

代码：

#include<iostream>
#include<cstdio>
#define N 3000010
using namespace std;
long long f[N],s[N],p(10007),x,y;
int n,t;
inline long long power(long long a,long long b){
    long long ans;
    for(ans=1;b;b>>=1,a=a*a%p) if(b&1) ans=ans*a%p;
    return ans;
}
inline long long cal(long long n,long long m){
    long long s1(1),s2(1);
    if (n<m) return 0;
    if (m>n-m) m=n-m;
    for (long long i=0;i<m;i++){(s1*=n-i)%=p;(s2*=i+1)%=p;}
    return s1*power(s2,p-2)%p;
}
inline long long lucas(long long n,long long m){
    if (m==0) return 1;
    else return cal(n%p,m%p)*lucas(n/p,m/p)%p;
}
int main(){
  scanf("%d",&t);
  for (int i=1;i<=t;i++){scanf("%lld%lld",&x,&y);printf("%lld\n",lucas(x,y));}
}