高精度乘法

高精度乘法
1）高精度乘高精度的朴素算法

传入参数约定：传入参数均为string类型，返回值为string类型
算法思想：倒置相乘，然后统一处理进位，再还原。
算法复杂度：o(n^2).

#include<iostream>
#include<cstring>
#include<algorithm>
using namespace std;
const int L=110;
string mul(string a,string b)//高精度乘法a,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;
}