pascal高精度运算( 乘法)

《pascal高精度运算( 乘法)》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第五讲:高精度运算乘法篇高精度运算涉及到的问题：1、数据的输入。字符串2、数据的存储。数组3、数据的运算：进位和借位。4、结果的输出：小数点的位置、处理多于的0等。整数i乘多精度数组a设x为当前位乘积和进位x:=x+a[j]*i;a[j]:=xmod10;x:=xdiv10精确计算n的阶乘n!（7

2、之间，j为数组变量a的下标，取值范围在1至100之间，X存放来自低位的进位数。Xj1009998…321constmax=100;n=20;vara:array[1..max]of0..9;i,j,k;x:integer;begink:=1;a[k]:=1;{a=1}fori:=2tondo{a*2*3….*n}beginx:=0;{进位初始化}forj:=1dokdo{a=a*i}beginx:=x+a[j]*i;a[j]:=xmod10;x:=xdiv10end;whilex>0do{处理最高位的进位}begink:=k+1;a[k]:=xmod10;x:=xdiv10ende

3、nd;writeln;fori:=kdowento1write(a[i]){输出a}end.三、高精度乘法1、高精度乘单精度（整数），求a*b，a<=10200,b<=108vari,len:integer;s:string;b,m:longint;a:array[1..250]oflongint;beginreadln(s);{读入被乘数a}readln(b);{读入乘数b}len:=length(s);fori:=1tolendoa[i]:=ord(s[len-i+1])-48;fori:=1tolendo{每项乘b}a[i]:=a[i]*b;fori:=1tolendo{处理

4、进位}begina[i+1]:=a[i+1]+a[i]div10;a[i]:=a[i]mod10;end;fori:=lendownto1dowrite(a[i]);end.{multiply}m:=a[len+1];{处理最高进位}whilem<>0dobegin{处理最高位的进位}inc(len);a[len]:=mmod10;m:=mdiv10;end;ifa[len+1]>0theninc(len);？乘法运算a←a*b(a、b为numtype类型)1、积的位数为la+lb-1或者la+lb；2、如果暂且不考虑进位关系，则ai*bj应该累加在积的第j+i-1位上：x:=a[

5、i]*b[j]+xdiv10+c[i+j-1];c[i+j-1]:=xmod10;3、可以先乘、后处理进位2、高精度乘高精度2、高精度乘高精度vari,j,la,lb,len:integer;s1,s2:string;m:longint;a,b,c:array[1..250]ofinteger;beginreadln(s1);la:=length(s1);fori:=1toladoa[i]:=ord(s1[la-i+1])-48;readln(s2);lb:=length(s2);fori:=1tolbdob[i]:=ord(s2[lb-i+1])-48;fori:=1tolado

6、forj:=1tolbdoc[i+j-1]:=c[i+j-1]+a[i]*b[j];len:=la+lb;fori:=1tolendobeginc[i+1]:=c[i+1]+c[i]div10;c[i]:=c[i]mod10;end;whilec[len]=0dodec(len);m:=c[len];whilem>0dobeginc[len]:=mmod10;m:=mdiv10;inc(len);end;fori:=len-1downto1dowrite(c[i])End.四、高精度除法{高精度a除以单精度b}typedata=array[0..300]ofinteger;vara

7、,c:data;s1,s2:string;len,k:integer;i:integer;b,d:longint;procedureinit;beginreadln(s1);readln(b);fillchar(a,sizeof(a),0);len:=length(s1);fori:=1tolendoa[i]:=ord(s1[len-i+1])-48;end;proceduredevide;{c:=adivb;d:=amodb}vari:integer;beginfi