高精度计算-安阳一中.DOC

安阳一中信息学奥赛辅导资料 第 PAGE 4 页 共 NUMPAGES 5 页 高精度计算(2) 教学目标 熟练掌握高精度乘除计算 重点难点分析 高精度乘法的实现； 教具或课件 使用多媒体演示文稿 主要教学过程 引入新课 本节主要介绍高精度乘法的相关内容。 教学过程设计 【例3】 高精度乘法。从键盘读入两个正整数，求它们的积。 分析：（1）、乘法运算a←a*c(a为高精度类型,c为字节型) 按照乘法规则,从a的第l位开始逐位与c相乘。在第i位乘法运算中(1≤i≤la),a的i位与c的乘积必须加上i-1位的进位(i-1位的乘积除以10的整商),然后规整积的i-l位(取i-1位的乘积对10的余数)。 procedure multiply(var a:numtype; c:byte); var i:byte； begin a[1] ←a[l]*c； {第1位初始化} for i←2 to la do {逐位相乘} begin a[i] ←a[i]*c； a[i] ←a[i]+a[i-l] div 10； a[i-1] ←a[i-l] mod 10 end:{for} while a[1a]=10 do {积的最高位进位} begin la←la+1; a[la] ←a[la-1] div 10; a[la-1] ←a[la-1]mod 10; end; {while} end;{multiply} （2）、乘法运算c←a*b(a,b为高精度类型,) 类似加法，可以用竖式求乘法。在做乘法运算时，同样也有进位，同时对每一位进乘法运算时，必须进行错位相加，如图3, 图4。 分析C 数组下标的变化规律，可以写出如下关系式：C i = C i +C i +…由此可见，C i跟A[i]*B[j]乘积有关，跟上次的进位有关，还跟原C i的值有关，分析下标规律，有 x:= A[i]*B[j]+ x DIV 10+ C[i+j-1]; C[i+j-1] := x mod 10; 类似，高精度乘法的参考程序： program exam3; const max=200; var a,b,c:array[1..max] of 0..9; n1,n2:string; lena,lenb,lenc,i,j,x:integer; begin write(’Input multiplier:’); readln(n1); write(’Input multiplicand:’); readln(n2); lena:=length(n1); lenb:=length(n2); for i:=1 to lena do a[lena-i+1]:=ord(n1[i])-ord(’0’); for i:=1 to lenb do b[lenb-i+1]:=ord(n2[i])-ord(’0’); for i:=1 to lena do begin x:=0; for j:=1 to lenb do begin {对乘数的每一位进行处理} x := a[i]*b[j] + x div 10 + c[i+j-1]; {当前乘积+上次乘积进位+原数} c[i+j-1] := x mod 10; end; c[i+j]:= x div 10; {进位} end; lenc:=i+j; while (c[lenc]=0) and (lenc1) do dec(lenc); for i:=lenc downto 1 do write(c[i]); writeln end. 【例4】 高精度除法。从键盘读入两个正整数，求它们的商（做整除）。 分析：做除法时，每一次上商的值都在０～９，每次求得的余数连接以后的若干位得到新的被除数，继续做除法。因此，在做高精度除法时，要涉及到乘法运算和减法运算，还有移位处理。当然，为了程序简洁，可以避免高精度乘法，用0～9 次循环减法取代得到商的值。这里，我们讨论一下高精度数除以单精度数的结果，采取的方法是按位相除法。 program exam4; const max=200; var a,c:array[1..max] of 0..9; x,b:longint; n1,n2:string; lena:integer; code,i,j:integer; begin write(’Input dividend:’); readln(n1); write(’Input divisor:’); readln(n2); lena:=le