数学知识学习课件.pptVIP

计算a和b的最大公约数gcd(a,b) 根据欧几里德定理 〔证明〕 关键是证明gcd(a，b)与gcd(b，a mod b)可互相整除。设ｄ＝gcd(a，b)，将a mod b 化成a与b的线性组合 a mod b＝a－a div b*b。由于d能整除a和b，因此d一定能整除a与b的线性组合，即d能整除(a mod b)。又因为ｄ能整除b和(a mod b)，显然d能整数gcd(b，a mod b)。同理可证gcd(b，a mod b)能整除gcd(a，b)。由于gcd(a，b)与gcd(b，a mod b)可互相整除，蕴涵着gcd(a，b)＝gcd(b，a mod b)，因此欧几里得公式成立。 计算n的最大互质数 首先必须明确一个基本事实 gcd(N-1,N)=1，或者说gcd(N-p,N)≤p 我们一般可根据n的奇偶性分析，在[1…n]中找出与n互质的最大数。 传球游戏 N个人围圈玩传球游戏，开始时第一个人拿着球，每个人把球传给左手的第K个人。满足1≤K≤N/2。求K的最大值，使得第一个人重新拿到球之前，每个人都拿过球。 输入：数串n（3 ≤ N ≤ 10200) 输出：K的最大值 规律 如果人数n和间隔k互质，则无论从哪个位置开始传球，都能够使得每一个人都拿到球。 数学模型 数学语言描述 给定一个整数N，要求一个尽可能大的整数K，满足1≤K≤N/2，且gcd(K, N) = 1。 基本事实 gcd(N-1, N)=1，或者说gcd(N-p, N)≤p 在N为奇数的情况下 ∵ K*2=N-1 ，gcd(K*2, N) = 1 ∴可以得到gcd(K, N)=1 ∵ K的上限是[N/2] ∴直接输出[N/2]即可。 在N为偶数的情况下 第一种情况:N是偶数且为4的倍数 ∵N/2是偶数，gcd（ N/2 ，N）≠1, ∴ K= N/2不可能符合要求； ∵ N/2-1是奇数， gcd(N/2-1, N/2)=1。∴K= N/2-1满足题意。 第二种情况: N是偶数但不是4的倍数 ∵N/2-1也是偶数， ∴ K= N/2-1不可能是正确答案; ∵ N/2-2是奇数。 根据刚才的介绍，有gcd(N/2-2, N/2)≤2， 但是N/2是奇数，所以gcd(N/2-2, N/2)=1。 最后同样可以得到gcd(N/2-2, N)=1。 ∴ K= N/2-2满足题意。 注意:n和answer为高精度数组 if n=2 then begin answer:=1; exit; end; answer := n div 2; //若n div 2是奇数，则返回n div 2; if n是偶数 then begin answer := answer-1; //n是2的倍数，则减1； if answer是偶数 then answer:=answer–1; //n不是4的倍数，再减1 end； 输出answer; 计算方程ax+by=c的初始解x0、y0和gcd(a,b) 根据欧几里德定理 bx’+(a-a div b*b)y’=c ay’+b(x’-a div b*y’)=c x=y’ y=x’-a div b*y’ 一直递归下去,直至a=0、c mod b=0为止。此时，gcd(a,b)=b, x=0, y=c div b 计算和返回方程ax + by =c的x0 , y0、a和b的最大公约数gcd(a,b)。 function get_root(A,B,C:int64; var x,y,gcdnum:int64):boolean; {计算和返回方程ax + by =c的x0 , y0、a和b的最大公约数gcdnum } var tx , ty : int64; begin if A = 0 then if C mod B = 0 then begin gcdnum:=B; x:=0; y:=C div B; get_root := true; end else get_root := false else if get_root(B, A mod B, C , tx , ty , gcdnum) then begin x:= ty; y := tx - (A div B) * ty; get_root := true; end else get_root := false; end;{ get_root } 若无解，则失败返回；否则计算 计算和返回方