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

题目描述

给你三个整数 a,b,pa,b,p，求 a^b \bmod pabmodp。

输入格式

输入只有一行三个整数，分别代表 a,b,pa,b,p。

输出格式

输出一行一个字符串 a^b mod p=s，其中 a,b,pa,b,p 分别为题目给定的值， ss 为运算结果。

输入输出样例

输入 #1复制

2 10 9

输出 #1复制

2^10 mod 9=7

说明/提示

样例解释

2^{10} = 1024210=1024，1024 \bmod 9 = 71024mod9=7。

数据规模与约定

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

对于这个题，首先，我们先了解一个式子

(a*b)%p=(a%p*b%p)%p

 通过这个式子，我们可以将一个数模P不断拆分

一.递归

第一种方式，构造出这样的式子，将数不断分小，最容易想到的便是递归的写法。

#include<iostream>
using namespace std;
typedef long long ll;
ll a = 0;
ll b = 0;
ll p = 0;
ll dg(ll b) {
	if (b <= 0) {
		return 1;
	}
	ll ans = dg(b / 2) % p;
	if(b&1){
		return (ans * ans * a) % p;
	}
	else {
		return (ans * ans) % p;
	}
}
int main() {
	cin >> a >> b >> p;
	cout << a << "^" << b << " mod " << p << "=" << dg(b);
	return 0;
}

代码详解：首先，递归的内容为dg(b/2),也就是说，每一次递归我都将数分为一半

例如：2^10%9=(2^5%9*2^5%9)%p,递归将2^5传下去，这样一直传到b=0，a^0=1,返回一，接着，一步步返回，一步步计算。

    ll ans = dg(b / 2) % p;是我们需要返回的一部分，但是，b这个数并不一定是2的倍数，那么这时候，我们就要多一步处理，2^5%9=(2^2%9*2^2%9)%p这样显然不对，于是，我们可以这样来写，当我们的b为奇数的时候，我们可以将其分解为(2)^5=(2)^4*(2)^4*2,提了一个a出来，因此，在返回值上，如果此时的b为偶数，那么说明下一个递归返回的值是没有问题的，不需要处理，而如果b为奇数，那么传下去的b是有问题的，就要在返回的时候乘上一个底数。

二.快速幂

首先，我们来写一个快速幂。

ll qp(ll a, ll b)
{
    ll base = a;
    ll ans = 1;
    while (b > 0) {
        if (b & 1) {
            ans *= base;
        }
        base *= base;
        b >>= 1;
    }
    return ans;
}

快速幂，本质上与递归类似，依旧是拆分，不过拆分的方式不一样，递归写法是直接对b去一半，通过其他步骤解决b不能等分的问题，而快速幂则是将b看做一个二进制数，通过将二进制拆分为加法，在将其变成乘法，再套用公式。

例如：

ll qp(ll a, ll b)
{
    ll base = a;
    ll ans = 1;
    while (b > 0) {
        base *= base;
    }
}

首先，这一个while循环中，设置了一个base变量，它的目的是将2的n次方全部遍历一遍，base*=base，就是a->a^2->a^4->a^8......即a^2^0->a^2^1->a^2^3......这样就得到了一个个我们需要的基数。

ll qp(ll a, ll b)
{
    ll base = a;
    ll ans = 1;
    while (b > 0) {
        base *= base;
        b >>= 1;
    }
    return ans;
}

看看b中是否有这个数，b>>=1,即二进制的数字向右移一位如10的二进制1010，第一次我们会发现，它没有2^0,向后移变成0101，此时有了2^1。

ll qp(ll a, ll b)
{
    ll base = a;
    ll ans = 1;
    while (b > 0) {
        if (b & 1) {
            ans *= base;
        }
        base *= base;
        b >>= 1;
    }
    return ans;
}

在b右移的过程中，我们会发现我们只需要含1的base，一次每一次都要判断是不是有这个1，如果有，就乘上这个base，达到最终下面的结果

而我们此题目所要求的是(a*b)%p=(a%p*b%p)%p

所以，我们只需要在每次得到新base和得到新ans之后加上%p,就可以了

ll qp(ll a, ll b, ll p)
{
    ll base = a;
    ll ans = 1;
    while (b > 0) {
        if (b & 1) {
            ans *= base;
            ans %= p;
        }
        base *= base;
        base %= p;
        b >>= 1;
    }
    return ans;
}

值得注意的是，这样%p的方式并不完全符合我们的公式，它这样%p的个数最终会超过我们公式中%p的数量，但模的p更多并不会影响答案。例如：7%9%9%9%9......最终的结果都是7，不会有改变。