C++算法之高精度计算

互联网的猫

已于 2024-02-02 17:24:22 修改

阅读量2.5k

点赞数 41

分类专栏： C++STL库及细节用法详解高精度算法专题文章标签： c++ 算法

于 2023-12-14 17:57:24 首次发布

本文链接：https://blog.csdn.net/m0_72674633/article/details/134991719

版权

C++STL库及细节用法详解同时被 2 个专栏收录

20 篇文章 2 订阅

订阅专栏

高精度算法专题

7 篇文章 0 订阅

订阅专栏

高精度算法介绍

在C/C++中，我们经常会遇到限定数据范围的情况，我们先来看一下常用的int和long long两种数据类型的范围。

C++标准规定：

int占一个机器字长，在32位系统中int占32位，即4个字节，所以int的范围是[- $2^{31}$ , $2^{31}-1$ ],为 $10^{9}$ 数量级。long long的范围则是[ $-2^{63}$ , $2^{63}-1$ ],为 $10^{18}$ 数量级。如果超过该数量级，就需要用到高精度算法。

高精度算法应用

高精度加法

1.输入两个整数a,b,输出它们的和(a,b<= $10^{9}$ ）

直接使用int/long long进行设变量计算。

2.输入两个整数a,b，输入它们的和(a,b<= $10^{500}$ ）

分析:int: $-2^{31}$ ~ $2^{31}-1$ ,数量级为 $10^{9}$ ，long long: $-2^{63}$ ~ $2^{63}-1$ ,数量级为 $10^{18}$ ，解决方案：用数组模拟高精度。

算法核心：c[i]=a[i]+b[i],c[i+1]=c[i]/10,c[i]=c[i]%10;

总结下来就是对于位数较大的整数的加法，使用数组来存储整数，数组的长度等于整数的位数，按照字符串的形式输入，先统计长度，然后倒置过来利用数组进位来模拟加法，最终将数组最大位置的前导0判断去除，倒叙输出即为最终的结果。

#include<iostream>
#include<cstring>
#include<cstdlib>
using namespace std;
char s1[505],s2[505];
int a[505],b[505],c[505];
int main(){
	int la,lb,lc;
	cin>>s1>>s2;//用数组模拟加法，数组输入
	la=strlen(s1);
	lb=strlen(s2);
	for(int i=0;i<la;i++){//倒置过来
		a[la-i]=s1[i]-'0';
	}
	for(int i=0;i<lb;i++){
		b[lb-i]=s2[i]-'0';
	}
	lc=max(la,lb)+1;
	for(int i=1;i<=lc;i++){//使用数组进位模拟加法
		c[i]+=a[i]+b[i];
		c[i+1]+=c[i]/10;
		c[i]%=10;
	}
	if(c[lc]==0&&lc>1){//去掉前导0
		lc--;
	}
	for(int i=lc;i>0;i--){//倒序输出
		cout<<c[i];
	}
	return 0;
}

高精度减法

要点：1.如果a<b，则需要交换a与b，最后计算完以后添加一个负号即可。2.如果a[i]<b[i]，需要高位借一当十使用。

总结来说，就是在加法的基础上最开始加上了一个比较的操作，将大的数做除数，同时记录是否交换，解决了正负号的问题，其次使用数组模拟减法，从高位借1的操作。

#include<iostream>
#include<cstring>
#include<algorithm>
#include<string>
#define maxn 10100
using namespace std;
int a[maxn],b[maxn],c[maxn],flag;
string m,n,s;
bool compare(string s1,string s2){
	if(s1.length()!=s2.length()){
		return s1.length()>s2.length();
	}
	for(int i=s1.length()-1;i>=0;i--){
		if(s1[i]!=s2[i]){
			return s1[i]>s2[i];
		}
	}
	return true;
}
int main(){
	cin>>m>>n;
	if(!compare(m,n)){
		flag=1;
		s=m;
		m=n;
		n=s;
	}
	int len=max(m.length(),n.length());
	for(int i=m.length()-1,j=1;i>=0;i--,j++){
		a[j]=m[i]-'0';
	}
	for(int i=n.length()-1,j=1;i>=0;i--,j++){
		b[j]=n[i]-'0';
	}
	for(int i=1;i<=len;i++){
		if(a[i]<b[i]){
			a[i+1]--;
			a[i]+=10;
		}
		c[i]=a[i]-b[i];
	}
	while(c[len]==0&&len>1){
		len--;
	}
	if(flag==1){
		cout<<"-";
	}
	for(int i=len;i>0;i--){
		cout<<c[i];
	}
	return 0;
}

高精度乘法

通过图示分析，我们发现规律，a[1]*b[1]对应c[1]位置，a[2]*b[1]对应c[2]位置，a[3]*b[1]对应c[3]位置，a[4]*b[1]对应c[4]位置，a[1]*b[2]对应c[2]，a[2]*b[2]对应c[3]位置，a[3]*b[2]对应c[4]位置，a[4]*b[2]对应c[5]位置。

容易得出a[i]*b[j]对应c[i+j-1]位置，然后考虑进位，我们可以得到c[i+j-1]的关系式：c[i+j-1]+=a[i]*b[j];c[i+j]=c[i+j-1]/10;c[i+j-1]%=10;

总结来说除了运算的操作采用上面的关系式，其他的操作和加法的是相似的。

#include<iostream>
#include<cstring>
using namespace std;
char s1[2005],s2[2005];
int a[2005],b[2005],c[4005];
int main(){
	int la,lb,lc;
	cin>>s1>>s2;
	la=strlen(s1);
	lb=strlen(s2);
	for(int i=0;i<la;i++){
		a[la-i]=s1[i]-'0';
	} 
	for(int i=0;i<lb;i++){
		b[lb-i]=s2[i]-'0';
	} 
	lc=la+lb;
	for(int i=1;i<=la;i++){
		for(int j=1;j<=lb;j++){
			c[i+j-1]+=a[i]*b[j];
			c[i+j]+=c[i+j-1]/10;
			c[i+j-1]=c[i+j-1]%=10;
		}
	}
	while(c[lc]==0&&lc>1){
		lc--;
	}
	for(int i=lc;i>0;i--){
		cout<<c[i];
	}
	return 0;
}

高精度除法

高精度除以低精度

也就是大数除以小数，使用逐位试商法（也就是不断求商和求余数的过程）

#include<iostream>
#include<cstring>
using namespace std;
char s1[5005];//存储大数 
long long b;//存储小数
long long a[5005],c[5005],la,lc,x; 
int main(){
	cin>>s1>>b;
	la=strlen(s1);//被除数的长度 
	for(int i=1;i<=la;i++){
		a[i]=s1[i-1]-'0';
	}
	for(int i=1;i<=la;i++){
		c[i]=(x*10+a[i])/b;//求商 
		x=(x*10+a[i])%b;//求余数 
	}
	lc=1;
	while(c[lc]==0&&lc<la){
		lc++;
	} 
	for(int i=lc;i<=la;i++){
		cout<<c[i];
	} 
	return 0;
}

高精度除以高精度

方法：减法模拟除法

#include<iostream>
#include<cstring>
using namespace std;
char s1[305],s2[305];
int a[305],b[305],c[305],tmp[305];//a为被除数，b为除数,c为商
void init(int *x){
	char s[305];
	cin>>s;
	x[0]=strlen(s);
	for(int i=0;i<x[0];i++){
		x[x[0]-i]=s[i]-'0';//将字符串转化为数字，并且倒叙存储 
	}
} 
void print(int a[]){
	if(a[0]==0){
		cout<<0;
		return;
	}
	for(int i=a[0];i>0;i--){
		cout<<a[i];
	}
	return;
}
int compare(int a[],int b[]){//返回1表示a>b,返回0表示a=b，返回-1表示a<b 
 
	if(a[0]>b[0]){
		return 1;
	}
	if(a[0]<b[0]){
		return -1;
	}
	for(int i=a[0];i>0;i--){
		if(a[i]>b[i]){
			return 1;
		}
		if(a[i]<b[i]){
			return -1;
		}
	}
	return 0;
}
void minu(int a[],int b[]){
	for(int i=1;i<=a[0];i++){
		if(a[i]<b[i]){
			a[i+1]--;
			a[i]+=10;
		}
		a[i]=a[i]-b[i];
	}
	while(a[a[0]]==0&&a[0]>0){
		a[0]--;
	}//删除a的前导0，修正a的位数 
}
void numcpy(int p[],int q[],int n){//将除数移动相应的位数 
	for(int i=1;i<=p[0];i++){
		q[i+n-1]=p[i];
	}
	q[0]=p[0]+n-1;//移动后数组的位数 
}
int main(){
	init(a);//输入a 
	init(b);//输入b
	c[0]=a[0]-b[0]+1;
	for(int i=c[0];i>=1;i--){
		memset(tmp,0,sizeof(tmp));
		numcpy(b,tmp,i);//除数移动相应的位数 
		while(compare(a,tmp)>=0){//当被除数大于除数时执行循环 
			c[i]++;
			minu(a,tmp);
		}
	} 
	while(c[c[0]]==0&&c[0]>1){
		c[0]--;
	}
	print(c);//商 
	cout<<endl;
	print(a);//余数 
	return 0;
}