用C语言实现快速求平方根倒数算法

以下是用C语言实现快速求平方根倒数算法的代码：

#include <stdio.h>
#include <math.h>

double sqrt_reciprocal(double x) {
    union {
        int i;
        double d;
    } u;
    u.d = x;
    u.i = 0x5fe6ec85e7de30da - (u.i >> 1);
    return u.d;
}

int main() {
    double x = 2.0;
    double y = sqrt_reciprocal(x);
    printf("sqrt(%.1f) = %.5f\n", x, sqrt(x));
    printf("sqrt_reciprocal(%.1f) = %.5f\n", x, y);
    printf("sqrt(%.1f) * sqrt_reciprocal(%.1f) = %.5f\n", x, x, sqrt(x) * y);
    return 0;
}

该算法使用IEEE 754浮点数格式的位操作，通过对浮点数的二进制表示进行修改，来实现快速求平方根倒数的计算。