【C++学习笔记】 模拟竖除法

最近写了一道有关除法的题目UVa202，该题需要用到模拟竖除法的思想。之前没了解过，于是花了一点时间学习了这个方法。

所谓模拟竖除法，就是模拟我们小学做除法时用到的除法竖式。由于计算机的浮点数的精度限制，我们可以用这种方式求出任意精度的小数。

第一种情况是除数是一个低精度浮点数。这种情况比较简单，只需要将被除数的每一位存储到一个int类型的数组中。接着用竖除法一步步计算，将每一步的计算结果存到一个新的int数组中，并用一个int类型变量保存商。

举例：

数组num[100]保存被除数的每一位，为了方便起见，可以用num[0]保存被除数的长度。

int a为被除数，int b为商。

假设输入被除数:112352；除数：9；

那么num[0] == 6,num[1] == 1,num[2] == 1,...,num[6] = 2;

第一次除法：1/9不够除。

第二次除法：往后借一位，11/9，此时b=2，结果11/9 = 1保存到新的数组中。

.........

代码：

#include<iostream>
using namespace std;
int main() {
	int num[256],result[256], a, b = 0, index = 0; 
	string snum;
	cin >> snum >> a;
	num[0] = snum.length();

	for (int i = 1; i <= num[0]; i++) {
		num[i] = snum[i - 1] - '0';
	}

	for (int i = 1; i <= num[0]; i++) {
		b = b * 10 + num[i];
		if (b >= a) {
			result[i] = b / a;
			b = b % a;
		}
		else {
			result[i] = 0;
		}
	}

	int flag = false;
	for (int i = 1; i <= num[0]; i++) {
		if (result[i] != 0 && flag == false) {
			flag = true;
		}
		if(flag)
			cout << result[i];
	}
	return 0;
}

该代码只有正整数部分相除，小数部分思想相同。

第二中情况是高精度被除数除以一个高精度除数，这种情况下除数也需要用一个数组保存，因此会复杂很多。在每次除以除数时，需要先比较这一步中的被除数和除数的大小关系，如果被除数大于除数，那么每次减去除数，知道被除数小于除数。同时，在高精度除以高精度中，被除数会直接用来减去除数，因此我们需要一个辅助数组来为除数补‘0’，使得被除数的位数和除数的位数相同。

#include<iostream>
using namespace std;

int cmp_arr(int a[],int b[]) {
	for (int i = 1; i <= b[0]; i++) {
		if (a[i] > b[i]) return 1;
		else if(a[i] < b[i])return -1;
	}
	return 0;
}

void copy_arr(int b[], int tmpb[], int dest) {
	for (int i = 1; i <= b[0]; i++) {
		tmpb[i + dest] = b[i];
	}
	tmpb[0] = b[0] + dest;
}

void sub_arr(int a[], int b[]) {
	for (int i = b[0]; i >= 1; i--) {
		if (a[i] < b[i]) {
			a[i - 1] --;
			a[i] += 10;
		}
		a[i] -= b[i];
	}
}

void print_arr(int c[]) {
	bool flag = false;
	for (int i = 1; i <= c[0]; i++) {
		if (!flag && c[i] == 0) continue;
		flag = true;
		cout << c[i];
	}
	if (!flag) cout << '0';
}

int main() {
	int a[256], b[256], c[256];
	string stra, strb;
	cin >> stra >> strb;
	a[0] = stra.length();
	b[0] = strb.length();
	for (int i = 1; i <= a[0]; i++)
		a[i] = stra[i - 1] - '0';
	for (int i = 1; i <= b[0]; i++)
		b[i] = strb[i - 1] - '0';

	memset(c, 0, sizeof(c));

	if (cmp_arr(a, b) >= 0) c[0] = a[0] - b[0] + 1;
	else c[0] = a[0] - b[0];

	int tmpb[256] = { 0 };
	for (int i = 1; i <= c[0]; i++) {
		memset(tmpb, 0, sizeof(tmpb));
		copy_arr(b, tmpb, i - 1);
		while (cmp_arr(a, tmpb) >= 1) {
			sub_arr(a, tmpb);
			c[i]++;
		}
	}

	print_arr(c);
	return 0;
}