高精度计算是一种程序设计的算法。由于中央处理器的字长限制,如32位CPU中一个整数最大只能取值4,294,967,295(=232-1)。因此在进行更大范围的数值计算中,往往要采取模拟手段。通常通过分离字符的方法通过数字数组进行输入、通过数组倒序输出、通过模拟竖式计算进行计算。一般而言,主要模拟的是按位运算,可以用不同的进位制达成不同的目的。
有许多程序库支持高精度计算,最著名的是GNU多重精度运算库。另外,Java,Python和Pascal也有本地的高精度运算支持。
高精度计算的一个常见应用是公开密钥加密,这些算法经常对长度上百位的整数进行运算。高精度计算的另一个应用是在需要没有人为限制位数和没有算术溢出的情况下使用。在检查固定精度计算的结果以及确定公式中系数的精确值或近似值时,高精度计算也很有用。比如,在高斯求积中,我们需要确定的值。
高精度加法是信息学的一种重要算法。这种算法使用多个存储单位进行计算,因此它的计算范围超过一般使用一个存储单位的算法。也是一些信息学竞赛的常考题目。
以358934760892734899+38960302975237462为例:
1、计算结果的位数
358934760892734899共18位
38960302975237462 共17位
故结果不会超过19位。
2、将要计算的数字分割成多段,按照顺序排列(这里以0-32767作为每一存储单位存储的数的限制):
(为提高空间利用效率,可以一个存储单位存储多位数。)
3、将两数相加。
4、输出结果。
从高位到低位依次输出。除最高位以外,其他低位上不足4位的要在前面补上0。
pascal:
var a,b,c:array of integer; n:string; lena,lenb,lenc,i,x:integer; begin readln(n); lena:=length(n); for i:=1 to lena do a:=ord(n)-ord('0'); readln(n); lenb:=length(n); for i:=1 to lenb do b:=ord(n)-ord('0'); i:=1; x:=0; while (i<=lena) or(i<=lenb) do begin c:=a+b+x; x := c div 10; c := c mod 10; i := i + 1; end; if x>0 then begin lenc:=i; c:=x; end else lenc:=i-1; for i:=lenc downto 1 do write(c); end.