Fourier–Motzkin消去法

Fourier–Motzkin消去法（英: Fourier-Motzkin elimination）とは、数理論理学および計算機科学において、一次不等式からなる一階述語論理式の限量子（∀や∃）を除去するアルゴリズム。(限量子消去法)（英: Quantifier elimination）の1つ。1826年にジョゼフ・フーリエが発見し、1918年に L. L. Dines が再発見し、1936年に Theodore Motzkin が再々発見した^[1]。

アルゴリズム

以下の手順を繰り返し行い、限量子を除去していく^[2]。変数の定義域は実数もしくは有理数。

1. 以下の置き換えを行う。
   1. 等式や不等式に ￢ がかかる物は反転させる。
   2. a ≧ b は b ≦ a へ
   3. a ＞ b は b ＜ a へ
   4. a ＝ b は (a ≦ b) ∧ (b ≦ a) へ
   5. ∀x. P(x) は ￢ ∃x. ￢ P(x) へ
2. 下記簡略化は随時行う。
   1. 常に成り立つ もしくは 常に不成立な不等式は真や偽に置き換える。
   2. 真や偽が ∧ や ∨ や ￢ に付く場合は適切に論理式を簡略化する。
   3. ∃x. P(x) の形式において、P(x) が x を使用してなければ、∃x. を取り除く。
3. ∃x. P(x) の形式で、P(x) に限量子が出てこない物を探し、P(x) を選言標準形に変換する。
4. ∃x. (P(x) ∨ Q(x)) ⇔ (∃x. P(x)) ∨ (∃x. Q(x)) の変換を行う。
5. Q が x を使用していない場合、∃x. (P(x) ∧ Q) ⇔ (∃x. P(x)) ∧ Q の変換を行う。
6. 下記公式を使い、限量子を除去する。

   ${\displaystyle \exists x.\left(\bigwedge _{h}c_{h}\leq a_{h}x\right)\land \left(\bigwedge _{i}c_{i}<a_{i}x\right)\land \left(\bigwedge _{j}b_{j}x\leq d_{j}\right)\land \left(\bigwedge _{k}b_{k}x<d_{k}\right)}$ 

   ${\displaystyle \Leftrightarrow \left(\bigwedge _{h,j}b_{j}c_{h}\leq a_{h}d_{j}\right)\land \left(\bigwedge _{h,k}b_{k}c_{h}<a_{h}d_{k}\right)\land \left(\bigwedge _{i,j}b_{j}c_{i}<a_{i}d_{j}\right)\land \left(\bigwedge _{i,k}b_{k}c_{i}<a_{i}d_{k}\right)}$ 

上記は一般的な形式だが、より分かりやすく、2つの不等式の場合に書き下すと以下の通り。

${\displaystyle (\exists x.c\leq ax\land bx\leq d)\Leftrightarrow bc\leq ad}$ 

${\displaystyle (\exists x.c<ax\land bx\leq d)\Leftrightarrow bc<ad}$ 

${\displaystyle (\exists x.c\leq ax\land bx<d)\Leftrightarrow bc<ad}$ 

${\displaystyle (\exists x.c<ax\land bx<d)\Leftrightarrow bc<ad}$ 

具体例

${\displaystyle \forall i.(1\leq i\leq n\to i\neq n)}$  に対して、Fourier–Motzkin消去法を使用し、簡略化する。

${\displaystyle \forall i.(1\leq i\leq n\to i\neq n)}$  ${\displaystyle \Leftrightarrow \forall i.((\lnot (1\leq i\leq n))\lor (i\neq n))}$  ${\displaystyle \Leftrightarrow \lnot \exists i.\lnot ((\lnot (1\leq i\leq n))\lor (i\neq n))}$  ${\displaystyle \Leftrightarrow \lnot \exists i.(1\leq i\leq n)\land (i=n)}$  ${\displaystyle \Leftrightarrow \lnot \exists i.(1\leq i)\land (i\leq n)\land (i=n)}$  ${\displaystyle \Leftrightarrow \lnot \exists i.(1\leq i)\land (i\leq n)\land (i\leq n)\land (n\leq i)}$  ${\displaystyle \Leftrightarrow \lnot \exists i.(1\leq i)\land (i\leq n)\land (n\leq i)}$  ${\displaystyle \Leftrightarrow \lnot \exists i.(1\leq i)\land (n\leq i)\land (i\leq n)}$  ここで公式を使用する。 ${\displaystyle \Leftrightarrow \lnot \exists i.(1\leq n)\land (n\leq n)}$  ${\displaystyle \Leftrightarrow \lnot \exists i.(1\leq n)}$  ${\displaystyle \Leftrightarrow \lnot (1\leq n)}$  ${\displaystyle \Leftrightarrow 1>n}$  ${\displaystyle \Leftrightarrow n<1}$  

整数への拡張

変数の定義域を整数にする場合の拡張を William Pugh が1992年に発表している^[3]。オメガテストと名付けた判定条件を追加している。

また、1972年に D. C. Cooper がFourier–Motzkin消去法の直接の拡張ではないが、整数変数に対する限量子消去法である(Cooper法)を発表している^[4]。

多項式の場合

一次式ではなく、より一般的な多項式の場合 George E. Collins が1975年に発表した Cylindrical Algebraic Decomposition (CAD) を使用することにより、限量子消去が出来る。

参照

1. ^ Satisfiability Checking Fourier–Motzkin Variable Elimination
2. ^ Arithmetic Decision Procedures: a simple introduction
3. ^ The Omega Test: a fast and practical integer programming algorithm for dependence analysis
4. ^ Theorem Proving in Arithmetic without Multiplication