高精度计算（代码加解析，洛谷p1601，p1303）除法待更新

目录

高精度加法

高精度减法

高精度乘法

---

高精度加法

我们知道在c++语言中任何数据类型都有一定的表示范围。当两个被加数很大时，正常加法不能得到精确解。在小学，我们做加法都采用竖式方法。那么我们也只需要按照加法进位的方式就能得到最终解。

  8 5 6
+ 2 5 5
-------
1 1 1 1

加法进位：

c[i] = a[i] + b[i];
if(c[i] >= 10){
    c[i] %= 10;
    c[i+1]++;
}

按位相加求和：

    int len = s1.length();
    int carry = 0;
    for(int i = len-1; i >= 0; i--){
        int tmp = s1[i] - '0' + s2[i] - '0' + carry;
        carry = tmp/10;
        tmp %= 10;
        ret = char(tmp + '0') + ret;
    }
    if(carry != 0) ret = char(carry + '0') + ret;

 完整代码：

#include<iostream>
#include<string>
using namespace std;

int main()
{
    string s1,s2,ret;
    cin >> s1 >> s2;
    int len1 = s1.length();
    int len2 = s2.length();
    // 不论s1，s2谁短，高位补0补到一样长
    if(len1<len2)
    {
        for(int i=1;i<=len2-len1;i++)
           s1="0"+s1;
    }
    else
    {
        for(int i=1;i<=len1-len2;i++)
           s2="0"+s2;
    }
    int len = s1.length();  //得到两个字符串的长度
    int carry = 0;  // 记录进位
    // i从字符串最后一位向前走，也就是从两个数的个位向高位走
    for(int i = len-1; i >= 0; i--){
        // 相当于c[i] = a[i] + b[i],如果c[i]>=10, c[i]%=10, c[i+1]++；
        int tmp = s1[i] - '0' + s2[i] - '0' + carry; //字符转换成数字，要-'0'
        carry = tmp/10;
        tmp %= 10;
        ret = char(tmp + '0') + ret; //把每个新位放到原字符串前面
    }
     //如果carry位不等于0，则表明两个数的最高位相加还有进位
    if(carry != 0) ret = char(carry + '0') + ret;  
    cout << ret;
}

高精度减法

类似加法，也可以用竖式求解。需要注意的是，被减数必须比减数大，同时需要处理借位。

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int main()
{
    int a[256],b[256],c[256],lena,lenb,lenc,i;
    char n[256],n1[256],n2[256];
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(c,0,sizeof(c));
    printf("Input minuend:");	gets(n1);
    printf("Input subtrahend:");	gets(n2);
    if(strlen(n1) < strlen(n2) || (strlen(n1) == strlen(n2) && strcmp(n1,n2) < 0)){
    	strcpy(n,n1);
    	strcpy(n1,n2);
    	strcpy(n2,n);
    	cout << "-"; 	//因为交换了减数和被减数，结果为负数 
	}
	lena = strlen(n1);
	lenb = strlen(n2);
	for(i = 0; i <= lena-1; i++) a[lena-i] = int(n1[i] - '0');	//被减数放入a数组 
	for(i = 0; i <= lenb-1; i++) b[lenb-i] = int(n2[i] - '0');	//减数放入b数组
	i= 1;
	while(i <= lena || i <= lenb){
		if(a[i] < b[i]){
			a[i]+=10;	//高位借1 
			a[i+1]--;
		}
		c[i] = a[i] - b[i];	//对应位相减
		i++; 
	} 
	lenc = i;
	while(c[lenc] == 0 && lenc > 1) lenc--;	//最高位0不输出
	for(i = lenc; i >= 1; i--) cout << c[i];
	cout << endl;
	return 0; 
}

高精度乘法

类似加法，可以用竖式求乘法。在做乘法运算时，同样也有进位。同时对每一位进行乘法运算时，必须进行错位相加。

     8  5  6
*       2  5 
----------------
  4  2  8  0
1 7  1  2
----------------
2 1  4  0  0

#include<iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int main()
{
    char a1[101],b1[101];
    int a[101],b[101],c[10001],lena,lenb,lenc,i,j,x;
    memset(a,0,sizeof(a));
    memset(b,0,sizeof(b));
    memset(c,0,sizeof(c));
    scanf("%s",a1);
    scanf("%s",b1);
    lena = strlen(a1); lenb = strlen(b1);
    for(i = 0; i <= lena-1; i++) a[lena-i] = int(a1[i] - 48);	
	for(i = 0; i <= lenb-1; i++) b[lenb-i] = int(b1[i] - 48);	
	for(i = 1; i <= lena; i++){
		x = 0;
		for(j = 1; j <= lenb; j++){
			c[i+j-1] += a[i]*b[j] + x;	//当前乘机+上次乘机进位+原数 
			x = c[i+j-1] / 10;
			c[i+j-1] %= 10;
		}
		c[i+lenb] = x;	//进位 
	} 
	lenc = lena + lenb;
	while(c[lenc] == 0 && lenc > 1)	lenc--;	//删除前导0 
	for(i = lenc; i > 0; i--) cout << c[i];
	cout << endl; 
	
	return 0; 
}

还有一种思路，把每一位的乘积和加起来先不做处理。最后再去处理计算进位。

        7    8    9
*            2    3
----------------------
       21   24    27
  14   16   18
----------------------
  14   37   42    27

 1 (1) 8 (4) 1 (4) 4 (2) 7