高精度乘法

#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
char a[2010],b[2010];
int x[2010],y[2010],ans[2050];
int main()
{
	cin>>a>>b;
	int len1=strlen(a),len2=strlen(b);
	for(int i=0;i<len1;i++) x[i+1]=a[len1-i-1]-'0';
	for(int i=0;i<len2;i++) y[i+1]=b[len2-i-1]-'0';
	for(int i=1;i<=len1;i++){
		for(int j=1;j<=len2;j++)
		  ans[i+j-1]+=x[i]*y[j];
	}
	int i;
	for(i=1;i<=len1+len2-1;i++)
	  ans[i+1]+=ans[i]/10,ans[i]%=10;
	while(ans[i]) ans[i+1]+=ans[i]/10,ans[i]%=10,++i;
	while(!ans[i]&&i) --i;
	if(!i) i=1; 
	for(int j=i;j;--j)
	  cout<<ans[j];
	return 0;
}