使用Boost cpp_int计算sqrt平方后与原数不等的问题求助

问题描述 #

使用Boost库的cpp_int类型处理大整数时，对整数x计算sqrt(x)后再平方，结果与原数x不相等，示例代码及输出如下：

示例代码 #

#include <iostream>
#include <boost/multiprecision/cpp_int.hpp>
using namespace boost::multiprecision;
using namespace std;

int main() {
   boost::multiprecision::cpp_int num = 7459874565236544789;
   
   boost::multiprecision::cpp_int num_sqrt = sqrt(num);
   boost::multiprecision::cpp_int num_sqrt_square = num_sqrt * num_sqrt;

   cout <<  "\n num: " << num <<
            "\n num_sqrt: " << num_sqrt <<
            "\n num_sqrt_square: " << num_sqrt_square <<
            "\n _num_sqrt_square: " << sqrt(num) * sqrt(num);
}

实际输出 #

num: 7459874565236544789
num_sqrt: 2731277094
num_sqrt_square: 7459874564209084836
_num_sqrt_square: 7459874564209084836

num != num_sqrt_square !!!

原因分析 #

Boost的cpp_int是纯整数类型，其中的sqrt()函数返回的是整数平方根——即向下取整后的最大整数n，满足n² ≤ x < (n+1)²。你的原数7459874565236544789并不是完全平方数，因此取整后的平方根平方自然小于原数。

解决方法 #

根据不同需求，可采用以下方式处理：

• 验证整数平方根有效性：计算整数平方根后，通过num_sqrt * num_sqrt == num判断原数是否为完全平方数，或通过(num_sqrt + 1) * (num_sqrt + 1) > num确认整数平方根的正确性。
• 获取精确平方根（含小数）：如果需要保留平方根的完整精度，可使用Boost多精度库中的浮点类型（如cpp_dec_float），示例代码如下：

#include <iostream>
#include <boost/multiprecision/cpp_dec_float.hpp>
#include <boost/multiprecision/cpp_int.hpp>

using namespace boost::multiprecision;
using namespace std;

int main() {
    cpp_int num = 7459874565236544789;
    cpp_dec_float_50 num_sqrt = sqrt(cpp_dec_float_50(num));
    // 平方后转换回整数类型（注意精度范围）
    cpp_int restored_num = static_cast<cpp_int>(num_sqrt * num_sqrt);
    
    cout << "\n num: " << num <<
            "\n num_sqrt (float): " << num_sqrt <<
            "\n restored_num: " << restored_num;
}

• 整数场景下的差值验证：若仅需确认整数平方根与原数的关系，可计算两者的差值：

cpp_int diff = num - num_sqrt * num_sqrt;
cout << "\n 原数与平方数的差值: " << diff; // 输出1027459953，即原数比平方数大该值

内容的提问来源于stack exchange，提问作者zkoria