other_____多精度数值的处理.pdfVIP

多精度数值的处理 朱全民 Free pascal中整数类型  Shortint : -128..127  SmallInt (Integer) :-32768 .. 32767 ,2B  Longint:-2147483648 .. 2147483647,4B  Int64 : -9223372036854775808 .. 9223372036854775807 , 8B  Byte : 0..255, 1B  Word : 0..65535, 2B  Cardinal (Longword) : 0..4294967295,4B  QWord: 0 ..18446744073709551615 ,8B  Boolean (ByteBool) : false .. true ,1B  wordbool: false .. true , 2B  LongBool : false .. true , 4B  Char : 0..255,1B 多精度数值处理需解决的问题 1.数据存储  用数组或字符串 2.计算结果位数确定  用数学方法确定。利用对数函数（log）。 3.进位与借位处理  加法：if a[i]10 then begin a[i]=a[i]-10 ； a[i+1]=a[i+1]+1  减法：if a[i]b[i] then begin a[i]=a[i]+10 ； a[I1]=a[i+1]-1  乘法：y= a[i]*b[i]+c; c= y div 10 ;a[i]=y mod 10 4.商与余数处理  除法：商=A DIV B, 余数=A MOD B 多精度数值相加  如果我们用数组A 、B分别存储加数和被加数，用数组C存储结果。  C[i]:= A[i]+B[i];  if C[i]10 then begin C[i]:= C[i] mod 10; C[i+1]:= C[i+1]+1 end;  因此，算法描述如下：  procedure add(a,b;var c);  { a,b,c都为数组，a存储被加数，b存储加数，c存储结果}  var i,x:integer;  begin  i:=1  while (i=a数组长度0) or(i=b数组的长度) do begin  x := a[i] + b[i] + x div 10; {第i位相加并加上次的进位}  c[i] := x mod 10; {存储第i位的值}  i := i + 1 {位置指针变量}  end  end; 多精度数值相减  类似加法，可以用竖式求减法。在做减法运算时，需要注意的是：被减 数必须比减数大，同时需要处理借位。  因此，可以写出如下关系式  if a[i]b[i] then begin a[i+1]:=a[i+1]-1;a[i]:=a[i]+10 end  c[i]:=a[i]-b[i]  因此，算法描述如下：  procedure minu (a,b;var c);  { a,b,c都为数组，c存储结果}  var i,x:integer;  Begin  处理被减数和减数, a存储被加数，b存储加数  i:=1  while (i=a数组长度0) or(i=b数组的长度) do begin  x := a[i] - b[i] + 10+ x; {不考虑大小问题，先往高位借10}  c[i] := x mod 10; {存储第i位的值}  x := x div 10 - 1; {将高位借掉的1减去}  i := i + 1 {位置指针变量}  end  end; 多精度数值相乘  C i跟A[i]*B[j]乘积有关，跟上次的进位有关，还跟原C i的值有关，分析下 标规律，有  x:= A[i]*B[j]+ x div 10+ C[i+j-1]