算法 - 高精度

所有代码都是手敲的，有bug欢迎提出。虽然没人看
本篇文章包括高精度加、减、乘、除。高精度加减的两个数都是高精度，高精度乘除是高精度 ×或÷ 低精度。

一、预备知识——C++与STL入门

（一）string

做标记的是这个本文会用到的：

（二）vector

使用vector是因为有size()函数，遍历很方便，而普通的数组获得其长度较麻烦。

（三）引用符&

在后续的代码中，引用的使用目的是提高代码效率。
这里先记住一点：当大型对象被传递给函数时，使用引用参数可使参数传递效率得到提高。关于引用和指针，请参考我的另一篇文章：

二、技巧

（一）如何存储

以加法为例，如果数组存储的是正序，那么最高位出现进位时，所有元素均需要向后移一位。如果数组存储的逆序，只需要在末尾补上即可。比如123，存储在数组中的顺序是321。输入输出的对应代码如下

string a, b;
vector<int> A, B;

cin>>a>>b;
for(int i = a.size() - 1; i >= 0; i --) A.push_back(a[i] - '0');
for(int i = b.size() - 1; i >= 0; i --) B.push_back(b[i] - '0');

auto C = add(A, B);

for(int i = C.size() - 1; i >= 0; i --) printf("%d", C[i]);

（二）去前导零

在减法和除法中，我们经常会遇到前面有前导0的结果，附上解决代码：

while(C.size() > 1 && C.back() == 0) C.pop_back();

三、思路和板子

之所以前面铺垫那么多，是想让这一部分的重点更明确，那就是明确每一种运算的本质——每一位是怎么来的。
一个共有的特点是：根据人的计算方式，从低位到高位计算，由于数组按逆序存储，所以也应从低位到高位处理。
其实一圈下来减法是最难的，突然有了计组的补码减法的感觉

（一）高精度加法

设进位为t，则：

• 该位 = ( A[i] + B[i] + t ) % 10
• 下一位的进位 = ( A[i] + B[i] + t ) / 10

//C = A + B, A >= 0, B >= 0
vector<int> add(vector<int> &A, vector<int> &B){
	vector<int> C;
	int t = 0; //进位
	for(int i = 0; i < A.size() || i < B.size(); i ++){
		if(i < A.size()) t += A[i];
		if(i < B.size()) t += B[i];
		C.push_back(t % 10);
		t /= 10;
	}
	//判断是否有进位
	if(t) C.push_back(1);
	
	return C;	
}

（二）高精度减法

设借位为t，则：

• 该位
  (1) 若A[i] - B[i] - t >= 0，则该位 = A[i] - B[i] - t
  (2) 若A[i] - B[i] - t < 0，则该位 = ( A[i] - B[i] - t + 10 ) % 10
  两者都可以用 ( A[i] - B[i] - t + 10 ) % 10 的形式表示
• 下一位的借位
  (1) 若A[i] - B[i] - t >= 0， t = 0
  (2) 若A[i] - B[i] - t < 0, t = 1

//C = A - B, A >= B, A >= 0, B >= 0
vector<int> sub(vector<int> &A, vector<int> &B){
	vector<int> C;
	int t = 0; //借位
	for(int i = 0; i < A.size(); i ++){
		t = A[i] - t;
		if(i < B.size()) t -= B[i];
		C.push_back((t + 10) % 10); 
		if(t < 0) t = 1;
		else t = 0;		
	}
	//已经保证了A >= B，因此不用再判断t
	//去前导零
	while(C.size() > 1 && C.back() == 0) C.pop_back();

	return C;
}

这个板子需要保证A >= B，因此需要一个cmp函数进行判断

bool cmp(vector<int> A, vector<int> B){
	if(A.size() != B.size())
		return A.size() > B.size();
	for(int i = A.size() - 1; i >= 0; i --){
		if(A[i] != B[i])
			return A[i] > B[i];
	}
	return true; //相等
}

//main
if(cmp(A, B)) auto C = sub(A, B);
else{
	auto C = sub(B, A);
	printf("-");
}
/* 
输出C
*/

（三）高精度乘法

高精度乘法与人的手算过程有区别。举个例子：A = 123，B = 12，每次相乘是乘B，而不是一位一位。
C0 = ( 3 * 12 + 0 ) % 10 = 6
t1 = ( 3 * 12 + 0 ) / 10 = 3

C1 = ( 2 * 12 + t1 ) % 10 = 7
t2 = ( 2 * 12 + t1 ) / 10 = 2

C2 = ( 1 * 12 + t2 ) % 10 = 4
t3 = ( 1 * 12 + t2 ) / 10 = 1

有进位：
C3 = t3= 1
结果：1476
抽象概括为：

• 当前位 = ( A[i] * B + t ) % 10
• 下一位的进位 = ( A[i] * B + t ) / 10

板子如下：

//C = A * b, A >= 0, b >= 0
vector<int> mul(vector<int> &A, int b){
	vector<int> C;
	int t = 0;
	for(int i = 0; i < A.size(); i ++){
		t += A[i] * b;
		C.push_back(t % 10);
		t /= 10;
	}
	//判断是否有进位
	if(t) C.push_back(t);
	//去前导零
	while(C.size() > 1 && C.back() == 0) C.pop_back();
	
	return C;
}

高精度除法

除法的过程遵循人类的计算方式。但注意除法是从高位开始计算。
没有改变存储方式是为了和前三种统一。
r表示余数

• 当前位 = ( A[i] + r * 10 ) / b
• 余数 = ( A[i] + r * 10 ) % b

// A/B = C ... r, A >= 0, b > 0
vector<int> div(vector<int> &A, int b, int &r){
	vector<int> C;
	r = 0;
	for(i = A.size() - 1; i >= 0; i --){
		r = r * 10 + A[i];
		C.push_back(r / b);
		r %= b;
	}
	// <algorithm>
	reverse(C.begin(), C.end());
	// 去前导零
	while(C.size() > 1 && C.back() == 0) C.pop_back();
	return C;
}

附上完整代码。

//加法
#include<iostream>
#include<vector>

using namespace std;

vector<int> add(vector<int>&A, vector<int>&B){
    vector<int> C;
    int t = 0; //进位
    for(int i = 0; i < A.size() || i < B.size(); i++){
        if(i < A.size()) t += A[i];
        if(i < B.size()) t += B[i];
        C.push_back(t % 10);
        t /= 10;
    }
    if(t) C.push_back(1);
    return C;
}

int main(){
    string a, b;
    vector<int> A, B;

    cin>>a>>b; //string类型不能用scanf和printf
    for(int i = a.size() - 1; i >= 0; i --) A.push_back(a[i] - '0');
    for(int i = b.size() - 1; i >= 0; i --) B.push_back(b[i] - '0');

    auto C = add(A, B);

    for(int i = C.size() - 1; i >= 0; i --) printf("%d", C[i]);

    return 0;
}

//减法
#include<iostream>
#include<vector>

using namespace std;

bool cmp(vector<int> &A, vector<int> &B){
    //高精度比较
    if(A.size() != B.size()) 
        return A.size() > B.size();
    for(int i = A.size() - 1; i >= 0; i --)
        if(A[i] != B[i])
            return A[i] > B[i];
    return true;
}

vector<int> sub(vector<int> &A, vector<int> &B){
    vector<int> C;
    int t = 0; //借位
    for(int i = 0; i < A.size(); i ++){
        t = A[i] - t;
        if(i < B.size()) t -= B[i];
        C.push_back((t + 10) % 10);
        if(t < 0) t = 1;
        else t = 0;
    }

    //去前导零
    while(C.size() > 1 && C.back() == 0) C.pop_back();
    return C;
}

int main(){
    string a, b;
    vector<int> A, B;

    cin >> a >> b;
    for(int i = a.size() - 1; i >= 0; i --) A.push_back(a[i] - '0');
    for(int i = b.size() - 1; i >= 0; i --) B.push_back(b[i] - '0');

    if(cmp(A, B)){
        auto C = sub(A, B);
        for(int i = C.size() - 1; i >= 0; i --) printf("%d", C[i]);
    }
    else{
        auto C = sub(B, A);
        printf("-");
        for(int i = C.size() - 1; i >= 0; i --) printf("%d", C[i]);
    }


    return 0;
}

//乘法
#include<iostream>
#include<vector>

using namespace std;

vector<int> mul(vector<int> &A, int b){
    vector<int> C;
    int t = 0; //进位
    for(int i = 0; i < A.size(); i ++){
        t = A[i] * b + t;
        C.push_back(t % 10);
        t /= 10;
    }
    if(t) C.push_back(t);
    
    //去前导零
    while(C.size() > 1 && C.back() == 0) C.pop_back();
    return C;
}

int main(){
    string a;
    int b;
    vector<int> A;
    
    cin >> a >> b;
    for(int i = a.size() - 1; i >= 0; i --) A.push_back(a[i] - '0');
    
    auto C = mul(A, b);
    
    for(int i = C.size() - 1; i >= 0; i --) printf("%d", C[i]);
    
    return 0;
    
}

//除法
#include<iostream>
#include<vector>
#include<algorithm>

using namespace std;

vector<int> div(vector<int> &A, int b, int &r){
    vector<int> C;
    r = 0; //余数
    
    for(int i = A.size() - 1; i >= 0; i--){
        r = r * 10 + A[i];
        C.push_back(r / b);
        r %= b;
    }
    
    reverse(C.begin(), C.end());
    while(C.size() > 1 && C.back() == 0) C.pop_back();
    
    return C;
}

int main(){
    string a;
    int b, r;
    cin >> a >> b;
    
    vector<int> A;
    for(int i = a.size() - 1; i >= 0; i --) A.push_back(a[i] - '0'); 
    
    auto C = div(A, b, r);
    
    for(int i = C.size() - 1; i >= 0; i --) printf("%d", C[i]);
    
    cout << endl << r;
    
    return 0;
}