高精度加减乘除阶乘计算的算法及代码.pdf

高精度计算 主要的方法是利用数组模拟计算比如： 高精度加法 12345678910111213 + 1111111111111111111 开两个数组存储： a[]={3,1,2,1,1,1,0,1,9,8,7,6,5,4,3,2,1}; b[]={1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1}; 两个数组分别把数值倒存,在一位一位的加,每位加后判断是否大于 10,在进位(注:如果太大的数 值,可以考虑4 位一存哦.) 注意下面的a1,b1,c1 为 数组的长度 if a1b1 then c1:=a1 else c1:=b1; for i:=1 to c1 do begin c[i]:=a[i]+b [i]+c[i]; c[i+1]:=c[i+1]+c[i]div 10; c[i]:=c[i] mod 10; end; while c[c1+1]0 do begin inc(c1); c[c1+1]:=c[c1+1]+c[c1]div 10; c[c1]:=c[c1] mod 10; end; 四位一存: var s:ansistring ; a,b:array [1..11000]of int64; x:array [1..51000]of int64; c:array [1..35000]of int64; i,j,a1,b1,c1,d,x1:longint; begin while not eof do begin fillchar(c,sizeof(c),0); readln (s); x1:=length(s); for i:=1 to x1 do x[i]:=ord (s[x1-i+1])-48; for i:=1 to (x1 div 5) do a[i]:=x[5*i]*10000+x[5*i-1]*1000+x[5*i-2]*100+x[5*i-3]*10+x[5*i-4]; a1:=x1 div 5; if x1 mod 50 then begin d:=0; for i:=1 to x1 mod 5 do d:=d*10+x[x1-i+1]; a1:=a1+1; a[a1]:=d; end; readln (s); fillchar(x,sizeof(x),0); x1:=length(s); for i:=1 to x1 do x[i]:=ord (s[x1-i+1])-48; for i:=1 to (x1 div 5) do b [i]:=x[5*i]*10000+x[5*i-1]*1000+x[5*i-2]*100+x[5*i-3]*10+x[5*i-4]; b1:=x1 div 5; if x1 mod 50 then begin d:=0; for i:=1 to x1 mod 5 do d:=d*10+x[x1-i+1]; b1:=b1+1; b [b1]:=d; end; for i:=1 to a1 do for j:=1 to b1 do c[i+j-1]:=c[i+j-1]+a[i]*b [j]; c1:=a1+b1-1; for i:=1 to c1 do if c[i]99999 then begin c[i+1]:=c[i+1]+c[i] div 100000; c[i]:=c[i] mod 100000; end; while c[c1+1]0 do