快速幂

题目1:$a^b\%p$

输入3个数字a, b, p 求a^b%p
数据范围很大，最大到10¹⁸

这个题目核心的就是用到了递推来防止溢出
数字b可以表示成$b={\sum }_{k=0}^{k=n}{d}_k\ast 2^k$,其中$d_k$是二进制表示中第k位的值，$d_k$只有两个取值，0或者1。那么$$a^b=a^{{\sum }_{k=0}^{k=n}{d}_k\ast 2^k}=\prod _{k=0}^{k=n}{a}^{d_k\ast 2^k}=a^{d_0\ast 2^0}\ast a^{d_1\ast 2^1}\ast a^{d_2\ast 2^2}\ast ...\ast a^{d_n\ast 2^n}$$

需要用到的递推公式如下 $$a^{2^{k+1}}=a^{2\ast 2^k}={(a^{2^k})}^2$$这样就可以通过$a^{2^k}$一步求出$a^{2^{k+1}}$(求个平方即可)

因为$d_k$只有0或者1两种可能，当取0的时候对最终的结果没有任何影响。
利用这个递推公式就可以逐项求出$$a^{2^0},a^{2^1},a^{2^2}...,a^{2^n}$$
$a^{2^0}$就是a，其他项依次递推即可。

结算结果的时候，我们只关系$d_k$是1的情况，如果是0，那相当于在结果上乘了一个$a^{0\ast 2^k}=a^0=1$可以直接跳过。

代码如下

#include<iostream>
int a, b, p;

using namespace std;

int main() {
    cin >> a >> b >> p;
    int res = 1 % p;
    while ( b ) {
        if ( b & 1 ) res = res * 1ll * a % p;
        a = a * 1ll * a % p;
        b >>= 1;
    }
    cout << res << endl;
    return 0;
}

题目2:$a\ast b\%p$

输入3个数字a, b, p 求a*b%p
思路和上一题一样，只不过这次是
$$a\ast b=a\ast \sum _{k=0}^{k=n}{d}_k\ast 2^k=\sum _{k=0}^{k=n}a\ast d_k\ast 2^k=a\ast d_0\ast 2^0+a\ast d_1\ast 2^1+a\ast d_2\ast 2^2+...+a\ast d_n\ast 2^n$$
递推式 $$a\ast 2^{k+1}=2\ast (a\ast 2^k)$$

#include<iostream>

using namespace std;
typedef unsigned long long ull;
ull a, b, p;

int main() {
    cin >> a >> b >> p;
    ull res = 0;
    while ( b ) {
        if ( b & 1 ) res = (res + a) % p;
        a = a * 2 % p;
        b >>= 1;
    }
    cout << res << endl;
    return 0;
}

代码中变量a存的是递推式算出来的值，即$(a\ast 2^k)\%p$(这里a是题目给的数字)