【算法】高精度算法（加减乘除）

为啥有高精度算法

A+B难道不是直接cout<<a+b;吗
不不不，如果你遇到了爆炸范围，就不能用int或者long long了

高精度算法-加法

高精度算法就是开数组，一个一个位数加。

算法核心：

c[i]+=a[i]+b[i]; // 和小学加法一样，按位相加。
c[i+1]=c[i]/10;  // 保存进位。例如，如果c[i]是12，那么要进1，c[i+1]这个时候就是进位1。
c[i]=c[i]%10; // 进位之后剩下的数字。c[i]原来是12，但是要进位，所以c[i]上留下的数字只有2。

计算A+Bproblem

#include <bits/stdc++.h>
using namespace std;
char s1[505],s2[505];
int a[505],b[505],c[505]; 
int main()
{
	int la,lb,lc;
	int i;
	cin>>s1>>s2;
	
	la=strlen(s1); // 求s1的长度
	lb=strlen(s2); // 求s2的长度
	 
	for(i=0;i<la;i++)
	{
		a[la-i]=s1[i]-'0';
	}
	for(i=0;i<lb;i++)
	{
		b[lb-i]=s2[i]-'0';
	}
	// 将字符转化为数字，并且将字符转置便于计算
	
	lc=max(la,lb)+1;
	// 加法结果的长度应该是两个数中最长的那个再加一
	
	for(i=1;i<=lc;i++)
	{
		c[i]+=a[i]+b[i];
		c[i+1]=c[i]/10;
		c[i]=c[i]%10;
	} 
	if(c[lc]==0&&lc>0) lc--; // 删除前导0
	for(i=lc;i>0;i--)
	{
		cout<<c[i];
	}
	return 0;
}

易错点
（1）前导0的意思是：如果计算出012这样的东西，要删除0。

if(c[lc]==0&&lc>0) lc--;

为什么要加上lc>0：如果你输入0+0，那么计算结果是0。如果不加上lc>0，那么程序会把结果的0也删掉。

（2）将字符转化为数字：

s1[i]-'0'

将字符转置便于计算：我们可以知道，例如输入一个数17839，那么一开始1是s[0]，7是s[1]，以此类推。但是为了方便计算（因为高精度算法按位相加是从最后一个位开始加的），我们需要把最后面一个数变成第一位，即s[0]，然后倒数第二位变成s[1]，以此类推。

for(i=0;i<la;i++)
{
	a[la-i]=s1[i]-'0';
}

高精度算法-减法

同理，开数组。

算法核心：
（1）如果 a<b ，那么ab交换（这是为了方便计算，只要输出结果的时候加个负号就好了）

（2）如果 a[i]<b[i] ， 需要高位借一当十用：
a[i]<b[i]，那么向a[i+1]借一，a[i+1]就需要减去一（因为被借掉了），然后a[i]这个位就可以加10，计算出来的结果c[i]就是a[i]-b[i]

if(a[i]<b[i])
{
	a[i+1]--;
	a[i]+=10;
} 
c[i]=a[i]-b[i];

计算A-Bproblem

#include <bits/stdc++.h>
using namespace std;
char s1[10090],s2[10090],s3[10090];
int a[10090],b[10090],c[10090];
bool judge(char s1[],char s2[])
{
	int u=strlen(s1),v=strlen(s2);
	if(u!=v) return u>v;
	/*
	如果u>v,那么肯定有s1>s2，返回ture
	如果u==v,那么就进行下面的比较
	*/
	for(int i=0;i<u;i++)
	{
		if(s1[i]!=s2[i]) return s1[i]>s2[i];
	} 
	return true;
} 
int main()
{
	int flag,la,lb,lc;
	int i;
	cin>>s1>>s2;
	
	if(!judge(s1,s2)) // 比较s1和s2的大小
	{
		flag=1;
		// 交换s1和s2，需要一个中介s3，这里因为都是字符数组，所以要用strcpy（）函数
		strcpy(s3,s1); 
		strcpy(s1,s2);
		strcpy(s2,s3);
	}
	la=strlen(s1); // 求s1的长度
	lb=strlen(s2); // 求s2的长度
	 
	for(i=0;i<la;i++)
	{
		a[la-i]=s1[i]-'0';
	}
	for(i=0;i<lb;i++)
	{
		b[lb-i]=s2[i]-'0';
	}
	// 将字符转化为数字，并且将字符转置便于计算
	
	lc=max(la,lb);
	// 减法结果的长度应该是两个数中最长的那个
	
	for(i=1;i<=lc;i++)
	{
		if(a[i]<b[i])
		{
			a[i+1]--;
			a[i]+=10;
		} 
		c[i]=a[i]-b[i];
	} 
	while(c[lc]==0&&lc>1) lc--;
	// 这里一定不能用if！不然会有一个测试点过不了我也不知道为啥
	if(flag==1) cout<<"-";
	for(i=lc;i>0;i--)
	{
		cout<<c[i];
	}
	return 0;
}   	

高精度算法-乘法

核心代码：

for(i=1;i<=la;i++)
{
	for(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]%=10;
	}
} 

A*B问题：

#include <bits/stdc++.h>
using namespace std;
char s1[2001],s2[2001],s3[2001];
int a[4000000],b[4000000],c[4000000];
int main()
{
	int la,lb,lc;
	int i,j;
	cin>>s1>>s2;
	
	la=strlen(s1); // 求s1的长度
	lb=strlen(s2); // 求s2的长度
	 
	for(i=0;i<la;i++)
	{
		a[la-i]=s1[i]-'0';
	}
	for(i=0;i<lb;i++)
	{
		b[lb-i]=s2[i]-'0';
	}
	// 将字符转化为数字，并且将字符转置便于计算
	
	lc=la+lb;
	// 减法结果的长度应该是两个数中最长的那个
	
	for(i=1;i<=la;i++)
	{
		for(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]%=10;
		}
	} 
	while(c[lc]==0&&lc>1) lc--;
	for(i=lc;i>0;i--)
	{
		cout<<c[i];
	}
	return 0;
}
   	

高精度算法-除法

高精度除以低精度

核心代码：

for(i=1;i<=la;i++)
{
	c[i]=(x*10+a[i])/b;
	x=(x*10+a[i])%b;
}

A/B problem：

#include <bits/stdc++.h>
using namespace std;
char s1[5005]; 
long long a[5005],c[5005],b,x;
int main()
{
	int la,lb,lc;
	int i;
	cin>>s1>>b;
	
	la=strlen(s1); // 求s1的长度
	for(i=1;i<=la;i++)
	{
		a[i]=s1[i-1]-'0'; // 不用转置，被除数一位位放入a数组 
	}
	for(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(i=lc;i<=la;i++)
	{
		cout<<c[i];
	}
	return 0;
}
   

高精度除以高精度

完整代码：

#include <bits/stdc++.h>
using namespace std;
char s1[305],s2[305]; 
int a[305],b[305],c[305],tmp[305];
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;
}

void numcpy(int p[],int q[],int n)//将p数组整体移动n-1位到q数组中 
{
	for(int i=1;i<=p[0];i++) q[i+n-1]=p[i];
	q[0]=p[0]+n-1; 
} 

int compare(int a[],int b[])
{
	if(a[0]>b[0]) return 1;//a的位数大于b的位数则a比b大
	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]=a[i]+10;
		}
		a[i]=a[i]-b[i];
	}
	while(a[a[0]]==0&&a[0]>0) a[0]--;//删除a的前导0，修正a的位数 
} 

int main()
{
	init(a);
	init(b);
	c[0]=a[0]-b[0]+1;
	// c[0]放商的位数，a[0]放被除数的位数，b[0] 放除数的位数 
	
	for(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]>0) c[0]--;
	print(c);//商 
	cout<<endl;
	print(a);//余数 
	return 0;
}