lucas定理裸题。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define maxn 20050
#define mod 10007
using namespace std;
long long t,n,m,inv1[maxn],inv2[maxn];
long long f_pow(long long x,long long y)
{
long long ans=,base=x;
while (y)
{
if (y&) ans=(ans*base)%mod;
base=(base*base)%mod;
y>>=;
}
return ans%mod;
}
void get_table()
{
inv1[]=;inv2[]=;
for (long long i=;i<=;i++)
{
inv1[i]=inv1[i-]*i%mod;
inv2[i]=f_pow(inv1[i],mod-);
}
}
long long comb(long long n,long long m)
{
return inv1[n]*inv2[m]%mod*inv2[n-m]%mod;
}
long long lucas(long long n,long long m)
{
if (!m) return ;
return comb(n%mod,m%mod)*lucas(n/mod,m/mod)%mod;
}
int main()
{
get_table();
scanf("%lld",&t);
for (long long i=;i<=t;i++)
{
scanf("%lld%lld",&n,&m);
printf("%lld\n",lucas(n,m));
}
return ;
}

最新文章

  1. LINUX下编译源码时所需提前安装的常用依赖包列表
  2. CSRF防御之ASP.NET MVC
  3. MFC Initinstance中DoModal()返回-1
  4. 【转】基于注解的SpirngMVC简单介绍
  5. Android Log图文详解
  6. Android Study ING
  7. Springmvc构造RESTful详细讲解
  8. Highcharts20151130
  9. Android屏幕保持唤醒状态
  10. 获取input标签的所有属性
  11. Windows phone 之Socket
  12. Matlab绘制三维图形以及提示框
  13. iOS8 用AVAudioPlayer播放音乐(Swift)
  14. chrome console 命令简记
  15. F#中的自定义隐式转换
  16. golang 栈操作
  17. 性能优化——Android图片压缩与优化的几种方式
  18. Day8 Python基础之遗漏知识点(六)
  19. lua 立即执行函数
  20. 2017-2018-2 20165206 实验二《Java面向对象程序设计》实验报告

热门文章

  1. nginx配置 php 单入口
  2. Compound Interest Calculator1.0
  3. alpha发布(技术随笔)
  4. hdu----(3068)最长回文(manacher)
  5. Spring配置文件解析--集合注入方法
  6. Excepion
  7. 小心buffer的拼接问题 --转
  8. K2工作流的使用
  9. async 和await
  10. VBA读取可选择文件夹下可选择txt文件内容