c/c++ 大整数乘法

运算思路:

其做法和小学学的有一点不一样。以 147 X 35 为例,这里把 147 视为
高精度 bign(自定义的结构下面代码会提到) 类型,而 35 视为 int 类型, 并且在下面的过程中, 始终将 35 作为一个整体看待。

  1. 7x35=245, 取个位数 5 作为该位结果, 高位部分 24 作为进位。
  2. 4 X 35 = 140, 加上进位 24, 得 164, 取个位数 4 为该位结果,高位部分 16 作为进位。
  3. 1 X 35 =35, 加上进位 16, 得 51, 取个位数 1 为该位结果,高位部分 5 作为进位。
  4. 没的乘了, 此时进位还不为0, 就把进位5直接作为结果的高位。
  5. 对某一步来说是这么一个步骤:取 bign 的某位与 int 型整体相乘, 再与进位相加, 所得结果的个位数作为该位结果, 高位部分作为新的进位。

代码实现和解析如下:

#include<bits/stdc++.h>
using namespace std;
struct bign {
	int d[1000];
	int len;
	bign() {
		memset(d,0,sizeof(d));
		len=0;
	}
};
bign change(char str[]) {//将整数转换为bign 
	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;
}
int compare(bign a,bign b) {//比较一下a和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;
	}
}
bign sub(bign a,bign b){//高精度a-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;//当前位加10 
		}
		c.d[c.len++]=a.d[i]-b.d[i];//减法结果为当前位结果 
	}
	while(c.len-1>=1&&c.d[c.len-1]==0){
		c.len--;//去除高位的0,同时至少保留一位最低位 
	}
	return c;
}
void print(bign a){//输出 
	for(int i=a.len-1;i>=0;i--) printf("%d",a.d[i]);
}
int main(){
	char str1[1000],str2[1000];
	scanf("%s%s",str1,str2);
	bign a=change(str1);
	bign b=change(str2);
	print(sub(a,b));
	return 0;
}

其中的bign储存方法可以连接:大整数储存和加法

特别注意:

如果a和b中存在负数, 需要先记录下其负号, 然后取它们的绝对值代入函数。

  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
以下是大整数乘法分治法的C代码: ```c #include <stdio.h> #include <stdlib.h> #include <string.h> // 用于将字符转换为数字 int char_to_int(char c) { return (int)(c - '0'); } // 用于将数字转换为字符 char int_to_char(int n) { return (char)(n + '0'); } // 用于将字符串表示的数字转换为整数数组 void string_to_array(char *str, int *a, int len) { int i, j; for (i = 0, j = len - 1; i < len; i++, j--) { a[i] = char_to_int(str[j]); } } // 用于将整数数组表示的数字转换为字符串 void array_to_string(int *a, int len, char *str) { int i, j; for (i = len - 1, j = 0; i >= 0; i--, j++) { str[j] = int_to_char(a[i]); } str[len] = '\0'; } // 用于向数组中添加零 void append_zero(int *a, int len, int n) { int i; for (i = len; i < len + n; i++) { a[i] = 0; } } // 用于删除数组中前导零 int remove_leading_zeros(int *a, int len) { int i; for (i = len - 1; i > 0; i--) { if (a[i] != 0) { break; } } return i + 1; } // 用于将两个整数数组相加 void add(int *a, int *b, int *c, int len) { int i, carry = 0; for (i = 0; i < len; i++) { c[i] = a[i] + b[i] + carry; carry = c[i] / 10; c[i] %= 10; } if (carry > 0) { c[i] = carry; } } // 用于将两个整数数组相减 void subtract(int *a, int *b, int *c, int len) { int i, borrow = 0; for (i = 0; i < len; i++) { c[i] = a[i] - b[i] - borrow; if (c[i] < 0) { c[i] += 10; borrow = 1; } else { borrow = 0; } } } // 用于将两个整数数组相乘 void multiply(int *a, int *b, int *c, int len) { if (len == 1) { c[0] = a[0] * b[0]; return; } int i, j, k; int *a1 = (int *)malloc(len / 2 * sizeof(int)); int *b1 = (int *)malloc(len / 2 * sizeof(int)); int *c1 = (int *)malloc(len * sizeof(int)); int *a2 = (int *)malloc(len / 2 * sizeof(int)); int *b2 = (int *)malloc(len / 2 * sizeof(int)); int *c2 = (int *)malloc(len * sizeof(int)); int *a1b1 = (int *)malloc(len * sizeof(int)); int *a2b2 = (int *)malloc(len * sizeof(int)); int *tmp = (int *)malloc(len * sizeof(int)); memset(a1, 0, len / 2 * sizeof(int)); memset(b1, 0, len / 2 * sizeof(int)); memset(c1, 0, len * sizeof(int)); memset(a2, 0, len / 2 * sizeof(int)); memset(b2, 0, len / 2 * sizeof(int)); memset(c2, 0, len * sizeof(int)); memset(a1b1, 0, len * sizeof(int)); memset(a2b2, 0, len * sizeof(int)); memset(tmp, 0, len * sizeof(int)); for (i = 0; i < len / 2; i++) { a1[i] = a[i]; b1[i] = b[i]; } for (i = len / 2; i < len; i++) { a2[i - len / 2] = a[i]; b2[i - len / 2] = b[i]; } multiply(a1, b1, c1, len / 2); multiply(a2, b2, c2, len / 2); add(a1, a2, tmp, len / 2); add(b1, b2, tmp + len / 2, len / 2); multiply(tmp, tmp + len / 2, a1b1, len / 2 + 1); subtract(a1b1, c1, a1b1, len * 2); subtract(a1b1, c2, a1b1, len * 2 - (len / 2)); append_zero(a1b1, len * 2 - (len / 2), len / 2); add(c1, a1b1, c1, len * 2); append_zero(c2, len, len / 2); add(c1, c2, c, len * 2); free(a1); free(b1); free(c1); free(a2); free(b2); free(c2); free(a1b1); free(a2b2); free(tmp); } int main() { char a[1001]; char b[1001]; int len1, len2, len, i; int *num1, *num2, *result; scanf("%s", a); scanf("%s", b); len1 = strlen(a); len2 = strlen(b); len = len1 > len2 ? len1 : len2; num1 = (int *)malloc(len * sizeof(int)); num2 = (int *)malloc(len * sizeof(int)); result = (int *)malloc(len * 2 * sizeof(int)); memset(num1, 0, len * sizeof(int)); memset(num2, 0, len * sizeof(int)); memset(result, 0, len * 2 * sizeof(int)); string_to_array(a, num1, len1); string_to_array(b, num2, len2); multiply(num1, num2, result, len); i = remove_leading_zeros(result, len * 2); if (i == 0) { printf("0\n"); } else { array_to_string(result, i, a); printf("%s\n", a); } free(num1); free(num2); free(result); return 0; } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值