高精度四则运算

b对一个整数来说，如果超过了当前数据类型所能表示的范围，那么就只能用数组去模拟输出，下面介绍高精度的储存、高精度比较大小和四则运算(只适用于正整数)。

一、高精度的储存

大整数的储存使用一组数组来储存的，数组的低位表示数的低位，比如整数123456789储存在数组中就是d[0] = 9, d[1] = 8 …d[8] = 1, 但是我们用字符串输入时，字符串的低位存的是整数的高位，我们在储存时需要反转一下，这是需要一个len来记录整数的位数，并和数组组成一个结构体。

struct bign{
	int d[1000];
	int len;
	bign(){
		memset(d, 0, sizeof(d));
		len = 0;
	}
};

// 将输入的字符串逆序存到数组中
bign reserve(char str[]){
	bign a;
	a.len = strlen(str);
	for(int i = 0; i < a.len; i++){
		a.d[i] = str[a.len - i -1] - '0';
	}
	return a;
}

二、高精度的比较

两个数比较大小，首先看整数的位数，位数大的整数大。如果位数一样，就从高位比到低位，直到出现不相等的情况，就可以得出两个数的大小，否则两个数相等。

//a > b返回1 ,a = b 返回0, a < b 返回-1 
int compare(bign a, bign b){
	if(a.len > b.len) return 1;
	else if(a.len < b.len) return -1;
	else{
		for(int i = a.len-1; i >= 0; i--){
			if(a.d[i] > b.d[i]) return 1;
			else if(a.d[i] < b.d[i]) return -1;
		}
	}
	return 0;
}

三、高精度加法

高精度加法其实就是模拟小学的加法计算，下面先模拟一下123+789的过程：$$\begin{array}{rccc}& 1& 2& 3\\ +& 7& 8& 9\\ & 9& 1& 2\end{array}$$
1.计算3+9，取个位2为结果，并产生高位的进位1。
2.计算2+8，加上来自低位的进位并取个位1为结果，产生向高位的进位1.
3.计算1+7，加上来自低位的进位产生结果9，无向高位进位，计算完毕。

代码实现如下：

bign add(bign a, bign b){
	bign c;
	int carry = 0;
	for(int i = 0; i < a.len || i < b.len; i++){
		int t = a.d[i] + b.d[i] + carry;//计算a, b和进位的和
		c.d[c.len++] = t%10; //结果为个位
		carry = t/10; //进位为十位
	}
	//计算到最后如果有进位，高一位的结果为该进位值
	if(carry != 0){
		c.d[c.len++] = carry;
	}
	return c;
}

四、高精度减法

同加法一样，也是模拟小学减法，下面模拟721-189的过程：
$$\begin{array}{rccc}& 7& 2& 1\\ -& 1& 8& 9\\ & 5& 3& 2\end{array}$$
1.计算1-9，不够减，向高位借1，本位加10，然后计算11-9，结果为2。
2.计算1-8，不够减，向高位借1，本位加10，然后计算11-8，结果为3。
3.计算6-1，结果为5。

代码实现如下：

bign sub(bign a, bign b){
	bign c;
	for(int i = 0; i < a.len || i < b.len; i++){
		if(a.d[i] < b.d[i]){
			a.d[i+1]--;
			a.d[i] += 10;
		}
		c.d[c.len++] = a.d[i] - b.d[i];
	}
	//去掉高位0
	while(c.len > 1 && c.d[c.len-1] == 0){
		c.len--;
	}
	return c; 
}

注：本代码只能计算大数减小数，如果a < b，则需要先交换两数，然后输出加上负号。

五、高精度乘法（高精度$\times$低精度）

同上，也是模拟小学乘法，不同的是高精度乘法实现的过程需要把低精度看成一个整体，而不是一位一位去乘，下面模拟$12\times 89$(12为低精度)的过程：
$$\begin{array}{rcccc}& & & 1& 2\\ \times & & & 8& 9\\ & & 1& 0& 8\\ & & 9& 6& \\ & 1& 0& 6& 8\end{array}$$
1.计算$9\times 12$, 取个位数8为结果，产生高位的进位10。
2.计算$8\times 12$并加上进位10, 取个位数6为结果，产生高位的进位10，没有乘数可乘，直接取进位为最高位，计算完毕。

代码实现如下：

bign mul(bign a, int b){
	bign c;
	int carry = 0;
	for(int i = 0; i < a.len; i++){
		int t = a.d[i] * b + carry;
		c.d[c.len++] = t % 10;
		carry = t / 10;
	}
	//乘法的进位可能不止一位，所以要用while
	while(carry != 0){
		c.d[c.len++] = carry % 10;
		carry /= 10;
	}
	return c;
}

六、高精度除法

也是模拟小学除法，同乘法一样，需要把除数当成一个整体，下面模拟$144÷12$：

1.计算$1÷12$，不够除，结果为0，余数为1。
2.余数乘10，计算$14÷12$，结果为1，余数为2。
3.余数乘10，计算$24÷12$，结果为2，余数为0。

bign div(bign a, int b){
	bign c;
	int r = 0;
	c.len = a.len; //商的位数与除数一一对应，不够除为0
	for(int i = a.len -1; i >=0; i--){
		int t = a.d[i] + r*10;
		c.d[i] = t / b;
		r = t % b;
	}
	while(c.len > 1 && c.d[c.len-1] == 0){
		c.len--;
	}
	c.d[c.len] = r; //利用多余数组储存余数
	return c;
}

以上就是高精度的全部内容，有问题欢迎私信我！

题目推荐：
P1009_阶乘之和(NOIP1998 普及组)