大整数加法与乘法(高精度算法)
一、大整数加法
原理:
当数据的位置太长时,我们品是剋有运用的short、 int 、 long long 型变量都无法满足要求,我们就需要用到大整数的加法运算思想了:
要点:
- 开一个数组,存储需要计算的数字,注意,数组的首位存的是数字的长度,其他位将数字的每一位倒过来存放
- 定义一个总和的数组,和上面的数组一样存
- 相加过程中,数字会出现进位的情况,进位有可能会导致综合的长度发生改变
例题
-
计算两个数字的和
#include<iostream> #include<cstring> using namespace std; int num1[1005], num2[1005] , sum[1005]; char s1[1005], s2[1005]; int main() { cin >> s1 >> s2; num1[0] = strlen(s1);//两个数组的首位存储字符串的长度 num2[0] = strlen(s2); //数字的每一位到这存储到num数组中 for (int i = 0, j = num1[0]; i < num1[0]; i++, j--) { num1[j] = s1[i] - '0'; } for (int i = 0, j = num2[0]; i < num2[0]; i++, j--) { num2[j] = s2[i] - '0'; } sum[0] = max(num1[0], num2[0]);//两个数相加,最小长度是较大长度的数字的长度 for (int i = 1; i <= sum[0]; i++) { sum[i] = num1[i] + num2[i]; //对应的位相加 //进位操作 if(sum[i] > 9) { sum[i + 1] += sum[i] / 10; sum[i] %= 10; if ( i == sum[0]) {//进位可能会改变数字的长度 sum[0]++; } } } for (int i = sum[0]; i > 0; i--) { cout << sum[i]; } cout << endl; return 0; }
-
euler13(具体清搜索网站“欧拉计划”)
#include<iostream>
#include<cstring>
using namespace std;
int num1[1005], num2[1005], sum[1100], t = 50;
char s1[1005], s2[1005];
int main() {
while(t--) {
cin >> s1 >> s2;
num1[0] = strlen(s1);
num2[0] = strlen(s2);
for (int i = 0, j = num1[0]; i < num1[0]; i++, j--) {
num1[j] = s1[i] - '0';
}
for (int i = 0, j = num2[0]; i < num2[0]; i++, j--) {
num2[j] = s2[i] - '0';
}
int k = max(num1[0], num2[0]);
sum[0] = max(sum[0], k);
for (int i = 1; i <= sum[0]; i++) {
sum[i] += num1[i] + num2[i];
if (sum[i] > 9) {
sum[i + 1] += sum[i] / 10;
sum[i] %= 10;
if ( i == sum[0] ) sum[0]++;
}
}
}
for (int i = sum[0]; i > 0; i--) {
cout << sum[i];
}
cout << endl;
return 0;
}
3.euler25
分析:
首先我们先用两个数据相加,比如 a, b, 将a加到b中,此时b变为了他们之和,再交换将b加到a中,就实现了斐波那契数列的相加,模拟一下:
一开始 : a = 1, b = 1, b+=a ,b= 2 -----> a+= b, a = 3, b= 2 -----> b+= a, b = 5 , a =3-----…就这样子循环下去就好了
#include<iostream>
using namespace std;
int func(int *n1, int *n2) {//n1是数组num[a], n2是数组num[b]
n2[0] = n1[0];
for (int i = 1; i <= n2[0]; i++) {
n2[i] += n1[i];
if (n2[i] > 9) {//进位
n2[i + 1] += n2[i] / 10;
n2[i] %= 10;
if (i == n2[0]) n2[0]++;
}
}
return n2[0] == 1000;//如果有1000位,返回1
}
int main() {
int num[2][1005] = {{1, 1}, {1, 1}};//num[0][0]代表长度,num[0][1]代表第一个数值是1
int a = 0, b = 1;
for (int i = 3; 1; i++) {
if (func(num[a], num[b])) {
cout << i << endl;
break;
}
swap(a, b);//num[a] --> num[b] ,变为num[b] --->num[a]
}
return 0;
}
二、大整数乘法
原理:
见下图:
要点:
- 数字储存到数组中位置: 比如 num1[i] * num2[j]—》 那么放在在ans[i + j -1]的位置,可以看上图模拟出来
- ans数组中乘积结果的初始长度 ,定义为乘积的最小长度 i + j - 1
代码
#include<iostream>
#include <cstring>
using namespace std;
char s1[1005], s2[1005];//要计算的字符串数据
int num1[1005], num2[1005], ans[2005]; //存两个乘数, 还有存相乘结果的数组
int main() {
cin >> s1 >> s2;
num1[0] = strlen(s1);
num2[0] = strlen(s2);
//数字倒过来存
for (int i = 1, j = num1[0] - 1; i <= num1[0]; i++, j--) {
num1[i] = s1[j] - '0';
}
for (int i = 1, j = num2[0] - 1; i <= num2[0]; i++, j--) {
num2[i] = s2[j] - '0';
}
ans[0] = num1[0] + num2[0] - 1;//ans中数字的长度(相乘之后的)
for (int i = 1; i <= num1[0]; i++) {//相乘结果存到ans中,顺序存下去
for (int j = 1; j <= num2[0]; j++) {
ans[i + j -1] += num1[i] * num2[j];
}
}
for (int i = 1; i <= ans[0]; i++) {//进位操作
if (ans[i] > 9) {
ans[i + 1] += ans[i] / 10;
ans[i] %= 10;
if (i == ans[0]) ans[0]++;
}
}
for (int i = ans[0]; i > 0; i--) {
cout << ans[i];
}
cout << endl;
return 0;
}
三、大整数减法
直接上代码:
#include<iostream>
#include <cstring>
using namespace std;
char s1[1005], s2[10055];
int num1[1005], num2[1005], ans[1005];
int main() {
cin >> s1 >> s2;
num1[0] = strlen(s1);
num2[0] = strlen(s2);
for (int i = 1, j = num1[0] - 1; i <= num1[0]; i++, j--) {
num1[i] = s1[j] - '0';
}
for (int i = 1, j = num2[0] - 1; i <= num2[0]; i++, j--) {
num2[i] = s2[j] - '0';
}
ans[0] = max(num1[0], num2[0]);
for (int i = 1; i <= ans[0]; i++) {
if (num1[i] < num2[i]) {//直接减,不够减,借1,即+10
num1[i] += 10;
ans[i] = num1[i] - num2[i];
num1[i + 1]--;//借了一位给上面,所以值减1
} else {
ans[i] = num1[i] - num2[i];
}
}
//下面步骤去前导0,因为计算结束后,可能会出现前导0的情况
int m = ans[0];
if (ans[m] == 0) ans[0]--;//去掉最高位的前导0
for (int i = ans[0]; i > 0; i--) {
cout << ans[i];
}
cout << endl;
return 0;
}
写于2021-07-08