高精度乘法入門詳解(高精乘高精)
阿新 • • 發佈:2019-01-07
高精度乘法。輸入兩個正整數,求它們的積。
【演算法分析】
類似加法,可以用豎式求乘法。在做乘法運算時,同樣也有進位,同時對每一位進行乘法運算時,必須進行錯位相加,如圖3、圖4。
分析c陣列下標的變化規律,可以寫出如下關係式:ci = c’i +c”i +…由此可見,c i跟a[i]*b[j]乘積有關,跟上次的進位有關,還跟原c i的值有關,分析下標規律,有c[i+j-1]= a[i]*b[j]+ x + c[i+j-1]; x=c[i+j-1]/10 ; c[i+j-1]%=10;
高精度乘法的參考程式:
#include<iostream>
#include<cstring>
#include<cstdio>
using namespace std;
int main()
{
char a1[100],b1[100];
int a[100],b[100],c[100],lena,lenb,lenc,i,j,x;
memset(a,0,sizeof(a));
memset(b,0,sizeof(b));
memset(c,0,sizeof(c));
gets(a1);gets(b1);
lena=strlen(a1);lenb=strlen(b1);
for (i=0;i<=lena-1;i++) a[lena-i]=a1[i]-48 ;
for (i=0;i<=lenb-1;i++) b[lenb-i]=b1[i]-48;
for (i=1;i<=lena;i++)
{
x=0; //用於存放進位
for (j=1;j<=lenb;j++) //對乘數的每一位進行處理
{
c[i+j-1]=a[i]*b[j]+x+c[i+j-1]; //當前乘積+上次乘積進位+原數
x=c[i+j-1]/10;
c[i+j-1] %= 10;
}
c[i+lenb]=x; //進位
}
lenc=lena+lenb;
while (c[lenc]==0&&lenc>1) //刪除前導0
lenc--;
for (i=lenc;i>=1;i--)
cout<<c[i];
cout<<endl;
return 0;
}