如何读取浮点字符串并转换为GMP的mpq_t类型（规避精度损失）

因为gmp_sscanf不支持直接将十六进制浮点格式读取为mpq_t，且通过mpf_t中转会丢失精度，最可靠的方式是手动解析十六进制浮点字符串的结构，直接转换为分数形式赋值给mpq_t，完全避免浮点数精度问题。

核心原理 #

十六进制浮点格式0xI.FpE的数值等价于：
$$(I + \frac{F}{16^{len(F)}}) \times 2^E$$
其中：

• I是十六进制整数部分
• F是十六进制小数部分，len(F)是其位数
• E是十进制指数（2的幂次）

将其转换为分数形式：
$$\frac{(I \times 16^{len(F)} + F) \times 2^E}{16{len(F)}}$$
之后对分子分母约分即可得到精确的mpq_t值。

实现代码 #

#include <gmp.h>
#include <string.h>
#include <stdlib.h>

// 将十六进制浮点字符串解析为mpq_t，成功返回0，失败返回-1
int parse_hex_float_to_mpq(mpq_t result, const char *str) {
    // 校验格式开头
    if (strncmp(str, "0x", 2) != 0) return -1;
    str += 2;

    char *dot_pos = strchr(str, '.');
    char *p_pos = strchr(str, 'p');
    if (!p_pos) return -1; // 缺少指数部分

    // 拆分整数、小数、指数部分
    char int_part[1024] = {0}, frac_part[1024] = {0}, exp_str[32] = {0};
    if (dot_pos && dot_pos < p_pos) {
        strncpy(int_part, str, dot_pos - str);
        strncpy(frac_part, dot_pos + 1, p_pos - dot_pos - 1);
    } else {
        strncpy(int_part, str, p_pos - str);
    }
    strncpy(exp_str, p_pos + 1, strlen(p_pos + 1));

    int exp = atoi(exp_str);
    mpz_t int_z, frac_z, pow16, numerator, denominator;
    mpz_inits(int_z, frac_z, pow16, numerator, denominator, NULL);

    // 解析十六进制部分为大整数
    strlen(int_part) ? mpz_set_str(int_z, int_part, 16) : mpz_set_ui(int_z, 0);
    strlen(frac_part) ? mpz_set_str(frac_z, frac_part, 16) : mpz_set_ui(frac_z, 0);

    // 计算16^len(F)
    size_t frac_len = strlen(frac_part);
    mpz_ui_pow_ui(pow16, 16, frac_len);

    // 构造分子：I*16^len(F) + F
    mpz_mul(numerator, int_z, pow16);
    mpz_add(numerator, numerator, frac_z);

    // 构造分母：16^len(F)
    mpz_set(denominator, pow16);

    // 处理指数部分
    if (exp > 0) {
        mpz_mul_2exp(numerator, numerator, exp); // 分子乘以2^E
    } else if (exp < 0) {
        mpz_mul_2exp(denominator, denominator, -exp); // 分母乘以2^|E|
    }

    // 赋值给mpq_t并约分
    mpq_set_num(result, numerator);
    mpq_set_den(result, denominator);
    mpq_canonicalize(result);

    mpz_clears(int_z, frac_z, pow16, numerator, denominator, NULL);
    return 0;
}

// 测试示例
int main() {
    mpq_t q;
    mpq_init(q);
    const char *test_str = "0x12f.3p3";

    if (parse_hex_float_to_mpq(q, test_str) == 0) {
        gmp_printf("解析结果: %Qd\n", q); // 输出应为4851/2，对应十进制2425.5
    }

    mpq_clear(q);
    return 0;
}

关键说明 #

• 全程使用mpz_t大整数运算，完全避免浮点精度丢失
• 手动拆分字符串格式，确保解析的准确性
• 最后通过mpq_canonicalize自动约分，得到最简分数形式

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