我的高精度乘法模板

来源:互联网 发布:淘宝2016全年交易额 编辑:程序博客网 时间:2024/06/07 12:45

O(n^2)复杂度的算法,以后会发布O(n*log(n))的基于快速傅里叶变换的高精乘

用字符串读取数据

倒序存储

模拟手工运算

最好自己用个比较小的数字按算法走一遍,然后就一目了然了。

#include<cstdio>#include<cstring>#include<cmath>#include<cstdlib>#include<iostream>#include<algorithm>#include<sstream>#include<vector>#include<map>#include<stack>#include<list>#include<set>#include<queue>#define LL long long#define lson l,m,rt<<1#define rson m+1,r,rt<<1 | 1using namespace std;const int L=1005;int n,m;string mul(string a,string b){    string s;    int na[L],nb[L],nc[L],La=a.size(),Lb=b.size();//na存储被乘数,nb存储乘数,nc存储积    fill(na,na+L,0);fill(nb,nb+L,0);fill(nc,nc+L,0);//将na,nb,nc都置为0    for(int i=La-1;i>=0;i--) na[La-i]=a[i]-'0';//将字符串表示的大整形数转成i整形数组表示的大整形数    for(int i=Lb-1;i>=0;i--) nb[Lb-i]=b[i]-'0';    for(int i=1;i<=La;i++)        for(int j=1;j<=Lb;j++)        nc[i+j-1]+=na[i]*nb[j];//a的第i位乘以b的第j位为积的第i+j-1位(先不考虑进位)    for(int i=1;i<=La+Lb;i++)        nc[i+1]+=nc[i]/10,nc[i]%=10;//统一处理进位    if(nc[La+Lb]) s+=nc[La+Lb]+'0';//判断第i+j位上的数字是不是0    for(int i=La+Lb-1;i>=1;i--)        s+=nc[i]+'0';//将整形数组转成字符串    return s;}int main(){    string a,b;    while(cin>>a>>b) cout<<mul(a,b)<<endl;    return 0;}


0 0
原创粉丝点击