heyongjun666/aicu-graph/math-symbolic-computation

import sympy as sp
x, n = sp.symbols('x n', real=True)   # 声明实数可帮助化简
y = sp.Function('y')
sp.simplify(expr)                     # 化简
sp.expand(expr)                       # 展开
sp.factor(expr)                       # 因式分解
sp.diff(expr, x)                      # 求导
sp.diff(expr, x, 2)                   # 高阶导
sp.integrate(expr, x)                 # 不定积分
sp.integrate(expr, (x, 0, 1))         # 定积分
sp.limit(expr, x, sp.oo)              # 极限
sp.series(expr, x, 0, 6)              # 级数展开（含 O 项）
sp.solve(sp.Eq(expr, 0), x)           # 解方程
sp.solveset(expr, x, sp.S.Reals)      # 集合式求解
sp.summation(expr, (n, 0, 10))        # 求和
sp.simplify(sp.diff(sp.integrate(expr, x), x))  # 微积分基本定理自检

import numpy as np
from mpmath import mp
mp.dps = 30                          # 高精度，避免浮点误判
np.sum(arr)                          # 求和
np.prod(arr)                         # 求积

# 有限情形穷举：报告搜索范围与命中结果
candidates = range(-100, 101)
hits = [k for k in candidates if cond(k)]