高精度乘法

#include 

#include 

#include 

#include 

 

#define MAX 10001

 

int bigchenfa(int *sum,int *a,int *b,int lsum,int la,int lb)

{

    int i,j;

    memset(sum,0,sizeof(sum));

    lsum = 0 ;

    for(i=1 ; i<= la ; i++) 

        for(j=1,lsum=i-1; j<= lb ; j++)

            sum[++lsum] += b[j] * a[i] ;

 

    for(i=1 ; i<= lsum ; i++)

    if (sum[i] >= 10)

    {

        if ( sum[lsum] >= 10)

            lsum ++ ;

        sum[i+1] += sum[i] / 10 ;

        sum[i] %= 10 ;

    }

 

    return lsum ;

}

int main(void)

{

 

    int la=0,lb=0,lsum=0;

    int i,j ;

    char sa[MAX],sb[MAX] ;

    while(2==scanf("%s%s",sa,sb))

    {

        int a[MAX]={0},b[MAX]={0},sum[MAX*2]={0} ;

        la = strlen(sa);

        lb = strlen(sb);

 

        for(i=1,j=la-1; i<= la ; i++,j--)

            a[i] = sa[j] - '0' ;

        for(i=1,j=lb-1; i<= lb ; i++,j--)

            b[i] = sb[j] - '0' ;

 

        lsum = bigchenfa(sum,a,b,lsum,la,lb) ;

 

        for(i=lsum ; i>= 1 ; i--)

            printf("%d",sum[i]);

 

        printf("\n");

    }

    return 0 ;

}