0x00 基本算法

0x01 位运算

前言

在我刚开始刷题的时候，觉得位运算很简单，只不过是一些与或非等等的一些操作罢了。现在看来，是我目光短浅，见识少的缘故，导致我无法窥得位运算的精妙之处。

位运算有四种基本操作：

此外，在$m$位二进制数中，通常称最低位为第0位，从右到左以此类推，最高位为第$m-1$位。

补码

以C++的 int 和 unsigned int 类型为例。

• 32位无符号整数 unsigned int：原码与补码相同。
• 32位有符号整数 int：最高位为符号位，0表示非负数，1表示负数。最高位为0时，直接可以当作32位二进制数，最高位为1时，需要将其按位取反并+1，得到负数的值。

移位运算

左移

二进制下每一位数字同时向左移动，低位填充0，高位越界后舍弃。
$$1\ll n=2^n,n\ll 1=2n$$

算术右移

二进制补码表示下每一位数字同时向右移动，高位以符号位填充，低位越界后舍弃。
$$n\gg 1=\lfloor \frac{n}{2.0}\rfloor$$

逻辑右移

在二进制补码表示下每一位数字同时向右移动，高位以0填充，地位越界后舍弃。

例题

1、$a^b$

求a的b次方对p取模的值，其中 $1\le a,b,p\le 10^9$。

思路：我们知道，每一个正整数都可以通过二进制数来表示。假设b在二进制表示下有k位，因此，$b$可以表示为：
$$b=c_{k-1}{2}^{k-1}+c_{k-2}{2}^{k-2}+...+c_0{2}^0$$
其中，$c_i(0\le i\le k)$表示$b$在二进制表示下第$i$位的数字，取值为0或者1。

因此：
$$a^b=a^{c_{k-1}\ast 2^{k-1}}\ast a^{c_{k-2}\ast 2^{k-2}}\ast ...\ast a^{c_0\ast 2^0}$$
因为$k=\lceil lo{g}_2(b+1)\rceil$，所以上式乘积项最多为$\lceil lo{g}_2(b+1)\rceil$个。又因为：
$$a^{2^i}=(a^{2^{i-1}}{)}^2$$
所以我们可以通过k次递推求出每个乘积项。

代码：

int power(int a, int b, int p)
{
    int ans = 1 % p;
    for(; b > 0; b >>= 1){
        // 如果该位为1
        if(b & 1){
            ans = (long long)ans * a % p;  // 通过long long强制转换，防止越界
        }
        a = (long long)a * a % p;
    }
    return ans;
}

算法时间复杂度：$O(lo{g}_{2}b)$。

如果$1\le a,b,c\le 10^{18}$，那么 long long 的强制类型转换不再可行，并且不存在一个可供强制类型转换的128位整型类型。因此，我们需要一些特殊的处理办法。

2、64位整数乘法

求a乘b对p取模的值，其中$1\le a,b,c\le 10^{18}$。

思路：采用快速幂的思想。把整数 b 用二进制表示，即$b=c_{k-1}{2}^{k-1}+c_{k-2}{2}^{k-2}+...+c_0{2}^0$，那么$a\ast b=c_{k-1}\ast a\ast 2^{k-1}+c_{k-2}\ast a\ast 2^{k-2}+...+c_0\ast a\ast 2^0$。

因为$a\ast 2^i=(a\ast 2^{i-1})\ast 2$，如果已经求出$a\ast 2^{i-1}$%p，那么预算过程中每一步的结果都不会超过$2\ast 10^{18}$，仍然在64位整数 long long 的表示范围内。所以可以通过 k 次递推求出每个乘积项。

代码：

long long mul_64(long long a, long long b, long long p)
{
    long long ans = 0;
    for(; b > 0; b >>= 1){
        if(b & 1){
            ans = (ans + a) % p;
        }
        a = a * 2 % p;
    }
}

时间复杂度：$O(lo{g}_{2}b)$