高精度计算地若干问题.pptVIP

高精度计算的若干问题 目 录： 一、高精度计算的若干问题 1. 高精度乘法与分治算法 2. 牛顿迭代法 3. 高精度除法 4. 高精度开方 5. π值计算 1. 高精度乘法与分治算法 高精度计算的核心运算是高精度乘法。加减法的计算复杂度是O(n) 阶的，不会有阶上的改进，除法、开方等运算可转换为乘法。 两个高精度数的乘法一般是O(n2) 阶的，采用快速Fourier变换(FFT)或快速数论变换(FNT)可降至O(n*logn) 阶,其算法较复杂。这里介绍一种分治算法，可将计算复杂度降至O(n1.58) 阶。 设 A, B 是两个 N(N=2n) 位数，且 A=a110N/2+a2 , B=b110N/2+b2 则：AB=(a110N/2+a2 )(b110N/2+b2) =a1b110N+(a1b2+a2 b1)10N/2+a2b2 =p10N+(r – p – q)10N/2+q 其中，p=a1b1, q=a2b2, r=(a1+a2 )(b1+b2) 这样，计算两个N位数的乘法 AB，可转换为计算３个N/2位数的乘法(p, q, r)。 　设Tn为 N=2n 时所做的１位乘法的次数，则 Tn =3Tn-1 =32Tn-2 =…= 3nT0 =3n, 由于　N=2n，n=log2N, Tn=3logN=Nlog3≈ N1.58