C++ 大数乘法

最新推荐文章于 2024-01-28 14:31:58 发布
最新推荐文章于 2024-01-28 14:31:58 发布
阅读量1.9k 收藏 4
点赞数
分类专栏： C++笔记簿 文章标签： c++
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_45464666/article/details/125476279
版权

C++ 大数乘法 分别包括 普通大数乘法、分治法大数乘法、改进的分治法大数乘法。

#include<iostream>
#include<iomanip>
#include<random>
#include<ctime>
#include<cmath>
using namespace std;

int randint(int a, int b);
int* rand_num(int n);
void show_num(int* c, int n);
void multiply_normal(int* c1, int* c2, int n, int flag);
void multiply_n(int* c1, int* c2, int n, int flag);
int bit_number(int* c, int n);
int* high_num_2(int* c, int n);
int* low_num_2(int* c, int n);
void add_list(int* c1, int* c2, int n, int m);
int* multiply_partition(int* c1, int* c2, int n);
int* x_y(int c1, int c2);
void multiply_p(int* c1, int* c2, int n, int f);
void sub_list(int* c1, int* c2, int n, int m);

int* multiply_partition_enhanced(int* c1, int* c2, int n);
void multiply_p_e(int* c1, int* c2, int n, int f);

int main()
{

	cout << "输入数字位数 n = ";
	int n;
	cin >> n;
	cout << endl;
	//生成两个随机 n 位数
	int* c1 = NULL;
	int* c2 = NULL;
	c1 = rand_num(n);
	c2 = rand_num(n);
	//show_num(c1, n);
	//show_num(c2, n);
	//cout <<"正确乘积： " << bit_number(c1, n) * bit_number(c2, n) << endl;
	//int a1[] = { 1, 2, 3, 4, 5, 6 };
	//int a2[] = { 3,6,8 };
	//show_num(a1, 6);
	//show_num(a2, 3);
	//sub_list(a1, a2, 3, 1);
	//show_num(a1, 6);
	cout << "是否显示乘积结果？ 是：1； 否：0  ";
	int flag = 0;
	cout<< "flag = ";
	cin >> flag;
	cout << endl;

	multiply_n(c1, c2, n, flag);
	multiply_p(c1, c2, n, flag);
	multiply_p_e(c1, c2, n, flag);

	delete[]c1;
	delete[]c2;
	return 0;
}



int randint(int a, int b)
{
	random_device rd;
	default_random_engine eng(rd());
	uniform_int_distribution<int> distr1(a, b);//均匀分布的方式产生随机数
	return distr1(eng);
}

int* rand_num(int n)
{
	int* c = new int[n];
	for (int i = 0; i < n - 1; i++)
	{
		c[i] = randint(0, 9);
	}
	c[n - 1] = randint(1, 9);
	return c;
}

void show_num(int* c, int n)
{
	/*cout << "生成的数字为：";*/
	if (c[n - 1] != 0)
	{
		cout << c[n - 1];
	}
	for (int j = n - 2; j >= 0; j--)
	{
		cout << c[j];
	}
	cout << endl;
}

void multiply_normal(int* c1, int* c2, int n, int flag)
{
	clock_t start = clock();      //获取当前系统时间

	int* c = new int[2 * n];
	int x = 0;
	for (int i = 0; i < 2 * n; i++)
	{
		c[i] = 0;
	}
	for (int i = 0; i < n; i++)
	{
		for (int j = 0; j < n; j++)
		{
			x = c[i + j] + c1[i] * c2[j];
			c[i + j] = x % 10;
			c[i + j + 1] += x / 10;
		}
	}

	clock_t end = clock();
	double programTimes = ((double)end - start) / CLOCKS_PER_SEC;
	if (flag)
	{
		cout << "Result." << setw(26) << "普通大数乘法乘积：";
		show_num(c, 2 * n);
	}
	cout << "Time." << setw(28) << "普通大数乘法时间：" << programTimes << " 秒" << endl << endl;
	delete[]c;
}

void multiply_n(int* c1, int* c2, int n, int flag)
{
	multiply_normal(c1, c2, n, flag);
}

int bit_number(int* c, int n)
{
	int x = 0;
	for (int i = 0; i < n; i++)
	{
		x += c[i] * pow(10, i);
	}
	return x;
}


int* x_y(int c1, int c2)
{
	int x = 0;
	int* c = new int[2];
	x = c1 * c2;
	c[0] = x % 10;
	c[1] = x / 10;
	return c;
}

int* high_num_2(int* c, int n)
{
	int* k = new int[n / 2];
	for (int i = 0; i < n / 2; i++)
	{
		k[i] = c[i + n / 2];
	}
	return k;
}

int* low_num_2(int* c, int n)
{
	int* k = new int[n / 2];
	for (int i = 0; i < n / 2; i++)
	{
		k[i] = c[i];
	}
	return k;
}

void add_list(int* c1, int* c2, int n, int m) //n为c2的长度，m为c2的进位
{
	int x = 0;
	for (int i = 0; i < n; i++)
	{
		x = c1[i + m] + c2[i];
		c1[i + m] = x % 10;
		if (x / 10 != 0)
			c1[i + m + 1] += x / 10;
	}
}

int* multiply_partition(int* c1, int* c2, int n)
{
	if (n == 1)
	{
		return x_y(c1[0], c2[0]);
	}
	int* x = NULL;
	int* y = NULL;
	if (n % 2 != 0)
	{
		x = new int[n + 1];
		y = new int[n + 1];
		for (int i = 0; i < n; i++)
		{
			x[i] = c1[i];
			y[i] = c2[i];
		}
		x[n] = 0;
		y[n] = 0;
		n += 1;
	}
	else
	{
		x = new int[n];
		y = new int[n];
		for (int i = 0; i < n; i++)
		{
			x[i] = c1[i];
			y[i] = c2[i];
		}
	}
	int* a1 = high_num_2(x, n);
	int* a0 = low_num_2(x, n);
	int* b1 = high_num_2(y, n);
	int* b0 = low_num_2(y, n);
	int* p = multiply_partition(a1, b1, n / 2);
	int n_p = p[n - 1] == 0 ? n - 1 : n;
	int* q = multiply_partition(a1, b0, n / 2);
	int n_q = q[n - 1] == 0 ? n - 1 : n;
	int* r = multiply_partition(a0, b1, n / 2);
	int n_r = r[n - 1] == 0 ? n - 1 : n;
	int* s = multiply_partition(a0, b0, n / 2);
	int n_s = s[n - 1] == 0 ? n - 1 : n;

	delete[]x;
	delete[]y;
	delete[]a1;
	delete[]a0;
	delete[]b1;
	delete[]b0;

	int* c = new int[2 * n];
	for (int i = 0; i < 2 * n; i++)
	{
		c[i] = 0;
	}
	add_list(c, p, n_p, n);
	add_list(c, r, n_r, n / 2);
	add_list(c, q, n_q, n / 2);
	add_list(c, s, n_s, 0);

	delete[]p;
	delete[]q;
	delete[]r;
	delete[]s;

	return c;
}

void multiply_p(int* c1, int* c2, int n, int flag)
{
	clock_t start = clock();      //获取当前系统时间
	int* c = multiply_partition(c1, c2, n);
	clock_t end = clock();
	double programTimes = ((double)end - start) / CLOCKS_PER_SEC;

	if (flag)
	{
		cout << "Result." << setw(26) << "分治法大数乘法乘积：";
		show_num(c, 2 * n);
	}
	cout << "Time." << setw(28) << "分治法大数乘法时间：" << programTimes << " 秒" << endl << endl;
	delete[]c;
}


void sub_list(int* c1, int* c2, int n, int m)
{
	for (int i = 0; i < n; i++)
	{
		if (c1[i + m] < c2[i])
		{
			c1[i + m] = c1[i + m] + 10 - c2[i];
			c1[i + m + 1] -= 1;
		}
		else
			c1[i + m] -= c2[i];
	}
	if (c1[n + m] < 0)
	{
		for (int i = 0;;i++)
		{
			c1[n + m + i] += 10;
			c1[n + m + i + 1] -= 1;
			if (c1[n + m + i + 1] >= 0)
				break;
		}
	}

}

int* multiply_partition_enhanced(int* c1, int* c2, int n)
{
	if (n == 1)
	{
		return x_y(c1[0], c2[0]);
	}
	int* x = NULL;
	int* y = NULL;
	if (n % 2 != 0)
	{
		x = new int[n + 1];
		y = new int[n + 1];
		for (int i = 0; i < n; i++)
		{
			x[i] = c1[i];
			y[i] = c2[i];
		}
		x[n] = 0;
		y[n] = 0;
		n += 1;
	}
	else
	{
		x = new int[n];
		y = new int[n];
		for (int i = 0; i < n; i++)
		{
			x[i] = c1[i];
			y[i] = c2[i];
		}
	}
	int* a1 = high_num_2(x, n);
	int* a0 = low_num_2(x, n);
	int* b1 = high_num_2(y, n);
	int* b0 = low_num_2(y, n);

	int* p = multiply_partition_enhanced(a1, b1, n / 2);
	int n_p = p[n - 1] == 0 ? n - 1 : n;
	int* a2 = new int[n / 2 + 1];
	int* b2 = new int[n / 2 + 1];
	for (int i = 0; i < n / 2 + 1; i++)
	{
		a2[i] = 0;
		b2[i] = 0;
	}
	add_list(a2, a1, n / 2, 0);
	add_list(a2, a0, n / 2, 0);
	add_list(b2, b1, n / 2, 0);
	add_list(b2, b0, n / 2, 0);

	int m = n / 2;
	if (a2[n / 2] != 0 || b2[n / 2] != 0)
	{
		m = n / 2 + 1;
	}
	int* q = multiply_partition_enhanced(a2, b2, m);	// 如果 n = 2, m = 1 or 2  66*77		12*14 = 168 nq = 3		n / 2 = 1
	int n_q = q[2 * m - 1] == 0 ? 2 * m - 1 : 2 * m;
	int* s = multiply_partition_enhanced(a0, b0, n / 2);
	int n_s = s[n - 1] == 0 ? n - 1 : n;
	delete[]x;
	delete[]y;
	delete[]a1;
	delete[]a0;
	delete[]b1;
	delete[]b0;
	delete[]a2;
	delete[]b2;

	int* c = new int[2 * n + 1];
	for (int i = 0; i < 2 * n + 1; i++)
	{
		c[i] = 0;
	}
	add_list(c, p, n_p, n);
	add_list(c, q, n_q, n / 2);
    add_list(c, s, n_s, 0);
	sub_list(c, p, n_p, n / 2);
	sub_list(c, s, n_s, n / 2);

	delete[]p;
	delete[]q;
	delete[]s;

	return c;
}

void multiply_p_e(int* c1, int* c2, int n, int flag)
{
	clock_t start = clock();      //获取当前系统时间
	int* c = multiply_partition_enhanced(c1, c2, n);
	clock_t end = clock();
	double programTimes = ((double)end - start) / CLOCKS_PER_SEC;

	if (flag)
	{
		cout << "Result." << setw(26) << "改进的分治法大数乘法乘积：";
		show_num(c, 2 * n);
	}
	cout << "Time." << setw(28) << "改进的分治法大数乘法时间：" << programTimes << " 秒" << endl << endl;
	delete[]c;
}

高山爱做梦
关注
  • 0
    点赞
  • 踩
  • 4
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
C++实现大数乘法
weixin_45703960的博客
02-12 1108
#include<stdio.h> #include<iostream> #include<string> #include<math.h> using namespace std; class A{ public: string plus(string a,string b){ string c; int m=0; int n=0; int i=0; int x=0; ...
C++大数乘法
Vioshine的博客
12-19 869
C++大数乘法
大整数乘法—C++实现
ksgwmzgd的博客
04-05 820
大整数乘法—C++实现 自己通过用二维数组表来改写此算法。 引用博客文章链接：https://blog.csdn.net/qq_44731019/article/details/109367811?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522164913780516782184653108%2522%252C%2522scm%2522%253A%252220140713.130102334…%2522%257D&request_id
C++大数乘法程序
最新发布
OverTheSunn的博客
01-28 187
【代码】C++大数乘法程序。
C++大整数运算
热门推荐
kedixa 的专栏
02-22 1万+
C++大整数运算系列文章讲述了实现大整数四则运算的基本思路和方法，对算法做简单的介绍并采用C++代码做进一步说明。本文并不试图写成一篇教程，而是对实现过程进行整理和概括。
C++实现大数乘法算法代码
09-03
在C++编程中，处理大数乘法是一个常见的任务，特别是在计算、加密和算法实现等领域。这段代码展示了如何使用C++实现一个简单的手动大数乘法算法，它基于学校教的列式乘法规则（也称为长乘法）。以下是这个算法的详细...
C++ 大数乘法 一般算法
05-03
大数乘法的一般算法！ 这个是关于C语言中大数乘法的一般算法。我愿意共享这个资源。。
大数乘法c++
11-05
大数乘法的c++编程实现，输入两个大数，每行输入一个，就可以得出结构
C++使用string的大数乘法运算（3）
08-25
C++使用string的大数乘法运算 本文主要介绍了使用C++中的string实现大数乘法运算的方法，并提供了详细的代码实现和解释。 一、大数乘法的概念 大数乘法是指对两个大于long long int的数字进行乘法运算。由于C++中...
double相乘
秋秋的博客
06-06 9693
1.报错程序 public class Multiply { public static Double multiply(Double a, Double b) { return a * b } } 2.分析造成问题的原因 直接执行报错 3.解决方案 public class Multiply { public static Double m...
大整数乘法全解（绝对可运行C++）
11-08
大整数乘法全解（绝对可运行C++），源代码，加说明！
大数浮点数幂运算(c++实现)
04-22
自己写得大数浮点数幂运算(c++实现)，系poj acm 的problem:1001的实现
C++大整数相乘的一个简单方法
weixin_43700556的博客
04-02 939
就是模拟算数乘法的一个方法 先由第一个数字的第一个数分别乘以第二个数的从个位数开始的个位数，进行适当的进位就行了 #include<iostream> #include<string> #include<algorithm> using namespace std; int main() { string a; getline(cin,a); ...
超大数字乘法问题（c++）
keyboard not found! Please press F1 to continue.
02-24 1166
问题描述：    用C++实现2个大数相乘。算法思路：    用指定len的2个整型数组代表2个大数，设乘数为A，被乘数为B，保存结果的数组为buffer，可得buffer的元素总数为A.len + B.len。逆序循环遍历A（从个位到首位），在每次循环体中，取A的一位数与B的每一位数相乘，并且每次乘积都保存到buffer指定位置上（当前拿的A的哪一位，结果的保存就从buffer的哪一位开始保存）...
C语言double数据加减,double类型数据的加减乘除
weixin_42297562的博客
05-23 2876
import java.math.BigDecimal;/*** 由于Java的简单类型不能够精确的对浮点数进行运算，这个工具类提供精 确的浮点数运算，包括加减乘除和四舍五入。*/public class DoubleUtils{private static final int DEF_DIV_SCALE = 2;/*** @Description 两个Double数相加** @param d1*...
double类型的乘法问题
weixin_30823683的博客
05-21 4917
遇到了一个问题,两个double相乘,结果不正确,好长,很奇怪 参考了BigDecimal类的加减乘除(解决double计算精度问题)和Java double 相乘的小问题记录 终于解决了问题,下面提供一下代码 1 public static void main(String[] args){ 2 BigDecimal a1 = new BigDecimal(454.5...
大数的加减乘除运算(C++)
xiayang1023的专栏
11-21 4496
自认为算法比较渣，最近经常遇到一些大数问题，虽然原理
帮我生成 大数乘法 的代码c++
06-06
以下是一个简单的实现大数乘法的 C++ 代码： ```c++ #include #include #include #include using namespace std; vector<int> multiply(vector<int>& num1, vector<int>& num2) { int n1 = num1.size(), n2 ...

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值