高精度算法（+-×÷）

众所周知，加减乘除是日常生活经常需要用到的四种运算符。但是学过c的朋友们都知道，在c里面最大的整数数据类型unsigned long long最多也只能表达20多位数，也就是说如果出现了更大的数字，则无法使用常规的+=×÷来解决，这个时候就需要高精度算法来解决大数问题。

所谓的高精度算法，核心思想是一种返璞归真的做法，就是利用小学刚学数学的时候排竖式的方法，进行计算。所以具体步骤就是将两个大数存放在两个数组里面，每一个下标存放其中一个数字，之后再在两个数组之间进行一位对应一位的计算。

下面是高精度加法运算的函数

string add(string sa, string sb) {
	int lena = sa.size();
	int lenb = sb.size();
	int ma = max(lena, lenb);
	int a[1005] = { 0 }, b[1005] = { 0 }, c[1005] = { 0 };//a数组和b数组分别是两个数字，c数组用来存放结果
	for (int i = 0; i < lena; i++) {
		a[i] = sa[lena - i - 1] - '0';
	}//先将字符串倒过来存放，因为加法计算是从低位到高位
	for (int i = 0; i < lenb; i++) {
		b[i] = sb[lenb - i - 1] - '0';
	}//同上
	for (int i = 0; i < ma; i++) {
		c[i] += a[i] + b[i];
		if (c[i] > 9) {
			c[i + 1] += c[i] / 10;
			c[i] = c[i] % 10;
		}
	}
	if (c[ma] != 0)
		ma++;
	string s = "";
	for (int i = ma - 1; i >= 0; i--) {
		s += char(c[i] + '0');
	}
	return s;
}//高精度加法

高精度减法的运算方式和高精度加法类似，但是需要对借位思想有一点理解

string substract(string sa, string sb) {//前面一个数字大于等于后面一个数字
	int lena = sa.size();
	int lenb = sb.size();
	int lenc = 0;
	int a[1005] = { 0 }, b[1005] = { 0 }, c[1005] = { 0 };
	for (int i = 0; i < lena; i++) {
		a[i] = sa[lena - i - 1] - '0';
	}
	for (int i = 0; i < lenb; i++) {
		b[i] = sb[lenb - i - 1] - '0';
	}
	for (int i = 0; i < lena; i++) {
		c[i] += a[i] - b[i];
		if (c[i] < 0) {
			c[i]+= 10;
			a[i + 1]--;
		}
		if (c[i] != 0)
			lenc = i + 1;
	}
	
	string s = "";
	for (int i = lenc - 1; i >= 0; i--) {
		s += char(c[i] + '0');
	}
	return s;
}//高精度减法

高精度乘法的核心思想是，让其中一个乘数每一次出一个数（另外一个数一次从低位向高位和这个数字相乘（不要忘记进位）并且需要在每一步的时候进行该位的相加，如果把所有都算出来之后再相加可能会时间复杂度过高影响性能）比方说，1234×123，先用3乘以1、2、3、4，然后利用123和1234的第几位（代码中的i和j来衡量应该加在结果c数组的哪一位）

string multiplication(string sa, string sb) {
	int lena = sa.size();
	int lenb = sb.size();
	int lenc = 0;

	int a[1005] = { 0 }, b[5000] = { 0 }, c[5000] = { 0 };
	string s1 = "";
	for (int i = 0; i < lena; i++) {
		a[i] = sa[lena - i - 1] - '0';
	}
	for (int i = 0; i < lenb; i++) {
		b[i] = sb[lenb - i - 1] - '0';
	}
	int x;
	for (int j = 0; j < lenb; j++) {
		x = 0;
		for (int i = 0; i < lena; i++) {
			c[i + j] = a[i] * b[j] + c[i + j] + x;
			x = c[i + j] / 10;
			c[i + j] %= 10;
		}
		c[lena + j] = x;
	}
	lenc = lena + lenb;
	if (c[lenc - 1] == 0)
		lenc = lenc - 1;
	for (int i = lenc - 1; i >= 0; i--) {
		s1 += char(c[i] + '0');
	}
	return s1;
}//高精度乘法

最后是高精度除法，高精度除法需要从高位向地位计算，这是与其它三个算法不相同的地方，具体步骤是先将mod（余数）设置为0，然后用（最高位数+10✖mod）除以被除数（sa），将所得的余数赋值给mod，进行循环。

void divide(int sa, string sb) {//(sb/sa)
	
	int lenb = sb.size();
	int lenc = 0;
	int mod = 0;
	int b[5000] = { 0 }, c[5000] = { 0 };
	string s1 = "";

	for (int i = 0; i < lenb; i++) {
		b[i] = sb[lenb - i - 1] - '0';
	}
	for (int i = lenb - 1; i >= 0; i--) {
		c[i] = (mod * 10 + b[i])/sa;
		mod = (mod * 10 + b[i])%sa;
	}
	lenc = lenb;
	while (lenc > 1 && c[lenc - 1] == 0) {
		lenc--;
	}
	for (int i = 0; i < lenc; i++) {
		cout << c[lenc - i - 1];
	}
	cout << "   " << mod << endl;
}

以上四段代码就是四项比较基本的高精度算法，其中我写的高精度除法只能用高精度除数除以低精度被除数。希望对uus有帮助。