高精度計算之大數乘法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;
}