2982: combination

Time Limit: 1 Sec Memory Limit: 128 MB
Submit: 664 Solved: 397
[Submit][Status][Discuss]

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

HINT

Source

組合數裸題，一般我們選用lucas定理實現

人話定義環節——

組合數：從n個數中取m個數，一共的方法數，記作C(n,m)

lucas定理：C(n,m)=C(n%mod,m%mod)*C(n/mod,m/mod)%mod

這裏直接給出實現代碼

#pragma GCC optimize("O2")
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<queue>
#include<stack>
#include<set>
#include<map>
#include<limits.h>
#include<ctime>
#define N 100001
typedef long long ll;
const int inf=0x3fffffff;
const int maxn=2017;
using namespace std;
inline ll read()
{
    ll f=1,x=0;char ch=getchar();
    while(ch>‘9‘|ch<‘0‘)
    {
        if(ch==‘-‘)
        f=-1;
        ch=getchar();
    }
    while(ch<=‘9‘&&ch>=‘0‘)
    {
        x=(x<<3)+(x<<1)+ch-‘0‘;
        ch=getchar();
    }
    return f*x;
}
const int mod=10007;
ll fac[mod<<1],inv[mod<<1];//fac函數存儲n的階乘，inv數組存放逆元的階乘 
void init()
{
	fac[0]=fac[1]=1,inv[0]=inv[1]=1;
	for(int i=2;i<mod;i++)fac[i]=fac[i-1]*i%mod;//預處理階乘 
	for(int i=2;i<mod;i++)inv[i]=(mod-mod/i)*inv[mod%i]%mod;//遞推求逆元 
	for(int i=2;i<mod;i++)inv[i]=inv[i-1]*inv[i]%mod;//預處理逆元的階乘 
}
ll lucas(ll n,ll m)
{
	if(m>n)return 0;
	if(n<=mod&&m<=mod)return fac[n]*inv[m]*inv[n-m]%mod;
	return lucas(n%mod,m%mod)*lucas(n/mod,m/mod)%mod;
}
int main()
{
	int t=read();
	init();
	while(t--)
	{
		ll n=read(),m=read();
		printf("%d\n",lucas(n,m));
	}
}

bzoj2982: combination(lucas定理板子)