四 大数除法。
大数除法,应该算是四则运算里面最难的一种了。不同于一般的模拟,除法操作步数模仿手工除法,而是利用减法操作实现的。
其基本思想是反复做除法,看从被除数里面最多能减去多少个除数,商就是多少。
逐个减显然太慢,要判断一次最多能减少多少个整的10的n次方。
以7546除23为例。
先减去23的100倍,就是2300,可以减3次,余下646。 此时商就是300;
然后646减去23的10倍,就是230,可以减2次,余下186。此时商就是320;
然后186减去23,可以减8次,此时商就是328.
根据这个思想,不难写出下面的代码。
还是那句话,可能算法效率不是很高。但是常规解题思路一般就是这样了。
如果以后有能力,有时间了。 我会试着去优化。
ps:大数系列学习资源来自 <c程序设计竞赛实训教程>一书和一些大牛的博客。
- #include<stdio.h>
- #include<string.h>
- #include<stdlib.h>
- #define MaxLen 200
-
-
-
-
- int SubStract( int *p1, int *p2, int len1, int len2 )
- {
- int i;
- if( len1 < len2 )
- return -1;
- if( len1 == len2 )
- {
- for( i=len1-1; i>=0; i-- )
- {
- if( p1[i] > p2[i] )
- break;
- else if( p1[i] < p2[i] )
- return -1;
- }
- }
- for( i=0; i<=len1-1; i++ )
- {
- p1[i] -= p2[i];
- if( p1[i] < 0 )
- {
- p1[i] += 10;
- p1[i+1]--;
- }
- }
- for( i=len1-1; i>=0; i-- )
- if( p1[i] )
- return (i+1);
- return 0;
- }
- int main()
- {
- int n, k, i, j;
- int len1, len2;
- int nTimes;
- int nTemp;
- int num_a[MaxLen];
- int num_b[MaxLen];
- int num_c[MaxLen];
- char str1[MaxLen + 1];
- char str2[MaxLen + 1];
-
- scanf("%d",&n);
- while ( n-->0 )
- {
- scanf("%s", str1);
- scanf("%s", str2);
-
- for ( i=0; i<MaxLen; i++ )
- {
- num_a[i] = 0;
- num_b[i] = 0;
- num_c[i] = 0;
- }
-
- len1 = strlen(str1);
- len2 = strlen(str2);
-
- for( j=0, i=len1-1; i>=0; j++, i-- )
- num_a[j] = str1[i] - '0';
- for( j=0, i=len2-1; i>=0; j++, i-- )
- num_b[j] = str2[i] - '0';
-
- if( len1 < len2 )
- {
- printf("0\n");
- continue;
- }
- nTimes = len1 - len2;
- for ( i=len1-1; i>=0; i-- )
- {
- if ( i>=nTimes )
- num_b[i] = num_b[i-nTimes];
- else
- num_b[i] = 0;
- }
- len2 = len1;
- for( j=0; j<=nTimes; j++ )
- {
- while((nTemp = SubStract(num_a,num_b + j,len1,len2 - j)) >= 0)
- {
- len1 = nTemp;
- num_c[nTimes-j]++;
- }
- }
-
-
- for( i=MaxLen-1; num_c[i]==0 && i>=0; i-- );
- if( i>=0 )
- for( ; i>=0; i-- )
- printf("%d", num_c[i]);
- else
- printf("0");
- printf("\n");
- }
- return 0;
- }