P1226 【模板】快速幂||取余运算

最新推荐文章于 2024-10-17 21:05:59 发布

静谧幽蓝_

最新推荐文章于 2024-10-17 21:05:59 发布

阅读量290

点赞数

分类专栏：洛谷题库题目洛谷【模板】文章标签： c++

本文链接：https://blog.csdn.net/ssl_ztc/article/details/120111894

版权

洛谷题库题目同时被 2 个专栏收录

36 篇文章 3 订阅

订阅专栏

洛谷【模板】

15 篇文章 1 订阅

订阅专栏

快速链接

原题链接
题目大意
数据范围
解题思路
上代码

原题链接

P1226
题目类型： ${\color{orange} 普及-}$
AC记录：Accepted

题目大意

给你三个数 $a, b, p$ ，求出 $a^b\ \%\ p$ 的结果。
输出格式为a^b mod p=s，其中 $a, b, p$ 如上， $s$ 为计算的结果。
$\mathbf{Sample}$ $\mathbf{Input}$

2 10 9

$\mathbf{Sample}$ $\mathbf{Output}$

2^10 mod 9=7

$\mathbf{Hint\&Explain}$
$2^{10}=1024，1024\mod 9=7$

数据范围

对于 $100\%$ 的数据，保证 $0\le a,b<2^{31}，a+b>0，2\le p <2^{31}$ 。

解题思路

看看数据，直接做肯定会溢出，正确的做法是用二进制来做。

我们可以把 $b$ 化为二进制，将原本的 $a^b$ 看成若干个互不相同的 $a$ 的 $2^k$ 次方相乘得到的结果，如：
$\ \ \ \ \ \ \ \ \ \ \ \ (10)_{10}=(1010)_2\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ \ \\ \ \ \ \ \ \ \ \ \ \ \ \ \ \ a^{10}= ( 1\times a^8 ) \times ( 0\times a^4 ) \times ( 1\times a^2 ) \times ( 0\times a^1 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \\ =a^8 \times a^2 \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$

由于每一笔运算的结果都有一个 $p$ 在限制运算结果的上限，所以在这里的每一笔运算都不会溢出。

但这里还有一个问题， $b$ 的最大值为 $2^{31}-1$ ，如果手动存储 $2^k$ 得到的结果，再运算 $a^{2^k}$ 的值还是拿不到满分，因为最大 $a^k$ 的值为 $a^{2^{30}}$ ，需要运行 $2^{30}$ 次，这个时间肯定会超时，所以这里有一个小技巧：
$\text{我们已知：}a^b \times a^c=a^{b+c} \\ \text{那么}\left( a^{2^i} \right)^2=\left( a^{2^i} \right) \times \left( a^{2^i} \right)=a^{2^i+2^i}=a^{2^i\times2}=a^{2^{i+1}}$
由上面的公式，我们可以只存储 $a$ 的值，当进入到下一位的时候，我们就将 $a$ 设为 $a\times a \ \% \ p$ ，这样就可以解决了。

最后，祝大家早日
请添加图片描述

上代码

#include<iostream>

using namespace std;

typedef unsigned long long ull;
typedef long long ll;

void POWER_AND_PRINT(ull a,ull b,ull p)
{
	ull initA=a,initB=b;
	
	ull ans=1%p;
	
	a%=p;
	
	for(; b; b>>=1)
	{
		if(b&1) ans=(ull)ans*a%p;
		
		a=(ull)a*a%p;
	}
		
	cout<<initA<<"^"<<initB<<" mod "<<p<<"="<<ans<<endl;
}

int main()
{
    ull n,m,k;
    
    cin>>n>>m>>k;

    POWER_AND_PRINT(n,m,k);

    return 0;
}