2025-11-20T17:49:15.132482

Complete Reduction for Derivatives in a Primitive Tower

Du, Gao, Li et al.
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.
academic

Complete Reduction for Derivatives in a Primitive Tower

基本信息

  • 论文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是域在其常数子域上的线性算子。该约简使我们能够将元素ff分解为导数和余项ϕ(f)\phi(f)的和。ϕ\phi的一个直接应用是ff在域内可积当且仅当ϕ(f)=0\phi(f) = 0。本文算法化地提出了原始塔中导数的完全约简。原始塔的典型例子是由(多重)对数函数和对数积分生成的微分域。利用余项和留数,我们提供了原始塔中元素具有初等积分的充要条件,并讨论了如何为某些特殊原始塔中的非D-finite函数构造telescoper。

研究背景与动机

问题背景

  1. 符号积分中的核心问题:在符号计算中,确定一个函数是否具有初等形式的积分是一个基本问题。对于超越Liouville函数,这个问题通常通过单项式扩张来描述。
  2. 完全约简的重要性:完全约简是一种线性算子,能够将微分域中的任意元素分解为导数部分和"最小"的余项。这种分解对于:
    • 判断函数的域内可积性
    • 基于约简的创造性telescoping
    • 有限项积分(求和)
  3. 现有方法的局限性
    • 加法分解(additive decomposition)不总是线性映射,缺乏理论和实践上的便利性
    • 现有的完全约简主要针对超指数函数、代数函数、D-finite函数等特定类型
    • 对于原始塔(primitive tower)这一重要类别缺乏系统的完全约简算法

研究动机

本文的研究动机源于:

  1. 理论需求:为原始塔建立完whole的完全约简理论框架
  2. 算法需求:开发高效的算法来计算原始塔中的完全约简
  3. 应用需求:将完全约简应用于初等积分计算和telescoper构造

核心贡献

  1. 建立了原始塔中导数完全约简的算法框架:提出了系统的三步方法来构造完全约简
  2. 开发了关键的辅助算法:包括辅助约简(AuxiliaryReduction)、基构造(Basis)和投影(Projection)算法
  3. 提供了初等积分的充要条件:基于余项和留数给出了原始塔中元素具有初等积分的判定准则
  4. 扩展了telescoper构造方法:为某些非D-finite函数提供了telescoper存在性的充分条件
  5. 实现了高效的算法:实验表明该方法在多数情况下优于现有方法

方法详解

任务定义

给定原始塔F0F1FnF_0 \subset F_1 \subset \cdots \subset F_n,其中Fi=Fi1(ti)F_i = F_{i-1}(t_i)tit_iFi1F_{i-1}上的原始单项式,目标是构造完全约简ϕ:FnFn\phi: F_n \to F_n使得:

  • 对任意fFnf \in F_n,存在唯一的gFng \in F_nrim(ϕ)r \in \text{im}(\phi)使得f=g+rf = g' + r
  • ϕ(f)=r\phi(f) = rff的余项
  • ker(ϕ)=Fn\ker(\phi) = F_n'(所有导数的集合)

核心算法架构

1. 三步构造方法

对于原始单项式扩张F(t)F(t),算法分三步:

步骤1:定义辅助子空间 定义A=im(ϕ)CC[t]\mathcal{A} = \text{im}(\phi) \otimes_C C[t]作为F[t]F[t]'F[t]F[t]中的辅助子空间,其中ϕ:FF\phi: F \to FFF上已有的完全约简。

步骤2:确定交集的基 构造F[t]AF[t]' \cap \mathcal{A}CC-基{v0,v1,v2,}\{v_0, v_1, v_2, \ldots\},其中:

  • v0=ϕ(t)v_0 = \phi(t')
  • vi=ϕ(t)tiMi,0(ti)v_i = \phi(t')t^i - M_{i,0}(t^i)(对i1i \geq 1

步骤3:固定补空间 通过有效基技术确定A\mathcal{A}F[t]F[t]中关于F[t]F[t]'的补空间Aθ\mathcal{A}_\theta

2. 关键算法组件

算法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]'θ\theta-补空间。

3. 技术创新点

引理3.6的关键结果:证明了{v0,v1,}\{v_0, v_1, \ldots\}构成F[t]AF[t]' \cap \mathcal{A}CC-基,其中每个viv_i的次数为ii且首项系数为ϕ(t)\phi(t')

定理3.13的主要结果F(t)=F(t)AθStF(t) = F(t)' \oplus \mathcal{A}_\theta \oplus S_t 其中StS_t是简单元素的集合,Aθ\mathcal{A}_\thetaθ\theta-补空间。

算法复杂度分析

  • 算法3.10将R-对的计算数量从O(d2)O(d^2)优化到O(d)O(d)(其中dd是多项式次数)
  • 通过递归构造,整个原始塔的完全约简可以高效计算

实验设置

测试环境

  • 计算平台:Intel Core i9 3.6GHz,16GB内存
  • 软件环境:Maple 2021
  • 对比系统:本文方法(CR)、AddDecompInField算法(AD)、Maple的int函数

测试数据

实验使用原始塔Q(x)(t1,t2,t3)\mathbb{Q}(x)(t_1, t_2, t_3),其中:

  • t1=log(x)t_1 = \log(x)
  • t2=log(x+1)t_2 = \log(x+1)
  • t3=log(t1)t_3 = \log(t_1)

测试了四组不同复杂度的被积函数,每组包含不同次数的多项式导数。

评价指标

  • 计算时间:三种方法的平均运行时间
  • 成功率:能否返回正确积分结果
  • 适用范围:处理不同复杂度问题的能力

实验结果

主要结果

性能对比表格

第一组(稠密有理函数系数)

次数CR(秒)AD(秒)int(秒)
11.420.961.15
28.3210.424.52
337.0147.3623.30
4122.55149.0253.43
51085.04>3600166.27

第二组(线性多项式系数)

次数CR(秒)AD(秒)int(秒)
60.901.233.83
82.094.2917.46
107.0512.3131.61
1212.5631.0866.22
1430.3557.67144.70
1662.11170.70322.19

关键发现

  1. CR方法总体优于AD方法,在大多数测试用例中表现更好
  2. 相比Maple的int函数,CR在复杂度较高时表现优异,但在简单情况下略慢
  3. 稳定性更好:CR和AD都能处理int函数无法处理的某些积分问题
  4. 算法组件分析:HermiteReduce和AuxiliaryReduction是最耗时的部分,Projection相对高效

案例分析

例4.5:对于函数 f=((x1)2t1+x)t23+x(x1)t1x2(x1)t22f = \frac{((x-1)^2 t_1 + x)t_2^3 + x(x-1)t_1}{x^2(x-1)t_2^2} CR成功找到其积分,而Maple和Mathematica都未能给出初等形式的结果。

例5.4:展示了完整的初等积分计算过程,包括余项分析和留数计算。

相关工作

主要研究方向

  1. 完全约简理论:已有针对超指数函数5、代数函数12,15、D-finite函数6,13,35的完全约简
  2. 加法分解:在reduction-based创造性telescoping中的应用2-4,24
  3. 符号积分:超越Liouville函数的初等积分算法8,17,26,28,34

本文的创新性

  • 首次系统化:为原始塔建立完整的完全约简理论
  • 避免复杂的情况分析:相比其他扩张类型,原始单项式的处理更为简洁
  • 双重技术结合:将积分分部法与参数Risch方程求解相结合

结论与讨论

主要结论

  1. 理论贡献:建立了原始塔中导数完全约简的完整理论框架
  2. 算法贡献:提供了高效的算法实现,在多数情况下优于现有方法
  3. 应用价值:成功应用于初等积分计算和telescoper构造

局限性

  1. 适用范围限制:主要针对原始塔,对其他类型的超越扩张需要进一步研究
  2. 计算复杂度:对于高次多项式,计算时间仍然较长
  3. 实现优化空间:HermiteReduce等基础算法仍有优化空间

未来方向

  1. 扩展到其他扩张类型:将方法推广到超指数扩张等其他情况
  2. 算法优化:进一步提高计算效率,特别是对于高维情况
  3. 理论深化:探索更一般的Liouville扩张中的完全约简

深度评价

优点

  1. 理论严谨性:数学推导严密,定理证明完整
  2. 算法创新性:三步构造方法具有原创性,避免了复杂的情况分析
  3. 实用价值高:解决了符号积分中的重要问题,具有直接应用价值
  4. 实验充分性:提供了详细的性能对比和案例分析

不足

  1. 写作密度高:技术内容密集,对读者的数学背景要求较高
  2. 算法复杂度分析不够详细:缺乏理论的复杂度分析
  3. 实验范围有限:主要在三层原始塔上测试,更高维情况的表现未知

影响力

  1. 学术贡献:为符号计算领域提供了重要的理论工具
  2. 实用价值:可直接应用于计算机代数系统的改进
  3. 可复现性:提供了详细的算法描述和实验数据

适用场景

  • 计算机代数系统中的符号积分模块
  • 涉及对数函数和对数积分的数学计算
  • 需要判断函数可积性的理论研究
  • 创造性telescoping的预处理步骤

参考文献

论文引用了36篇相关文献,涵盖了符号积分、完全约简、创造性telescoping等相关领域的重要工作,为本研究提供了坚实的理论基础。