高精度计算

这里使用最简单的方式处理高精度计算问题：使用数组保存高精度数的每一位，模仿四则运算竖式规则依次操作每一位，完成高精度计算。

1 高精度加法：模仿加法，对应的每一位相加，可以分别处理每一次产生的进位也可以最后集中处理。很简单，你懂得。

2 高精度乘法：模仿乘法，因数1的第i位与因数2的第j位相乘的结果加到积的第i+j位上。最后集中处理进位问题。

3 高精度除法：模仿除法，其实是做减法。除法的每次商是几，我们是凭观察看出来的。计算机"观察“不出来，只能一次一次的减，每减一次，商的对应位就增一。我们做商是从高位做到低位直到个位时结束。其实是将除数扩大到与被除数同一个数量级上去除（也就是去减）。模仿这个小学时学会的计算方法，我们就可以教会计算机做除法了。下面是除法的代码。

char sz1[MAX_LENGTH+10]; char sz2[MAX_LENGTH+10]; int an1[MAX_LENGTH+10]; int an2[MAX_LENGTH+10]; int ans[MAX_LENGTH+10]; int Substract(int *p1,int *p2, int nLen1, int nLen2){ if(nLen1<nLen2){ return -1; } int nLen=0; int i,j; if(nLen1==nLen2){ for(i=nLen1-1;i>=0;i--){ if(p1[i]>p2[i]){ break; } if(p1[i]<p2[i]){ return -1; } } } for(i=0;i<nLen1;i++){ p1[i] -= p2[i]; if(p1[i]<0){ p1[i] += 10; p1[i+1]--; } } for(i=nLen1-1;i>=0;i--){ if(p1[i]){ return i+1; } } return 0; } void Divide(){ memset(an1,0,sizeof(an1)); memset(an2,0,sizeof(an2)); memset(ans,0,sizeof(ans)); puts("Input 2 numbers"); gets(sz1); gets(sz2); int nLen1=strlen(sz1); int nLen2=strlen(sz2); int i,j; for(i=nLen1-1,j=0; i>=0; i--,j++){ an1[j]=sz1[i]-'0'; } for(i=nLen2-1,j=0; i>=0; i--,j++){ an2[j]=sz2[i]-'0'; } nLen1=Substract(an1,an2,nLen1,nLen2); if(nLen1<0){ printf("0/n"); return ; // divide is done. } if(nLen1==0){ printf("1/n"); return ; } ans[0]++; //an1-an2 has done. int nTimes=nLen1-nLen2; if(nTimes<0){ // Remaintation < an2 goto OutputResult; } if(nTimes>0){ //Length(Remaintation) == Length(an2)+nTimes for(i=nLen2+nTimes-1;i>=0;i--){ if(i>=nTimes){ an2[i]=an2[i-nTimes]; } else{ an2[i]=0; } } } nLen2=nLen1; // after for-clause two lengths have become same. int nTmp; for(i=0;i<=nTimes;i++){ while((nTmp=Substract(an1,an2+i,nLen1,nLen2-i))>=0){ nLen1=nTmp; ans[nTimes-i]++; } } OutputResult: for(i=0;i<=nTimes;i++){ if(ans[i]>=10){ ans[i+1] += ans[i]/10; ans[i] %= 10; } } if(ans[nTimes+1]){ printf("%d",ans[nTimes+1]); } if(ans[nTimes]){ printf("%d",ans[nTimes]); } for(i=nTimes-1;i>=0;i--){ printf("%d",ans[i]); } }  

ps：

1 memset()函数将一片连续的内存区域填充为某一指定内容，wmemset()函数可以设置宽字符。原型如下

 

void *memset(void *dest, int c, size_t count);

wchar_t *wmemset( wchar_t *dest, wchar_t c, size_t count);

dest ：Pointer to destination.

c ：Character to set.

count ：Number of characters.

2 strlen()自然是获取字符串长度了。