利用字符实现大数除法运算(纯c实现)

那年还是语言刚入门的我，不知纠结了几多次，才完成了能支持大数的算术四则运算，现在看起来别是一番滋味哈。

算术四则运算中就数除法的实现最为繁杂，基本思路按照除法笔算的步骤进行，计算商时一位一位的进行，当计算某一位商时则不断地与除数相减，得到的相减次数作为本位商的结果，实现代码如下，代码可读性不高，效率有待优化，英雄莫见笑。

/* * author:Zhanhang * * qq:273711460 * E-mail:zhanxinhang@gmail.com * HomePage:http://blog.csdn.net/zhanxinhang */ #include<stdio.h> #include<string.h> #include<malloc.h> //清多余的零 void clz(char r[]) { if(r[0]=='-')r++; //使r指针向前移一位； int len = strlen(r); int i=0,k=0; if(len<2)return; if(r[i]=='0') { for(;r[k]=='0'&&r[k+1]!='.';k++); for(;i<len-k;i++) { r[i]=r[i+k]; } r[i]='\0'; } } //清除数可不要的小数点 void clDot(char r[]) { int len = strlen(r); if(r[len-1]=='.') r[len-1]=0; } //为除计算之用的减法 void dsub(char a[],char b[]) { int i,alen,blen; char r,t,jiewei=0; char *pb ; alen = strlen(a); blen = strlen(b); pb = (char *)calloc(alen,sizeof(char)); //使pb与a位数对齐// memset(pb,'0',alen); t = alen; for(i=blen-1;i>=0;i--) pb[--t]=b[i]; for(i=alen-1;i>=0;i--) //相减运算 { t = a[i]-jiewei; if(t<pb[i]) jiewei=1; else jiewei = 0; r = (t+jiewei*10)-pb[i]; a[i] = r+'0'; } a[alen]='\0'; clz(a); free(pb); } //比较两数的大小 int cmp(const char a[],const char b[]) { int alen = strlen(a); int blen = strlen(b); if(alen>blen) return 1; else if(alen == blen) { if(strcmp(a,b)>0) return 1; if(strcmp(a,b)==0) return 0; } else return -1; } //除法计算 void div(const char a[],const char b[],char result[]) { bool isNegative = false; char *op1,*pa,*pb,*pr; int up,alen,blen,adotp,bdotp,i,k,dotp,t,t1,j,quo_size; /////////////判定符号/////////////// //如果为异号 if((a[0] == '-'||b[0] == '-')&&a[0] != b[0]) result[0] = '-',isNegative = true; //去除负号 if(a[0] == '-')a++; if(b[0] == '-')b++; /////////////////////////////////// alen = strlen(a)-1; //减去一位小数点 blen = strlen(b)-1; ///////获取被除数小数点移位后的位置////////// adotp = strchr(a,'.')-a; bdotp = strchr(b,'.')-b; //计算商小数点位置 dotp = adotp+blen-bdotp; if(isNegative)dotp++; //////////准备数据///////////// op1 = (char *)calloc(alen+blen+1,sizeof(char)); pa = (char *)calloc(alen+blen+1,sizeof(char)); pb = (char *)calloc(blen+1,sizeof(char)); pr = (char *)calloc(alen+blen+1,sizeof(char)); for(i = 0,t=0; i<=alen; i++) { if(a[i]!='.') pa[t++] = a[i]; } for(;t<dotp-1;t++) { pa[t] = '0'; } for(;t<blen;t++) { pa[t] = '0'; } pa[t] = '\0'; for(i = 0,t=0; i<=blen; i++) { if(b[i]!='.') pb[t++] = b[i]; } pb[t] = '\0'; clz(pa); clz(pb); ////////取得被除数的高位数op1,且op1大于被除数b////////// strncpy(op1,pa,strlen(pb)); if(strcmp(op1,pb)<0) { strncpy(op1,pa,strlen(pb)+1); } /////计算////// j = k = strlen(op1); t1=0; quo_size = strlen(pa)+1-k; //获取商的长度 while(t1<quo_size) { up = 0; t = cmp(op1,pb); while(t>=0) { dsub(op1,pb); t = cmp(op1,pb); up++; } pr[t1++] = up+'0'; op1[strlen(op1)]=pa[j++]; clz(op1); } quo_size+=50;//加50精度 while(t1<quo_size&&(cmp(op1,(char *)"0")>0)) { up = 0; op1[strlen(op1)]='0'; t = cmp(op1,pb); while(t>=0) { dsub(op1,pb); t = cmp(op1,pb); up++; } pr[t1++] = up+'0'; } ////////////////////////////// if(isNegative)t=1; else t=0; for(i=0;i<=t1;i++)//复制结果并给商加上小数点 { if(t==dotp) result[t++] = '.'; result[t++]=pr[i]; } clz(result); clDot(result); free(op1); free(pa); free(pb); } main() { char result[500]={0}; div("-9999.","3.",result); printf("%s\n",result); div("10000000000000000000000000.","5.",result); printf("%s\n",result); div("12345678900.","12345678900.",result); printf("%s\n",result); div("1.","3.",result); printf("%s\n",result); div("-444.","-2.",result); printf("%s\n",result); div("6660.","-3330.",result); printf("%s\n",result); div("1.","3333333333333333333333.",result); //bug printf("%s\n",result); }
结果：（每一行对应于main中的每一次调用div所得到的结果）

======= welcome to my HomePage(http://blog.csdn.net/zhanxinhang) to have a communication =======