7第七讲高精度计算.pptVIP

高精度计算 北大李文新《在线程序设计….》第七章 基本类型的表示范围 大整数加法 大整数减法 大整数乘法 基本类型的表示范围 int 类型 32位： 能表示的十进制数不超过10位 double 类型 64位 有效数字15-16位 如果需要100位甚至1000位整数，或要求精确到小数点后100位，如何解决？ 用数组存储，模拟手工计算 我们一般称这种基本类型无法表示的整数为大整数 7.1 例题：大整数加法 问题描述 求两个不超过200 位的非负整数的和。 输入数据 有两行，每行是一个不超过200 位的非负整，没有多余的前导0。 输出要求 一行，即相加后的结果。结果里不能有多余的前导0，即如果结果是342，那么就不能输出为0342。 输入样例 22222222222222222222 33333333333333333333 输出样例 55555555555555555555 解题思路 首先要解决的就是存储200 位整数的问题。显然，任何C/C++固有类型的变量都无法保存它。最直观的想法是可以用一个字符串来保存它。用str[200]来保存一个200 位的整数，让str[0]存放个位数，str[1]存放十位数，str[2]存放百位数…… 那么如何实现两个大整数相加呢？方法很简单，就是模拟小学生列竖式做加法，从个位开始逐位相加，超过或达到10 则进位。 大整数加法_实现过程 void add(char str_a[N], char str_b[N], char str_c[N]) { （1）把str_a,str_b的各位转变成对应整数，个位从下标0开始存入a, b （2）用变量len_a, len_b分别存储a,b中最高位的下标 （3）从低位到高位，对应位相加，左对齐存入c, 注意进位 （4）把c转化为字符串str_c } （1）把str_a,str_b的各位转变成对应整数，个位从下标0开始存入a[], b[] (1)将str_a复制给a; 将str_b复制给b; (2)将a,b逆转 (3)将a,b的每个元素都由字符变成对应整数 a[i]是字符，a[i]-48或a[i]=‘0’为其对应整数 （1）把str_a,str_b的各位转变成对应整数，个位从下标0开始存入a[], b[] (1) strcpy(a,str_a); strcpy(b,str_b); (2) strrev(a); strrev(b); (3)for(i=0;str[i]!=‘\0’;i++) a[i]=a[i]-’0’; for(i=0;str[i]!=‘\0’;i++) b[i]=b[i]-’0’; （2）用变量len_a, len_b分别存储a,b中最高位的下标 Len_a=strlen(str_a); Len_b=strlen(str_b); (3) 从低位到高位，对应位相加，左对齐存入c, 注意进位 4 5 6 7 9 3 8 1 6 5 4 3 2 1 0 6 5 4 3 2 1 0 carry 0 8 0 4 1 9 0 3 1 1 (3) 从低位到高位，对应位相加，左对齐存入c, 注意进位 (1)求len_a, len_b的最大值存入len (2)进位carry初始化为0 （3）for(i=0;ilen;i++)//从低位到高位 { 对应位相加加进位，存入sum 个位存入c[i]; 进位存入carry; } （4）若最高位有进位，向高位扩展一位，存入进位 （5）为下一步逆转做准备，在末尾加上字符串结束标志 (3) 从低位到高位，对应位相加，左对齐存入c, 注意进位 len=max(a_len,b_len); carry=0; for(i=0;ilen;i++) //对应位相加，注意进位 { sum=a[i]+b[i]+carry; c[i]=sum%10; carry=sum/10; } if(carry==1) //若最高位有进位 { c[i]=carry; i++; } len=i; for(i=0;ilen;i++) //各元素中的数字再转化为字符 c[i]=c[i]+0; c[len]=\0; //或c[len]=0; 字符串结束标志 } 课堂练习 一直某系统用基本类型int最大能表达的整数是32767， 现在要求1+2+…+n的和，n不大于32767，请编程实现 大整数减法 void sub(char str_a[N], char str_b[N], char str_c