[論文レビュー] The method of characteristics, and "problem 89" of Graham, Knuth and Patashnik

Herbert S. Wilf|ArXiv.org|Jun 30, 2004

Advanced Mathematical Identities参考文献 3被引用数 10

ひとこと要約

この論文は、線形偏微分方程式を解くための特性曲線法を、DNA配列照合問題と、グレアマー、クランク、パタシュニクの一般化された二項係数に関する問題89に適用する。後者について、6つのパラメータのうち1つが消える場合の正確な閉形式解を導出し、この条件下で関連する多項式がすべて実数解しかも持たないことを証明し、行和の母関数および漸近的解析を提供する。

ABSTRACT

We apply the method of characteristics for the solution of pde's to two combinatorial problems. The first is finding an explicit form for a distribution that arises in bio-informatics. The second is a question raised by Graham, Knuth and Patashnik abiout a sequence of generalized binomial coefficients. We find an exact formula, which factors in an interesting way, in the case where one of the six parameters of the problem vanishes. We also show that the associated polynomial sequence has real zeros only, provided that one parameter vanishes, and the other five are nonnegative.

研究の動機と目的

• DNA配列照合に由来する再帰関係を母関数と偏微分方程式技術を用いて解くこと。
• グレアマー、クランク、パタシュニクの一般化された二項係数を定義する2変数線形再帰に関する問題89を解くこと。
• 6つのパラメータのうち1つが消える場合の一般化された二項係数の正確な閉形式表現を導出すること。
• 関連する多項式列がすべて実数解しかも持たない条件を確立すること。
• 特定のパrameter制約下での一般化されたパスカルの三角形における行和の漸近的成長を分析すること。

提案手法

• DNA照合問題の再帰関係を、母関数 $ Q(x,t) $ に対する1階線形偏微分方程式として形式化する。
• 特性曲線法を用いて偏微分方程式 $ (x-1)\frac{\partial Q}{\partial x} + (1-t)\frac{\partial Q}{\partial t} = Q $ を解き、$ Q(x,t) = \frac{1 - (1 - (1-t)(1-x))^n}{1-t} $ を得る。
• 母関数を展開して、$ Q_k(x) $ の3通りの同等な表現を導出し、特にコンactな二項係数和表現を含む。
• ゴスパーのアルゴリズムを用いて、$ Q_k(x) $ の異なる表現の同等性を確認し、特に $ Q_k(x) = (1-x)^{k+1} \sum_{j=0}^{n-k-1} \binom{j+k}{j} x^j $ の形を裏付ける。
• 一般化された二項係数問題に関して、母関数 $ u(x,y) $ を導出し、それを偏微分方程式に変換し、特性曲線法を用いて一般解を任意関数 $ f $ の形で得る。
• 初期条件を課して $ f $ を特定し、閉形式母関数を導出し、その後係数を抽出する。

実験結果

リサーチクエスチョン

• RQ1 DNA配列照合に由来する列 $ Q_k(x) $ の正確な閉形式表現は何か？
• RQ2 特性曲線法は、母関数に対して偏微分方程式を生成する組合せ的再帰に効果的に適用可能か？
• RQ3 一般化された二項係数再帰が実数解のみをもつ多項式を生成するパrameterの条件は何か？
• RQ4 $ \alpha' = 0 $ のとき、行和 $ \rho(n) = \sum_k {n \atopwithdelims| k} $ の漸近的挙動は何か？
• RQ5 再帰式のパラメータの1つが消えると、母関数の構造はどのように変化するか？

主な発見

• $ Q_k(x) $ の正確な閉形式表現として $ Q_k(x) = (1-x)^{k+1} \sum_{j=0}^{n-k-1} \binom{j+k}{j} x^j $ が得られ、これは簡潔でかつ列の累積分布の性質を明らかにする。
• $ k \geq n $ のとき $ Q_k(x) $ は恒等的に0となり、固定された $ k $ と $ |x| < 1 $ に対して $ \lim_{n \to \infty} Q_k(x) = 1 $ であるため、退化分布への収束を示す。
• $ \beta' = 0 $ で他のすべてのパラメータが非負のとき、多項式 $ \phi_n(x) = \sum_k {n \atopwithdelims| k} x^k $ はすべて実数かつ非正の解しかも持たない。これは帰納法とロルの定理を用いて証明される。
• 行和 $ \rho(n) $ の母関数は $ F(x) = \sum_n \rho(n) \frac{x^n}{n!} = \frac{(1 - \alpha x)^{-\gamma/\alpha}}{\left(1 + \frac{\beta'}{\beta}(1 - (1 - \alpha x)^{-\beta/\alpha})\right)^{1 + \gamma'/(\beta')}} $ として導出され、$ \alpha' = 0 $ のとき有効である。
• 母関数の解は超幾何関数の形で表現され、係数抽出により、二項係数と有理関数の累乗を含む二重和の公式が $ \rho(n) $ に対して得られる。
• 特性曲線法は、1つのパラメータが消える特殊ケースにおいて、問題89の完全な解を効果的に得ることができたが、一般解は依然として扱いにくい。