[論文レビュー] Computation of the highest coefficients of weighted Ehrhart quasi-polynomials of rational polyhedra

Velleda Baldoni, Nicole Berline+3 more|arXiv (Cornell University)|Nov 7, 2010

Advanced Combinatorial Mathematics参考文献 20被引用数 4

ひとこと要約

この論文は、格子点に多項式重みが付与された有理数多面体の重み付きEhrhart準多項式の最高次係数を、局所的近似に基づく効率的なアルゴリズムで計算する手法を提示する。Barvinokの手法を拡張し、重み付きの場合に適応させるもので、生成関数と段階多項式表現を用いて、これらの係数を閉形式で計算する。LattE macchiatoなどの最先端ソフトウェアと比較しても、基本的なMaple実装でも競争力のある性能を示している。

ABSTRACT

This article concerns the computational problem of counting the lattice points inside convex polytopes, when each point must be counted with a weight associated to it. We describe an efficient algorithm for computing the highest degree coefficients of the weighted Ehrhart quasi-polynomial for a rational simple polytope in varying dimension, when the weights of the lattice points are given by a polynomial function h. Our technique is based on a refinement of an algorithm of A. Barvinok [Computing the Ehrhart quasi-polynomial of a rational simplex, Math. Comp. 75 (2006), pp. 1449--1466] in the unweighted case (i.e., h = 1). In contrast to Barvinok's method, our method is local, obtains an approximation on the level of generating functions, handles the general weighted case, and provides the coefficients in closed form as step polynomials of the dilation. To demonstrate the practicality of our approach we report on computational experiments which show even our simple implementation can compete with state of the art software.

研究の動機と目的

• 組合せ論、最適化、統計に中心的な問題である、有理数多面体における格子点の数え上げに伴う計算課題に対処すること。
• 各格子点に多項式関数 $ h(x) $ を用いて値が割り当てられる一般の重み付きケースに、Barvinokの非重み付きEhrhart係数計算を拡張すること。
• すべての係数を計算するNP困難性を避けるために、重み付きEhrhart準多項式の最高次係数を効率的に計算するアルゴリズムを開発すること。
• 希少な係数を、拡張パrameter $ n $ の段階多項式として閉形式で得る方法を提供し、実用的かつスケーラブルな計算を可能にすること。
• LattE macchiatoなどの最先端ツールと同等の性能を示す計算実験を通じて、手法の実用性を示すこと。

提案手法

• 生成関数を用いた重み付き格子点数の重み付き和の近似を実行する、Barvinokのアルゴリズムの改良版を用い、多項式重み $ h(x) $ を処理する。
• 単体なアフィン錐上で局所的近似戦略を適用し、頂点付近の生成関数の振る舞いを捉えるパッチ関数を用いて近似を行う。
• 中間の生成関数を導入して問題を扱いやすい部分に分解し、統合と体積に基づく公式を用いて最高次項を計算する。
• 主な革新点は、周期的係数 $ E_m(n \bmod q) $ を段階多項式として表現することで、最高次係数の閉形式での記号的計算を可能にすること。
• 高次元の単体に対して効率を高めるために、Barvinokの分解をユニモジュラー錐に分割するプライマルバージョンを活用する。
• Mapleで実装され、LattE macchiatoとの比較を通じて検証されており、ランダムな格子単体をテストケースとして用いている。

実験結果

リサーチクエスチョン

• RQ1 有理数多面体に多項式重みが付与された場合、重み付きEhrhart準多項式の最高次係数を、有理数多面体に対して効率的に計算できるか？
• RQ2 Barvinokの非重み付きアルゴリズムを、一般の多項式重み関数 $ h(x) $ を扱えるように拡張する方法は何か？ その際、計算効率を保つことができるか？
• RQ3 単体錐上での生成関数の局所的近似が、重み付きEhrhart準多項式の最高次係数をどれほど正確に再現できるか？
• RQ4 得られた係数を、拡張パrameter $ n $ の段階多項式として閉形式で表現できるか？ これにより記号的・アルゴリズム的処理が可能になるか？
• RQ5 新しいアルゴリズムの性能は、LattE macchiatoなどの既存の最先端ソフトウェアと比較して、全Ehrhart多項式を計算する際の性能に比べてどの程度優れているか？

主な発見

• 固定された $ k_0 $ に対して、次元が入力に含まれる場合でも、有理数多面体に対して、重み付きEhrhart準多項式の上位 $ k_0 $ 個の係数を多項式時間で計算可能である。
• 最高次係数 $ E_{d+M} $ は、正確に積分 $ \int_{\mathfrak{p}} h(x) \, dx $ に等しく、既知の幾何的解釈と整合していることを確認した。
• 計算実験では、著しい高速化が達成された：次元11では、上位3つの係数を1681CPU秒で計算したが、LattE macchiatoのプライマル法による全Ehrhart多項式計算は30分以上を要した。
• 解釈型のMapleで実装された本手法は、次元5以上ではLattE macchiatoのプライマル法を上回った。これは、コンパイル言語での最適化によるさらなる性能向上の可能性を示している。
• Barvinok分解のカットオフ（例：行列式1000で停止）を適用した場合、性能優位性が顕著に現れ、LattE macchiatoのPrimal 1000バージョンでも同様の傾向が確認された。
• 結果として、本アルゴリズムは理論的にも妥当であり、実用的にも競争力があることが確認された。特に、全Ehrhart計算が非現実的になる高次元単体において顕著な優位性を示した。