[C++]高精度算法

其他 2018-11-25 09:50:58 阅读次数: 0
版权声明：Copyright © 2018-2018 Takaway 复习用，转载注明一下就行了~ https://blog.csdn.net/LynlinBoy/article/details/84146881

目录

高精度加法    高精度减法    高精度乘法    高精度除法    高精度阶乘

高精度加法

用程序来模拟竖式加法即可，注意在输出的时候除去多余的前导零。

程序代码

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int MAX_SIZE = 3333;
inline void High_precision_Addition(int *a, int *b, int *c) {
	c[0] = 1; int x = 0;
	while (c[0] <= a[0] || c[0] <= b[0]) {
		c[c[0]] = a[c[0]] + b[c[0]] + x;
		x = c[c[0]] / 10;
		c[c[0]] %= 10;
		++c[0];	
	};
	c[c[0]] = x;
	while (c[c[0]] == 0 && c[0] > 1) --c[0];	
};
int main(int argc, char *argv[]) {
	char str1[MAX_SIZE], str2[MAX_SIZE];
	cin >> str1 >> str2;
	int a[MAX_SIZE], b[MAX_SIZE], c[MAX_SIZE];
	memset(a, 0, sizeof a); memset(b, 0, sizeof b); memset(c, 0, sizeof c);
	a[0] = strlen(str1); b[0] = strlen(str2);
	for (int i = 0; i < a[0]; ++i) a[a[0] - i] = str1[i] - '0'; 
	for (int i = 0; i < b[0]; ++i) b[b[0] - i] = str2[i] - '0';
	High_precision_Addition(a, b, c);
	for (int i = c[0]; i >= 1; --i) printf("%d", c[i]); cout.put('\n'); //和
	return 0;	
};

高精度减法

用程序来模拟竖式减法即可，与高精度加法类似，不过得考虑被减数与减数的大小关系还得考虑借位，所以代码相对于高精度加法更为复杂。

程序代码

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int MAX_SIZE = 3333;
inline bool Is_swap(char *__x_, char *__y_) {
	char Mid_Element[MAX_SIZE];
	if ((strlen(__x_) < strlen(__y_)) || (strlen(__x_) == strlen(__y_) && strcmp(__x_, __y_) < 0)) {
		strcpy(Mid_Element, __x_);
		strcpy(__x_, __y_);
		strcpy(__y_, Mid_Element);
		return true;
	};	
	return false;
};
inline void High_precision_Subtraction(int *a, int *b, int *c) {
	 c[0] = 1;
	while (c[0] <= a[0] || c[0] <= b[0]) {
		if (a[c[0]] < b[c[0]]) {
			a[c[0]] += 10;
			--a[c[0] + 1];	
		};
		c[c[0]]	= a[c[0]] - b[c[0]];
		++c[0];
	};
	while (c[c[0]] == 0 && c[0] > 1) --c[0];
};
int main(int argc, char *argv[]) {
	char str1[MAX_SIZE], str2[MAX_SIZE];
	cin >> str1 >> str2; bool flag = Is_swap(str1, str2);
	int a[MAX_SIZE], b[MAX_SIZE], c[MAX_SIZE];
	memset(a, 0, sizeof a); memset(b, 0, sizeof b); memset(c, 0, sizeof c);
	a[0] = strlen(str1); b[0] = strlen(str2);
	for (int i = 0; i < a[0]; ++i) a[a[0] - i] = str1[i] - '0'; 
	for (int i = 0; i < b[0]; ++i) b[b[0] - i] = str2[i] - '0';
	High_precision_Subtraction(a, b, c);
	if (flag == true) cout.put('-');
	for (int i = c[0]; i >= 1; --i) printf("%d", c[i]); cout.put('\n'); //差
	return 0;	
};

高精度乘法

用程序模拟竖式乘法即可，在做乘法运算时，同样也有进位，同时对每一位进行乘法运算时，必须进行错位相加。

程序代码

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int MAX_SIZE = 3333;
inline void High_precision_Multiplication(int *a, int *b, int *c) {
	int x;
	for (int i = 1; i <= a[0]; ++i) {
		x = 0;
		for (int j = 1; j <= b[0]; ++j) {
			c[i + j - 1] = a[i] * b[j] + x + c[i + j - 1];
			x = c[i + j - 1] / 10;
			c[i + j - 1] %= 10;	
		};
		c[i + b[0]] = x;	
	};
	c[0] = a[0] + b[0];
	while (c[c[0]] == 0 && c[0] > 1) --c[0];
};
int main(int argc, char *argv[]) {
	char str1[MAX_SIZE], str2[MAX_SIZE];
	cin >> str1 >> str2;
	int a[MAX_SIZE], b[MAX_SIZE], c[MAX_SIZE];
	memset(a, 0, sizeof a); memset(b, 0, sizeof b); memset(c, 0, sizeof c);
	a[0] = strlen(str1); b[0] = strlen(str2); 
	for (int i = 0; i < a[0]; ++i) a[a[0] - i] = str1[i] - '0'; 
	for (int i = 0; i < b[0]; ++i) b[b[0] - i] = str2[i] - '0';
	High_precision_Multiplication(a, b, c);
	for (int i = c[0]; i >= 1; --i) cout << c[i]; cout.put('\n'); //积
	return 0;
};

高精度除法

高精度除法一般有高精度除以低精度与高精度除以高精度，前者代码较为简单，后者则作用范围更加广阔。高精度除以低精度与前三种高精度算法类似，都是按位运算；而高精度除以高精度则是以减法来模拟除法。

程序代码1

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int MAX_SIZE = 3333;
inline int High_precision_Division(int *a, int b, int *c) {
	int x;
	for (int i = 1; i <= a[0]; ++i) {
		c[i] = (x * 10 + a[i]) / b;
		x = (x * 10 + a[i]) % b;	
	};
	c[0] = 1;
	while (c[c[0]] == 0 && c[0] < a[0]) ++c[0];
    return x;
};
int main(int argc, char *argv[]) {
	char str[MAX_SIZE]; int b; 
	cin >> str >> b;
	int a[MAX_SIZE], c[MAX_SIZE];
	memset(a, 0, sizeof a); memset(c, 0, sizeof c);
	a[0] = strlen(str); 
	for (int i = 0; i < a[0]; ++i) a[i + 1] = str[i] - '0'; 
	int x = High_precision_Division(a, b, c);
	for (int i = c[0]; i <= a[0]; ++i) cout << c[i]; cout.put('\n'); //商
	cout << x << endl; //余数
	return 0;	
};

程序代码2

#include <iostream>
#include <cstring>
#include <cstdio>
#include <string>
using namespace std;
const int MAX_SIZE = 3333;
inline void Input(int *radix) {
	string str; cin >> str; radix[0] = str.length();
	for (int i = 0; i < radix[0]; ++i) radix[radix[0] - i] = str[i] - '0';	
};
inline void Print(const int *radix) {
    if (radix[0] == 0) {cout << 0 << endl; return;};
	for (int i = radix[0]; i >= 1; --i) printf("%d", radix[i]);
	cout.put('\n'); return;	   	   	
};
inline int Compare(const int *__x_, const int *__y_) {
	if (__x_[0] > __y_[0]) return 1;
	if (__x_[0] < __y_[0]) return -1;
	for (int i = __x_[0]; i >= 1; --i) {
		if (__x_[i] == __y_[i]) continue;
		else if (__x_[i] > __y_[i]) return 1;
		else if (__x_[i] < __y_[i]) return -1;	
	};
	return 0;   	
};
inline void Subtraction(int *__x_, int *__y_) {
	int sign = Compare(__x_, __y_);
	if (sign == 0) {__x_[0] = 0; return;}; 
	if (sign == 1) {	
		for (int i = 1; i <= __x_[0]; ++i) {
			if (__x_[i] < __y_[i]) {
				__x_[i] += 10; --__x_[i + 1];	
			};
			__x_[i] -= __y_[i];	
		};
		while (__x_[__x_[0]] == 0 && __x_[0] > 0) --__x_[0]; return; 	
	};
};
inline void Numcpy(const int *radix, int *pending, int pos) {
	for (int i = 1; i <= radix[0]; ++i) pending[i + pos - 1] = radix[i];
	pending[0] = radix[0] + pos - 1; 
};
inline void High_precision_Division(int *a, int *b, int *c) {
	int sign = Compare(a, b);
	if (sign == -1) {c[0] = 1; c[1] = 0; return;}; 
	if (sign == 1 || sign == 0) { 
		int Temp[MAX_SIZE]; c[0] = a[0] - b[0] + 1;
		for (int i = c[0]; i >= 1; --i) {
			memset(Temp, 0, sizeof Temp); Numcpy(b, Temp, i);
			while (Compare(a, Temp) != -1) {++c[i]; Subtraction(a, Temp);};	
		};
		while (c[c[0]] == 0 && c[0] > 0) --c[0]; return;	
	};
};
int main(int argc, char *argv[]) {
	int a[MAX_SIZE], b[MAX_SIZE], c[MAX_SIZE];
	memset(a, 0, sizeof a); memset(b, 0, sizeof b); memset(c, 0, sizeof c);
	Input(a); Input(b); High_precision_Division(a, b, c);
	Print(c); Print(a); //商与余数
	return 0;
};

高精度阶乘

高精度阶乘的本质是什么？实际上就是一个高精度乘法，所以code就很好写了呀~不过还得做一个小优化来应付几千甚至几万的阶乘，要不然按位算的话就妥妥的TLE了呀。

程序代码

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
const int MAX_SIZE = 3333;
inline void High_precision_Factorial(const int N, int *c) {
	c[1] = 1; int x, temp, index = 2;
	for (int i = 2; i <= N; ++i) {
		x = 0;
		for (int j = 1; j <= index - 1; ++j) { //小优化
			temp = c[j] * i + x;
			c[j] = temp % 10;
			x = temp / 10;
		};	
		while (x != 0) {
			c[index++] = x % 10;
			x /= 10;	
		};
	};
	for (int i = MAX_SIZE; i >= 1; --i) if (c[i] != 0) {c[0] = i; break;}; 
};
int main(int argc, char *argv[]) {
	int N; cin >> N;
	int c[MAX_SIZE]; memset(c, 0, sizeof c);
	High_precision_Factorial(N, c);
	for (int i = c[0]; i >= 1; --i) printf("%d", c[i]); cout.put('\n');
	return 0;	
};

猜你喜欢

转载自blog.csdn.net/LynlinBoy/article/details/84146881
今日推荐
面壁智能发布 Eurux-8x22B 开源大模型 —— 堪称「理科状元」
开源日报 | 谷歌扶持鸿蒙上位；开源Rabbit R1；Docker加持的安卓手机；微软的焦虑和野心；海尔电器把开放平台关了
中国码农的“35岁魔咒”
蘭雅 CorelDRAW 插件 2024.5.1 国际劳动节版，免费下载
Arc Browser for Windows 1.0 正式 GA
90后程序员开发视频搬运软件、不到一年获利超 700 万，结局很刑！
周排行
【转】spring中对控制反转和依赖注入的理解
tms webcore 安装和使用
java程序员进阶相关书籍
SpringMVC接受请求参数、
如何保存训练好的机器学习模型
MyEclipse、Eclipse设置项目JDK的三个地方
商超行业微信小程序开发定制一般多少钱 （行业技术人员解读）
Markdown编辑器语言——30分钟入门到到精通
Linux系统下MongoDB的简单安装与基本操作
Power Strings
每日归档
更多
2024-05-07(14)
2024-05-06(40)
2024-05-05(0)
2024-05-04(7)
2024-05-03(19)
2024-05-02(0)
2024-05-01(4)
2024-04-30(1)
2024-04-29(40)
2024-04-28(0)

代码天地 © 2018 codetd.com

境外服务器   最新电视剧2022