高精度计算（附高精度加法程序）

高精度计算

有些数值计算要求精度高，希望计算的位数可达几十位甚至几百位，虽然计算机的计算精度也算较高了，但因受到硬件的限制，往往达不到实际问题所要求的精度，利用计算机处理高精度计算应注意以下问题：

<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]

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
#include
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;
}