数论基础及应用课件.pptVIP

数论基础及应用 数论是研究数的性质的学科是一门古老而充满现代魅力的数学学科。数论基本上可分为初等数论、解析数论、代数数论等几个较大的分支。 在古代，我国对数论的研究曾有过辉煌的成就, 如孙子定理（国外文献一般称为中国剩余定理）、商高定理(勾股数)、圆周率的计算等等。在现代，我国一些著名的数学家，如华罗庚、王元、陈景润、潘承洞、丁夏畦等都在数论领域做出了一些举世公认的重要成果。 1．欧几里德转辗相除法 利用gcd(a,b)=gcd(b,a mod b)求a，b的最大公约数： Function gcd(a,b:longint):longint;beginif b=0 then gcd:=a Else gcd:=gcd(b,a mod b);?end; 思考：如何把上述算法写成迭代形式？ 2．扩展的欧几里德算法 如果gcd(a,b)=d，一定存在整数x和y满足gcd(a,b)=ax+by。 求d及满足gcd(a,b)=ax+by的整数对（x，y）的算法 function exgcd(a,b:longint;var ,y:longint):longint;vart:longint;beginif b=0 then??begin? exgcd:=a;? x:=1;? y:=0;?end 求d及满足gcd(a,b)=ax+by的整数对（x，y）的算法（续） else?begin?exgcd:=exgcd(b,a mod b,x,y);?t:=x;?x:=y;?y:=t-(a div b)*y;?end; end; 应用1：求解二元一次不定方程ax+by=c整数解 解二元一次不定方程 ax+by=c ① 其中a,b,c都是整数，所求的解(x,y)也是整数 关于方程①的可解性，有下面的两个重要的结论： （1）设gcd(a,b)表示整数a,b的最大公约数。方程①有解的充分必要条件是gcd(a,b)|c。（记号“x|y”表示x能整除y，即存在整数k，使y=kx）。 （2）如果(x0,y0)是方程①的一组解，则对任何整数t,(x0+bt,y0-at)也都是方程①的解。 下面我们讨论具体求解的方法。 为了避免计算中对负数和0的讨论，我们假定a0，b0，并且a=b。 假定方程①有解，即系数满足：gcd(a,b)|c，这时，c’=c/gcd(a,b)一定是个整数。我们先讨论下面的方程： ax+by= gcd(a,b) ② 根据上述扩展的欧几里德算法，一定存在整数x0和y0满足ax+by =gcd(a,b)。 显然，如果(x0,y0)是方程②的一组解，则(c’x0,c’y0)也是方程①的一组解，即 a(c’x0)+b(c’y0)=(c’f)=c。 求二元一次不定方程ax+by=c一组整数解（x0,y0）的算法 procedure equation(a,b,c:longint;var x0,y0:longint);var d,x,y:longint;begin?d:=exgcd(a,b,x,y);{参见扩展的欧几里德算法}?if c mod d0 then?begin? writeln(no answer);? halt;?end else?begin? x0:=x*(c div d);? y0:=y*(c div d);?end;end; 递推法求二元一次不定方程ax+by=c一组整数解（x0,y0） 我们先讨论下面的方程： ax+by=f ②其中f=gcd(a,b) 例如 方程 107x+73y=1 ③ 其中a=107,b=73,f=1我们用类似于求最大公约数的辗转相除的方法求这个解。利用辗转相除，可以得到： 107=73*1+34 (1) 73=34*2+5 (2) 34=5*6+4 (3) 5=4*1+1 (4) 4=1*4 (5) 递推法求二元一次不定方程ax+by=c一组整数解（x0,y0）(续） 为了消去(3)中的”4”，令 (3)*1-(4): 34=5*7-1 (6) 为了消去(2)中的”5”，令 (2)*7-(6): 73*7=34*15+1 (7) 为了消去(1)中的”34”，令(1)*15-(7): 107*15=73*22-1, 即：107*(-15)+73*22=1,于是，③的一组解为（-15，22）。 下面归纳一般的算法：