(不)恢复余数法

(不)恢复余数法

(不)恢复余数法是一种用于二进制除法的算法。

算法原理

与数学中竖式除法类似,只是人为应用竖式除法时,可以直观的看出每一步余数与除数的大小,而(不)恢复余数法是计算机计算除法的一中方法,每一步都需要计算比较后才能得出结果。

计算过程

被除数 MnM_nMn,除数 D,商 Q = 0,余 R = 0,n 为 M 的整数位数,M[i] 为 M 的第 i 位数。

恢复余数法

  • 循环条件:i = n - 1; i >= 0; i-- 或者 Q 达到指定精度;
  • 循环体:
    • R = R << 1 + M[i] - D;
    • 如果 R >= 0, 则 Q = Q << 1 + 1;
    • 如果 R < 0, 则 Q = Q << 1 + 0, R = R + D(恢复余数);

不恢复余数法

  • 循环条件:i = n - 1; i >= 0; i-- 或者 Q 达到指定精度;
  • 循环体:
    • 如果 R >= 0, 则 R = R << 1 + M[n - 1] - D, 否则 R = R << 1 + M[n - 1] + D;
    • 如果 R >= 0, 则 Q = Q << 1 + 1;
    • 如果 R < 0, 则 Q = Q << 1 + 0;
  • 纠正余数: 如果 R < 0, 则 R = R + D;

两者区别

对于任意数 A, B, 总是存在关系:A×2±B=(A±B)×2±BA × 2 ± B = (A ± B) × 2 ± B

### 补码恢复余数在小数运算中的应用 #### 1. 补码恢复余数概述 补码恢复余数是一种用于定点数除法的有效算法,特别适用于二进制小数的运算。该方通过利用补码来简化硬件设计并提高效率。符号位单独处理,而数值部分则使用绝对值的补码形式进行操作[^4]。 #### 2. 小数运算的具体实现过程 对于小数而言,假设给定两个有符号的小数X和Y (其中|X|<|Y|),目标是求解 X/Y 的商Q以及最终余数R: - 设置初始条件:令被除数为[X]补,即X的补码表示; - 设定计数器i=0, 并准备存储每一位商的结果Qi。 - **循环执行以下步骤直到完成n次迭代(n取决于所需精度)**: - 判断当前余数是否为正: - 若为正,则说明够减,此时应将商位置1(Qi←1), 同时更新新的余数为旧余数减去除数[Y]补; - 反之当余数为负时表示够减,这时应该把当前位置设成零(Qi←0),并将之前产生的错误纠正回来——也就是加上除数[-Y]补以恢复正常状态下的余数。 - 左移一位新形成的余数以便于下一轮比较,并继续上述流程直至达到预定次数为止。 - **调整结果** - 对于最后一次可能存在的负余数情况做特殊处理,必要时再次添加除数使其变为正值。 - 整理所有得到的部分商构成完整的商值Q,并给出剩余未满整数倍关系的终态余数R。 ```python def complement_division(dividend, divisor, n_bits): """ Perform division using the non-restoring method with two's complement representation. :param dividend: Dividend as an integer value represented by its binary string form ('0b...') :param divisor: Divisor as an integer value represented by its binary string form ('0b...') :param n_bits: Number of bits used to represent numbers including sign bit :return: Tuple containing quotient and remainder both in their original signed decimal forms """ # Convert inputs from strings into integers while preserving signs abs_x = int(''.join(['0' if c == '0' else '1' for c in bin(abs(int(dividend)))[2:].zfill(n_bits)]).lstrip('0') or '0', 2) abs_y = int(''.join(['0' if c == '0' else '1' for c in bin(abs(int(divisor)))[2:].zfill(n_bits)]).lstrip('0') or '0', 2) # Initialize variables needed during computation q = [''] * n_bits # Quotient digits will be stored here one at a time starting MSB first r = format(abs_x, f'#0{n_bits}b')[2:] # Start off with full width absolute version of input number for i in range(0, n_bits): temp_r = int(r, base=2) << 1 if ((temp_r >= abs_y)): new_rem = '{:#0{width}b}'.format(temp_r-abs_y, width=n_bits+2)[2:] q[i], r = "1", new_rem elif((temp_r < abs_y)): new_rem = '{:#0{width}b}'.format(temp_r+abs_y, width=n_bits+2)[2:] q[i], r = "0", new_rem final_quotient = ''.join(q[::-1]) last_remainder = int(r,base=2)>>1 return (-int(final_quotient, 2) if '-' in str(dividend)+str(divisor) else int(final_quotient, 2)),last_remainder ``` 此函数实现了基于补码的恢复余数来进行两数相除的操作,返回商及其对应的余数。注意输入参数应当是以带有前缀`0b...`的形式传递过来的二进制字符串表达式[^5]。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

伸手不FUN

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值