如何使用SymPy内置函数生成无理数非十进制数字序列

SymPy没有单独封装“生成无理数任意进制数字序列”的专用独立函数，但内置的高精度求值接口原生支持多进制输出，完全可以满足需求，不需要自行实现核心的数值计算、进制转换逻辑。

针对圆周率前1024位二进制序列的实现代码 #

核心调用SymPy常量pi的evalf()方法，该方法支持通过base参数直接指定求值的进制，通过prec参数指定对应进制下的有效位数，全程无自定义算法逻辑，同版本SymPy在任意环境下输出完全一致，满足公开复现要求：

from sympy import pi

# 二进制下多取2位精度：覆盖整数位占用、规避末尾舍入误差影响
pi_bin_value = pi.evalf(prec=1026, base=2)
# 提取小数部分的二进制位，截取前1024位转为列表结构
bit_sequence = [
    int(bit) 
    for bit in str(pi_bin_value).split('.')[1][:1024]
]

执行后bit_sequence就是存储了0/1值的类数组列表，可直接用于组合优化测试用例生成。

通用说明 #

• 该方法对所有SymPy内置的无理数常量（自然常数E、黄金分割率GoldenRatio等）、任意目标进制（8进制、16进制等）均适用，仅需替换对应常量、调整base参数和截取长度即可。
• 精度设置时多留2-3位冗余量即可，避免浮点数舍入规则导致末尾位偏差，冗余位最终会被截断，不会影响目标序列的正确性。
• 检索不到对应函数的原因：多进制高精度求值属于SymPy通用数值求值模块的内置功能，没有单独拆分为面向“进制序列生成”场景的专用API，因此按场景关键词检索很难直接定位到evalf()的base参数用法。

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