高精度相乘（大数乘大数）

原理：倒置相乘，然后统一处理进位，再还原

时间复杂度：O（n2）

//C
#include<stdio.h>
#include<string.h>
void mul(char a[],char b[])
{
	int i,j,la,lb,lc;
	int na[100]={0},nb[100]={0},nc[100]={0};
	la = strlen(a);
	lb = strlen(b);
	for(i = 0; i < la; i ++)
		na[la-i] = a[i] - '0';
	for(i = 0; i < lb; i ++)
		nb[lb-i] = b[i] - '0';
	for(i = 1; i <= la; i ++)
		for(j = 1; j <= lb; j ++)
			nc[i+j-1] += na[i]*nb[j];
	for(i = 1; i <= la+lb; i ++)
	{
		nc[i+1] += nc[i]/10;
		nc[i] %= 10;
	}
	if(nc[la+lb])
		printf("%d",nc[la+lb]);
	for(i = la+lb-1; i >= 1; i --)
		printf("%d",nc[i]);
}
int main()
{
	char s1[100],s2[100];
	while(scanf("%s %s",s1,s2) != EOF)
	{
		mul(s1,s2);
		printf("\n");
	}
	return 0;
}

//C++ 
#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';
    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;
}