AcWing - 求组合数 III(lucas&逆元)
发布日期:2021-07-01 00:21:42
浏览次数:4
分类:技术文章
本文共 1177 字,大约阅读时间需要 3 分钟。
题目链接:
时/空限制:1s / 64MB题目描述
给定n组询问,每组询问给定三个整数a,b,p,其中p是质数,请你输出的值。
输入格式
第一行包含整数n。
接下来n行,每行包含一组a,b,p。
输出格式
共n行,每行输出一个询问的解。
数据范围
1≤n≤20,
1≤b≤a≤10^18, 1≤p≤10^5,
输入样例
3
5 3 7 3 1 5 6 4 13
输出样例
3
3 2
解题思路
题意:求的值。
思路:- 100000组, 1 <= b <= a <= 2000, 递推 -> O(N^2).
- 10000组, 1 <= b <= a <= 10^5, 预处理 -> O(NlogN).
- 20组, 1 <= b <= a <= 10^18, p <= 10^5, 卢卡斯定理(lucas)
Accepted Code:
/* * @Author: lzyws739307453 * @Language: C++ */#includeusing namespace std;typedef long long ll;int Fast_Power(int a, int b, int p) { int res = 1; while (b) { if (b & 1) res = 1ll * res * a % p; a = 1ll * a * a % p; b >>= 1; } return res;}int Com_Num(int a, int b, int p) { int res = 1; for (int i = 1, j = a; i <= b; i++, j--) res = 1ll * res * j % p * Fast_Power(i, p - 2, p) % p; return res;}int lucas(ll a, ll b, int p) { if (a < p && b < p) return Com_Num(a, b, p); return 1ll * Com_Num(a % p, b % p, p) * lucas(a / p, b / p, p) % p;}int main() { int t; scanf("%d", &t); while (t--) { int p; ll a, b; scanf("%lld%lld%d", &a, &b, &p); printf("%d\n", lucas(a, b, p)); } return 0;}
转载地址:https://lzyws739307453.blog.csdn.net/article/details/99713778 如侵犯您的版权,请留言回复原文章的地址,我们会给您删除此文章,给您带来不便请您谅解!
发表评论
最新留言
网站不错 人气很旺了 加油
[***.192.178.218]2024年04月20日 06时11分37秒
关于作者
喝酒易醉,品茶养心,人生如梦,品茶悟道,何以解忧?唯有杜康!
-- 愿君每日到此一游!
推荐文章
软件质量的8个特性
2021-07-03
2021年不可错过的17种JS优化技巧(一)
2021-07-03
在 Vue 中用 Axios 异步请求API
2021-07-03
MySQL进阶查询(SELECT 语句高级用法)
2021-07-03
Mysql 之主从复制
2021-07-03
【NLP学习笔记】中文分词(Word Segmentation,WS)
2021-07-03
对于时间复杂度的通俗理解
2021-07-03
如何输入多组数据并输出每组数据的和?
2021-07-03
行阶梯型矩阵
2021-07-03
MATLAB指定路径保存图片方法
2021-07-03
JAVA学习笔记6 - 数组
2021-07-03
【学习笔记】Android Activity
2021-07-03
location区段
2021-07-03
linux内存的寻址方式
2021-07-03
how2heap-double free
2021-07-03
how2heap-fastbin_dup_consolidate
2021-07-03
tf keras SimpleRNN源码解析
2021-07-03
MyBatisPlus简单入门(SpringBoot)
2021-07-03
攻防世界web进阶区web2详解
2021-07-03