高精度乘法

#include<stdio.h>
#include<string.h>
const int maxn=1000;
int mul(char a[],char b[],int c[])
{
    int la=strlen(a);
    int lb=strlen(b);
    int lc=la+lb-1;
    int ta[maxn],tb[maxn];
    for(int i=0;i<la;++i)
      ta[i]=a[la-i-1]-'0';
    for(int i=0;i<lb;++i)
      tb[i]=b[lb-i-1]-'0';
    for(int i=0;i<=lc;++i)
      c[i]=0;
    for(int i=0;i<la;++i)
      for(int j=0;j<lb;++j)
        c[i+j]=c[i+j]+ta[i]*tb[j];
    int cnt=0;
    for(int i=0;i<=lc;++i)
    {
      int s=cnt+c[i];
      c[i]=s%10;
      cnt=s/10;
    }
    while(lc>0&&!c[lc])
    lc--;
    return lc;
}
int main()
{
    char a[maxn],b[maxn];
    int c[maxn*2];
    while(~scanf("%s %s",a,b))
    {
    int lc=mul(a,b,c);
    for(int i=lc;i>=0;--i)
      printf("%d",c[i]);
    puts("");
    }
}