高精度除法的实现

高精度除法与高精度加法的定义、前置过程都是大致相同的，如果想了解具体内容，可以移步至我的这篇博客：高精度加法计算的实现

在这里就不再详细讲解，只讲解主体过程qwq

主体过程

高精度除法的原理和小学学习的竖式除法是一样的。

概括来说，假如被除数长度为，除数长度为，为了减少冗余运算，我们从从后往前开始计算，将被除数与除数相对应的每一位相（整）除，实际上这一步可以看作一个逐次减法的过程，然后存进商的对应位置上，再将余数乘并放进下一位。

用高精度计算，先除百位，将减去一次后变为，小于，所以将存入百位，将存入十位；

再除十位，将减去三次后变为，小于，所以将存入十位，将存入个位；

最后除个位，将减去八次后变为，小于，所以将存入个位，将存入余数数组。

其实，高精度除法按理来说不需要反转存储，正序存储会更方便，但大部分题目，如果需要高精度除法去做，那么很有可能也需要其他的高精度计算，为了统一，我们还是使用反转存储。

接下来，我们这里实现一个函数，它判断了被除数以下标为最低位，是否可以再减去除数而保持非负。这个函数分为三部分：

被除数剩余的部分比除数长，这个情况下最多多出 1 位，函数返回真。
如第一步判断为假，就说明被除数与除数一样长，那我们就从高位到低位，逐位比较：如果被除数当前位比除数当前位大，函数返回真；反之，函数返回假。
如第二步也判断为假，就说明被除数与除数相等，相等的情形下也是可行的，函数返回真。

下面给出高精度除法的代码：

cpp 复制代码

bool big(int a[],int b[],int low,int L){
	if(a[low+L]!=0) return 1;
	for(int i=L-1;i>=0;--i){
		if(a[low+i]>b[i]) 
			return 1;
		if(a[low+i]<b[i])
			return 0;
	}
	return 1;
}
void div(int a[],int b[],int c[],int d[]){
	clear(c);
	clear(d);
	int la,lb;
	for(la=L-1;la>0;la--){
		if(a[la-1]!=0)
			break;
	}
	for(lb=L-1;lb>0;lb--){
		if(b[lb-1]!=0)
			break;
	}
	if(lb==0) return;
	for(int i=0;i<la;i++) d[i]=a[i];
	for(int i=la-lb;i>=0;i--){
		while(big(d,b,i,lb)){
			for(int j=0;j<lb;j++){
				d[i+j]-=b[j];
				if(d[i+j]<0){
					d[i+j+1]-=1;
					d[i+j]+=10;
				}
			}
			c[i]++;
		}
	}
}

高精度计算器（总结）

到这里，我们的高精度计算就全部完成了。

下面给出高精度计算器的代码：

cpp 复制代码

const int L=10000;
string s;
int a[L],b[L],c[L],d[L];
void clear(int a[]){
	for(int i=0;i<L;i++)
		a[i]=0;
}
void read(int a[]){
	cin>>s;
	int L=s.size();
	for(int i=0;i<L;i++)
		a[i]=s[L-1-i]-'0';
}
void print(int a[]){
	int i;
	for(i=L-1;i>=1;i--){
		if(a[i]!=0)
			break;
	}
	for(;i>=0;i--)
		cout<<a[i];
	cout<<endl;
}
void add(int a[],int b[],int c[]){
	clear(c);
	for(int i=0;i<L-1;++i){
		c[i]+=a[i]+b[i];
		if(c[i]>=10){
			c[i+1]+=1;
			c[i]-=10;
		}
	}
}
void sub(int a[],int b[],int c[]){
	clear(c);
	for(int i=0;i<L-1;++i){
		c[i]+=a[i]-b[i];
		if(c[i]<0){
			c[i+1]-=1;
			c[i]+=10;
		}
	}
}
void mul(int a[],int b[],int c[]){
	clear(c);
	for(int i=0;i<L-1;i++){
		for(int j=0;j<=i;j++)
			c[i]+=a[j]*b[i-j];
		if(c[i]>=10){
			c[i+1]+=c[i]/10;
			c[i]%=10;
		}
	}
}
bool big(int a[],int b[],int low,int L){
	if(a[low+L]!=0) return 1;
	for(int i=L-1;i>=0;--i){
		if(a[low+i]>b[i]) 
			return 1;
		if(a[low+i]<b[i])
			return 0;
	}
	return 1;
}
void div(int a[],int b[],int c[],int d[]){
	clear(c);
	clear(d);
	int la,lb;
	for(la=L-1;la>0;la--){
		if(a[la-1]!=0)
			break;
	}
	for(lb=L-1;lb>0;lb--){
		if(b[lb-1]!=0)
			break;
	}
	if(lb==0) return;
	for(int i=0;i<la;i++) d[i]=a[i];
	for(int i=la-lb;i>=0;i--){
		while(big(d,b,i,lb)){
			for(int j=0;j<lb;j++){
				d[i+j]-=b[j];
				if(d[i+j]<0){
					d[i+j+1]-=1;
					d[i+j]+=10;
				}
			}
			c[i]++;
		}
	}
}

每周六更新一篇文章，内容一般是自己总结的经验或是在其他网站上整理的优质内容

点个赞，关注一下呗~