快速幂、取整取模

快速幂

例子：求 $a^{b}modm$的值。

如$a^7$ 7 二进制为（111）$a^7$ = $a^4$ * $a^2$ * $a^1$
$a^{18}$ 18 二进制为（10010）$a^{18}$ = $a^{16}$ * $a^2$
a的任意次幂都可以由多个不同的a的$2^k$次幂相乘得到，可以通过对a不断平方并把结果赋值给a，依次得到 $a^{2^0}$ （a）， $a^{2^1}$（$a^2$ = a * a）， $a^{2^2}$（$a^4$ = $a^2$ * $a^2$）， …… $a^{2^k}$，即通过不断地进行 a = a * a；

而$a^b$ 𝑚𝑜𝑑 𝑚 = ($a^k$ 𝑚𝑜𝑑 𝑚 )* ($a^t$ 𝑚𝑜𝑑 𝑚 )（k + t = b），因此求 $a^b$ 𝑚𝑜𝑑 𝑚 的值可由以下代码实现

int quick_pow (int a, int b)
{
	int res = 1;
    while (b) 
    {
        if (b&1) 
            res = res * a % mod;
        a = a*a % mod;
        b >>= 1;
    }
    return res;
}

以b = 5为例来看上面代码，5二进制为（101），“&”为按位与运算
b = 101，b&1 = 1； res = a; a = a*a;
b = 10，b&1 = 0; （res = res * a % mod;不执行），res = a；a = $a^2$ * $a^2$
b = 1，b&1 = 1；res = a * $a^4$ ；a = $a^4$ * $a^4$
b = 0，退出循环
res的值即为$a^5$ 。

应用

快速幂可以用来计算斐波那契数列。
已知斐波那契数列 $F_{}n$ = $F_{}n-1$ + $F_{n-2}$ ​​(n>=3),$F_1$ = 1,$F_{}2$ ​​= 1
由斐波那契数列定义可得

只需利用快速幂计算$${\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]}^{n-1}$$再与$$\left[\begin{array}{c}f(2)\\ f(1)\end{array}\right]$$ 相乘即可得出f(n)。

关于取整的经典题目

${\sum }_{k=1}^n\lfloor \frac{n}{k}\rfloor$，n <= $10^9$

容易发现使得$\lfloor \frac{n}{k}\rfloor$每种取值的k一定是连续的一段区间。

例如，使10/k值为1 的 k 在连续区间[6，10]。
因此，只要知道每种不同取值k的区间长度，把不同n/k的值乘以区间长度，再相加，即可得到结果。
左端点 l 从1开始，易知，每次新区间的左端点等于上一个区间的右端点加1
x = n / l，右端点 r 是满足n / r = x的最大正整数， n 除以r 向下取整为x，所以n >= x * r，只有当 r = n / x 向下取整时，r才为最大。
所以r = n / x = n / (n / l)；

int count (int n)
{
	int ans = 0;
	for (int l = 1, r = 1; l <= n; l = r + 1)
	{
		r = n / (n / l);
		ans += n / l * (r - l + 1);
	}
	return ans;
}

取整与取模的转换

𝑛 𝑚𝑜𝑑 𝑚 = 𝑛 −𝑚 $\lfloor \frac{n}{m}\rfloor$

求${\sum }_{k=1}^{n}n\mathrm{m}\mathrm{o}\mathrm{d}k$，n <= $10^9$

${\sum }_{k=1}^{n}n\mathrm{m}\mathrm{o}\mathrm{d}k$ = ${\sum }_{k=1}^{n}n-k\ast \lfloor \frac{n}{k}\rfloor$ =n * n - ${\sum }_{k=1}^{n}k\ast \lfloor \frac{n}{k}\rfloor$

# include <iostream>
using namespace  std;

int count (int n)
{
	int ans = 0;
	for (int l = 1, r = 1; l <= n; l = r + 1)
	{
		r = n / (n / l);
		ans += (l+r)*(r-l+1)/2*(n/l);
	}
	return ans;
}

int main ()
{
    int n;
    cin >> n;
    cout <<n*n - count(n);
    return 0;
}

对每个区间[l，r]，$\lfloor \frac{n}{k}\rfloor$ = n / l，是个固定值

${\sum }_{k=l}^{r}k\ast \lfloor \frac{n}{k}\rfloor$ = $\lfloor \frac{n}{k}\rfloor$ * ${\sum }_{k=l}^{r}k$ = $\lfloor \frac{n}{k}\rfloor$ * (l+r) * (r-l+1)/2 = (n/l) * (l+r)*(r-l+1)/2

由此，可以求出每个区间${\sum }_{k=l}^{r}k\ast \lfloor \frac{n}{k}\rfloor$的值，从而得到${\sum }_{k=1}^{n}k\ast \lfloor \frac{n}{k}\rfloor$的值。

由 ${\sum }_{k=1}^{n}n\mathrm{m}\mathrm{o}\mathrm{d}k$ = n * n - ${\sum }_{k=1}^{n}k\ast \lfloor \frac{n}{k}\rfloor$可以求出最终结果