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

超大数字超高精度几何运算 　　【摘 要】针对计算机精度位数的限制，按照位运算原理，创意设计加、减、乘、除和乘方的大数五则运算新算法。给出具体的C语言程序实现代码，可以精确计算加数、乘数、被减数、减数、被除数、除数、底数在十进制位数(含小数、负数)1000位以内，指数在十进制位数(正整数)8位以内，和数、积数、差数、商数、乘方等结果(含小数、负数)在十进制位数10000位以内的运算。 　　【关键词】大数 算法原理 编程 　　 　　1 研究背景 　　 　　华罗庚说:“宇宙之大，粒子之微，火箭之速，化工之巧，地球之变，生物之谜，日用之繁，无处不用数学。” 　　在生活中，计算是一件很常见的事。在大多数时间，我们可以利用手中的计算器或者用计算机操作系统自带的计算器进行一些十分精确的计算，但这些计算器计算的数字却不能够超过32位。实际上，“宇宙之大，粒子之微”，有许许多多的方面，需要我们使用非常巨大或者极其微小的数据进行计算，像天文学、物理学、地理学、生物学、力学等。此外，在有限元、流体力学计算、天气和气候预报、基因工程、蛋白质和分子分析模拟、石油储量分析等应用领域，高精度计算也是很常用的。当然，在科学和数学方面，则更加会频繁地使用高精度的、大数字的高速加减乘除和乘方运算。在历届的数学竞赛里，也总会出现一些较大数字的运算，如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