[논문 리뷰] Sparse Gröbner Bases: the Unmixed Case

Jean‐Charles Faugère, Pierre-Jean Spaenlehauer외 1명|arXiv (Cornell University)|2014. 02. 28.

Polynomial and algebraic computation참고 문헌 31인용 수 17

한 줄 요약

이 논문은 희소 다항식계의 기호적 해법을 가속화하기 위해 반군 대수에서의 희소 그로브너 기저를 도입한다. 특히 모든 다항식이 동일한 뉴턴 다면체를 공유하는 균일한 경우에 초점을 맞춘다. 뉴턴 다면체의 조합적 성질을 활용하는 F5 및 FGLM 알고리즘의 희소 변형을 제안하며, 실험에서 100배 이상의 가속도를 달성하고, 다면체의 부피, 정규성, 에흐라르 다항식에 따라 복잡도 한계를 제시한다.

ABSTRACT

Toric (or sparse) elimination theory is a framework developped during the last decades to exploit monomial structures in systems of Laurent polynomials. Roughly speaking, this amounts to computing in a \emph{semigroup algebra}, \emph{i.e.} an algebra generated by a subset of Laurent monomials. In order to solve symbolically sparse systems, we introduce \emph{sparse Gröbner bases}, an analog of classical Gröbner bases for semigroup algebras, and we propose sparse variants of the $F_5$ and FGLM algorithms to compute them. Our prototype "proof-of-concept" implementation shows large speed-ups (more than 100 for some examples) compared to optimized (classical) Gröbner bases software. Moreover, in the case where the generating subset of monomials corresponds to the points with integer coordinates in a normal lattice polytope $\mathcal P\subset\mathbb R^n$ and under regularity assumptions, we prove complexity bounds which depend on the combinatorial properties of $\mathcal P$. These bounds yield new estimates on the complexity of solving $0$-dim systems where all polynomials share the same Newton polytope (\emph{unmixed case}). For instance, we generalize the bound $\min(n_1,n_2)+1$ on the maximal degree in a Gröbner basis of a $0$-dim. bilinear system with blocks of variables of sizes $(n_1,n_2)$ to the multilinear case: $\sum n_i - \max(n_i)+1$. We also propose a variant of Fröberg's conjecture which allows us to estimate the complexity of solving overdetermined sparse systems.

연구 동기 및 목표

• 희소 다항식계의 효율적인 알고리즘 개발을 위해 반군 대수를 통해 단항식의 구조를 활용한다.
• 밀도가 낮은 단항식 지지 집합을 가진 시스템에 적용할 때 전통적인 그로브너 기저 방법의 계산적 병목 현상을 해결한다.
• 뉴턴 다면체의 내재된 조합적 불변량에 따라 희소 그로브너 기저 계산의 복잡도 한계를 수립한다.
• 희소 그로브너 기저 이론을 활용해 이항계에 대한 알려진 차수 한계를 다항계로 일반화한다.
• 과잉 결정된 희소 시스템에서의 복잡도를 추정하기 위해 프로베르그의 추측을 수정한 형태를 제안한다.

제안 방법

• 유한 부분집합 $M \subset \mathbb{Z}^n$ 에 의해 생성되는 아핀 반군 $S$ 를 갖는 반군 대수 $k[S]$ 에서 희소 그로브너 기저를 정의한다.
• $(\alpha,1)$ 이 $\alpha \in M$ 인 경우에 대해 생성되는 동차 반군 $S^{(h)}$ 를 갖는 동차 대수 $k[S^{(h)}]$ 에서 희소 MatrixF5 알고리즘의 변형을 도입하여 희소 그로브너 기저를 계산한다.
• 모노미얼 사상에 의해 희소 그로브너 기저를 유리 함수의 단변수 표현으로 변환하기 위한 FGLM 알고리즘의 희소 변형을 개발한다.
• 감소된 희소 그로브너 기저 내에서의 최대 희소 차수를 제한하는 '워치독 차수'를 복잡도 지표로 활용한다.
• 뉴턴 다면체 $\mathscr{P}$ 의 조합적 불변량, 즉 정규화된 부피 $\operatorname{vol}(\mathscr{P})$, 캐스텔누오보-무디드 정규성 $\operatorname{reg}(k[S^{(h)}])$, 에흐라르 다항식 $\operatorname{HP}_{\mathscr{P}}(\ell)$ 을 활용하여 복잡도 추정을 유도한다.
• 정규성 가정 하에 $\mathscr{P} \cap \mathbb{Z}^n$ 에서의 지지 집합을 갖는 시스템에 알고리즘을 적용하여 정확성과 효율성을 보장한다.

실험 결과

연구 질문

• RQ1 뉴턴 다면체의 구조를 활용함으로써 희소 다항식계에 대한 그로브너 기저 계산을 어떻게 가속화할 수 있는가?
• RQ2 정규 래티스 다면체 $\mathscr{P}$ 의 조합적 불변량에 따라 희소 그로브너 기저 계산의 복잡도 한계는 무엇인가?
• RQ3 기존의 이항계에 대한 차수 한계 $\min(n_1,n_2)+1$ 은 희소 그로브너 기저 이론을 통해 다항계로 일반화될 수 있는가?
• RQ4 실제 적용에서 희소 F5 및 FGLM 알고리즘의 변형이 기존 알고리즘보다 얼마나 뛰어나게 성능을 발휘하는가?
• RQ5 워치독 차수를 기반으로 한 수정된 프로베르그의 추측은 과잉 결정된 희소 시스템의 복잡도를 예측할 수 있는가?

주요 결과

• 희소-MatrixF5 및 희소-FGLM 알고리즘은 실제 적용에서 기존 F5에 비해 100배 이상의 가속도를 달성하며, 일부 예제에서는 FGb 대비 19,000배 이상의 향상을 보였다.
• 뉴턴 다면체 $\mathscr{P}$ 를 갖는 균일한 경우, 워치독 차수는 몫환의 Hilbert 급수에서 첫 번째 음수 계수에 의해 제한되며, 이는 이항계의 $\min(n_1,n_2)+1$ 를 다항계에서 $\sum n_i - \max(n_i) + 1$ 으로 일반화한다.
• 희소 그로브너 기저 계산의 복잡도 한계는 정규화된 부피 $\operatorname{vol}(\mathscr{P})$, 캐스텔누오보-무디드 정규성 $\operatorname{reg}(k[S^{(h)}]) = n+1 - \ell$, 그리고 에흐라르 다항식 $\operatorname{HP}_{\mathscr{P}}(\ell)$ 에 따라 결정된다.
• 제안된 프로베르그의 추측의 변형은 $\mathbb{C}[\mathscr{P}]$ 내의 준정규 수열이 Hilbert 급수 $\operatorname{HS}_{\mathbb{C}[\mathscr{P}]/\langle f_1,\ldots,f_m \rangle}(t) = \left[ \operatorname{HS}_{\mathscr{P}}(t) \prod_{i=1}^m (1 - t^{d_i}) \right]_+$ 를 갖는다고 예측하며, 워치독 차수는 첫 번째 0 계수의 인덱스로 제한된다.
• 희소 단항식 지지 집합을 갖는 소수항 다항식계에서는, 비정규 반군으로 인해 복잡도 분석이 적용되지 않더라도 희소 알고리즘이 정확하게 작동하며 최대 16,800배의 가속도를 달성한다.
• 희소 FGLM을 통해 확보된 유리 함수의 단변수 표현은 해를 다항식 $Q(T)$ 의 근으로 표현하며, 단항식 $X^\alpha$ 는 유리 함수 $Q_\alpha(T)$ 로 매핑되며, 이는 효율적인 해 복원을 가능하게 한다.