AcWing - 求组合数 III(lucas&逆元)
发布日期:2021-07-01 00:21:42 浏览次数:3 分类:技术文章

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

题目链接:

时/空限制:1s / 64MB

题目描述

给定n组询问,每组询问给定三个整数a,b,p,其中p是质数,请你输出C_a^b\ mod\ 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

解题思路

题意:C_a^b\ mod\ p的值。

思路:

  • 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)
    C_n^m\ \%\ p=C_{n/p}^{m/p}*C_{n\%p}^{m\%p}\ \%\ p

Accepted Code:

/*  * @Author: lzyws739307453  * @Language: C++  */#include 
using 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 如侵犯您的版权,请留言回复原文章的地址,我们会给您删除此文章,给您带来不便请您谅解!

上一篇:AcWing - 求组合数 IV(分解质因数)
下一篇:AcWing - 求组合数 II(预处理&逆元)

发表评论

最新留言

初次前来,多多关照!
[***.217.46.12]2024年04月15日 04时56分21秒