bzoj——2982: combination

最新推荐文章于 2022-04-04 11:42:22 发布

weixin_33755557

最新推荐文章于 2022-04-04 11:42:22 发布

阅读量51

点赞数

原文链接：http://www.cnblogs.com/z360/p/7359869.html

版权

2982: combination

Time Limit: 1 Sec Memory Limit: 128 MB
Submit: 611 Solved: 368
[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

思路：

卢卡斯定理裸题

代码：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#define N 10008
#define ll long long 
#define mod 10007
using namespace std;
ll t,n,m,ans,f[N];
ll read()
{
    ll x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1; ch=getchar();}
    while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();}
    return x*f;
}
ll qpow(ll n,ll k)
{
    ll res=1;
    while(k)
    {
        if(k&1) res=res*n%mod;
        n=n*n%mod; k>>=1;
    }return res;
}
ll c(ll n,ll m)
{
    if(m>n) return 0;
    return f[n]*qpow(f[n-m]*f[m],mod-2)%mod;
}
ll lus(ll n,ll m)
{
    if(m==0) return 1;
    return c(n%mod,m%mod)*lus(n/mod,m/mod)%mod;
}
int main()
{
    t=read();f[0]=1;
    for(int i=1;i<=mod;i++)
     f[i]=f[i-1]*i%mod;
    while(t--)
    {
        n=read(),m=read();
        ans=lus(n,m);
        printf("%lld\n",ans);
    }
    return 0;
}