高精度计算
有些数值计算要求精度高,希望计算的位数可达几十位甚至几百位,虽然计算机的计算精度也算较高了,但因受到硬件的限制,往往达不到实际问题所要求的精度,利用计算机处理高精度计算应注意以下问题:
<1>数据的接收与存贮
①利用字符串读入数据
a[i]=s[s.length()-i]-'0';
②直接用循环加数组方法输入数据
注意:倒序存储
<2>确定位数
位数即为字符串长度
<3>进位,借位处理
①加法进位
c[i]=a[i]+b[i];
if(c[i]>=10){c[i]%=10;c[i+1]++;}
②减法借位
if(a[i]<b[i]){a[i+1]--;a[i]+=10;}
c[i]=a[i]-b[i];
③乘法进位
c[i+j-1]=a[i]*b[i]+x+c[i+j-1];
x=c[i+j-1]/10;
c[i+j-1]%=10;
<4>商和余数的求法:视被除数和除数的位数情况进行处理
下附高精度加法C++程序
#include<cstdio> #include<cstring> using namespace std; int main() { char a1[100],b1[100]; int a[100],b[100],c[100],lena,lenb,lenc,i,x; memset(a,0,sizeof(a)); memset(b,0,sizeof(b)); memset(c,0,sizeof(c)); printf("加数1:"); gets(a1); printf("加数2:"); gets(b1); lena=strlen(a1); lenb=strlen(b1); for (i=0;i<=lena-1;i++) a[lena-i]=a1[i]-48; for (i=0;i<=lenb-1;i++) b[lenb-i]=b1[i]-48; lenc=1; x=0; while(lenc<=lena||lenc<=lenb) { c[lenc]=a[lenc]+b[lenc]+x; x=c[lenc]/10; c[lenc]%=10; lenc++; } c[lenc]=x; if (c[lenc]==0) lenc--; printf("和: "); for (i=lenc;i>=1;i--) { printf("%d",c[i]); } printf("\n"); return 0; }