使用Pari/GP提取实数各位数字的更优实现方法咨询

原代码的核心问题 #

• 用strprintf("%g", n)做字符串转换是最大隐患：%g默认仅输出6位有效数字，超过精度的数值会直接截断，极大/极小数还会自动转为科学计数法格式，完全没法得到正确的数字序列
• 遍历查找小数点的逻辑冗余，用字符串替换即可直接去掉小数点

方案1：纯数值运算实现（推荐，无边界问题） #

real_to_int(n) = {
    \\ 整数直接返回
    if(type(n) == "t_INT" || n == floor(n), return(floor(n)));
    my(dec = frac(n), scale = 1);
    \\ 计算把小数部分转成整数需要乘的10的幂次
    while(frac(dec * scale) != 0, scale *= 10);
    my(res = n * scale);
    \\ 移除尾部零
    return(res / 10^valuation(res, 10));
}

优势：完全避开字符串转换的格式问题，只要输入在Pari/GP的设置精度范围内，不管多大/多小的实数都能正确处理，不会出现截断、格式识别错误的问题。

方案2：优化字符串处理实现（速度更快） #

real_to_int(n) = {
    if(n == floor(n), return(floor(n)));
    \\ 直接用convert转字符串，不会有%g的截断问题
    my(s = strreplace(convert(n, str), ".", ""));
    my(res = eval(s));
    return(res / 10^valuation(res, 10));
}

注意：如果输入是超出convert默认输出范围的极大/极小数，会输出科学计数法格式导致转换错误，仅适合处理常规精度的普通实数。

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