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(記号計算)発表会議 : ISSAC'25(国際記号計算・代数計算シンポジウム)論文リンク : https://arxiv.org/abs/2510.13456 微分体における導数の完全約簡φは、体がその定数部分体上の線形作用素である。この約簡により、要素fを導数と剰余項φ(f)の和に分解することが可能になる。φの直接的な応用は、f が体内で積分可能であることと当且つつφ(f) = 0であることが同値であることである。本論文は、原始塔における導数の完全約簡をアルゴリズム的に提示する。原始塔の典型的な例は、(多重)対数関数と対数積分によって生成される微分体である。剰余項と留数を利用して、原始塔における要素が初等積分を持つための必要十分条件を提供し、特定の原始塔における非D-finite関数に対するtelescoperの構成方法について論じる。
記号積分における核心的問題 :記号計算において、関数が初等形式の積分を持つかどうかを判定することは基本的な問題である。超越Liouville関数に対して、この問題は通常、単項式拡張によって記述される。完全約簡の重要性 :完全約簡は線形作用素であり、微分体における任意の要素を導数部分と「最小」の剰余項に分解することができる。この分解は以下の点で重要である:関数の体内可積性の判定 約簡に基づく創造的なtelescoping 有限項積分(和) 既存方法の限界 :加法分解(additive decomposition)は常に線形写像ではなく、理論的および実践的な利便性に欠ける 既存の完全約簡は主に超指数関数、代数関数、D-finite関数など特定の型に対するものである 重要なカテゴリーである原始塔(primitive tower)に対する体系的な完全約簡アルゴリズムが欠けている 本研究の動機は以下に由来する:
理論的必要性 :原始塔における完全約簡の完全な理論的枠組みを確立するアルゴリズム的必要性 :原始塔における完全約簡を計算するための効率的なアルゴリズムを開発する応用的必要性 :完全約簡を初等積分計算とtelescoper構成に応用する原始塔における導数の完全約簡のアルゴリズム的枠組みを確立 :完全約簡を構成するための体系的な3段階法を提案重要な補助アルゴリズムを開発 :補助約簡(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 ) に対して、アルゴリズムは3段階に分かれる:
段階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:補空間の固定
有効基技術を通じて、F [ t ] ′ F[t]' F [ t ] ′ に関する A \mathcal{A} A の 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 ) 異なる複雑度の被積分関数の4つのグループをテスト。各グループは異なる次数の多項式導数を含む。
計算時間 :3つの方法の平均実行時間成功率 :正しい積分結果を返すことができるかどうか適用範囲 :異なる複雑度の問題を処理する能力第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
第2グループ(線形多項式係数) :
次数 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拡張における完全約簡の探索理論的厳密性 :数学的導出は厳密であり、定理証明は完全であるアルゴリズムの革新性 :3段階構成法は独創的であり、複雑な場合分析を回避している実用価値が高い :記号積分における重要な問題を解決し、直接的な応用価値を持つ実験の充分性 :詳細な性能比較とケース分析を提供執筆密度が高い :技術的内容が密集しており、読者の数学的背景に対する要求が高いアルゴリズム複雑度分析が不十分 :理論的複雑度分析が欠けている実験範囲が限定的 :主に3層の原始塔でテストされており、より高次元の場合の性能は不明学術的貢献 :記号計算分野に重要な理論的ツールを提供実用価値 :計算機代数システムの改善に直接応用可能再現性 :詳細なアルゴリズム記述と実験データを提供計算機代数システムにおける記号積分モジュール 対数関数と対数積分を含む数学計算 関数の可積性を判定する必要がある理論研究 創造的telescopingの前処理ステップ 論文は36篇の関連文献を引用しており、記号積分、完全約簡、創造的telescopingなど関連分野の重要な研究をカバーしており、本研究に堅実な理論的基礎を提供している。