CINTA第二次选做作业2

第三章第九题

题目： 定义 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 都可以。
算法： 首先由公式可得：${\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]}^n={\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]}^{n-1}\ast \left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]=\left[\begin{array}{cc}F(n)& F(n-1)\\ F(n-1)& F(n-2)\end{array}\right]\ast \left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]$
$=\left[\begin{array}{cc}F(n)+F(n-1)& F(n)\\ F(n-1)+F(n-2)& F(n-1)\end{array}\right]$
由此可得，矩阵每进行一次乘法，结果的第一列的各行数据就是原来两列的两行数据之和，第二列就是原来的第一列，因此仅需进行加法和赋值即可。
答案： C++代码如下

int F1(int n)
{
	int f1 = 1, f2 = 1, f3 = 1, f4 = 0;
	n -= 1;//矩阵的n次方会进行n-1次相乘
	while (n--)
	{
		int tempf1 = f1, tempf3 = f3;//第一列
		f1 = f1 + f2;
		f3 = f3 + f4;
		f2 = tempf1;
		f4 = tempf3;
	}
	return f2;
}

第四章第十四题

题目： 给定任意正整数$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$即题目得证。