CINTA——2选做

CINTA——2选做

第三章第九题

计算斐波那契数列$F(n)$，就是解题中已给出的矩阵的n次幂，那么怎样才可以快速的解出来呢？

当然，最简单的算法就是对矩阵进行循环乘法了，时间复杂度为$O(n)$，但是，如果我们采用教材中3.3节所给出的模指数运算的算法包含的中的快速幂算法，就可以将时间复杂度降为$O(logn)$

那么对于模指数运算的算法$x^y(\mathrm{mod}n)$来说，这里的$x$就相当于矩阵
$$\begin{array}{c}\left[\begin{array}{cc}1& 1\\ 1& 0\end{array}\right]\end{array}$$
这里的$y$就相当于n,然后我们只需要将书中算法的整数$x$，换为矩阵就好了，而且矩阵已知，我们的参数列表只需要给出即可，不过矩阵乘法需要自己写代码来实现，下面给出递归版本的代码↓

//递归快速幂来求解斐波那契数列
int qpow(inta[][2],int n){
    //a是二维数组代表题中的矩阵
    //注意，这里的*要重载为满足矩阵乘法运算规则
    if(n==0)return 0;
    int temp=qpow(a,n/2);
    else{
        if(n%2==0){
            return temp*temp;//若为避免使用temp*temp，可以直接使用return(a*a,n2)
            }//阶为偶数
        else{
            if(n%2==1){//阶为奇数
                return temp*temp*a;//如果是模指数算法，则在上面两个如return后面加上%p即可
            }
        }
    }
}

注意，教材中的是模指数算法，但是求解本题需要的是快速幂算法，区别就在于模指数算法多了模运算。

第四章第十四题

$(1)$因为p是素数，所以p的因子只有1和p，根据题中函数的定义，则$F(p)=\varphi (1)+\varphi (p)=1+p-1=$得证。

$(2)$同上面小题的思路，$p^2$只有$1,p,p^2$因子，所以$F(p^2)=\varphi (1)=\varphi (p)+\varphi (p^2)=1+p-1+p^2-p=p^2$得证。

$(3)$思路同第二小问，$p^k$的因子有$1,p^k,p,p^{k-1},p^{k-2},p^2\cdots$,则$F(p^k)=\varphi (1)+\varphi (p)+\varphi (p^2)+\cdots +\varphi (p^k)=1+p-1+p^2-p+\cdots +p^k-p^{k-1}=p^k$得证。

$(4)$对任意素数$p,q$，$pq$的因子有$1,p,q,pq$。$F(pq)=\varphi (1)+\varphi (p)+\varphi (q)+\varphi (pq)=1+p-1+q-1+(p-1)(q-1)=p+q-1+pq-p-q+1=pq$得证。

$(5)$如果$m,n$互为素数，那么它们的因子也要互为素数才可以。故假设m的因子是$a_1,a_2,\cdots a_x$，n的因子为$b_1,b_2,\cdots ,b_y$。

那么$m\cdot n$相乘可表示为因子相乘，即$a_1{b}_1,a_1{b}_2\cdots a_1{b}_y,a_2{b}_1,a_2{b}_2\cdots a_2{b}_y,\cdots a_x{b}_1,a_x{b}_2\cdots a_x{b}_y$

那么函数$F(mn)=\varphi (a_1{b}_1)+\varphi (a_1{b}_2)+\cdots +\varphi (a_1{b}_1)+\varphi (a_2{b}_2)+\cdots +\varphi (a_x{b}_1)+\varphi (a_x{b}_2)+\cdots +\varphi (a_x{b}_y)$根据欧拉函数公式，可继续有
$$\varphi (a_1)(\varphi (b_1)\cdots \varphi (b_y))+\varphi (a_2)(\varphi (b_1)\varphi (b_1)\cdots \varphi (b_y))+\cdots +\varphi (a_x)(\varphi (b_1)\cdots \varphi (b_y))$$
$$=(\varphi (a_1)\cdots \varphi (a_x))(\varphi (b_1)\cdots \varphi (b_y))=F(m)F(n)$$
得证。

$(6)$根据提示：任意整数n都是素数的乘积，即$n=p_0^{k_0}{p}_1^{k_1}\cdots p_t^{k_t}$，则$F(n)=F(p_0^{k_0}{p}_1^{k_1}\cdots p_t^{k_t})$，根据$(5)$的结论可知上式等于$F(p_0^{k_0})F(p_1^{k_1})\cdots F(p_t^{k_t})$，再由$(3)$中的结论可知上式可继续等于$p_0^{k_0}\cdots p_t^{k_t}=n$得证。