大数乘法程序

┏━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┓

┃作为ＣＳＤＮ博客第一篇文章，一时不知道写什么好。　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　┃

┃拿最近写的一个大数乘法的程序，是上ＰＫＵ的ＡＣＭ遇到的一道题。其实很简单的，就开头的题目。做完之后想了想决定扩┃

┃展成一个适用性强一点的程序，因为Ｃ＋＋语言本就没有直接计算大整数的函数之类的可以调用，于是就写成了下面的代码。┃

┃除了可以计算整数，小数也可以同样进行计算，暂时只能算正数。　　　　　　　　　　　　　　　　　　　　　　　　　　┃

┃ＰＳ：算法虽然很笨，也没有精力去改进，不过计算精度和准确性还是有保证的。　　　　　　　　　　　　　　　　　　　┃

┗━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━┛

 

#include <stdio.h> #include <memory.h> #include <string.h> void BigDecAdd(int *add1,int *add2,int nLengthAdd1=1,int nLengthAdd2=1); //大整数的加法函数,做乘法时需要用 void Add(int *add1,int add2); //小整数数加法函数,大整数加法内调用 int BigDecMul(int *nMul_1, int *nMul_2, int *result, int nLengthMul_1=1,int nLengthMul_2=1); //大整数乘法计算 int BigDecMul(char *nMul_1,char *nMul_2,char *result); //真正直接调用的函数 int main() { //理论上可以计算的大整数长度无上限,main函数只是一个使用实例,假设大整数的长度不超过90,结果长度不超过180 char strMul1[91]; char strMul2[91]; char strResult[181]; //C++中无数据类型可以直接存大整数,所以用字符串存. scanf("%90s",strMul1); scanf("%90s",strMul2); memset(strResult,0,181); //将输入的两个相乘的大整数传给int BigDecMul(char *nMul_1,char *nMul_2,char *result)即可 //计算结果存储在strResult中 BigDecMul(strMul1,strMul2,strResult); return 0; } int BigDecMul(char *nMul_1,char *nMul_2,char *result) { int strlen1,strlen2; int floatLength1,floatLength2; int nLen1,nLen2; int i,j; bool zeroEnd; int length1,length2; int temp; strlen1=(int)strlen(nMul_1); strlen2=(int)strlen(nMul_2); nLen1=(strlen1%5>0)?(strlen1/5+1):(strlen1/5); nLen2=(strlen2%5>0)?(strlen2/5+1):(strlen2/5); int *mul1=new int[nLen1+1]; int *mul2=new int[nLen2+1]; int *nResult; zeroEnd=true; length1=strlen1; temp=1; floatLength1=-1; memset(mul1,0,sizeof(int)*(nLen1+1)); for(i=strlen1-1,j=0;i>=0;i--) { if(nMul_1[i]=='0'&&zeroEnd) //去0 length1--; else { zeroEnd=false; if(nMul_1[i]=='.') floatLength1=length1-i-1; else if(nMul_1[i]>='0'||nMul_1[i]<='9') { mul1[j]=mul1[j]+(nMul_1[i]-'0')*temp; temp*=10; if(temp>=100000) { temp=1; j++; } } } } nLen1=mul1[j]!=0?j+1:j; zeroEnd=true; length2=strlen2; temp=1; floatLength2=-1; memset(mul2,0,sizeof(int)*(nLen2+1)); for(i=strlen2-1,j=0;i>=0;i--) { if(nMul_2[i]=='0'&&zeroEnd) length2--; else { zeroEnd=false; if(nMul_2[i]=='.') floatLength2=length2-i-1; else { mul2[j]=mul2[j]+(nMul_2[i]-'0')*temp; temp*=10; if(temp>=100000) { temp=1; j++; } } } } nLen2=mul2[j]!=0?j+1:j; nResult=new int[nLen1+nLen2]; memset(nResult,0,nLen1+nLen2); int nLenR=BigDecMul(mul1,mul2,nResult,nLen1,nLen2); char *strtemp=new char[strlen1+strlen2+1]; int lentemp; sprintf(strtemp,"%d",nResult[nLenR]); lentemp=(int)strlen(strtemp); for(i=nLenR-1;i>=0;i--) { sprintf(strtemp+lentemp,"%05d",nResult[i]); lentemp+=5; } temp=0; int floatLength=0; if(floatLength1<0) { temp+=strlen1-length1; } else { floatLength+=floatLength1; } if(floatLength2<0) { temp+=strlen2-length2; } else { floatLength+=floatLength2; } if(floatLength<temp) { memset(strtemp+lentemp,'0',temp-floatLength); lentemp+=(temp-floatLength); strtemp[lentemp]=0; } else if(floatLength>temp) { temp=floatLength-temp; if(temp<lentemp) { memcpy(strtemp+lentemp-temp+1,strtemp+lentemp-temp,temp); strtemp[lentemp-temp]='.'; strtemp[lentemp+1]=0; lentemp+=1; } else { memcpy(strtemp+temp-lentemp+2,strtemp,lentemp); strtemp[0]='0'; strtemp[1]='.'; memset(strtemp+2,'0',temp-lentemp); strtemp[temp+2]=0; lentemp+=(temp-lentemp+2); } } memcpy(result,strtemp,lentemp); result[lentemp]=0; delete mul1; delete mul2; //delete nResult; //delete strtemp; return 0; } int BigDecMul(int *nMul_1, int *nMul_2, int *result, int nLengthMul_1,int nLengthMul_2) { int i,j; __int64 tempresult,temp1,temp2; int *a=new int[nLengthMul_1+nLengthMul_2]; memset(result,0,(nLengthMul_1+nLengthMul_2)*sizeof(int)); for(i=0;i<nLengthMul_2;i++) { memset(a,0,(nLengthMul_1+nLengthMul_2)*sizeof(int)); for(j=0;j<nLengthMul_1;j++) { temp1=nMul_1[j]; temp2=nMul_2[i]; tempresult=temp1*temp2; temp2=tempresult%100000; Add(a+j,(int)temp2); Add(a+j+1,(int)((tempresult-temp2)/100000)); } BigDecAdd(result+i,a,nLengthMul_1+nLengthMul_2,j); } for(i=nLengthMul_1+nLengthMul_2-1;i>=0;i--) { if(result[i]!=0) return i; } delete a; return 0; } void Add(int *add1, int add2) { add1[0]+=add2; if(add1[0]>=100000) { add1[0]-=100000; Add(add1+1,1); } } void BigDecAdd(int *add1,int *add2,int nLengthAdd1,int nLengthAdd2) { int nBound=nLengthAdd1>nLengthAdd2?nLengthAdd1:nLengthAdd2; for(int i=0;i<nBound;i++) Add(add1+i,add2[i]); }