高精度计算之大数乘法C/C++语言实现

高精度计算之大数乘法C/C++语言实现

问题分析

当我们平时做乘法时,遵循的其实是从个位起,乘以另一个数的数位,等到Result,对Result整除10得到进位数res,除于10得到当前位保留数C[i].
下一位的Result则等于新的两个个位数相乘+进位数res+C[i],再重复上一过程.
因此只需将输入的两个大数倒序,再实现上述算法过程即可.

C语言实现

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#define n 201

unsigned int a[n],b[n],c[n*n]={0},lena,lenb;
char reversea[n],reverseb[n];
int t;

void Reverse()
{
	int i;
	for (i=1;i<=lena;i++) a[i]=reversea[lena-i]-48;
	for (i=1;i<=lena;i++) b[i]=reverseb[lenb-i]-48;
	return;
}

void Multiply()
{
	int result,res,i,j;
	for (i=1;i<=lenb;i++)
	{
		res=0;
		for (j=1;j<=lena;j++)
		{
			t=i+j;
			result=a[j]*b[i]+res+c[t-1];
			res=result/10;
			c[t-1]=result%10;
		}
		while(res)
		{
			c[t++]=res%10;
			res/=10;
		}
	}
	while (c[t]==0) t--;
	return;
}

void print(int c[])
{
	int i;
	for (i=t;i>=1;i--) printf("%d",c[i]);
	printf("\n");
	return;
}

int main(int argc, char *argv[]) {
	gets(reversea);
	gets(reverseb);
	lena=strlen(reversea);
	lenb=strlen(reverseb);
	Reverse();
	Multiply();
	print(c);
	return 0;
}

C++实现

#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
int main( ) 
{
    int a[250]={0},b[250]={0},c[40000]={0},la,lb;
    int result=0,res=0,t=1; 
    char a1[201],b1[201];
    gets(a1);
    gets(b1);
    la=strlen(a1);
    lb=strlen(b1);
    for(int i=0;i<=la-1;i++) 
    {
        a[la-i]=a1[i]-48;   
    }
    for(int i=0;i<=lb-1;i++)   
    {
        b[lb-i]=b1[i]-48;
    }
    for(int i=1;i<=lb;i++)
    {
        res=0;
        for(int j=1;j<=la;j++)
        {
            t=j+i;
            result=a[j]*b[i]+res+c[t-1];  
            res=result/10;     
            c[t-1]=result%10;  
        }
        while(res) 
        {
            c[t++]=res%10;  
            res=res/10;
        }
    }
    while(c[t]==0) 
    {
        t--;
    }
    for(int i=t;i>=1;i--) 
    {
        cout<<c[i];
    }
    return 0;
}