高精度算法总结

加法：其算法原理其实和列竖式进行加法一样

用三个数组分别存两个加数以及结果，比如 123+987，a[]={3,2,1}，b[]={7,8,9}（因为加法是从个位数开始加，所以顺序要反转），

c[1]=a[1]+b[1]+c[1]=3+7+0=10;

c[1+1]=c[2]=c[1]/10=1;   c[1]=c[1]%10=0;

c[2]=a[2]+b[2]+c[2]=2+8+1=11;

c[2+1]=c[3]=c[2]/10=1;   c[2]=c[2]%10=1;

c[3]=a[3]+b[3]+c[3]=1+9+1=11;

c[3+1]=c[4]=c[3]/10=1;   c[3]=c[3]%10=1;

所以结果就为数组c倒序输出1110

for (int i = 0; i < ans; i++) {
        c[i] += a[i] + b[i];
        c[i + 1] = c[i] / 10;
        c[i] %= 10;
    }

 其中ans为两个数组中的最大长度，最后输出结果的时候要判断一下结果的位数是与最大长度的加数一样还是要比其大一位。

全部代码：

#include<iostream>
#include<string>
using namespace std;

const int N = 510;
int a[N], b[N], c[N];
int main()
{
    string str1;
    string str2;
    cin >> str1;
    cin >> str2;
    for (int i = 0; i < str1.size(); i++)
        a[str1.size() - 1 - i] = str1[i] - '0';
    for (int i = 0; i < str2.size(); i++)
        b[str2.size() - 1 - i] = str2[i] - '0';
    int ans = max(str1.size(), str2.size());
    for (int i = 0; i < ans; i++) {
        c[i] += a[i] + b[i];
        c[i + 1] = c[i] / 10;
        c[i] %= 10;
    }
    ans += 1;
    if (c[ans - 1] == 0 && ans > 1)
        ans -= 1;
    for (int i = 0; i < ans; i++)
        cout << c[ans - 1 - i];
    return 0;
}

 乘法：也是与竖式乘法一样

 将每一列的数加起来，进位和加法一样。

每一列的数怎么来的呢：比如右数第三列，上面乘数是456，下面是234，从上面的左数第一个，下面的右数第一个开始相乘，4*4=16，5*3=15，6*2=12。

代码：

for (int i = 0; i < num; i++)
	{
		for (int j = 0; j <= i; j++)c[i] += x[j] * y[i - j];
		if (c[i] >= 10)
		{
			c[i + 1] += c[i] / 10;
			c[i] %= 10;
		}
	}

数组x，y为两个乘数，数组c存储结果，if语句中的为进位，与加法一样。num为结果长度=乘数的长度相加再减1。

全部代码:

#include<iostream>
#include<string>
using namespace std;

const int N = 3000;
int x[N], y[N];
int c[N];
string a, b;
int num;
void mults()
{
	num = a.size() + b.size() - 1;
	for (int i = 0; i < num; i++)
	{
		for (int j = 0; j <= i; j++)c[i] += x[j] * y[i - j];
		if (c[i] >= 10)
		{
			c[i + 1] += c[i] / 10;
			c[i] %= 10;
		}
	}
	//去掉多余的前导零 
	for (int i = num; i > 0; i--)
	{
		if (c[i] == 0)num--;
		else break;
	}
}
int main()
{
	cin >> a >> b;
	memset(x, 0, sizeof(x));
	memset(y, 0, sizeof(y));
	for (int i = a.size() - 1; i >= 0; i--)	x[a.size() - 1 - i] = a[i] - '0';

	for (int i = b.size() - 1; i >= 0; i--) y[b.size() - 1 - i] = b[i] - '0';

	mults();

	for (int i = num; i >= 0; i--)	 printf("%d", c[i]);
	return 0;

}