pascal语言 高精度计算

《pascal语言 高精度计算》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第一章高精度计算利用计算机进行数值计算，有时会遇到这样的问题：有些计算要求精度高，希望计算的数的位数可达几十位甚至几百位，虽然计算机的计算精度也算较高了，但因受到硬件的限制，往往达不到实际问题所要求的精度。我们可以利用程序设计的方法去实现这样的高精度计算。介绍常用的几种高精度计算的方法。高精度计算中需要处理好以下几个问题：(1)数据的接收方法和存贮方法数据的接收和存贮：当输入的数很长时，可采用字符串方式输入，这样可输入数字很长的数，利用字符串函数和操作运算，将每一位数取出，存入数组中。另一种方法是直接用循环加数组方法输入数据。(2)高精度数位数的确定位数的确定：接收时往往

2、是用字符串的，所以它的位数就等于字符串的长度。(3)进位,借位处理加法进位：C[i]:=A[i]+B[i];ifC[i]>=10thenbeginC[i]:=C[i]mod10;C[i+1]:=C[i+1]+1end;减法借位：ifa[i]

3、例1】输入两个正整数，求它们的和。【分析】输入两个数到两个变量中，然后用赋值语句求它们的和，输出。但是，我们知道，在pascal语言中任何数据类型都有一定的表示范围。而当两个被加数很大时，上述算法显然不能求出精确解，因此我们需要寻求另外一种方法。在读小学时，我们做加法都采用竖式方法，如图1。这样，我们方便写出两个整数相加的算法。856+2551111图1A3A2A1+B3B2B1C4C3C2C1图2如果我们用数组A、B分别存储加数和被加数，用数组C存储结果。则上例有A[1]=6，A[2]=5，A[3]=8，B[1]=5，B[2]=5，B[3]=2，C[4]=1，C[3]=

4、1，C[2]=1，C[1]=1，两数相加如图2所示。因此，算法描述如下：procedureadd(a,b;varc);//a,b,c都为数组，a，b，c分别存储被加数、加数、结果vari,x:integer;begini:=1；x:=0;//x是进位while(i<=a数组长度)or(i<=b数组的长度)dobeginc[i]:=a[i]+b[i]+x;//第i位相加并加上次的进位x:=c[i]div10;//向高位进位c[i]:=c[i]mod10;//存储第i位的值i:=i+1//位置指针变量end通常，读入的两个整数用可用字符串来存储，程序设计如下：programe

5、xam1;constmax=200;vara,b,c:array[1..max]of0..9;n:string;lena,lenb,lenc,i,x:integer;beginwrite('Inputaugend:');readln(n);//输入加数lena:=length(n);//加数放入a数组fori:=1tolenadoa[lena-i+1]:=ord(n[i])-ord('0');write('Inputaddend:');readln(n);//输入被加数lenb:=length(n);//被加数放入b数组fori:=1tolenbdob[lenb-i+1]

6、:=ord(n[i])-ord('0');i:=1;x:=0;while(i<=lena)or(i<=lenb)dobeginc[i]:=a[i]+b[i]+x;//两数相加，然后加前次进位x:=c[i]div10;//向高位进位c[i]:=c[i]mod10;//保存第i位的值i:=i+1end;ifx>0then//处理最高进位beginlenc:=i;c[i]:=x;endelselenc:=i-1;fori:=lencdownto1dowrite(c[i]);//输出结果writelnend.【例2】高精度减法。输入两个正整数，求它们的差。【算法分析】类似加法，可

7、以用竖式求减法。在做减法运算时，需要注意的是：被减数必须比减数大，同时需要处理借位。高精度减法的参考程序：programexam2;constmax=200;vara,b,c:array[1..max]of0..9;n,n1,n2:string;lena,lenb,lenc,i,x:integer;beginwrite('Inputminuend:');readln(n1);//输入被减数write('Inputsubtrahend:');readln(n2);//输入减数if(length(n1)