CINTA第一次作业

第一章

迭代版本的乘法

int mul(int a, int b)
{
    int res = 0;
    while(b)
    {
        if(b & 1) res += a;
        a <<= 1;
        b >>= 1;
    }
    return res;
}

除法算法的证明

除法算法：
对任意给定的整数$a$和$b$，其中$b>0$，存在唯一的整数对$q$（商）和 $r$（余数）使得，
$$a=qb+r$$且 $0\le r<b$。

由于存在性的证明在书上已有提及，故我们只需要证明：
①$0\le r<b$
②$q$和$r$的唯一性
①假设$r\ge b$，则：
对于集合S： S = { a − b k : k ∈ Z 且 a − b k ≥ 0 } S=\left \{ a-bk:k\in \Z且a-bk\geq 0 \right \} S={a−bk:k∈Z且a−bk≥0}
$r$是其中的最小元且满足$r=a-qb$.
由于$r\ge b$，则：
$$r^{\prime }=a-(q+1)b\ge 0$$
因而$r^{\prime }\in S$且$r^{\prime }<r$，这与$r$是最小元矛盾.
故假设不成立，原命题成立.
②由$q$的唯一性不难知道$r=a-qb$同样具有唯一性,
故只需证明$q$的唯一性.
不妨假设存在${q_1}^{},{r_1}^{},{q_2}^{},{r_2}^{}\in Z,(q1\ne q2,r1\ne r2)$使：
$$a={q_1}^{}b+{r_1}^{}={q_2}^{}b+{r_2}^{}$$
则：
$$\left|{r_1}^{}-{r_2}^{}\right|=\left|{q_1}^{}b-{q_2}^{}b\right|\ge b$$
这显然与$0\le {r_1}^{},{r_2}^{}<b$矛盾,从而假设不成立.
故$q_1=q_2$，$q$的唯一性得证.

第二章

迭代版本的gcd

int gcd(int a,int b)
{
    while(b)
    {
        int t = a % b;
        a = b;
        b = t;
    }
    return a;
}

迭代版本的egcd

//求gcd(a,b)的Bezout系数x和y,即xa+yb=gcd(a,b)
int egcd(int a, int b, int *x, int *y)
{
    int r1 = 1, s1 = 0, r2 = 0, s2 = 1;
    while(b)
    {
        int q = a / b, t = a % b, tr = r1, ts = s1;
        r1 = r2, s1 = s2, a = b;
        r2 = tr - q * r2, s2 = ts - q * s2, b = t;
    }
    *x = r1, *y = s1;
    return a;
}

递归版本的egcd

//从最深处(即b=0)开始递归计算
//递归从bx+(a%b)y=d得到ax'+by'=d
//其中，x'=y, y'=x-(a/b)*y
int exgcd(int a, int b, int *x, int *y)
{
    if(!b)
    {
        *x = 1, *y = 0;
        return a;
    }
    int d = exgcd(b, a % b, x, y);

	//完成x'和y'的赋值
    int t = *x;
    *x = *y;
    *y = t - a / b * (*x);
    return d;
}

欧拉函数

int phi(int n)
{
	int res = 0;
	for(int i = 1; i < n; i ++ )
	{
		if(gcd(i, n) == 1) ++ res;
	}
	return res;
}

课后习题6

题目：假设$g^a\equiv 1(\mathrm{mod}m)$且 $g^b\equiv 1(\mathrm{mod}m)$，请证明:
$$g^{gcd(a,b)}\equiv 1(\mathrm{mod}m)$$
法一：无穷递降法
详细见链接

法二：利用同余运算的性质
先介绍一个引理(其证明思路放在法二的最后)

$g^a\equiv 1(\mathrm{mod}m)\Rightarrow g^{ka}\equiv 1(\mathrm{mod}m),k\in \mathbb{Z}$

由Bezout定理有：
$$\exists r,s\in \mathbb{Z},使ar+bs=gcd(a,b)$$
由上述引理不难得到：
$$g^a\equiv 1(\mathrm{mod}m)\Rightarrow g^{ar}\equiv 1(\mathrm{mod}m)$$
$$g^b\equiv 1(\mathrm{mod}m)\Rightarrow g^{bs}\equiv 1(\mathrm{mod}m)$$
从而有：
$$g^{ar+bs}\equiv g^{ar}\cdot g^{bs}\equiv 1(\mathrm{mod}m)$$
即：
$$g^{gcd(a,b)}\equiv 1(\mathrm{mod}m)$$

---

引理的证明：
①当$k\ge 0$时：
$$g^{ka}\equiv \begin{array}{c}\underbrace{g^a\cdot g^a\cdots g^a}\\ k\end{array}\equiv 1(\mathrm{mod}m)$$
②当$k<0$时：
由于$g^a\equiv 1(\mathrm{mod}m)$，故：$g^{-a}\equiv 1(\mathrm{mod}m)$
进而有，
$$g^{ka}\equiv (g^{-a}{)}^{-k}\equiv \begin{array}{c}\underbrace{(g^{-a})\cdot (g^{-a})\cdots (g^{-a})}\\ -k\end{array}\equiv 1(\mathrm{mod}m)$$
综上所述，引理成立。

---

法三：从循环群得到的一点小启发：
构造如下集合S:
S = { a : g a ≡ 1 ( m o d m ) , a > 0 } S=\left \{a:g^a\equiv1\pmod m ,a>0\right \} S={a:ga≡1(modm),a>0}
不难知道：S是全体自然数集合${\mathbb{N}}^{\ast }$的一个子集.
由良序原则知：S中一定存在最小元，记作$x$.
不妨假设$x\nmid a$，设$a=qx+r,q,r\in \mathbb{Z},0<r<x$，则：
$$g^a\equiv g^{qx+r}\equiv g^r\equiv 1(\mathrm{mod}m)$$
不难发现$r\in S$，但是$r<x$，这就与$x$是最小元矛盾.
故假设不成立，即$x\mid a$.
同理可证得$x\mid b$，进而有：
$$x\mid a,x\mid b\Rightarrow x\mid gcd(a,b)$$
不妨设$gcd(a,b)=kx,k\in \mathbb{Z}$，利用同余运算的性质有：
$$g^{gcd(a,b)}\equiv g^{kx}\equiv 1(\mathrm{mod}m)$$
证毕.

课后习题8

证明：如果$gcd(a,b)=d$，则 $gcd(a/d,b/d)=1$。
证明思路：利用Bezout定理即可。
$$gcd(a,b)=d\Rightarrow \exists x,y\in \mathbb{Z},xa+yb=d$$
由于$d=gcd(a,b)$，则：$d|a$，$d|b$
因而：
$$x\frac{a}{d}+y\frac{b}{d}=1$$
即：$$\exists x,y\in \mathbb{Z},x\frac{a}{d}+y\frac{b}{d}=1$$
由Bezout定理知：$$gcd(\frac{a}{d},\frac{b}{d})=1$$

---

事实上，对原命题稍加改动可以使该命题是充分且必要的，即：
$$gcd(\frac{a}{d},\frac{b}{d})=1,d>0\iff gcd(a,b)=d$$
$\Leftarrow$已证，我们证明$\Rightarrow$：
由Bezout定理有：
$$\exists x,y\in \mathbb{Z},x\frac{a}{d}+y\frac{b}{d}=1$$
即：
$$xa+yb=d$$
再次使用Bezout定理可知：$gcd(a,b)|d$
注意隐含条件$d|a,d|b$，从而有：$d|gcd(a,b)$
又$d>0$，显然有：$d=gcd(a,b)$