高精度加法的优化

改善高精度运算的效率

•以上接触到的高精度存储方法是用一个整型数组来表示一个很大的数,数组中的每一个数表示一位十进制数字。但这种方法的缺点是,如果十进制数的位数很多,则对应数组的长度会很长,并增加了高精度计算的时间。那么有什么方法可以改善高精度运算的效率呢?

 

•我们可以考虑用一个数记录2位数字、3位数字或更多位数字。理论上来说,数组中的每个数表示的数字越多,数组的长度就越短,程序运行的时间也就越短。但是,我们还需考虑到计算机中的一个数的取值范围,必须保证它们在运算过程中不会越界。在权衡了两方面的情况后得出:如果只有加减法运算，那么可以用一个longint记录9位数字（十亿进制），如果有加减乘除运算，那么用一个longint记录4位数字（万进制）是最佳方案。我们以万进制为例，来看一下优化的方法：

 

var
 s1,s2:string;
procedure plus(s1,s2:string);
var
 a,b,c:array[1..500]of longint;
 i,j1,j2,l,k,p,code:longint;
begin
 fillchar(a,sizeof(a),0);
 fillchar(b,sizeof(b),0);
 fillchar(c,sizeof(c),0);
 j1:=0; l:=length(s1);
 while l>0 do
 begin
  inc(j1);
  val(copy(s1,l-3,4),a[j1],code);
  s1:=copy(s1,1,l-4); l:=length(s1);
 end;
 j2:=0; l:=length(s2);
 while l>0 do
 begin
  inc(j2);
  val(copy(s2,l-3,4),b[j2],code);
  s2:=copy(s2,1,l-4); l:=length(s2);
 end;
 if j1>j2 then p:=j1 else p:=j2;
 k:=0;
for i:=1 to p do
 begin
  c[i]:=a[i]+b[i]+k;
  k:=c[i] div 10000;
  c[i]:=c[i] mod 10000;
 end;
 if k=1 then
 begin
  inc(p); inc(c[p]);
 end;
 write(c[p]);
 for i:=p-1 downto 1 do
 begin
  if c[i]<10 then write('000',c[i]) else
  if c[i]<100 then write('00',c[i]) else
  if c[i]<1000 then write('0',c[i]) else
  write(c[i]);
 end;
 writeln;
end;
begin
 readln(s1);
 readln(s2);
 plus(s1,s2);
end.