高精度运算

前言

高精度算法用于解决long long也解决不了的计算。

高精度的存储是把每一位单独存储，且是倒序存储，数组 c[0] 是这个数的个位，c[1] 是这个数的十位，以此类推。

一、高精度乘法

1.1 高精度 x 低精度

（1）首先，令进位为 tmp = 0；
（2）让c 每位的数值乘以低精度数d，并存储在当前位上。
（3）从最低位开始，当前位的数值加上进位数值tmp赋值给进位tmp，进位tmp模10结果就是当前位的结果。
（4）进位tmp对10取余结果就是新的进位tmp。
（5）重复（3）（4）步骤，直到完成位置从0到 lc-1 的计算。
（6）剩余的进位tmp不为0说明还能进位，用while循环即可完成存储。

int c[1000], lc;	//c用来存储高精度，lc是c的长度
void muti(int d)
{
	for (int i = 0; i < lc; i++)
		c[i] *= d;
	int tmp = 0;	//进位
	for (int i = 0; i < lc; i++)	//位置从0到lc-1的计算
	{
		tmp += c[i];
		c[i] = tmp % 10;
		tmp /= 10;
	}
	while (tmp != 0)	//超出原来长度的数（进位），继续存储
	{
		c[lc] = tmp % 10;
		tmp /= 10;
		lc++;
	}
    
	for (int i = lc - 1; i >= 0; i--)	//输出
		cout << c[i];
}

1.2 高精度 x 高精度

对于高精度 a[100] = {2, 1, 3}，对应十进制个位是2，十位是1，百位是3；
对于高精度 b[100] = {5, 2, 4}，对应十进制个位是5，十位是2，百位是4；
用c[100]存储。

通过模拟乘法可以知道，下标和相等的两数的乘积是存储在同一列，c[i + j] = a[i] * b[j]。

所以最终结果是每一列相加，就是i+j这一列所有数相加。所以可以用c[i+j] += a[i] * b[j]。

int a[1000], b[1000], c[1000];
int la, lb, lc;
void muti()
{
	for (int i = 0; i < la; i++)
		for (int j = 0; j < lb; j++)
			c[i + j] += a[i] * b[j];
	for (int i = 0; i < la + lb; i++)
		if (c[i] >= 10) {
			c[i + 1] += c[i] / 10;
			c[i] %= 10;
		}
	lc = la + lb;
	while (lc > 1 && c[lc - 1] == 0)	lc--;	//去除前面多余的0

	for (int i = lc-1; i >= 0; i--)	//输出
		cout << c[i];
}

二、高精度除法

2.1 高精度 / 低精度

利用竖式计算，c = a / d。

（1）首先，令余数为 x = 0；
（2）让余数乘以10，和高位的第一个位置的数相加，然后除以d，得到此位置的商;
（3）求出此时的余数x；
（4）接着往高位的下一位走，重复（2）（3）步骤，直到最后一位；

int a[1000], c[1000], la, lc;	//a是原数，c是除后的结果
void div(int d)
{
	int x = 0;	//余数
	for (int i = la - 1; i >= 0; i--)
	{
		c[i] = (x * 10 + a[i]) / d;
		x = (x * 10 + a[i]) % d;
	}
    lc = la;
	while (lc > 1 && c[lc - 1] == 0)	lc--;	//去除前面多余的0
    
    for (int i = lc-1; i >= 0; i--)	//输出
		cout << c[i];
}