超大数字超高精度的几何运算.docVIP

超大数字超高精度的几何运算 2??1000?-2??999?2??998?。虽然通过一些很巧妙的方法，可以计算出这道题目的答案应该是2，但这无法用直接计算的方法去验证结果。想用家用电脑验证，却因为数字太大超出计算范围而无法计算，或者只能得到一些近似的结果。 　　本程序按照位运算原理，创意设计加、减、乘、除和乘方的大数五则运算新算法。可以精确计算加数、乘数、被减数、减数、被除数、除数、底数在十进制位数(含小数、负数)1000位以内，指数在十进制位数(正整数)8位以内，和数、积数、差数、商数、乘方等结果(含小数、负数)在十进制位数10000位以内的运算。使用了C语言中的malloc函数，与普通数组进行对比，可改动答案和操作数的最大位数，且可避免因某个位数较大而出现的“未知错误”。因此，这种方法使得程序更具通用性、稳定性与可扩充性，保证程序正常运行。 　　2 算法原理 　　根据四则运算的原理和位运算原理设计了加减乘除和乘方的新算法，这些运算的算法介绍如下。 　　2.1 乘法算法原理 　　将两个乘数的每一位转换为单个的数字，存放于数组中。两数从左边第一位开始编号分别是:b(0)，b(1)，b(2)，b(3)……b(S1-1)；c(0)，c(1)，c(2)，c(3)……c(S2-1)，若b(I1)位(从左到右)与c(I2)位(从左到右)相乘的数将增加到答案的a(I3)位(从右到左，与正常思路相反)，则可以得到I3=(S1-I1)+(S2-I2)-2，只要再通过进位就可以很轻松得到结果。(见图1) 　　乘法示例：16*32 　　（1）分解到数组1 6 * 3 2(S1=2,S2=2),将答案数组初值赋零； 　　（2）I1 = 0 TO S1-1 ,I2 = 0 TO S2-1分别计算I3 = (S1-I1)+(S2-I2)-2，可以得到: 　　I1 = 0,I2 = 0 I3 = 2第二位(从右到左，与正常思路相反)增加1*3=3；第二位为3 　　I1 = 0,I2 = 1 I3 = 1第一位(从右到左，与正常思路相反)增加1*2=2；第一位为2 　　I1 = 1,I2 = 0 I3 = 1第一位(从右到左，与正常思路相反)增加6*3=18；第一位为20 　　I1 = 1,I2 = 1 I3 = 0第零位(从右到左，与正常思路相反)增加6*2=12；第零位为12 　　（3）进位 K=0 I = 0 TO R(R为答案的最大个数，下同)分别计算第I位(从右到左,与正常思路相反)上一位进的数和应该进下一位的数: 　　I=0第零位(从右到左，与正常思路相反)为12+K=12应进的位为K=FLOOR(12/10)=1第零位进位后为12MOD10=2 　　I=1第一位(从右到左，与正常思路相反)为20+K=21应进的位为K=FLOOR(21/10)=2第一位进位后为21MOD10=1 　　I=2第二位(从右到左，与正常思路相反)为3+K=5应进的位为K=FLOOR(5/10)=0第二位进位后为5MOD10=5 　　I=3第三位(从右到左，与正常思路相反)为0+K=0应进的位为K=FLOOR(0/10)=0第三位进位后为0MOD10=0 　　…… 　　I=R第R位(从右到左,与正常思路相反)为0+K=0应进的位为K=FLOOR(0/10)=0第R位进位后为0MOD10=0 　　最后得到的结果为:2 1 5 0 0 0…… 0 　　（4）打印，从结果的最后一位开始寻找不为0的数,输出这之后的所有数，本题为512。如有小数，小数位数为两数小数位数的总和。（见图2) 　　2.2 乘方算法原理 　　和乘法的思想一样，x^y=x*x*x*……*x(y个x相乘)，运用乘法原理加以实现。 　　乘方示例：16^3 　　（1）分解到数组1 6和数字3,将答案数组初值赋零，计数器K = 0 　　（2）将1 6存入答案数组,K + 1 　　（3）计算1 6 * 1 6并将结果存入答案数组，答案数组为：2 5 6,K + 1 　　（4）计算2 5 6 * 1 6并将结果存入答案数组，答案数组为：4 0 9 6,K + 1 　　（5）K = 3，结束。 　　（6）打印，从最后一位开始寻找不为0的数，输出这之后的所有数，本题为4096。小数位数为底数小数位数乘指数。 　　2.3 加法算法原理 　　首先，需要将b、c两个加数的小数点对齐(先计算两数的位数S1,S2，然后计算两数的小数位数L1,L2)，然后，就可以通过每一位相加，并进位的方法得出答案。不失一般性，设l1≤l2,i=s1-l1-s2+l2,i≥0: 　　加法示例：156.78+78.156 　　（1）分解到数组 1 5 6 7 8 + 7 8 1 5