第三章 循环结构练习

---

❥(^_-)第一题 $a\wedge b$ 【acwing 89】

求 a 的 b 次方对 p 取模的值。

输入格式： 三个整数 $a,b,p$ ,在同一行用空格隔开。
输出格式：输出一个整数，表示$a\stackrel{ˆ}{}b\; mod\; p$的值。

数据范围：$0\le a,b\le 10^9$ ; $1\le p\le 10^9$

输入样例1：
$3\; 2\; 7$
输出样例1：
$2$

输入样例2：
$5201314\; 0\; 1$
输出样例2：
$0$

输入样例3：
$131452020\; 5201314\; 1314$
输出样例3：
$486$

---

思路: 本题需要使用快速幂的思想。下面简述一下快速幂是如何进行计算的，假设我们现在需要求解$2^{31}$的结果。计算过程如下所示：
$$\begin{array}{rl}{2}^{31}& =1\times 2^{31}写成k\times a^b的形式\\ & =1\times (2\times 2^{30})因为上一步的b为奇数，所以把a^b拆分成a\times a^{b-1}的形式\\ & =2\times 2^{30}将上一步拆分出来的a合并到k中，即k=k\times a\\ & =2\times (2^2{)}^{15}因为上一步的b为偶数,所以此处变形为(a^2{)}^{\frac{b}{2}}的形式\\ & =2\times 4^{15}合并上一步的操作，即a=a^2,b=\frac{b}{2}\\ & =2\times (4\times 4^{14})上一步的b是奇数，所以拆分一个a出来，指数减1\\ & =8\times 4^{14}\\ & =8\times (4^2{)}^7因为上一步指数部分b是偶数，所以底数翻倍，指数除以2\\ & =8\times 16^7\\ & =8\times (16\times 16^6)因为上一步指数部分b是奇数，所以拆分一个a出来，指数减1\\ & =128\times 16^6\\ & =128\times (16^2{)}^3因为上一步指数部分b是偶数，所以底数翻倍，指数除以2\\ & =128\times 256^3\\ & =128\times (256\times 256^2)因为上一步指数部分b是奇数，所拆分一个a出来，指数减1\\ & =32768\times 256^2\\ & =32768\times (256^2{)}^1因为上一步指数b是偶数，所以底数翻倍，指数除以2\\ & =32768\times 65536^1\\ & =32768\times 65536\times 65536^0因为上一步指数部分为奇数，所以拆分一个a出来，指数减1\\ & =2147483648\times 65536^0指数部分为0，结束操作。\\ & =2147483648\times 1\\ & =2147483648\end{array}$$
从上面过程中可以看出如果我们正常的对$2^{31}$进行计算的话循环一共需要31次，而上述的快速幂算法只需要循环7次（除去上述啰嗦的一些算式）即可以解决。可以看出快速幂算法是比普通的循环要快。那么下面给出快速幂算法最原始的代码：

long long fast_pow(long long a, long long b)
{
	long long ans  =  1;
	while(b != 0) // 判断指数部分是否为0
	{
		if(b % 2 == 0)  // 若指数是偶数，底数翻倍，指数除以2
		{
			a = a * a;
			b /= 2;
		}
		else  // 指数部分为奇数的情况，分离一个a出去，指数部分减1
		{
			ans = ans * a;
			b = b - 1;
		}
	}
	return ans;
}

对于上述判断 $b\%2!=0$ 以及 $b/=2$的时候，我们可以使用位运算来操作。$b\%2$的操作和$b\&1$的操作是等效的，而$b/=2$的操作和$b>>1$的操作也是等效的。另外当$b$为奇数的时候，$b-1$其实等价于$b>>1$，向右位移一位去掉二进制中最后一位的1。位运算可以加快运算速度，所以上述代码也可以改成（知道过程就可以自己写出来，自己能够写出来就行）：

long long fast_pow(long long a, long long b)
{
	long long ans  =  1;
	while(b) // 判断指数部分是否为0
	{
		if(b & 1)  // 指数部分是否为奇数
		{
			
			ans *= a;
			b >>= 1;
		}
		else 
		{
			a *= a;
			b >>= 1;
		}
	}
	return ans;
}

这里面发现可以合并的代码现在合并在一起如下：

long long fast_pow(long long a, long long b)
{
	long long ans  =  1;
	while(b)
	{
		if(b & 1) 	ans *= a;
		else 	a *= a;
		b >>= 1;
	}
	return ans;
}

在这里指数奇数的情况可以和偶数的情况同时计算，代码如下（最好记住，记不住的话尽量把过程记住）：

long long fast_pow(long long a, long long b)
{
	long long ans  =  1;
	while(b)
	{
		if(b & 1) 	ans *= a; // 指数是否为奇数，奇数执行，执行完之后立即执行偶数的情况
		a *= a; 
		b >>= 1;
	}
	return ans;
}

当然这里也有一个for循环的版本：

long long fast_pow(long long a, long long b)
{
	long long ans  =  1;
	for(;b;b>>=1)
	{
		if(b & 1) 	ans *= a;
		a *= a;
	}
	return ans;
}

对于这道题来说还有一个问题，就是求出题的答案。假设给定两个数，一个为$a$，另一个数为$b$，给定$a$和$b$的表达式如下：
$$\begin{array}{rl}& a=mp+i\\ & b=np+j\end{array}$$
式子中的$m,n,i,j$均为整数，且满足$0\le i,j\le p$，求$ab\%p$的结果。现在进行一个推导：
$$\begin{array}{rl}ab& =(mp+i)(np+j)\\ & =[mnp+mj+ni]p+ij\end{array}$$
对于$ab=(mp+i)(np+j)$来说，肯定能够被$p$整除。而对于$ij$来说，若$0\le ij<p$，则余数为$ij$，若$ij\ge p$，则余数为$ij\%p$。综上所述，可以得出如下结论：
$$\begin{array}{r}ab\%p=[(a\%p)(b\%p)]\%p\end{array}$$
接下来难度提高一点，求解$a^b\%p$。根据上述的条件，我们可以得到一下结论：
$$\begin{array}{rl}{a}^b\%p& =\begin{array}{c}(\underbrace{a\times a\times \cdots \times a})\\ b\end{array}\%p\\ & =\begin{array}{c}[\underbrace{(a\%p)(a\%p)\cdots (a\%p)}]\\ b\end{array}\%p\end{array}$$
若多了个系数$k$，则：
$$\begin{array}{r}(k\times a^b)=\begin{array}{c}[(k\%p)\underbrace{(a\%p)(a\%p)\cdots (a\%p)}]\\ b\end{array}\%p\end{array}$$
根据这个式子可以推导出快速幂模板如下：
$算法模板一$：

long long fast_pow(long long a, long long b, long long p)
{
	long long ans  =  1;
	while(b)
	{
		if(b & 1) 	ans = ans * a % p;  // 系数k%p
		a = a * a % p;  // b个a%p
		b >>= 1;
	}
	return ans % p;  // 当b直接等于0的时候循环无法解决，所以此处直接对ans取模
}

$算法模板二$：

long long fast_pow(long long a, long long b, long long p)
{
	long long ans  =  1;
	for(;b;b>>=1)
	{
		if(b & 1) ans = ans * a % p;
		a = a * a % p;
	}
	return ans % p;
}

上面的两个 模板的区别仅在于是用$for$循环还是是用$while$循环。这两个模板需要熟练掌握其中一个模板，以此来应对快速幂算法的题。

---

本题题解代码如下：

#include <stdio.h>

int main()
{
	long long a, b, p, ans = 1;
	scanf("%lld%lld%lld", &a, &b, &p);
	for(; b; b>>=1)
	{
		if(b&1) ans = ans * a % p;
		a = a * a % p;
	}
	ans %= p;  // 特殊处理当b直接输入0的情况，其他情况此处ans还是ans。
	printf("%lld", ans);
	return 0;
}

---

❥(^_-)第二题 $64$位整型乘法 【acwing 90】

求 a 乘 b 对 p 取模的值。

输入格式： 第一行输入整数$a$，第二行输入整数$b$,第三行输入整数$c$。
输出格式：输出一个整数，表示$ab\; mod\; p$的值。

数据范围：$0\le a,b,p\le 10^{18}$

输入样例1：
$$\begin{gathered} 3 \\ 2 \\ 7 \end{gathered}$$
输出样例1：
$6$

输入样例2：
$$\begin{gathered} 5201314 \\ 520 \\ 1314 \end{gathered}$$
输出样例2：
$868$

输入样例3：
$$\begin{gathered} 131452013145201314 \\ 520131452013145201 \\ 770880131452013920 \end{gathered}$$
输出样例3：
$101736777287696674$

---

思路:
首先这道题我们可以看成是$a$个$b$或者是$b$个$a$相乘再求余数，则可以得出如下表达式：
$$\begin{array}{rl}ab\%p& =\begin{array}{c}\underbrace{(a\times a\times \cdots \times a)}\\ b\end{array}\%p\\ & =\begin{array}{c}[\underbrace{(a\%p)(a\%p)\cdots (a\%p)}]\\ b\end{array}\%p\end{array}$$

从上述推导过程中可知只需要使用循环运算$b$次即可。可以得出如下代码：

long long mul(long long a, long long b, long long p)
{
	long long ans = 1;
	for(; b; b -= 1)
		ans = ans * a  % p;
	return ans;
}

但是在前面曾经提到过一秒可以运算的次数，而此处的$b$的范围是$0\le b\le 10^{18}$，所以运算$b$次一定会TLE（Time Limit Exceeded，时间超时）。所以换一个思路，根据上一个题推导出来的式子：
$$\begin{array}{r}ab\%p=[(a\%p)(b\%p)]\%p\end{array}$$
此时可以得出如下计算的代码：

long long mul(long long a, long long b, long long p)
{
	a %= p;
	b %= p;
	return (a * b) % p;
}

但是很遗憾，这个式子也不能帮我们解决这个问题。因为$0\le a,b,p\le 10^{18}$，所以此处的$(a\%p)(b\%p)$可能会超数据类型范围。此时我们需要考虑一个新的算法。接下来我们来进行推导。假设$b$的二进制有$k$位且满足如下式子：
$$\begin{array}{r}b=c_{k-1}\times 2^{k-1}+c_{k-2}\times 2^{k-2}+\cdots +c_1\times 2^1+c_0\times 2^0\end{array}$$
则
$$\begin{array}{r}ab=c_{k-1}\times a\times 2^{k-1}+c_{k-2}\times a\times 2^{k-2}+\cdots +c_1\times a\times 2^1+c_0\times a\times 2^0\end{array}$$
又
$$\begin{array}{r}a\times 2^i=(a\times 2^{i-1})\times 2\end{array}$$
所以可以得出如下结论:
$$\begin{array}{r}(a\times 2^i)\%p=\{[(a\times 2^{i-1})\%p]\times 2\}\%p\end{array}$$
上述结论看不懂可参考$ab\%p=[(a\%p)(b\%p)]\%p$，此处默认$b\%p=2$。故可以得出来以下代码模板：

long long mul(long long a, long long b, long long p)
{
	long long ans = 0;
	for(; b; b >>= 1)
	{
		if(b&1) ans = (ans + a) % p;  // b二进制最后以为是否为1，保持ans答案在数据范围中。
		a = a * 2 % p;	// 执行结论操作
	}
	return ans;
}

---

代码：

#include <stdio.h>

int main()
{
	long long a, b, p, ans = 0;
	scanf("%lld%lld%lld", &a, &b, &p);
	for(; b; b >>= 1)
	{
		if(b&1) ans = (ans + a) % p;
		a = a * 2 % p;
	}
	printf("%lld", ans);
	return 0;
}

---

❥(^_-)第三题 斐波那契数列

$Fibonacci$数列的递推公式为：$F_n=F_{n-1}+F_{n-2}$，其$F_1=F_2=1$。当$n$比较大，$F_n$也非常大，现在我们想知道，$F_n$除以$5201314$的余数是多少。

输入格式：输入包含一个整数$n$。
输出格式：输出一行，包含一个整数，表示$F_n$除以$5201314的余数$。
数据范围：$1\le n\le 5201314$。

输入样例1：
$520$
输出样例1：
$97115$

输入样例2：
$1314$
输出样例2：
$1136042$

输入样例3：
$5201314$
输出样例3：
$1906391$

---

思路： 本题的思路其实和前面的两道题大同小异，由于斐波那契数列后面的数越来越大，所以可以对每一项都进行取余操作即可。

---

代码：

#include <stdio.h>

int main()
{
	int a = 1, b = 1, n, ans = 1;
	scanf("%d", &n);
	for(int i = 3; i <= n; i ++)
	{
		ans = (a + b) % 5201314;
		a = b;
		b = ans;
	}
	printf("%d", ans);
	return 0;
}