修复Python中求解y*x²=定值的整数x、y的代码问题

原代码的核心问题在于：使用sympy.solve求解带整数约束的方程时，返回的是参数化通解而非具体的整数解集合，导致后续的过滤逻辑完全失效——sol[x]是符号变量而非具体整数，sympy.simplify(sol[x]).is_integer永远返回False，最终得不到任何有效解。

修复方案1：枚举约数直接计算 #

既然方程是y = value / x²，那么x必须是value的整数约数，且x²能整除value。直接枚举所有可能的x值并验证即可：

import sympy

def find_integer_values(value):
    integer_solutions = []
    # 获取value绝对值的所有正约数，再扩展正负值
    positive_divisors = sympy.divisors(abs(value))
    all_x_candidates = [d for d in positive_divisors] + [-d for d in positive_divisors]
    
    for x_val in all_x_candidates:
        if x_val == 0:
            continue  # x=0时左边为0，仅当value=0时有解，此处单独处理
        x_squared = x_val ** 2
        if value % x_squared == 0:
            y_val = value // x_squared
            integer_solutions.append((x_val, y_val))
    
    # 去重（若value=0时会有重复，此处可选）
    return list(set(integer_solutions))

# 示例调用
given_value = 1920
integer_solutions = find_integer_values(given_value)

if integer_solutions:
    print("Integer solutions for x and y:", integer_solutions)
else:
    print("No integer solutions found.")

修复方案2：使用SymPy的丢番图方程求解器 #

SymPy内置的diophantine函数专门用于求解整数解问题，能直接返回所有有效整数解集合：

import sympy

def find_integer_values(value):
    x, y = sympy.symbols('x y', integer=True)
    equation = sympy.Eq(y * x**2, value)
    
    # 求解丢番图方程，返回整数解的集合
    solution_set = sympy.diophantine(equation)
    # 转换为易读的元组列表
    integer_solutions = [(sol[0], sol[1]) for sol in solution_set]
    
    return integer_solutions

# 示例调用
given_value = 1920
integer_solutions = find_integer_values(given_value)

if integer_solutions:
    print("Integer solutions for x and y:", integer_solutions)
else:
    print("No integer solutions found.")

原代码错误点说明 #

1. sympy.solve在处理带整数约束的多元方程时，不会枚举所有具体整数解，而是返回参数化的通解表达式（比如x为任意满足x²|value的整数，y=value/x²），因此sol[x]是符号变量而非具体数值。
2. 后续的过滤逻辑试图判断符号变量是否为整数，这完全不符合SymPy的符号计算逻辑，导致所有解被错误过滤。

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