博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
HDU 3037 Saving Beans(Lucas定理的直接应用)
阅读量:6048 次
发布时间:2019-06-20

本文共 943 字,大约阅读时间需要 3 分钟。

解题思路:

直接求C(n+m , m) % p , 由于n , m ,p都非常大,所以要用Lucas定理来解决大组合数取模的问题。

#include 
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#include
#define LL long long#define FOR(i, x, y) for(int i=x;i<=y;i++)using namespace std;const int maxn = 100000 + 10;LL Pow_mod(LL a, LL b, LL mod){ LL ret = 1; while(b) { if(b & 1) ret = ret * a % mod; a = a * a % mod; b >>= 1; } return ret;}LL fac[maxn];void init(LL p){ fac[0] = 1; FOR(i, 1, p) fac[i] = (fac[i-1] * i) % p;}LL Lucas(LL n, LL m, LL p){ LL ret = 1; while(n && m) { LL a = n % p, b = m % p; if(a < b) return 0; ret = (ret * fac[a] * Pow_mod(fac[b] * fac[a-b] % p , p - 2 , p)) % p; n /= p, m /= p; } return ret;}int main(){ int T; scanf("%d",&T); while(T--) { int n, m, p; scanf("%d%d%d", &n, &m, &p); init(p); printf("%d\n", (int)Lucas(n + m, m, p)); } return 0;}

转载地址:http://mkxex.baihongyu.com/

你可能感兴趣的文章
Vbs脚本编程简明教程之十四
查看>>
如何UDP/TCP端口是否通了
查看>>
pxe实现系统的自动化安装
查看>>
Redis高可用技术解决方案总结
查看>>
Scale Out Owncloud 高可用(2)
查看>>
何为敏捷
查看>>
HA集群之四:Corosync+Pacemaker+DRBD实现HA Mysql
查看>>
服务器定义
查看>>
我的友情链接
查看>>
分布式系统的面试题15
查看>>
个人代码库の创建快捷方式
查看>>
由strcat函数引发的C语言中数组和指针问题的思考
查看>>
无锁编程
查看>>
如何在loadrunner中做关联
查看>>
二叉树的六种遍历方法汇总(转)
查看>>
用wxpython制作可以用于 特征筛选gui程序
查看>>
【转载】 [你必须知道的.NET]目录导航
查看>>
数据存储小例
查看>>
C++中构造函数详解
查看>>
电商网站中添加商品到购物车功能模块2017.12.8
查看>>