大整数乘法c语言代码_[算法]大整数乘法

最新推荐文章于 2023-11-21 23:50:40 发布
最新推荐文章于 2023-11-21 23:50:40 发布
阅读量826 收藏
点赞数
文章标签： 大整数乘法c语言代码

利用二分法和递归计算任意长度整数相乘

以下复杂度分析有问题，在于 划分为 A12(n2)，这样才相当于移位；

程序中采用string直接+'0'的方式来*10

第一次的代码有漏洞，已更正

我们可以把规模n变成n/2和n/2（把以1位为单位规模为n的问题 变成 以n/2为单位的规模为2的问题），把规模m变成m/2和m/2（把以1位为单位规模为m的问题 变成 以m/2为单位的规模为2的问题），如此，原来的大整数相乘就变成了两个2位数相乘，只不过低位向高位的进制不再是10，而是和。更一般地，我们把整数A由规模n分为n1和n2，把整数B由规模m分为m1和m2，如下图：

651864cf16eff5ad7eb0baf08998c280.png

b994ed02bb9509943b87f50992bbbd0b.png

则A分为n1位的A1和n2位，B分为m1位的B1和m2位的B2，如下式所示：

以此类推，我们可以把A1、A2、B1、B2继续划分，直至最小单位。（这里在编程时需要用递归实现)

于是

注意用2^n做乘法相当于移位运算，需要O(n)时间。此公式中有4次乘法运算，3次加法运算，则

故T(n)=O(n^2)

而

A1B1 A2B2已经计算过，故乘法运算次数减少为3次

产生递推式

// 四次乘法，3次加法
#include <iostream>
#include <string>
#include <sstream>

using namespace std;

string multi(string A, string B); //计算大整数相乘
string Plus(string q, string w, string e, string r); //计算大整数相加
stringstream ss;

int main() {
      string A, B;
      while (cin >> A >> B) {
      cout << multi(A, B) << endl;
      }
      return 0;
}

string multi(string A, string B) {
      int len_A = A.length();
      int len_B = B.length();
      if (len_A == 1) {
      if (len_B == 1) { //最基本的情况：A和B都是一位数，把A、B从string转为int（我这里用的stringstream），然后相乘后转回为string型return回去。
      ss << A;
      int a;
      ss >> a;
      ss.clear();
      ss << B;
      int b;
      ss >> b;
      ss.clear();
      ss << b*a;
      string str_out;
      ss >> str_out;
      ss.clear();
      return str_out;
      }
      else {//A是个位数，B不是的情况下，按照分治的思想把B分开分别与A相乘。
      string B1, B2;
      B1 = B.substr(0, len_B / 2);
      B2 = B.substr(len_B / 2);
      string AB1 = multi(A, B1);
      string AB2 = multi(A, B2);
      cout<<AB1<<endl<<AB2<<endl;
      if (AB2.length() > B2.length()) {   //此时AB2最多比B2多出一位
      string str = AB2.substr(0, 1);
        /*ss << AB1;     漏洞,此时AB1可能已经超出int范围了
        int ab1;
        ss >> ab1;*/
        string tp0(AB1.length()-1,'0');
        str=tp0+str;
        string u0(AB1.length(),'0');
        AB1=Plus(AB1,str,u0,u0);
        return AB1 + AB2.substr(1);
      }
      else
        return AB1 + AB2;
        }
    }
      else {
      if (len_B == 1) {//B是个位数，A不是的情况与上述A是个位数B不是的情况相同。
      string A1, A2;
      A1 = A.substr(0, len_A / 2);
      A2 = A.substr(len_A / 2);
      string A1B = multi(A1, B);
      string A2B = multi(A2, B);
      if (A2B.length() > A2.length()) {
        string str = A2B.substr(0, 1);
        string tp0(A1B.length()-1,'0');
        str=tp0+str;
        string u0(A1B.length(),'0');
        A1B=Plus(A1B,str,u0,u0);
        return A1B + A2B.substr(1);
        }
      else {
        return A1B + A2B;
        }
        }
      else {//A和B都不是个位数，就按照上述方法分治就可以了，只是为了最后相加的时候方便，把返回的四个部分都用0凑成了位数相同的。
      string A1, A2, B1, B2;
      A1 = A.substr(0, len_A / 2);
      A2 = A.substr(len_A / 2);
      B1 = B.substr(0, len_B / 2);
      B2 = B.substr(len_B / 2);
      string part1_ = multi(A1, B1);
      string part1_0(A2.length()+B2.length(), '0');   //长度为A2.length()+B2.length()的0串
      part1_ = part1_ + part1_0;                      //A1B1*10(n2,m2)
      string part2_ = multi(A2, B2);
      string part2_00(part1_.length() - part2_.length(), '0');
      part2_ = part2_00 + part2_;                     //A2B2，高位补0
      string part3_ = multi(A1, B2);
      string part3_0(A2.length(), '0');
      part3_ = part3_ + part3_0;                    //A1B2*10(n2)
      string part3_00(part1_.length() - part3_.length(), '0');
      part3_ = part3_00 + part3_;
      string part4_ = multi(A2, B1);
      string part4_0(B2.length(), '0');
      part4_ = part4_ + part4_0;
      string part4_00(part1_.length() - part4_.length(), '0');
      part4_ = part4_00 + part4_;
      return Plus(part1_, part2_, part3_, part4_);   //未改进的第一种算法
    }
    }
}

string Plus(string q, string w, string e, string r) { //大整数相加，qwer位数相同
      int len_q = q.length();
      string y, out;
      int a, b, c, d;
      for (int i = 1; i <= len_q; i++) {  //逐位相加
    ss << q.substr(len_q - i, 1);
    ss >> a;
    ss.clear();
    ss << w.substr(len_q - i, 1);
    ss >> b;
    ss.clear();
    ss << e.substr(len_q - i, 1);
    ss >> c;
    ss.clear();
    ss << r.substr(len_q - i, 1);
    ss >> d;
    ss.clear();
    ss << a + b + c + d;
    ss >> y;
    ss.clear();
      if (i == 1)
      out = y;
      else if (out.length() > i - 1) {  //进一位
      ss << out.substr(0, 1);
      ss >> a;
      ss.clear();
      ss << y;
      ss >> b;
      ss.clear();
      ss << a + b;
      ss >> y;
      ss.clear();
      out = y + out.substr(1);
    }
      else {
      out = y + out;
    }
    }
	return out;
}

3858520f672963cd365a574a8f8458a8.png 
// 三次乘法。由于还没有实现大整数减法（或支持负数的加法），暂时还是伪代码
#include <iostream>
#include <string>
#include <sstream>

using namespace std;

string multi(string A, string B); //计算大整数相乘
string Plus(string q, string w, string e, string r); //计算大整数相加
stringstream ss;

int main() {
      string A, B;
      while (cin >> A >> B) {
      cout << multi(A, B) << endl;
      }
      return 0;
}

string multi(string A, string B) {
      int len_A = A.length();
      int len_B = B.length();
      if (len_A == 1) {
      if (len_B == 1) { //最基本的情况：A和B都是一位数，把A、B从string转为int（我这里用的stringstream），然后相乘后转回为string型return回去。
      ss << A;
      int a;
      ss >> a;
      ss.clear();
      ss << B;
      int b;
      ss >> b;
      ss.clear();
      ss << b*a;
      string str_out;
      ss >> str_out;
      ss.clear();
      return str_out;
      }
      else {//A是个位数，B不是的情况下，按照分治的思想把B分开分别与A相乘。
      string B1, B2;
      B1 = B.substr(0, len_B / 2);
      B2 = B.substr(len_B / 2);
      string AB1 = multi(A, B1);
      string AB2 = multi(A, B2);
      cout<<AB1<<endl<<AB2<<endl;
      if (AB2.length() > B2.length()) {   //此时AB2最多比B2多出一位
      string str = AB2.substr(0, 1);
        /*ss << AB1;     漏洞,此时AB1可能已经超出int范围了
        int ab1;
        ss >> ab1;*/
        string tp0(AB1.length()-1,'0');
        str=tp0+str;
        string u0(AB1.length(),'0');
        AB1=Plus(AB1,str,u0,u0);
        return AB1 + AB2.substr(1);
      }
      else
        return AB1 + AB2;
        }
    }
      else {
      if (len_B == 1) {//B是个位数，A不是的情况与上述A是个位数B不是的情况相同。
      string A1, A2;
      A1 = A.substr(0, len_A / 2);
      A2 = A.substr(len_A / 2);
      string A1B = multi(A1, B);
      string A2B = multi(A2, B);
      if (A2B.length() > A2.length()) {
        string str = A2B.substr(0, 1);
        string tp0(A1B.length()-1,'0');
        str=tp0+str;
        string u0(A1B.length(),'0');
        A1B=Plus(A1B,str,u0,u0);
        return A1B + A2B.substr(1);
        }
      else {
        return A1B + A2B;
        }
        }
      else {//A和B都不是个位数，就按照上述方法分治就可以了，只是为了最后相加的时候方便，把返回的四个部分都用0凑成了位数相同的。
      string A1, A2, B1, B2;
      A1 = A.substr(0, len_A / 2);
      A2 = A.substr(len_A / 2);
      B1 = B.substr(0, len_B / 2);
      B2 = B.substr(len_B / 2);
      string part1_ = multi(A1, B1);
	  part1_ = multi(2, part1_);
      string part1_0(A2.length()+B2.length(), '0');   //长度为A2.length()+B2.length()的0串
      part1_ = part1_ + part1_0;                      //2*A1B1*10(n2,m2)
      string part2_ = multi(A2, B2);
	  part2_ =multi(2,part2_);
      string part2_00(part1_.length() - part2_.length(), '0');
      part2_ = part2_00 + part2_;                     //2*A2B2，高位补0
	  
	  
      string part3_ = A1;
      string part3_0(A2.length(), '0');             //(A1*10(n2)-A2)
      part3_ = part3_ + part3_0;                    //A1*10(n2)  伪代码
	  part3_= part_3-A2                             //大整数减法（或支持负数的大整数加法），未实现
      
	  
	  
      string part4_ = B1;                           //(B2-B1*10(m2))
      string part4_0=(B2.length(), '0');
      part4_ = part4_ + part4_0;
	  part4_= B2-part4_;                        
	  
	  part5_= multi(part4_,part3_)
	  
      string part5_00(part1_.length() - part5_.length(), '0');
      part5_ = part5_00 + part5_;
	  string u0(part1_.length(), '0');
      return Plus(part1_, part2_, part5_, u0);   
    }
    }
}

string Plus(string q, string w, string e, string r) { //大整数相加，qwer位数相同
      int len_q = q.length();
      string y, out;
      int a, b, c, d;
      for (int i = 1; i <= len_q; i++) {  //逐位相加
    ss << q.substr(len_q - i, 1);
    ss >> a;
    ss.clear();
    ss << w.substr(len_q - i, 1);
    ss >> b;
    ss.clear();
    ss << e.substr(len_q - i, 1);
    ss >> c;
    ss.clear();
    ss << r.substr(len_q - i, 1);
    ss >> d;
    ss.clear();
    ss << a + b + c + d;
    ss >> y;
    ss.clear();
      if (i == 1)
      out = y;
      else if (out.length() > i - 1) {  //进一位
      ss << out.substr(0, 1);
      ss >> a;
      ss.clear();
      ss << y;
      ss >> b;
      ss.clear();
      ss << a + b;
      ss >> y;
      ss.clear();
      out = y + out.substr(1);
    }
      else {
      out = y + out;
    }
    }
	return out;
}

参考

https://blog.csdn.net/qq_36165148/article/details/81132525​blog.csdn.net
weixin_39755890
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
大整数乘法C语言版本
09-24
分治思想大整数乘法 基于王晓东版 算法设计与分析
大整数乘法---C语言实现
无言月梧桐的博客
06-13 5819
一.前言 如上 二.内容 在科学计算中，无法计算非常大的整数，于是就有大整数乘法问题。 （int存放4字节，float…已有的数据类型数据范围有限），硬件无法实现，只能靠软件（代码）实现了。 三.代码 #include<stdio.h> #include<string.h> int main(){ int i,j;//循环变量 char char1[256],char2[256];//存放两个大整数的字符 int int1[256],int2[256];//按低位到高位存
C语言之两个大数相乘
最新发布
m0_60963435的博客
11-21 979
​int main()//输入两个大数​。
C语言实现大整数乘法
热门推荐
selina8921的博客
03-30 2万+
转载自：点击打开链接乘法规律，一个数的第i位和另一个数的第j位相乘，一定会累加到结果的第i+j位，结果的数组一个数组元素存2位数，最后对结果处理进位，最后打印出来方法一见上面链接https://www.cnblogs.com/king-ding/p/bigIntegerMul.html方法二void IntMultiply(int a[], int b[], int c[], int ma, in...
C语言——大整数乘法
肥不溜秋大梅子
12-01 4350
Description   求两个不超过200位的非负整数的积。 Input   输入有两行，每行是一个不超过200位的非负整数，没有多余的前导0。 Output   一行，即相乘后的结果。结果里不能有多余的前导0，即如果结果是342，那么就不能输出为0342。 Sample Input   12345678900 98765432100 Sample Output   1...
大整数乘法 C语言。。。
zy
03-24 285
用字符数组接受大整数 然后根据两数相乘的基本法则进行运算 #include<stdio.h> #include<string.h> int c[200] = { 0 }; int len; void cheng(char a[], char b[]) { for(int i = strlen(a) - 1; i >= 0; i--) for (int j ...
大整数相乘算法c语言实现源代码
05-02
这个东西还是很不错的，我们正在学习算法，学完之后，对这个问题有了很深的理解，希望大家能够从中获益。
c语言大整数算法链表 加减乘
11-23
此算法用单链表实现，代码量有点多。功能基本上能实现。
大整数乘法
03-18
cout乘法，除法和求余只能完成一定范围内数的运算"; } void add(char a[],char b[],int x,int y,int f,int d) { char c[102]; int i; int j; int jin=0; int s=0; int m,ma; int xx,yy; if(d==0) { xx=0;yy=0; } ...
分治法求两个大整数相乘
03-12
分治法求两个大整数相乘C++实现。
大整数乘法(c语言+java+py)
m0_46625346的博客
10-09 867
1.乘法规律，一个数的第i位和另一个数的第j位相乘，一定会累加到结果的第i+j位，结果的数组一个数组元素存2位数，最后对结果处理进位，最后打印出来。 注：相乘后的位数k不会超过m+n。 乘法计算步骤如下： a数组存一个数，b数组存一个数，按位相乘后放入c数组中。 C语言代码：（时间复杂度为小于O(n²)，分治法，错位相加） #include<stdio.h> #include<math.h> #include<string.h> #define M 10000 char
大数乘法运算(C语言)
小黯
06-02 4403
大数运算 利用乘法法则，相乘，然后进位，取余 详细过程-如图 源代码 #include<stdio.h> #include<string.h> #define N 100 int main() { int i,j,L[N]; char m[N],n[N]; scanf("%s %s",m,n); int len1=strlen(m); int len2=strlen(n); memset(L,0,sizeof(L));//初始化为0 //第一步 for(i=0;
大数乘法（C语言实现）
萧海的博客
06-04 4375
实现过程分析： 我们回忆一下，在我们小时候刚接触多位数的乘法，我们的数学老师会教给我们一个方法，那就是“乘法的竖式计算”。在这里我们就采用该思想解决大数乘法的问题。 以下是我们经常进行乘法的竖式运算： 根据以上的竖式运算，我们实现过程总结如下： 1、先使用两个字符数组保存两个大数据； 2、用第一个数据的个位与第二个数据的所有位相乘，并将每一位的运算结果保存在暂存字符数组temp中，并进行进位调整，即如果该位的数值大于9，就将该数值的十位加到前一位，并将该位的个位保存在原位。 3、将temp与结果数组rst
C语言实现大数乘法(附完整源码)
希望我的博客，能帮上你解决学习中工作中所遇到的问题
04-23 157
C语言实现大数乘法(附完整源码)
用C语言实现动态数组的大整数乘法
FlushHip
04-28 1601
描述： 要求你用C语言实现两个数相乘，这两个数字的位数可能很大，而且还可能有不合法的输入。 具体做法： 首先，大整数乘法的主要部分要知道： for (i = 0; i < lenA; i++) for (j = 0; j < lenB; j++) ans[i + j] += (int)(pA[i] - '0') * (
超大整数乘法运算——C语言
hedafighter2011的专栏
04-10 4213
超大整数乘法运算——C语言         在计算机中，长整型(long int)变量的范围是 -2147483648 至 2147483647，因此若用长整型变量做乘法运算，乘积最多不能超过 10位数。即便用双精度型(double)变量，也仅能保证 16 位有效数字的精度。在某些需要更高精度的乘法运算的场合，需要用别的办法来实现乘法运算。          比较容易想
【C语言】大整数乘法代码
weixin_71027076的博客
10-12 29
具体的内容都写在注释里面啦，还有不理解的评论问我哦！
C语言大整数运算
Solidao的博客
04-15 2710
在上面的代码中，我们定义了一个sub()函数，该函数接受两个大整数a和b，并将它们相减后存储到结果数组res中。从最低位开始遍历两个大整数，依次相减，并将结果存储到结果数组中。在C语言中，通常使用数组来存储大整数，并通过编写相应的函数实现加、减、乘、除等运算。在上面的代码中，我们定义了一个add()函数，该函数接受两个大整数a和b，并将它们相加后存储到结果数组res中。总之，在进行大整数运算时需要注意细节，并且需要编写相应的函数来实现各种运算。即从高位开始逐步计算商和余数，并将结果存储到对应的数组中。
分治法大整数乘法c语言代码
09-14
分治法是一种算法设计策略，通过将问题分解成更小的子问题，再通过将解合并起来来解决原始问题。在大整数乘法问题中，可以使用分治法来提高算法效率。 下面是一个使用C语言编写的分治法大整数乘法的代码示例： #include <stdio.h> #include <math.h> long long int power(int a, int b) { if (b == 0) { return 1; } else { return a * power(a, b - 1); } } long long int multiply(long long int x, long long int y) { int n; if (x < 10 || y < 10) { return x * y; } n = fmax(log10(x) + 1, log10(y) + 1); int m = ceil(n / 2.0); long long int a = x / power(10, m); long long int b = x % power(10, m); long long int c = y / power(10, m); long long int d = y % power(10, m); long long int ac = multiply(a, c); long long int bd = multiply(b, d); long long int ad_bc = multiply(a + b, c + d) - ac - bd; return ac * power(10, 2 * m) + ad_bc * power(10, m) + bd; } int main() { long long int x, y; printf("请输入两个整数："); scanf("%lld %lld", &x, &y); long long int result = multiply(x, y); printf("两数乘积为：%lld\n", result); return 0; } 在这个代码示例中，使用了递归来实现分治法的思想。首先判断输入的两个数是否小于10，如果是，则直接返回乘积。如果不是，则计算出这两个数的位数n，然后取n的一半向上取整得到m。分别将两个数的高位和低位等分为a、b、c、d四部分。 接着，使用递归调用multiply函数来计算ac、bd和ad_bc三个结果，最后将这三个结果按照指定的规则相加得到最终的乘积。 通过使用分治法，可以显著提高大整数乘法的效率。因为将问题分解成更小的子问题，每个子问题处理的数据规模更小，计算量相对减少，从而提高了算法的效率。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值