[BZOJ2982] combination

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 <cstdio>
#include <cstring>
#include <cmath>
#include <queue>
#include <algorithm>
using namespace std;
int t;
int n,m;
int p=10007;
int jie[20000],yuan[20000];
void init()
{
    jie[0]=1;
    for(int i=1;i<p;i++) jie[i]=jie[i-1]*i%p;
    yuan[0]=yuan[1]=1;
    for(int i=2;i<p;i++) yuan[i]=(p-p/i)*yuan[p%i]%p; 
    for(int i=1;i<p;i++) yuan[i]=yuan[i]*yuan[i-1]%p;
}
int lucas(int n,int m)
{
    if(n<m) return 0;
    if(n<p&&m<p) return jie[n]%p*yuan[m]%p*yuan[n-m]%p;
    return lucas(n/p,m/p)*lucas(n%p,m%p)%p;
}
int main()
{
    scanf("%d",&t);
    init();
    while(t--) 
    {
        scanf("%d%d",&n,&m);
        printf("%d\n",lucas(n,m));
    }
    return 0;
}

发表评论

邮箱地址不会被公开。 必填项已用*标注