CINTA第二次选做作业

第三章第九题

**题目：**定义 Fibonacci 数列如下：$F(0)=0，F(1)=1$，且对于 $n\ge 2，F(n)=F(n-1)+F(n-2)$。所以，该数列是：0, 1, 1, 2, 3, 5, 8, 13, 21, . . . 。如何能快速地求出 $F(n)$呢？
$${\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]}^n=\left[\begin{array}{cc}F(n+1)& F(n)\\ F(n)& F(n-1)\end{array}\right]$$
虽然，看上去该算法需要一次矩阵的指数运算，但是借助快速指数运算的方法，这 里可以产生一个快速求解$F(n)$ 的算法。请给出算法，并编程实现，C 语言或者 Python 都可以。
**答案：**C++代码如下
先声明一个二阶矩阵的结构体，并进行二阶矩阵乘法的运算符重载

//二阶矩阵结构体
struct Matrix
{
	long long Data[2][2];
};
//二阶矩阵乘法运算符重载
Matrix operator*(Matrix A, Matrix B)
{
	Matrix C;
	C.Data[0][0] = A.Data[0][0] * B.Data[0][0] + A.Data[0][1] * B.Data[1][0];
	C.Data[1][0] = A.Data[1][0] * B.Data[0][0] + A.Data[1][1] * B.Data[1][0];
	C.Data[0][1] = A.Data[0][0] * B.Data[0][1] + A.Data[0][1] * B.Data[1][1];
	C.Data[1][1] = A.Data[1][0] * B.Data[0][1] + A.Data[1][1] * B.Data[1][1];
	return C;
}

矩阵快速指数运算的递归版本：

Matrix F(int n)//递归
{
	Matrix A;
	A.Data[0][0] = 1;
	A.Data[0][1] = 1;
	A.Data[1][0] = 1;
	A.Data[1][1] = 0;
	if (n == 1)
		return A;
	Matrix temp = F(n / 2);
	if (n % 2 == 0)
		return temp * temp;
	else if (n % 2 == 1)
		return A * temp * temp;
}

矩阵快速运算的迭代版本

Matrix F2(int n)//迭代
{
	Matrix A;
	A.Data[0][0] = 1;
	A.Data[0][1] = 1;
	A.Data[1][0] = 1;
	A.Data[1][1] = 0;
	Matrix res;//置为单位矩阵
	res.Data[0][0] = 1;
	res.Data[0][1] = 0;
	res.Data[1][0] = 0;
	res.Data[1][1] = 1;
	while (n != 0)
	{
		if (n & 1)
			res = res * A;
		n >>= 1;
		A = A * A;
	}
	return res;
}

第四章第十四题

**题目：**给定任意正整数$n，n$ 的所有因子分别记为 $d_0,d_1,\cdots ,d_r，$其中包括 $1$和$n$。记函数：
$$F(n)=\varphi (d_0)+\varphi (d_1)+\cdot \cdot \cdot +\varphi (d_r)$$
现在需要证明 $F(n)=n$。为完成这个任务，请依次完成以下小任务：
$(a)$. 证明，对任意素数 $p，F(p)=p$。
$(b)$. 证明，对任意素数$p，F(p^2)=p^2$。
$(c)$. 证明，对任意素数 $p，F(p^k)=p^k$，$k$ 是任意正整数。
$(d)$. 证明，对任意素数$p和q$，$F(pq)=pq$。
$(e)$. 证明函数 F 是积性函数，即对任意的正整数 $m和n$，如果 $gcd(m,n)=1$，则$F(mn)=F(m)F(n)$。[提示：既然 m 与 n 互素，则它们之间就不存在大于一 的公因子，那么它们的因子相乘之后作用于欧拉函数会满足什么特性？]
$(f)$. 最后，完成证明$F(n)=n$的任务。[提示：任意整数 n 都是素数的乘积，即 $n=p_0^{k_0}{p}_1^{k_1}\cdot \cdot \cdot p_t^{k_t}$ 。]
解：
$(a).$对于任意素数$p$，有两个因子$1和p$，所以有$F(p)=\varphi (1)+\varphi (p)$，根据欧拉$Phi$函数性质有$\varphi (p)=p-1$，所以$F(p)=\varphi (1)+\varphi (p)=1+p-1=p$，题目得证。
$(b).$对于任意素数$p$，$p^2$有三个因子$1,p,p^2$
$\therefore F(p^2)=\varphi (1)+\varphi (p)+\varphi (p^2)$
根据欧拉$Phi$函数性质有$\varphi (p^k)=p^k-p^{k-1}$
$\therefore F(p^2)=\varphi (1)+\varphi (p)+\varphi (p^2)=1+p-1+p^2-p=p^2$即题目得证。
$(c).$对于任意素数$p,p^k$有因子$1,p,p^2\cdots p^k$，
$\therefore F(p^k)=\varphi (1)+\varphi (p)+\cdots +\varphi (p^k)$
$=1+p-1+p^2-p+\cdots +p^{k-1}-P^{k-2}+p^k-p^k-1$
$=P^k$，$\therefore$ 题目得证。
$(d).$对于任意素数$p和q$，$pq$有因子$1,p,q,pq$
$\therefore F(pq)=\varphi (1)+\varphi (p)+\varphi (q)+\varphi (pq)$
$=1+p-1+q-1+(p-1)(q-1)$
$=pq\therefore$题目得证
$(e).$对任意的正整数 $m和n$，如果 $gcd(m,n)=1$，即$m与n$互素，则$m$的因子与$n$的因子也都是互素的，设$m$的因子为$d_{m1}，d_{m2}\cdots d_{mx}，$包括$1和m$，共有$x$个，$n$的因子设为$d_{n1}，d_{n2}\cdots d_{ny}$，包括$1和n$，共有$y$个，
则$mn$的因子有$d_{m1}{d}_{n1}，d_{m1}{d}_{n2}\cdots d_{m1}{d}_{ny}$，$，d_{m2}{d}_{n1}，d_{m2}{d}_{n2}\cdots d_{m2}{d}_{ny}\cdots$ $\cdots ，d_{mx}{d}_{n1}，d_{mx}{d}_{n2}\cdots d_{mx}{d}_{ny}$共$xy$个因子。
$\therefore F(mn)=\varphi (d_{m1}{d}_{n1})+\varphi (d_{m1}{d}_{n2})+\cdots \varphi (d_{m1}{d}_{ny})$$+\varphi (d_{m2}{d}_{n1})+\varphi (d_{m2}{d}_{n2})\cdots +\varphi (d_{m2}{d}_{ny})+\cdots$$\cdots \varphi (d_{mx}{d}_{n1})+\varphi (d_{mx}{d}_{n2})\cdots \varphi (d_{mx}{d}_{ny})$
$=\varphi (d_{m1})\varphi (d_{n1})+\varphi (d_{m1})\varphi (d_{n2})+\cdots \varphi (d_{m1})\varphi (d_{ny})$$+\varphi (d_{m2})\varphi (d_{n1})+\varphi (d_{m2})\varphi (d_{n2})\cdots +\varphi (d_{m2})\varphi (d_{ny})+\cdots$$\cdots \varphi (d_{mx})\varphi (d_{n1})+\varphi (d_{mx})\varphi (d_{n2})\cdots \varphi (d_{mx})\varphi (d_{ny})$
$=\varphi (d_{m1})[\varphi (d_{n0})+\varphi (d_{n1})+\varphi (d_{n2})+\cdots +\varphi (d_{ny}]$$+\varphi (d_{m2})[\varphi (d_{n0})+\varphi (d_{n1})+\varphi (d_{n2})+\cdots +\varphi (d_{ny}]+\cdots \cdots$
$\varphi (d_{mx})[\varphi (d_{n0})+\varphi (d_{n1})+\varphi (d_{n2})+\cdots +\varphi (d_{ny}]$
$=[\varphi (d_{m0})+\varphi (d_{m1})+\varphi (d_{m2})+\cdots +\varphi (d_{mx}]\ast [\varphi (d_{n0})+\varphi (d_{n1})+\varphi (d_{n2})+\cdots +\varphi (d_{ny}]$
=$F(m)F(n)$，即题目得证。
$(f).$任意整数 n 都是素数的乘积，即 $n=p_0^{k_0}{p}_1^{k_1}\cdot \cdot \cdot p_t^{k_t}$ 。可得$F(n)=F(p_0^{k_0}{p}_1^{k_1}\cdot \cdot \cdot p_t^{k_t})$，又由$(e)$中的结论可得
$F(n)=F(p_0^{k_0}{p}_1^{k_1}\cdots p_t^{k_t})$
$=F(p_0^{k_0})F(p_1^{k_1})\cdots F(p_t^{k_t})$
又由$(c)$的结论$F(p^k)=p^k$可得
$F(n)=F(p_0^{k_0})F(p_1^{k_1})\cdots F(p_t^{k_t})=p_0^{k_0}{p}_1^{k_1}\cdots p_t^{k_t}=n$即题目得证。