本文转载至:http://blog.csdn.net/hitwhylz/article/details/9700935,并加以完善。
完善内容:增加了余数的输出。
大数除法,应该算是四则运算里面最难的一种了。不同于一般的模拟,除法操作步数模仿手工除法,而是利用减法操作实现的。
其基本思想是反复做除法,看从被除数里面最多能减去多少个除数,商就是多少。
逐个减显然太慢,要判断一次最多能减少多少个整的10的n次方。
以7546除23为例。
先减去23的100倍,就是2300,可以减3次,余下646。 此时商就是300;
然后646减去23的10倍,就是230,可以减2次,余下186。此时商就是320;
然后186减去23,可以减8次,此时商就是328.
根据这个思想,不难写出下面的代码。
还是那句话,可能算法效率不是很高。但是常规解题思路一般就是这样了。
如果以后有能力,有时间了。 我会试着去优化。
ps:大数系列学习资源来自 <c程序设计竞赛实训教程>一书和一些大牛的博客。
注意:程序不保留小数(只有商,没有余数),看了很多程序都是没有小数。
1 /*
2 本程序说明:
3
4 大数除法 http://blog.csdn.net/hitwhylz/article/details/9700935
5
6 */
7
8 #include<stdio.h>
9 #include<string.h>
10 #include<stdlib.h>
11 #define MaxLen 200
12 //函数SubStract功能:
13 //用长度为len1的大整数p1减去长度为len2的大整数p2
14 // 结果存在p1中,返回值代表结果的长度
15 //不够减 返回-1 正好够 返回0
16 int SubStract( int *p1, int *p2, int len1, int len2 )
17 {
18 int i;
19 if( len1 < len2 )
20 return -1;
21 if( len1 == len2 )
22 { //判断p1 > p2
23 for( i=len1-1; i>=0; i-- )
24 {
25 if( p1[i] > p2[i] ) //若大,则满足条件,可做减法
26 break;
27 else if( p1[i] < p2[i] ) //否则返回-1
28 return -1;
29 }
30 }
31 for( i=0; i<=len1-1; i++ ) //从低位开始做减法
32 {
33 p1[i] -= p2[i];
34 if( p1[i] < 0 ) //若p1<0,则需要借位
35 {
36 p1[i] += 10; //借1当10
37 p1[i+1]--; //高位减1
38 }
39 }
40 for( i=len1-1; i>=0; i-- ) //查找结果的最高位
41 if( p1[i] ) //最高位第一个不为0
42 return (i+1); //得到位数并返回
43 return 0; //两数相等的时候返回0
44 }
45 int main()
46 {
47 int n, k, i, j; //n:测试数据组数
48 int len1, len2; //大数位数
49 int nTimes; //两大数相差位数
50 int nTemp; //Subtract函数返回值
51 int num_a[MaxLen]; //被除数
52 int num_b[MaxLen]; //除数
53 int num_c[MaxLen]; //商
54 char str1[MaxLen + 1]; //读入的第一个大数
55 char str2[MaxLen + 1]; //读入的第二个大数
56
57 scanf("%d",&n);
58 while ( n-->0 )
59 {
60 scanf("%s", str1); //以字符串形式读入大数
61 scanf("%s", str2);
62
63 for ( i=0; i<MaxLen; i++ ) //初始化清零操作
64 {
65 num_a[i] = 0;
66 num_b[i] = 0;
67 num_c[i] = 0;
68 }
69
70 len1 = strlen(str1); //获得大数的位数
71 len2 = strlen(str2);
72
73 for( j=0, i=len1-1; i>=0; j++, i-- )
74 num_a[j] = str1[i] - '0'; //将字符串转换成对应的整数,颠倒存储
75 for( j=0, i=len2-1; i>=0; j++, i-- )
76 num_b[j] = str2[i] - '0';
77
78 if( len1 < len2 ) //如果被除数小于除数,结果为0
79 {
80 printf("0\n");
81 continue; //利用continue直接跳出本次循环。 进入下一组测试
82 }
83 nTimes = len1 - len2; //相差位数
84 for ( i=len1-1; i>=0; i-- ) //将除数扩大,使得除数和被除数位数相等
85 {
86 if ( i>=nTimes )
87 num_b[i] = num_b[i-nTimes];
88 else //低位置0
89 num_b[i] = 0;
90 }
91 len2 = len1;
92 for( j=0; j<=nTimes; j++ ) //重复调用,同时记录减成功的次数,即为商
93 {
94 while((nTemp = SubStract(num_a,num_b + j,len1,len2 - j)) >= 0)
95 {
96 len1 = nTemp; //结果长度
97 num_c[nTimes-j]++;//每成功减一次,将商的相应位加1
98 }
99 }
100
101 //输出商
102 for( i=MaxLen-1; num_c[i]==0 && i>=0; i-- );//跳过高位0
103 if( i>=0 )
104 for( ; i>=0; i-- )
105 printf("%d", num_c[i]);
106 else
107 printf("0");
108 printf("\n");
109
110
111 //此时的num_a存的就是余数(其实就是取模)
112 for( i=MaxLen-1; num_a[i]==0 && i>=0; i-- );//跳过高位0
113 if( i>=0 )
114 for( ; i>=0; i-- )
115 printf("%d", num_a[i]);
116 else
117 printf("0");
118 printf("\n");
119
120 }
121 return 0;
122 }
同类文章:
【模板小程序】 十进制大数相乘(正数、负数、0均可),包含合法性检查:http://www.cnblogs.com/xiaoxi666/p/7272255.html
【模板小程序】十进制大数相加(正整数版本+整数版本【正负0】),包含合法性检查:http://www.cnblogs.com/xiaoxi666/p/7258312.html
『注:本文来自博客园“小溪的博客”,若非声明均为原创内容,请勿用于商业用途,转载请注明出处http://www.cnblogs.com/xiaoxi666/』