大整数乘除法

 参考https://blog.youkuaiyun.com/Wang_1997/article/details/68241892 

       https://blog.youkuaiyun.com/Wchenchen0/article/details/80286915?utm_source=blogxgwz27

大整数乘法

在下面的例子程序中,用 unsigned an1[200]和 unsigned an2[200]分别存放两个乘数,用aResult[400]来存放积。计算的中间结果也都存在 aResult 中。

aResult 长度取 400 是因为两个200 位的数相乘,积最多会有 400 位。 an1[0], an2[0], aResult[0]都表示个位。计算的过程基本上和小学生列竖式做乘法相同。

为编程方便,并不急于处理进位,而将进位问题留待最后统一处理。现以 835×49 为例来说明程序的计算过程。

        先算 835×9。 5×9 得到 45 个 1, 3×9 得到 27 个 10, 8×9 得到 72 个 100。由于不急于处理进位,所以 835×9 算完后, aResult 如下:

 

 

        接下来算 4×5。此处 4×5 的结果代表 20 个 10,因此要 aResult[1]+=20,变为:

 

 

        再下来算 4×3。此处 4×3 的结果代表 12 个 100,因此要 aResult[2]+= 12,变为:

 

 

        最后算 4×8。此处 4×8 的结果代表 32 个 1000,因此要 aResult[3]+= 32,变为:

 

 

        乘法过程完毕。接下来从 aResult[0]开始向高位逐位处理进位问题。 aResult[0]留下 5,把 4 加到 aResult[1]上, aResult[1]变为 51 后,应留下 1,把 5 加到 aResult[2]上……最终使得 aResult 里的每个元素都是 1 位数,结果就算出来了:

 

 

       总结一个规律,即一个数的第 i 位和另一个数的第 j 位相乘所得的数,一定是要累加到结果的第 i+j 位上。这里 i, j 都是从右往左,从 0 开始数。
     下面是根据这位博主思路自己模拟的一个demo

#include<bits/stdc++.h>
using namespace std;
char array1[200];
char array2[200];
int sum[400];
f(string s1,string s2)
{
reverse(s1.begin(),s1.end());
reverse(s2.begin(),s2.end());
// cout << s1 << s2<< endl;
memset(array1,0,sizeof(s1));
memset(array2,0,sizeof(s2));
int l1=s1.length();
int l2=s2.length();
strcpy(array1,s1.c_str());
strcpy(array2,s2.c_str());
//cout << l1 << l2<<endl;
for(int i=0;i<l1;i++){
for(int j=0;j<l2;j++){
sum[i+j]+=((int)array1[i]-48)*((int)array2[j]-48);
}
}
int f=0;
while(true){

if(f<l1+l2){
int de=sum[f]/10;
sum[f]=sum[f]%10;
sum[f+1]+=de;
f++;
}
else{
break;
}

}

for(int i=l1+l2-2;i>=0;i--){
cout << sum[i];
}
}
int main()
{
string s1;
string s2;
cin >> s1 >> s2;
f(s1,s2);

}

大整数除法

#include<bits/stdc++.h>
#define Max 100005
using namespace std;
int main()
{
int T;
cin>>T;
while(T--)
{
string a;
char s[Max];
cin>>a;
int b;
cin>>b;
long long sum=0;
//long long c;
for(int i=0;a[i];i++)
{
sum=sum*10+a[i]-'0';
s[i]=sum/b+'0';
sum=sum%b;
}
int j=0;
while(s[j]=='0')
j++;
for(;j<a.size();j++)
cout<<s[j];//商
cout<<endl;
cout<<sum<<endl;//余数
}
return 0;
}
/*
for(i = 0; i < len1; i++){
b = b * 10 + num1[i] - '0';
num3[i] = b / d + '0';
b = b % d;
}
这个实际上就是笔算的原理,比如,128 / 12, 我们笔算的时候,
先计算 1 / 12, 此时的b = 1, num3[0] = 0, b % d = 1;
然后 b = 1 * 10 + 2 = 12, 所以相当于 12 / 12, num3[1] = 1, b % d = 0;
最后 b = 8, num3[2] = 0, b % d = 8;
所以商为10,余数为8,而余数也就是模,所以如果说求模,就必须理解好这两步的反复执行:
b = b * 10 + num1[i] - '0'; 这是b在与要看的下一位组成一个新的数字
b = b % d; 这可以理解为b除以d之后的余数,所以既然能求出商,自然最后循环的结果也就能求出
*/

转载于:https://www.cnblogs.com/henuliulei/p/9867127.html

以下是C语言实现大整数加减乘除法的示例代码: 1. 大整数加法 ```c #include<stdio.h> #include<string.h> #define MAXN 1010 int a[MAXN], b[MAXN], c[MAXN]; void add(int a[], int b[], int c[]) { int carry = 0; for (int i = 0; i < MAXN; i++) { c[i] = a[i] + b[i] + carry; if (c[i] >= 10) { carry = 1; c[i] -= 10; } else { carry = 0; } } } int main() { memset(a, 0, sizeof(a)); memset(b, 0, sizeof(b)); memset(c, 0, sizeof(c)); char str1[MAXN], str2[MAXN]; scanf("%s %s", str1, str2); int len1 = strlen(str1), len2 = strlen(str2); for (int i = 0; i < len1; i++) { a[len1 - i - 1] = str1[i] - '0'; } for (int i = 0; i < len2; i++) { b[len2 - i - 1] = str2[i] - '0'; } add(a, b, c); int i; for (i = MAXN - 1; i >= 0; i--) { if (c[i] != 0) { break; } } for (; i >= 0; i--) { printf("%d", c[i]); } return 0; } ``` 2. 大整数减法 ```c #include<stdio.h> #include<string.h> #define MAXN 1010 int a[MAXN], b[MAXN], c[MAXN]; void sub(int a[], int b[], int c[]) { for (int i = 0; i < MAXN; i++) { c[i] = a[i] - b[i]; if (c[i] < 0) { c[i] += 10; a[i + 1]--; } } } int main() { memset(a, 0, sizeof(a)); memset(b, 0, sizeof(b)); memset(c, 0, sizeof(c)); char str1[MAXN], str2[MAXN]; scanf("%s %s", str1, str2); int len1 = strlen(str1), len2 = strlen(str2); for (int i = 0; i < len1; i++) { a[len1 - i - 1] = str1[i] - '0'; } for (int i = 0; i < len2; i++) { b[len2 - i - 1] = str2[i] - '0'; } sub(a, b, c); int i; for (i = MAXN - 1; i >= 0; i--) { if (c[i] != 0) { break; } } for (; i >= 0; i--) { printf("%d", c[i]); } return 0; } ``` 3. 大整数乘法 ```c #include<stdio.h> #include<string.h> #define MAXN 1010 int a[MAXN], b[MAXN], c[MAXN]; void mul(int a[], int b[], int c[]) { memset(c, 0, sizeof(c)); for (int i = 0; i < MAXN; i++) { for (int j = 0; j < MAXN; j++) { c[i + j] += a[i] * b[j]; } } for (int i = 0; i < MAXN - 1; i++) { c[i + 1] += c[i] / 10; c[i] %= 10; } } int main() { memset(a, 0, sizeof(a)); memset(b, 0, sizeof(b)); memset(c, 0, sizeof(c)); char str1[MAXN], str2[MAXN]; scanf("%s %s", str1, str2); int len1 = strlen(str1), len2 = strlen(str2); for (int i = 0; i < len1; i++) { a[len1 - i - 1] = str1[i] - '0'; } for (int i = 0; i < len2; i++) { b[len2 - i - 1] = str2[i] - '0'; } mul(a, b, c); int i; for (i = MAXN - 1; i >= 0; i--) { if (c[i] != 0) { break; } } for (; i >= 0; i--) { printf("%d", c[i]); } return 0; } ``` 4. 大整数除法 ```c #include<stdio.h> #include<string.h> #define MAXN 1010 int a[MAXN], b[MAXN], c[MAXN], d[MAXN]; void div(int a[], int b[], int c[], int d[]) { memset(c, 0, sizeof(c)); memset(d, 0, sizeof(d)); int lena, lenb, len; for (lena = MAXN - 1; lena >= 0 && a[lena] == 0; lena--); for (lenb = MAXN - 1; lenb >= 0 && b[lenb] == 0; lenb--); if (lenb > lena) { len = 1; d[0] = 0; return; } len = lena - lenb + 1; for (int i = lena; i >= lenb; i--) { int t = 0; while (a[i] >= b[lenb - 1] && t < 10) { t++; int carry = 0; for (int j = 0; j < lenb; j++) { int tmp = a[i - j] - t * b[lenb - j - 1] - carry; if (tmp < 0) { tmp += 10; carry = 1; } else { carry = 0; } a[i - j] = tmp; } c[i - lenb + 1] = t; } } for (int i = len - 1; i >= 0; i--) { d[i] = a[i]; } } int main() { memset(a, 0, sizeof(a)); memset(b, 0, sizeof(b)); memset(c, 0, sizeof(c)); memset(d, 0, sizeof(d)); char str1[MAXN], str2[MAXN]; scanf("%s %s", str1, str2); int len1 = strlen(str1), len2 = strlen(str2); for (int i = 0; i < len1; i++) { a[len1 - i - 1] = str1[i] - '0'; } for (int i = 0; i < len2; i++) { b[len2 - i - 1] = str2[i] - '0'; } div(a, b, c, d); int i; for (i = MAXN - 1; i >= 0; i--) { if (c[i] != 0) { break; } } for (; i >= 0; i--) { printf("%d", c[i]); } printf("\n"); for (i = MAXN - 1; i >= 0; i--) { if (d[i] != 0) { break; } } for (; i >= 0; i--) { printf("%d", d[i]); } return 0; } ```
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值