1. 程式人生 > >高精度乘法入門詳解(高精乘高精)

高精度乘法入門詳解(高精乘高精)

高精度乘法。輸入兩個正整數,求它們的積。

【演算法分析】

類似加法,可以用豎式求乘法。在做乘法運算時,同樣也有進位,同時對每一位進行乘法運算時,必須進行錯位相加,如圖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; }