第二章课后题

CINTA作业二

---

7. 手动计算以下模 $m$ 下 $a$ 的乘法逆元。
(a) $m=11,a=5;$
$$\left[\begin{array}{ccc}1& 0& 11\\ 0& 1& 5\end{array}\right]$$

解：

$$\left[\begin{array}{ccc}0& 1& 5\\ 1& -2& 1\end{array}\right]\text{\hspace{1em}(1)}$$


$$\{\begin{array}{rlr}r& =& 1\\ s& =& -2\\ b& =& 1\end{array}\text{\hspace{1em}(2)}$$

乘法逆元为：-2

(b) $m=121,a=13;$
$$\left[\begin{array}{ccc}1& 0& 121\\ 0& 1& 13\end{array}\right]$$
解：
$$\left[\begin{array}{ccc}0& 1& 13\\ 1& -9& 4\end{array}\right]\text{\hspace{1em}(1)}$$


$$\left[\begin{array}{ccc}1& -9& 4\\ -3& 28& 1\end{array}\right]\text{\hspace{1em}(2)}$$


$$\{\begin{array}{rlr}r& =& -3\\ s& =& 28\\ b& =& 1\end{array}\text{\hspace{1em}(3)}$$

乘法逆元为：28

© $m=1021,a=131;$
$$\left[\begin{array}{ccc}1& 0& 1021\\ 0& 1& 131\end{array}\right]$$
解：
$$\left[\begin{array}{ccc}0& 1& 131\\ 1& -7& 104\end{array}\right]\text{\hspace{1em}(1)}$$


$$\left[\begin{array}{ccc}1& -7& 104\\ -1& 8& 27\end{array}\right]\text{\hspace{1em}(2)}$$


$$\left[\begin{array}{ccc}-1& 8& 27\\ 4& -31& 23\end{array}\right]\text{\hspace{1em}(3)}$$


$$\left[\begin{array}{ccc}4& -31& 23\\ -5& 39& 4\end{array}\right]\text{\hspace{1em}(4)}$$


$$\left[\begin{array}{ccc}-5& 39& 4\\ 29& -226& 3\end{array}\right]\text{\hspace{1em}(5)}$$


$$\left[\begin{array}{ccc}29& -226& 3\\ -34& 265& 1\end{array}\right]\text{\hspace{1em}(6)}$$


$$\begin{array}{c}\{\begin{array}{l}r=-34\\ s=265\\ b=1\end{array}\end{array}$$

乘法逆元为：265

8. 编写C语言程序完成模指数运算，即给定整数 $x,y$ 和 $m$ 为输入，计算并计算返回值 $x^y\mathrm{mod}m$。

#include<iostream>
using namespace std;
long long int powANDmod(int x,int y,long long m){
	if(y<0) y=~y+1;
	long long res=1;
	while(y!=0){
		if(y&1) res=(res*x)%m;
		x=(x*x);
		res%=m;			//自己想的：猜测此层面只需 res%=m;
//		x=(x*x)%m;		//bintou老师的解法
		y=y>>1;
	}
	return res;
}

int main(){
    int a,b;
    long long c;
    cout<<"Please input numbers (x,y,m): "
    cin>>a>>b>>c;
	long long int t=powANDmod(a,b,c);
	cout<<t<<endl;
	return 0;
}

9. 定义 Fibonacci 数列如下：$F(0)=0,F(1)=1$，且对于 $n⩽2,F(n)=F(n-1)+F(n-2)$。所以，该数列是：$0,1,1,2,3,5,8,13,21,..$。如何能快速地求出 $F(n)$ 呢？很幸运，我们有以下等式：

$$\begin{array}{c}{\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]\end{array}$$

虽然，看上去该算法需要一次矩阵的指数运算，但是借助快速指数运算的方法，这里可以产生一个快速求解 $F(n)$ 的算法。请给出算法，并编程实现，C 语言或者 Python 都可以。

#include<iostream>
using namespace std;
#define M 2

long long Fibo(int y){
	while(y<1){
		cerr<<"Error, this number is lower than 1."<<endl
		<<"Please input n again: ";
		cin>>y;
	}
	long long x[M][M]={{1,1},{1,0}},oth[M][M]={{1,1},{1,0}},temp[M][M];
	y--;				//这里多了一次偏移
	while(y!=0){
		if(y&1){
			//res*=x;	实际上没必要的矩阵乘法：因为有限个值（4个）可以直接算
			for(int i=0;i<M;i++){
				for(int j =0;j<M;j++){
					temp[i][j]=0;
					for(int k=0;k<M;k++)
						temp[i][j]+=oth[i][k]*x[k][j]; 
				}	
			}
			
			//进行oth值的更新
			for(int i=0;i<M;i++){
				for(int j=0;j<M;j++){
					oth[i][j]=temp[i][j];
				}
			}
		}
		
		//x*=x;
		for(int i=0;i<M;i++){
			for(int j=0;j<M;j++){
				temp[i][j]=0;
				for(int k=0;k<M;k++){
					temp[i][j]+=x[i][k]*x[k][j];
				}
			}
		}
		//进行x值的更新
		for(int i=0;i<M;i++){
			for(int j=0;j<M;j++){
				x[i][j]=temp[i][j];
			}
		}
		y=y>>1;
	}
	return oth[1][0];
}

int main(){
	for(int i=1;i<50;i++) cout<<Fibo(i)<<endl;
	return 0;
}

10. 给定互素的正整数 $c$ 和 $m$，请证明在 $mod$ $m$ 的意义上存在唯一确定的整数值 $c^{-1}$，它使得 $c{c}^{-1}\equiv 1$

存在性证明：
对于互素的两数：$gcd(c,m)=1$
那么我们由扩展欧几里得算法必定可以得到：$cr+ms=1$
而这代表着 $xr+ms=1$ 必定成立（无论 $c$ 与 $c^{-1}$ 是否相等）

设 $1\le c^{-1}<m$ 证明其唯一性：
假设存在另一数 $c_1$
$c{c}^{-1}\equiv c{c}_1\equiv 1(\mathrm{mod}m)$
那么：$m|(c^{-1}-c_1)$
所以在规定区间，仅有一个值

编程题：编写一个 Python 程序计算乘法逆元，即输入互素的正整数 $c$ 和 $m$，返回 $c$ ，使得 $c{c}^{-1}\equiv (\mathrm{mod}m)$。要求:只返回为正整数的 $c^{-1}$。

# 乘法逆元算法，要求带有c,m
def egcd(c, m):
    r0, s0 = 1, 0
    r1, s1 = 0, 1
    while m % c != 0:
        q = m // c
        c, m = m % c, c
        r0, r1 = r1, r0 - q * r1
        s0, s1 = s1, s0 - q * s1
    # return r1, s1, c
    if s1 <= 0:
        return s1 + m
    return s1


c = eval(input("c: "))
m = eval(input("m: "))

# li = list(egcd(c, m))

# print("The result:\n\t\t r={}, s={}, b={}".format(li[-3], li[-2], li[-1]))
print("The result is: ", egcd(c, m))