对于浮点数和双精度,快速乘法/除2(C/C++)

#include<iostream>
#include <cmath>

float fast_multiply(float a, float b) {
    int a_bits = *reinterpret_cast<int*>(&a);
    int b_bits = *reinterpret_cast<int*>(&b);
    int result_bits = (a_bits & 0x807fffff) + (b_bits & 0x807fffff);
    result_bits += 0x00800000;
    return *reinterpret_cast<float*>(&result_bits);
}

float fast_divide(float a, float b) {
    int a_bits = *reinterpret_cast<int*>(&a);
    int b_bits = *reinterpret_cast<int*>(&b);
    int result_bits = (a_bits & 0x807fffff) - (b_bits & 0x807fffff);
    result_bits += 0x00800000;
    return *reinterpret_cast<float*>(&result_bits);
}

int main() {
    float a = 1.23f;
    float b = 4.56f;
    float c = fast_multiply(a, b);
    float d = fast_divide(a, b);
    std::cout << "a * b = " << c << std::endl;
    std::cout << "a / b = " << d << std::endl;
    return 0;
}