高精度乘法同样是模拟题,可以通过模拟竖式乘法的方法进行,但是本文提供的进位方式并不完全相同,小伙伴们可以看看区别。
废话不多说,先上代码:
#include<iostream>
#include<cstring>
using namespace std;
int a[100000001],b[100000001],c[100000001];
int main()
{
string A,B;
cin>>A>>B;
int lena=A.size(),lenb=B.size(),lenc=max(lena,lenb);
for(int i=1;i<=lena;i++)
a[i]=A[lena-i]-'0';
for(int i=1;i<=lenb;i++)
b[i]=B[lenb-i]-'0';
for(int i=1;i<=lena;i++)
for(int j=1;j<=lenb;j++)
{
c[i+j-1]+=a[i]*b[j];//关键,a[i]*b[j]最终是加在答案的第i+j-1位
}
int len=lena+lenb;
for(int i=1;i<=len;i++)
{
c[i+1]+=c[i]/10;
c[i]%=10;
}//统一进行进位
while(!c[len])
{
len--;
}
for(int i=max(len,1);i>0;i--)//max能保证有输出,这样写比较简便,不用单独写出来
cout<<c[i];
cout<<endl;
return 0;
}
跟普通的竖式乘法不同的地方是在于,竖式乘法进行一次分步的数字相乘便需进行一次进位判断,而本文提供的方法是先把每位算出来再统一进位,比如12x16,个位相乘是2x6=12,c[i+j-1]=c[1]=12,而并非竖式乘法中的已经进位了的2。
这样的好处是可以优化效率,避免多次进位处理,计算机中取模的效率远低于加法跟乘法。(该方法的进位次数就是最终结果的位数-1)
表格表示
数a | 第4位 | 第3位 | 第2位 | 第1位 |
---|---|---|---|---|
数a | 3 | 2 | ||
数b | 4 | 5 | ||
a[i]*b[1] | 15 | 10 | ||
a[i]*b[2] | 12 | 8 | ||
中间产物 | 12 | 23 | 10 | |
结果 | 1 | 4 | 4 | 0 |
其中这个第i+j-1位要理解,就跟普通乘法一样。
这种方法的复杂度是O(n2),其中n是数字的位数。
但这种方法依然对于较大的数处理起来很慢,可以利用快速傅里叶变换的方式改进算法,将复杂度优化到O(nlogn)。
本人也在继续学习,共勉!