A complete reduction $Ï$ for derivatives in a differential field is a linear operator on the field over its constant subfield. The reduction enables us to decompose an element $f$ as the sum of a derivative and the remainder $Ï(f)$. A direct application of $Ï$ is that $f$ is in-field integrable if and only if $Ï(f) = 0.$
In this paper, we present a complete reduction for derivatives in a primitive tower algorithmically. Typical examples for primitive towers are differential fields generated by (poly-)logarithmic functions and logarithmic integrals. Using remainders and residues, we provide a necessary and sufficient condition for an element from a primitive tower to have an elementary integral, and discuss how to construct telescopers for non-D-finite functions in some special primitive towers.
论文ID : 2510.13456标题 : Complete Reduction for Derivatives in a Primitive Tower作者 : Hao Du (北京邮电大学), Yiman Gao (约翰内斯开普勒大学), Wenqiao Li (中科院数学机械化重点实验室), Ziming Li (中科院数学机械化重点实验室)分类 : cs.SC (Symbolic Computation)发表会议 : ISSAC'25 (International Symposium on Symbolic and Algebraic Computation)论文链接 : https://arxiv.org/abs/2510.13456 微分域中导数的完全约简ϕ \phi ϕ 是域在其常数子域上的线性算子。该约简使我们能够将元素f f f 分解为导数和余项ϕ ( f ) \phi(f) ϕ ( f ) 的和。ϕ \phi ϕ 的一个直接应用是f f f 在域内可积当且仅当ϕ ( f ) = 0 \phi(f) = 0 ϕ ( f ) = 0 。本文算法化地提出了原始塔中导数的完全约简。原始塔的典型例子是由(多重)对数函数和对数积分生成的微分域。利用余项和留数,我们提供了原始塔中元素具有初等积分的充要条件,并讨论了如何为某些特殊原始塔中的非D-finite函数构造telescoper。
符号积分中的核心问题 :在符号计算中,确定一个函数是否具有初等形式的积分是一个基本问题。对于超越Liouville函数,这个问题通常通过单项式扩张来描述。完全约简的重要性 :完全约简是一种线性算子,能够将微分域中的任意元素分解为导数部分和"最小"的余项。这种分解对于:判断函数的域内可积性 基于约简的创造性telescoping 有限项积分(求和) 现有方法的局限性 :加法分解(additive decomposition)不总是线性映射,缺乏理论和实践上的便利性 现有的完全约简主要针对超指数函数、代数函数、D-finite函数等特定类型 对于原始塔(primitive tower)这一重要类别缺乏系统的完全约简算法 本文的研究动机源于:
理论需求 :为原始塔建立完whole的完全约简理论框架算法需求 :开发高效的算法来计算原始塔中的完全约简应用需求 :将完全约简应用于初等积分计算和telescoper构造建立了原始塔中导数完全约简的算法框架 :提出了系统的三步方法来构造完全约简开发了关键的辅助算法 :包括辅助约简(AuxiliaryReduction)、基构造(Basis)和投影(Projection)算法提供了初等积分的充要条件 :基于余项和留数给出了原始塔中元素具有初等积分的判定准则扩展了telescoper构造方法 :为某些非D-finite函数提供了telescoper存在性的充分条件实现了高效的算法 :实验表明该方法在多数情况下优于现有方法给定原始塔F 0 ⊂ F 1 ⊂ ⋯ ⊂ F n F_0 \subset F_1 \subset \cdots \subset F_n F 0 ⊂ F 1 ⊂ ⋯ ⊂ F n ,其中F i = F i − 1 ( t i ) F_i = F_{i-1}(t_i) F i = F i − 1 ( t i ) 且t i t_i t i 是F i − 1 F_{i-1} F i − 1 上的原始单项式,目标是构造完全约简ϕ : F n → F n \phi: F_n \to F_n ϕ : F n → F n 使得:
对任意f ∈ F n f \in F_n f ∈ F n ,存在唯一的g ∈ F n g \in F_n g ∈ F n 和r ∈ im ( ϕ ) r \in \text{im}(\phi) r ∈ im ( ϕ ) 使得f = g ′ + r f = g' + r f = g ′ + r ϕ ( f ) = r \phi(f) = r ϕ ( f ) = r 是f f f 的余项ker ( ϕ ) = F n ′ \ker(\phi) = F_n' ker ( ϕ ) = F n ′ (所有导数的集合)对于原始单项式扩张F ( t ) F(t) F ( t ) ,算法分三步:
步骤1:定义辅助子空间
定义A = im ( ϕ ) ⊗ C C [ t ] \mathcal{A} = \text{im}(\phi) \otimes_C C[t] A = im ( ϕ ) ⊗ C C [ t ] 作为F [ t ] ′ F[t]' F [ t ] ′ 在F [ t ] F[t] F [ t ] 中的辅助子空间,其中ϕ : F → F \phi: F \to F ϕ : F → F 是F F F 上已有的完全约简。
步骤2:确定交集的基
构造F [ t ] ′ ∩ A F[t]' \cap \mathcal{A} F [ t ] ′ ∩ A 的C C C -基{ v 0 , v 1 , v 2 , … } \{v_0, v_1, v_2, \ldots\} { v 0 , v 1 , v 2 , … } ,其中:
v 0 = ϕ ( t ′ ) v_0 = \phi(t') v 0 = ϕ ( t ′ ) v i = ϕ ( t ′ ) t i − M i , 0 ( t i ) v_i = \phi(t')t^i - M_{i,0}(t^i) v i = ϕ ( t ′ ) t i − M i , 0 ( t i ) (对i ≥ 1 i \geq 1 i ≥ 1 )步骤3:固定补空间
通过有效基技术确定A \mathcal{A} A 在F [ t ] F[t] F [ t ] 中关于F [ t ] ′ F[t]' F [ t ] ′ 的补空间A θ \mathcal{A}_\theta A θ 。
算法3.4 (AuxiliaryReduction) :
输入:p ∈ F[t]
输出:(q,r) ∈ F[t] × A 使得 p = q' + r
1. 初始化 p̃ ← p, q ← 0, r ← 0
2. while p̃ ≠ 0 do
d ← deg(p̃), l ← lc(p̃)
计算l的R-对(g, φ(l))
q ← q + gt^d, r ← r + φ(l)t^d
p̃ ← p̃ - lt^d - (dgt')t^(d-1)
3. return (q,r)
算法3.12 (Projection) :
将辅助子空间中的元素投影到F [ t ] ′ F[t]' F [ t ] ′ 和θ \theta θ -补空间。
引理3.6的关键结果 :证明了{ v 0 , v 1 , … } \{v_0, v_1, \ldots\} { v 0 , v 1 , … } 构成F [ t ] ′ ∩ A F[t]' \cap \mathcal{A} F [ t ] ′ ∩ A 的C C C -基,其中每个v i v_i v i 的次数为i i i 且首项系数为ϕ ( t ′ ) \phi(t') ϕ ( t ′ ) 。
定理3.13的主要结果 :
F ( t ) = F ( t ) ′ ⊕ A θ ⊕ S t F(t) = F(t)' \oplus \mathcal{A}_\theta \oplus S_t F ( t ) = F ( t ) ′ ⊕ A θ ⊕ S t
其中S t S_t S t 是简单元素的集合,A θ \mathcal{A}_\theta A θ 是θ \theta θ -补空间。
算法3.10将R-对的计算数量从O ( d 2 ) O(d^2) O ( d 2 ) 优化到O ( d ) O(d) O ( d ) (其中d d d 是多项式次数) 通过递归构造,整个原始塔的完全约简可以高效计算 计算平台 :Intel Core i9 3.6GHz,16GB内存软件环境 :Maple 2021对比系统 :本文方法(CR)、AddDecompInField算法(AD)、Maple的int函数实验使用原始塔Q ( x ) ( t 1 , t 2 , t 3 ) \mathbb{Q}(x)(t_1, t_2, t_3) Q ( x ) ( t 1 , t 2 , t 3 ) ,其中:
t 1 = log ( x ) t_1 = \log(x) t 1 = log ( x ) t 2 = log ( x + 1 ) t_2 = \log(x+1) t 2 = log ( x + 1 ) t 3 = log ( t 1 ) t_3 = \log(t_1) t 3 = log ( t 1 ) 测试了四组不同复杂度的被积函数,每组包含不同次数的多项式导数。
计算时间 :三种方法的平均运行时间成功率 :能否返回正确积分结果适用范围 :处理不同复杂度问题的能力第一组(稠密有理函数系数) :
次数 CR(秒) AD(秒) int(秒) 1 1.42 0.96 1.15 2 8.32 10.42 4.52 3 37.01 47.36 23.30 4 122.55 149.02 53.43 5 1085.04 >3600 166.27
第二组(线性多项式系数) :
次数 CR(秒) AD(秒) int(秒) 6 0.90 1.23 3.83 8 2.09 4.29 17.46 10 7.05 12.31 31.61 12 12.56 31.08 66.22 14 30.35 57.67 144.70 16 62.11 170.70 322.19
CR方法总体优于AD方法 ,在大多数测试用例中表现更好相比Maple的int函数 ,CR在复杂度较高时表现优异,但在简单情况下略慢稳定性更好 :CR和AD都能处理int函数无法处理的某些积分问题算法组件分析 :HermiteReduce和AuxiliaryReduction是最耗时的部分,Projection相对高效例4.5 :对于函数
f = ( ( x − 1 ) 2 t 1 + x ) t 2 3 + x ( x − 1 ) t 1 x 2 ( x − 1 ) t 2 2 f = \frac{((x-1)^2 t_1 + x)t_2^3 + x(x-1)t_1}{x^2(x-1)t_2^2} f = x 2 ( x − 1 ) t 2 2 (( x − 1 ) 2 t 1 + x ) t 2 3 + x ( x − 1 ) t 1
CR成功找到其积分,而Maple和Mathematica都未能给出初等形式的结果。
例5.4 :展示了完整的初等积分计算过程,包括余项分析和留数计算。
完全约简理论 :已有针对超指数函数5 、代数函数12,15 、D-finite函数6,13,35 的完全约简加法分解 :在reduction-based创造性telescoping中的应用2-4,24 符号积分 :超越Liouville函数的初等积分算法8,17,26,28,34 首次系统化 :为原始塔建立完整的完全约简理论避免复杂的情况分析 :相比其他扩张类型,原始单项式的处理更为简洁双重技术结合 :将积分分部法与参数Risch方程求解相结合理论贡献 :建立了原始塔中导数完全约简的完整理论框架算法贡献 :提供了高效的算法实现,在多数情况下优于现有方法应用价值 :成功应用于初等积分计算和telescoper构造适用范围限制 :主要针对原始塔,对其他类型的超越扩张需要进一步研究计算复杂度 :对于高次多项式,计算时间仍然较长实现优化空间 :HermiteReduce等基础算法仍有优化空间扩展到其他扩张类型 :将方法推广到超指数扩张等其他情况算法优化 :进一步提高计算效率,特别是对于高维情况理论深化 :探索更一般的Liouville扩张中的完全约简理论严谨性 :数学推导严密,定理证明完整算法创新性 :三步构造方法具有原创性,避免了复杂的情况分析实用价值高 :解决了符号积分中的重要问题,具有直接应用价值实验充分性 :提供了详细的性能对比和案例分析写作密度高 :技术内容密集,对读者的数学背景要求较高算法复杂度分析不够详细 :缺乏理论的复杂度分析实验范围有限 :主要在三层原始塔上测试,更高维情况的表现未知学术贡献 :为符号计算领域提供了重要的理论工具实用价值 :可直接应用于计算机代数系统的改进可复现性 :提供了详细的算法描述和实验数据计算机代数系统中的符号积分模块 涉及对数函数和对数积分的数学计算 需要判断函数可积性的理论研究 创造性telescoping的预处理步骤 论文引用了36篇相关文献,涵盖了符号积分、完全约简、创造性telescoping等相关领域的重要工作,为本研究提供了坚实的理论基础。