大数(正整数)除法之多精除以单精
大数运算中会经常遇到多精除以单精的情况，有时也会把某些除法分解成多精除以单精加快运行速度，这里所说单精数是不超过cpu位长二分这一的数，32位机是65535，在此以千进制数为例讲解一下具体实现，
For i = 1 To 被除数数组长度
运算结果(i) = 被除数(i) \ 单精数
被除数(i + 1) = (被除数(i) Mod 单精数) * 1000 + 被除数(i + 1)
Next
用实例演示一下：123 456 789 125 /234
运算结果(1)=123 \ 234=0 '被除数(1) 中初始存放123，234为除数
被除数(2)=(123 mod 234) *1000+456=123456 ‘1000为千进制数的基，(123 mod 234) *1000的结果为123000，被除数(2)中初始存放456
运算结果(2)=123456\ 234=527 '被除数(2) 中存放123456，234为除数
被除数(3)=(123456 mod 234) *1000+789=138789 ‘1000为千进制数的基，(123456 mod 234) *1000的结果为138000，被除数(3)中初始存放789
运算结果(3)=138789\ 234=593 '被除数(3) 中存放138789，234为除数
被除数(4)=(138789 mod 234) *1000+125=027125 ‘1000为千进制数的基，(138789 mod 234) *1000的结果为027000，被除数(4)中初始存放125
运算结果(4)=027125\ 234=115 '被除数(4) 中存放138789，234为除数
被除数(5)=(027125 mod 234) *1000+0=215000 ‘1000为千进制数的基，(027125 mod 234) *1000的结果为215000 ，被除数(5)中初始存放0
此时运算结果数组中依次存放000 527 593 115 为该除法的整数商，若要计算到小数后某某位，可以通过扩充被除数长度实现。
这段代码在提升落叶高精度表达式计算器1.2的速度中发挥出关键作用。
下面是网上收集的一段多精除以单精c代码：
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int main()
{
char a;
int c;
int b,l;
while (~scanf("%s%d",&a,&b))
{
memset(c,0,sizeof(c));
l=strlen(a);
int ys=0;
int t=0,temp,i;
for (i=0;i<l;i++)
{
temp=a-'0'+10*ys;
if (temp>=b)
{
c=temp/b;
ys=temp % b ;
t++;
}
else ys=temp;
}
while (c==0) i--;
for (int j=i;j>=0;j--)
printf("%d",c);
if (ys) printf(" %d",ys);
printf("\n");
}
return 0;
}
页:
[1]