如何用SymPy求解线性不等式并计算最优整数a值？

首先，你遇到的核心问题有两个：一是如何正确判断代入数值后的不等式是否成立，二是在不等式不成立时快速计算满足条件的最大整数a。

问题分析 #

1. 不等式成立判断：SymPy代入数值后得到的是SymPy布尔对象，不能直接用if not ...判断，需要转换成Python原生布尔值。
2. 最优整数a计算：满足a ≤ n/(s+16)的最大整数就是n/(s+16)的向下取整值，用SymPy的floor函数即可直接得到。

修正后的代码示例 #

import sympy

num_attractors = 4
attractor_size = 64
network_size = 256
a, s, n = sympy.symbols('a s n')
expr = a * (s + 16) <= n

# 代入数值并转换为Python布尔值判断不等式是否成立
is_valid = bool(expr.subs([(a, num_attractors), (s, attractor_size), (n, network_size)]))

if not is_valid:
    # 计算满足a ≤ n/(s+16)的最大整数a
    optimal_a = sympy.floor(n / (s + 16)).subs([(s, attractor_size), (n, network_size)])
    print(f"最优整数a的值为: {optimal_a}")

关于solve和reduce_inequalities的说明 #

• solve(expr, a)会返回不等式的符号化解集，比如对于a*(s+16) <=n，输出是a ≤ n/(s+16)这样的区间表达式，并非具体数值解。
• reduce_inequalities(expr, a)同样是对不等式进行符号化简，输出符号形式的解集范围。如果你需要具体的整数数值，直接代入后取floor是更直接高效的方式。

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