数学ビジュアライズ

AI数学単発記事

OpenAIのAIが証明した「ほぼ n⁻² は成り立たない」:ハイルブロンの三角形問題の新しい下界

正方形に n 点を置いたとき、3点が作るいちばん小さい三角形の面積をどこまで大きくできるかを問うのがハイルブロンの三角形問題です。自明な上界とエルデシュの放物線の構成を証明つきで確かめ、ハイルブロンの予想が1982年に否定された経緯と、今回OpenAIのモデルが「ほぼ n⁻²」という見方を否定した定理 Δ(n) ≥ c n^(−2+η) を、点を動かせる図とともに解説します。

目次
  1. 正方形の中の「いちばん小さい三角形」
  2. 問題:三角形の最小面積を大きくする
  3. 上からの評価:どう置いても 程度まで
  4. 下からの評価:エルデシュの放物線
  5. ハイルブロンの予想と「ほぼ 」という見方
  6. 今回の定理
  7. 証明の考え方
  8. どこまで確かなのか
  9. 残された問題
  10. 参考文献
この記事の要点
問題
正方形の中に nn 個の点を置き、どの3点を選んでも三角形の面積がなるべく大きくなるようにします。いちばん小さい三角形の面積を、どこまで大きくできるでしょうか。この最大値を Δ(n)\Delta(n) と書きます。
これまで
エルデシュの構成から Δ(n)\Delta(n) は 1/n21/n^2 の定数倍以上です。ハイルブロンは Δ(n)\Delta(n) が 1/n21/n^2 の定数倍以下だと予想しましたが、1982年にコムローシュ・ピンツ・セメレディが log⁡n\log n 倍の差で否定しました。それでも「1/n21/n^2 とほとんど同じ大きさだろう」と考える研究者は多くいました。
今回
OpenAIのモデルが書いた論文が、十分大きいすべての nn で Δ(n)≥c1n−2+η\Delta(n)\ge c_1n^{-2+\eta} となる定数 c1,η>0c_1,\eta>0 があることを示しました。log⁡n\log n ではなく nn のべきの差なので、「ほとんど同じ」という見方は成り立たなくなります。
ただし
η\eta の値は具体的に書けますが、約 3×10−2363\times10^{-236} ととても小さい数です。Δ(n)\Delta(n) の本当の大きさは、今も分かっていません。
確かさ
論文は査読前のプレプリントです。Leanでの形式化があり、リポジトリの照合の設定に載っています。ただし形式化されているのは、点の個数 nn が限りなく大きくなる列についての形です。

この記事で読む論文

A power improvement in the Heilbronn triangle lower bound

OpenAI 2026年9月25日付 プレプリント(査読前) PDF 433KB

PDFは2026年10月9日にGitHubから取得した版のコピーで、Apache License 2.0(全文)に従って掲載しています。修正版が出た場合は、GitHubのほうが新しくなります。

定義・命題の下の「証明」は、開くとステップが並びます。もう少し詳しくは大学1〜2年程度の数学の補足、専門的な話は論文の中身に踏み込む補足です。どちらも開かずに読み進められます。

正方形の中の「いちばん小さい三角形」

2026年10月6日、OpenAIは未公開の社内モデルが書いた数学の論文をGitHubのリポジトリ openai/math で公開しました[1]。その191番が、この記事で読む「A power improvement in the Heilbronn triangle lower bound」です[2]。

扱っているのは、ハイルブロンの三角形問題と呼ばれる離散幾何の問題です。正方形に点をばらまくと、そのうち3点が作る三角形の中には、とても小さいものが必ず混じります。点の置き方を工夫して、いちばん小さい三角形をどこまで大きくできるかを問います。

問題文は中学生でも分かりますが、答えの大きさは70年以上たっても決まっていません。この記事では、まず問題を定義し、上からの評価と下からの評価を証明つきで確かめます。そのうえで、これまでの予想のどれが否定されていて、今回は何が新しく示されたのかを整理します。

この論文は、10月7日付の修正・取り下げには含まれていません[3]。リポジトリの更新は修正と取り下げの記録にまとめています。

問題:三角形の最小面積を大きくする

平面の3点 p,q,rp,q,r が作る三角形の面積は、2つのベクトル q−pq-p と r−pr-p を並べた行列式で書けます。

Area⁡(pqr)=12∣det⁡(q−p, r−p)∣=12∣(q1−p1)(r2−p2)−(q2−p2)(r1−p1)∣\begin{aligned} \operatorname{Area}(pqr)&=\frac12\left|\det(q-p,\ r-p)\right|\\ &=\frac12\left|(q_1-p_1)(r_2-p_2)\right.\\ &\qquad\left.-(q_2-p_2)(r_1-p_1)\right|\end{aligned}

3点が一直線に並ぶと、この値は 00 です。

定義1 最小三角形の面積とハイルブロンの関数 [2, §1]

単位正方形 [0,1]2[0,1]^2 の中の、3個以上の点からなる有限集合 PP について、3つの異なる点が作る三角形の面積の最小値を Δ(P)\Delta(P) と書く。

Δ(P)=min⁡{p,q,r}⊆PArea⁡(pqr)\Delta(P)=\min_{\{p,q,r\}\subseteq P}\operatorname{Area}(pqr)Δ(P)=min⁡{p,q,r}⊆PArea⁡(pqr)\begin{aligned}\Delta(P) &= \min_{\{p,q,r\}\subseteq P}\operatorname{Area}(pqr)\end{aligned}

n≥3n\ge3 について、nn 点の集合 P⊆[0,1]2P\subseteq[0,1]^2 全体での Δ(P)\Delta(P) の最大値を Δ(n)\Delta(n) と書く。

Δ(n)=max⁡P⊆[0,1]2, ∣P∣=nΔ(P)\Delta(n)=\max_{P\subseteq[0,1]^2,\ |P|=n}\Delta(P)Δ(n)=max⁡P⊆[0,1]2, ∣P∣=nΔ(P)\begin{aligned}\Delta(n) &= \max_{P\subseteq[0,1]^2,\ |P|=n}\Delta(P)\end{aligned}

Δ(P)\Delta(P) は「PP の中でいちばん小さい三角形の面積」です。3点が一直線に並ぶ組が1つでもあれば Δ(P)=0\Delta(P)=0 になります。Δ(n)\Delta(n) は、nn 点をうまく置いたときに、その最小面積をどこまで大きくできるかを表します。

最大値が存在することは、論文の序論で確かめられています[2, §1]。点の重なりを許して nn 点の組を考えると、その全体は有界な閉集合です。組を動かすと最小面積は連続的に変わるので、この全体の中に最小面積が最大になる組があります。

重なる点を含む組では、その点を2回使う三角形の面積が 00 なので、最小面積は 00 です。一方、放物線の上の nn 点のように最小面積が正になる配置があります。よって最大になる組は重なりのない nn 点で、これが Δ(n)\Delta(n) を与えます。

例2 四隅の4点と格子の9点

正方形の四隅 (0,0),(1,0),(1,1),(0,1)(0,0),(1,0),(1,1),(0,1) から3点を選ぶと、どれも直角をはさむ2辺が 11 の直角三角形になる。面積はすべて 1/21/2 なので、Δ(P)=1/2\Delta(P)=1/2 である。

一方、3×33\times3 の格子点 {1/6,1/2,5/6}2\{1/6,1/2,5/6\}^2 では、横一列の3点が一直線に並ぶので Δ(P)=0\Delta(P)=0 である。

格子のように規則正しく並べると、一直線に並ぶ3点ができてしまいます。かといって点がかたまると、近くの3点で小さい三角形ができます。下の図で、点を動かして確かめてみてください。

00
11
11

点の数 n=n=7 いちばん小さい三角形の面積 0.01020 n2×n^2\times{}面積 0.500

図1具体例単位正方形の中の点と、そのうちいちばん小さい三角形(橙)です。点はドラッグで動かせます。正方形の空いているところを押すと点が増え、「最後の点を消す」で減ります。3点が一直線に並ぶと、赤い線で示します。「エルデシュの7点」は次の節で扱う配置で、n2×n^2\times{}面積はちょうど 0.50.5 です。

点を1つ動かして最小の三角形を大きくすると、たいてい別の三角形が小さくなります。点が増えるほどこの兼ね合いは厳しくなり、最小面積は小さくなっていきます。問題は、nn を大きくしたときに Δ(n)\Delta(n) がどのくらいの速さで小さくなるかです。

上からの評価:どう置いても 1/n1/n 程度まで

まず、どんなに工夫しても最小面積は 1/n1/n 程度より大きくならないことを確かめます。使うのは鳩の巣原理だけです。

命題3 自明な上界

n≥3n\ge3 なら Δ(n)≤1n−2\Delta(n)\le\dfrac{1}{n-2} である。

証明
  1. m=⌊(n−1)/2⌋m=\lfloor(n-1)/2\rfloor とおく。n≥3n\ge3 なので m≥1m\ge1 で、n≥2m+1n\ge2m+1、2m≥n−22m\ge n-2 が成り立つ。
  2. 正方形を幅 1/m1/m の縦長の帯 mm 本に分ける。n≥2m+1n\ge2m+1 個の点を mm 本の帯に入れるので、3点以上を含む帯がある(鳩の巣原理)。
  3. 幅 ww、高さ 11 の長方形に入る三角形の面積は w/2w/2 以下である。

    3つの頂点のうち、xx 座標が真ん中のものを通る縦の線を引く。この線は三角形を2つの三角形に分け、どちらも縦の線分(長さ ℓ\ell)を共通の辺にもつ。

    線分は正方形の中にあるので ℓ≤1\ell\le1 である。2つの三角形の高さ(横方向の距離)を h1,h2h_1,h_2 とすると h1+h2≤wh_1+h_2\le w なので、面積の和は次のようになる。

    ℓh12+ℓh22=ℓ(h1+h2)2≤w2\frac{\ell h_1}{2}+\frac{\ell h_2}{2}=\frac{\ell(h_1+h_2)}{2}\le\frac{w}{2}ℓh12+ℓh22=ℓ(h1+h2)2≤w2\begin{aligned}\frac{\ell h_1}{2}+\frac{\ell h_2}{2} &= \frac{\ell(h_1+h_2)}{2}\\&\le \frac{w}{2}\end{aligned}
  4. 3点を含む帯の中の三角形の面積は 1/(2m)1/(2m) 以下なので、Δ(P)≤1/(2m)≤1/(n−2)\Delta(P)\le1/(2m)\le1/(n-2) である。PP は任意なので Δ(n)≤1/(n−2)\Delta(n)\le1/(n-2) が従う。

たとえば n=101n=101 なら m=50m=50 で、どの101点の置き方にも面積 1/1001/100 以下の三角形があります。

この評価は粗いものです。帯の中の3点がたまたま帯の幅いっぱいに広がっているとは限らないからです。この上界をどこまで下げられるかも、この問題で問われていることの一つです。

下からの評価:エルデシュの放物線

反対に、最小面積が大きい置き方を作れば、Δ(n)\Delta(n) を下から抑えられます。よく知られた構成は、素数 pp で割った余りの世界の放物線を使います。

この構成はハイルブロンの問題を扱ったロスの1951年の論文の付録に、エルデシュによるものとして記録されています[2, §1][4]。ザハロフの解説は「おそらくエルデシュによる」と書いています[5, §1]。

命題4 エルデシュの構成 [5, §1]

pp を奇素数とし、x=0,1,…,p−1x=0,1,\dots,p-1 について x2x^2 を pp で割った余りを x2 mod px^2\bmod p と書く。

Pp={(xp, x2 mod pp):(x2 mod pp)x=0,1,…,p−1}\begin{aligned} P_p&=\left\{\left(\frac{x}{p},\ \frac{x^2\bmod p}{p}\right):\right.\\ &\qquad\left.\vphantom{\left(\frac{x^2\bmod p}{p}\right)}x=0,1,\dots,p-1\right\} \end{aligned}

は単位正方形の中の pp 点で、Δ(Pp)≥12p2\Delta(P_p)\ge\dfrac{1}{2p^2} をみたす。とくに Δ(p)≥12p2\Delta(p)\ge\dfrac{1}{2p^2} である。

証明
  1. pp 倍した点 (x, yx)(x,\ y_x)(yx=x2 mod py_x=x^2\bmod p)で考える。面積は p2p^2 倍になるので、どの三角形も面積 1/21/2 以上であることを示せばよい。
  2. 3点 x=a,b,cx=a,b,c(互いに異なる)の行列式 D=(b−a)(yc−ya)−(yb−ya)(c−a)D=(b-a)(y_c-y_a)-(y_b-y_a)(c-a) は整数で、pp で割ると (b−a)(c−a)(c−b)(b-a)(c-a)(c-b) と同じ余りになる。

    yx≡x2(modp)y_x\equiv x^2\pmod p なので、次の合同式が成り立つ。

    D≡(b−a)(c2−a2)−(b2−a2)(c−a)=(b−a)(c−a)((c+a)−(b+a))=(b−a)(c−a)(c−b)(modp)\begin{aligned} D&\equiv(b-a)(c^2-a^2)\\ &\qquad-(b^2-a^2)(c-a)\\ &=(b-a)(c-a)\\ &\qquad\bigl((c+a)-(b+a)\bigr)\\ &=(b-a)(c-a)(c-b)\\ &\qquad\pmod p\end{aligned}

    これは3次のヴァンデルモンドの行列式と同じ式である。

  3. a,b,ca,b,c は 00 以上 p−1p-1 以下の異なる整数なので、b−a, c−a, c−bb-a,\ c-a,\ c-b はどれも pp の倍数でない。pp は素数なので積も pp の倍数でなく、D≠0D\ne0 である。
  4. DD は 00 でない整数なので ∣D∣≥1|D|\ge1 で、三角形の面積 ∣D∣/2|D|/2 は 1/21/2 以上である。

証明の要点は、「行列式が整数」と「pp で割った余りが 00 でない」の2つを組み合わせたことです。余りの世界では3点が一直線に並ばないので、行列式は 00 にならず、整数なので 11 以上になります。今回の論文も、この2つの組み合わせを大きく発展させた形になっています。

00
11
11

p=p=11:三角形は全部で 165 個。いちばん小さい面積は格子の目で 1/21/2 で、それが 6 個。単位正方形では 1/(2p2)=1/1/(2p^2)=1/242 = \,=\,0.00413

図2具体例命題4の点を p×pp\times p の格子に打った図です。点は見やすいようにマスの中央に置いています(平行移動しても面積は変わりません)。緑の線は放物線 y=x2y=x^2 を pp ずつ下にずらして正方形に折り込んだもので、点はすべてこの線の上にあります。橙はいちばん小さい三角形です。p=5,…,23p=5,\dots,23 のどれでも、いちばん小さい三角形の面積はマスの目でちょうど 1/21/2 です。

pp が素数でない nn についても、ベルトランの仮説を使えば同じ程度の評価が得られます。n≤p≤2nn\le p\le2n となる素数 pp をとり、PpP_p から nn 点を選びます。点を減らしても最小面積は小さくならないので、次の評価が成り立ちます。

Δ(n)≥Δ(Pp)≥12p2≥18n2\Delta(n)\ge\Delta(P_p)\ge\frac{1}{2p^2}\ge\frac{1}{8n^2}Δ(n)≥Δ(Pp)≥12p2≥18n2\begin{aligned}\Delta(n) &\ge \Delta(P_p)\\&\ge \frac{1}{2p^2}\\&\ge \frac{1}{8n^2}\end{aligned}

命題3と合わせると、Δ(n)\Delta(n) は次の範囲にあります。

18n2≤Δ(n)≤1n−2(n≥3)\frac{1}{8n^2}\le\Delta(n)\le\frac{1}{n-2}\qquad(n\ge3)18n2≤Δ(n)≤1n−2(n≥3)\begin{gathered}\frac{1}{8n^2}\le\Delta(n)\le\frac{1}{n-2}\\(n\ge3)\end{gathered}18n2≤Δ(n)≤1n−2(n≥3)\begin{aligned}\frac{1}{8n^2} &\le \Delta(n)\\&\le \frac{1}{n-2}\\[4pt]&(n\ge3)\end{aligned}

1/n21/n^2 と 1/n1/n のあいだには大きな開きがあります。本当の Δ(n)\Delta(n) はどちらに近いのか、というのがこの問題の中心です。

ランダムに置いてから消す

同じ 1/n21/n^2 の大きさは、点をでたらめに置く方法でも得られます。この方法も、ザハロフの解説ではエルデシュによるものと推測されています[5, §1]。今回の論文の最後の段階もこの形をしているので、ここで確かめておきます。

命題5 削除による構成

mm 個の点の集合 QQ と A>0A>0 について、面積が AA 以下の三角形が ZZ 個あるとする。Z≤m−nZ\le m-n なら、QQ から nn 個の点を選んで、どの三角形も面積が AA より大きいようにできる。

証明
  1. 面積が AA 以下の三角形のそれぞれから、頂点を1つずつ選んで消す。消す点は ZZ 個以下なので、m−Z≥nm-Z\ge n 個以上の点が残る。
  2. 残った点の3点が作る三角形は、もとの「面積 AA 以下の三角形」のどれでもない。それらは頂点を1つ失っているからである。よって残った三角形の面積はすべて AA より大きく、残りから nn 点を選べばよい。

単位正方形に一様にでたらめに置いた3点が、面積 AA 以下の三角形を作る確率は、ある定数 CC について CACA 以下です[5, §2]。期待値の線形性から、mm 点のうち面積 AA 以下の三角形の個数の期待値は CA(m3)CA\binom m3 以下になります。

m=2nm=2n、A=c/n2A=c/n^2(cc は小さい定数)とすると期待値は nn 以下になり、期待値以下になる置き方が実際にあります。それに命題5を使うと、Δ(n)≥c/n2\Delta(n)\ge c/n^2 が得られます。

00
11
11

2020 点のうち、面積 1/10001/1000 以下の三角形が 13 個。消した点は 6 個で、残り 14 点のいちばん小さい三角形の面積は 0.00172

図3具体例命題5を、ランダムに置いた 2020 点、A=1/1000A=1/1000 で試した図です。「小さい三角形」で面積 1/10001/1000 以下の三角形を赤で示し、「消した後」でそれぞれの三角形から頂点を1つ消した結果を示します。点線の丸が消した点、青い三角形が残った点のうちいちばん小さい三角形です。最初に表示している配置では、13個の小さい三角形をなくすために6点を消し、残った14点の最小面積は 0.001720.00172 です。
もう少し詳しく3点の三角形が小さくなる確率が CACA 以下である理由

ザハロフの解説の計算をなぞります[5, §2]。三角形の面積は、2点 pi,pjp_i,p_j の距離 dd と、3点目 pkp_k から直線 pipjp_ip_j までの距離 ee で de/2de/2 と書けます。

dd が tt 以上 2t2t 未満になる確率は t2t^2 の定数倍程度です。そのとき面積が AA 以下になるには、pkp_k が直線のまわりの幅 2A/t2A/t 程度の帯に入る必要があり、その確率は min⁡(A/t,1)\min(A/t,1) の定数倍程度です。

2点の距離は 2<2\sqrt2<2 以下なので、t=1,1/2,1/4,…t=1,1/2,1/4,\dots の区間 [t,2t)[t,2t) で距離の全範囲を覆えます。これらについて足し合わせると、次のように AA の定数倍で抑えられます。

Pr⁡[Area⁡(pipjpk)≤A]≲∑j≥0tj2min⁡ ⁣(Atj,1)≤∑j≥0tjA=2A(tj=2−j)\begin{aligned} &\Pr\bigl[\operatorname{Area}(p_ip_jp_k)\le A\bigr]\\ &\quad\lesssim\sum_{j\ge0}t_j^2\min\!\left(\frac A{t_j},1\right)\\ &\quad\le\sum_{j\ge0}t_jA=2A\\ &\qquad(t_j=2^{-j})\end{aligned}

ハイルブロンの予想と「ほぼ n−2n^{-2}」という見方

ここまでで c/n2≤Δ(n)≤C/nc/n^2\le\Delta(n)\le C/n が分かりました。ここからは、どの予想がいつ誰によって否定されたのかを、順に整理します。

ハイルブロンの予想(1982年に否定)

この問題を考えたのはハイルブロンで、1940年代の終わりごろのことです[5, §1]。ハイルブロンは、下からの評価 1/n21/n^2 が正しい大きさだと考えました。

予想6 ハイルブロン [2, §1][5, §1]

ある定数 CC があって、すべての n≥3n\ge3 について次が成り立つ。

Δ(n)≤Cn2\Delta(n)\le\frac{C}{n^2}Δ(n)≤Cn2\begin{aligned}\Delta(n) &\le \frac{C}{n^2}\end{aligned}

論文はこの予想の出典として、ロスの1951年の論文の198ページを挙げています[4]。

この予想は1982年に、コムローシュ、ピンツ、セメレディによって否定されました。ある定数 c>0c>0 について次が成り立つことを示したのです[6][5, §2]。

Δ(n)≥clog⁡nn2\Delta(n)\ge\frac{c\log n}{n^2}Δ(n)≥clog⁡nn2\begin{aligned}\Delta(n) &\ge \frac{c\log n}{n^2}\end{aligned}

log⁡n\log n はいくらでも大きくなるので、C/n2C/n^2 以下にはおさまりません。証明では、まずランダムに点を多めに置き、小さい三角形をハイパーグラフの辺とみなします。そのうえで、局所的にまばらなハイパーグラフについての強い独立集合の評価を使って、小さい三角形を含まない大きな部分集合を選び出します[2, §1][5, §2]。

「ほぼ n−2n^{-2}」という見方

元の予想は否定されましたが、外れたのは log⁡n\log n 倍だけです。そこで、1/n21/n^2 に「ほぼ」等しいという弱い形なら正しいだろう、と考えられてきました。

予想7 ほぼ n−2n^{-2} の上界 [2, 式 (1.1)][5, §3]

どの ε>0\varepsilon>0 についても定数 CεC_\varepsilon と n0(ε)n_0(\varepsilon) があって、次が成り立つ。

Δ(n)≤Cε n−2+ε(n≥n0(ε))\Delta(n)\le C_\varepsilon\,n^{-2+\varepsilon}\qquad(n\ge n_0(\varepsilon))Δ(n)≤Cε n−2+ε(n≥n0(ε))\begin{gathered}\Delta(n)\le C_\varepsilon\,n^{-2+\varepsilon}\\(n\ge n_0(\varepsilon))\end{gathered}Δ(n)≤Cε n−2+ε(n≥n0(ε))\begin{aligned}\Delta(n) &\le C_\varepsilon\,n^{-2+\varepsilon}\\[4pt]&(n\ge n_0(\varepsilon))\end{aligned}

log⁡n\log n は nεn^{\varepsilon} よりゆっくり増えるので、コムローシュ・ピンツ・セメレディの結果とは矛盾しません。ザハロフの2026年の解説は、セメレディを含む多くの研究者が、これが正しいと考えていると書いています[5, §3]。特定の一人の名前がついた予想ではないので、この記事では「ほぼ n−2n^{-2} という見方」とも呼びます。

今回の論文が否定したのは、この見方です[2, §1]。ハイルブロンの元の予想は、1982年の時点ですでに否定されていました。

上からの評価の進歩

上からの評価も、1/n1/n から少しずつ下がってきました。下の表は、論文の序論とザハロフの解説の表によります[2, §1][5, §6]。

年示した人上からの評価
—自明な評価(命題3)Δ(n)≤C n−1\Delta(n)\le C\,n^{-1}
1951ロス[4]Δ(n)=o(n−1)\Delta(n)=o(n^{-1})
1972シュミット[7]Δ(n)≤C n−1(log⁡n)−1/2\Delta(n)\le C\,n^{-1}(\log n)^{-1/2}
1972ロス[8,9]Δ(n)≤C n−1−μ\Delta(n)\le C\,n^{-1-\mu}(ザハロフの表では指数が約 −1.117-1.117)
1981コムローシュ、ピンツ、セメレディ[10]Δ(n)≤Cε n−8/7+ε\Delta(n)\le C_\varepsilon\,n^{-8/7+\varepsilon}
2023コーエン、ポホアタ、ザハロフ[11]十分大きい nn で Δ(n)≤n−8/7−1/2000\Delta(n)\le n^{-8/7-1/2000}
2025コーエン、ポホアタ、ザハロフ[12]Δ(n)≤Cε n−7/6+ε\Delta(n)\le C_\varepsilon\,n^{-7/6+\varepsilon}

Δ(n)\Delta(n) を nan^a の形で比べたときの指数 aa で、分かっていることを数直線に描くと次のようになります。

−1-1
自明な上界
−87-\tfrac87
KPS 1981
−76-\tfrac76
CPZ 2025
−2-2
エルデシュ・KPS 1982
今回:−2+η-2+\eta
本当の指数がありうる範囲(今回で −2-2 そのものは除かれた)
aa
Δ(n)\Delta(n) を nan^{a} と比べたときの指数 aa
図4概念図Δ(n)\Delta(n) を nan^a と比べたときの指数 aa を横軸にとった図です。青い印は上からの評価で、−1-1 から −7/6-7/6 まで左へ動いてきました。橙の印は下からの評価です。エルデシュの構成とコムローシュ・ピンツ・セメレディの評価はどちらも指数 −2-2 のところにあり(後者は log⁡n\log n 倍だけ大きい)、今回の結果は指数を −2+η-2+\eta に上げました。η\eta は約 3×10−2363\times10^{-236} なので、この図の縮尺では −2-2 と区別できません。

ザハロフは、指数 7/67/6 は今の手法の自然な限界で、予想される n−2+o(1)n^{-2+o(1)} とはまだ大きく離れていると書いていました[5, §6]。今回の結果で、その「予想される大きさ」のほうが成り立たないことになります。

今回の定理

定理8 ハイルブロンの問題の下界のべきの改良 [2, Theorem 1.1]

絶対定数 η,c1>0\eta,c_1>0 と整数 n0≥3n_0\ge3 があって、n≥n0n\ge n_0 をみたすすべての整数 nn について次が成り立つ。

Δ(n)≥c1 n−2+η\Delta(n)\ge c_1\,n^{-2+\eta}Δ(n)≥c1 n−2+η\begin{aligned}\Delta(n) &\ge c_1\,n^{-2+\eta}\end{aligned}

「絶対定数」は、nn にも何にもよらない決まった数という意味です。つまり、十分大きいどの nn についても、最小の三角形の面積が c1n−2+ηc_1n^{-2+\eta} 以上になる nn 点の置き方があります。

系9 ほぼ n−2n^{-2} の上界は成り立たない [2, §1]

予想7は成り立たない。つまり Δ(n)\Delta(n) は、ほぼ n−2n^{-2} の上界をもたない。

証明
  1. ε=η/2\varepsilon=\eta/2 とし、Δ(n)≤C n−2+η/2\Delta(n)\le C\,n^{-2+\eta/2}(n≥n0′n\ge n_0')となる定数 CC と n0′n_0' があったと仮定して矛盾を導く。
  2. 定理8と合わせると、十分大きい nn で c1n−2+η≤C n−2+η/2c_1n^{-2+\eta}\le C\,n^{-2+\eta/2}、つまり nη/2≤C/c1n^{\eta/2}\le C/c_1 となる。
  3. η>0\eta>0 なので nη/2n^{\eta/2} は n→∞n\to\infty でいくらでも大きくなり、これは成り立たない。

η\eta はどれくらいか

論文は η\eta を具体的な式で与えています[2, §8, 式 (8.7)]。使う数は次のとおりです[2, §3, 式 (3.1)]。

d=41M=(4d−1d)=(16341)T=(M3)k=T2+1η=245435k+16d=41\qquad M=\binom{4d-1}{d}=\binom{163}{41}\qquad T=\binom{M}{3}\qquad k=T^2+1\qquad \eta=\frac{2}{45435k+16}d=41M=(4d−1d)=(16341)T=(M3)k=T2+1η=245435k+16\begin{gathered}d=41\\M=\binom{4d-1}{d}=\binom{163}{41}\\T=\binom{M}{3}\\k=T^2+1\\\eta=\frac{2}{45435k+16}\end{gathered}d=41M=(4d−1d)=(16341)T=(M3)k=T2+1η=245435k+16\begin{aligned}d &= 41\\[4pt]M &= \binom{4d-1}{d}\\&= \binom{163}{41}\\[4pt]T &= \binom{M}{3}\\[4pt]k &= T^2+1\\[4pt]\eta &= \frac{2}{45435k+16}\end{aligned}

計算すると、MM は39桁、TT は116桁、kk は232桁の整数です。45435k+1645435k+16 は約 6.3×102356.3\times10^{235} なので、η\eta は約 3.2×10−2363.2\times10^{-236} になります。論文自身も、この指数は「固定されているが極端に小さい」と書いています[2, Abstract]。

どれくらい小さいかを、コムローシュ・ピンツ・セメレディの log⁡n\log n 倍と比べてみます。定数の違いを無視して nηn^{\eta} と log⁡n\log n を比べると、nηn^\eta のほうが大きくなるのは、log⁡n\log n がおよそ 1.7×102381.7\times10^{238} を超えてからです。

ですから、この定理は具体的な大きさの nn で良い配置を与えるものではありません。示したのは、n→∞n\to\infty での指数が −2-2 ではない、という漸近的な事実です。

証明の考え方

証明は論文の24ページの大部分を占め、この記事では再現しません。ここでは、エルデシュの構成の考え方「行列式が整数」と「余りが 00 でない」を、論文がどう発展させているかを見ます。

点を整数の3つ組で表す

論文はまず、点を整数の3つ組(列ベクトル)u=(u1,u2,u3)u=(u_1,u_2,u_3) で表します。NN を大きな整数として、次の箱の中の整数の点を使います[2, §1.1, §6]。

0≤u1,u2<NN≤u3<2Nπ(u)=(u1u3, u2u3)∈[0,1)20\le u_1,u_2<N\qquad N\le u_3<2N\qquad \pi(u)=\left(\frac{u_1}{u_3},\ \frac{u_2}{u_3}\right)\in[0,1)^20≤u1,u2<NN≤u3<2Nπ(u)=(u1u3, u2u3)∈[0,1)2\begin{gathered}0\le u_1,u_2<N\\N\le u_3<2N\\\pi(u)=\left(\frac{u_1}{u_3},\ \frac{u_2}{u_3}\right)\in[0,1)^2\end{gathered}0≤u1,u2<NN≤u3<2Nπ(u)=(u1u3, u2u3)∈[0,1)2\begin{aligned}0 &\le u_1,u_2\\&< N\\[4pt]N &\le u_3\\&< 2N\\[4pt]\pi(u) &= \left(\frac{u_1}{u_3},\ \frac{u_2}{u_3}\right)\in[0,1)^2\end{aligned}

3つの列を並べた 3×33\times3 の整数行列 AA の行列式で、対応する三角形の面積が書けます。

命題10 列ベクトルと三角形の面積 [2, 式 (1.2), (8.3)]

u(1),u(2),u(3)∈Z3u^{(1)},u^{(2)},u^{(3)}\in\mathbb{Z}^3 の第3成分が正のとき、これらを列にもつ行列を AA とすると、次が成り立つ。

Area⁡(π(u(1)) π(u(2)) π(u(3)))=∣det⁡A∣2 u3(1)u3(2)u3(3)\begin{aligned} &\operatorname{Area}\bigl(\pi(u^{(1)})\,\pi(u^{(2)})\,\pi(u^{(3)})\bigr)\\ &\quad=\frac{|\det A|}{2\,u^{(1)}_3u^{(2)}_3u^{(3)}_3}\end{aligned}
証明
  1. 各列 u(j)u^{(j)} をその第3成分 u3(j)u^{(j)}_3 で割ると、列は (xj,yj,1)(x_j,y_j,1) の形になる。ここで (xj,yj)=π(u(j))(x_j,y_j)=\pi(u^{(j)}) である。
  2. 列を正の数で割ると行列式も同じ数で割られるので、割った後の行列式は det⁡A/(u3(1)u3(2)u3(3))\det A/(u^{(1)}_3u^{(2)}_3u^{(3)}_3) である。
  3. 最後の行がすべて 11 の行列式は、三角形の面積の2倍に符号をつけたものである。

    第1列と第2列から第3列を引いても行列式は変わらない。そのあと第3行で展開すると、2次の行列式になる。

    det⁡(x1x2x3y1y2y3111)=det⁡(x1−x3x2−x3y1−y3y2−y3)\begin{aligned} &\det\begin{pmatrix}x_1&x_2&x_3\\y_1&y_2&y_3\\1&1&1\end{pmatrix}\\ &\qquad=\det\begin{pmatrix}x_1-x_3&x_2-x_3\\y_1-y_3&y_2-y_3\end{pmatrix}\end{aligned}

    右辺の絶対値の半分が、3点 (xj,yj)(x_j,y_j) の三角形の面積である。

例11 N=5N=5 の箱の3列

列 (1,2,5), (4,1,6), (3,4,7)(1,2,5),\ (4,1,6),\ (3,4,7) は N=5N=5 の箱に入り、点 (1/5,2/5), (2/3,1/6), (3/7,4/7)(1/5,2/5),\ (2/3,1/6),\ (3/7,4/7) を表す。

det⁡(143214567)=1⋅(7−24)−4⋅(14−20)=+3⋅(12−5)=28\begin{aligned} &\det\begin{pmatrix}1&4&3\\2&1&4\\5&6&7\end{pmatrix}\\ &\quad=1\cdot(7-24)-4\cdot(14-20)\\ &\quad\phantom{{}={}}+3\cdot(12-5)\\ &\quad=28\end{aligned} Area⁡=282⋅5⋅6⋅7=115\operatorname{Area}=\frac{28}{2\cdot5\cdot6\cdot7}=\frac{1}{15}Area⁡=282⋅5⋅6⋅7=115\begin{aligned}\operatorname{Area} &= \frac{28}{2\cdot5\cdot6\cdot7}\\&= \frac{1}{15}\end{aligned}

3点から直接計算しても、面積は 1/151/15 になる。

エルデシュの構成は、第3成分をすべて pp にした場合にあたります。点 (x/p, yx/p)(x/p,\ y_x/p) を列 (x, yx, p)(x,\ y_x,\ p) で表すと、det⁡A=p⋅D\det A=p\cdot D(DD は命題4の証明の行列式)となり、面積は p∣D∣/(2p3)=∣D∣/(2p2)p|D|/(2p^3)=|D|/(2p^2) です。

命題10から、どの3列についても ∣det⁡A∣>τ|\det A|>\tau なら、第3成分が 2N2N 未満なので、すべての三角形の面積が τ/(16N3)\tau/(16N^3) 以上になります[2, 式 (8.4)]。エルデシュの構成は「行列式が 00 でない、つまり 11 以上」でした。今回はこれを「行列式が大きな数 τ\tau より大きい」に強めます。

ラベルで行列式を小さくさせない

行列式を大きくする仕掛けが、論文の第3節です。各列に、元が rdr^d 個の有限体 Frd\mathbb{F}_{r^d} の元をラベル ξ\xi として付けます(rr は大きくしていく素数、d=41d=41)。

ラベルが異なる3列について、有限体の放物線 (1,ξ,ξ2)(1,\xi,\xi^2) の行列式(ヴァンデルモンドの行列式)は 00 になりません。ここはエルデシュの構成と同じ考え方で、論文もそう述べています[2, §1.1]。

この「00 でない」という情報を、体のノルムで Fr\mathbb{F}_r の元に移し、さらに BB 進法の1つの桁に書き込みます。BB は大きな素数で、h=Bkh=B^k で割った余りを考えます。その結果、次の補題が得られます。

補題12 行列式の障害 [2, Lemma 3.3]

第3節の作り方で作った3つの列(成分は h=Bkh=B^k で割った余り)を並べた行列を CC とする。3つのラベルが互いに異なるなら、det⁡C\det C の hh で割った余りは、−τ-\tau 以上 τ\tau 以下の整数で表せない。ここで τ=⌊Bk−1/2⌋\tau=\lfloor B^{k-1}/2\rfloor である。

証明の骨組みは次のとおりです。各成分を BB 進法の桁 civc_{iv} で ci=∑vcivBvc_i=\sum_v c_{iv}B^v と書き、行列式を BB の多項式とみて展開します。桁の選び方を工夫して、Bk−1B^{k-1} の係数 γk−1\gamma_{k-1} が、ラベルの体のノルムと rr で割った余りで一致するようにしてあります。

ラベルが異なれば体のノルムは 00 でないので、γk−1≠0\gamma_{k-1}\ne0、つまり ∣γk−1∣≥1|\gamma_{k-1}|\ge1 です。一方、低い桁からの寄与は Bk−1B^{k-1} の 3/503/50 倍以下におさえられています。そのため、行列式の代表 SCS_C は絶対値が (47/50)Bk−1(47/50)B^{k-1} 以上、h/2h/2 未満の範囲に入り、τ\tau 以下の範囲から外れます[2, Lemma 3.3 の証明]。

00
τ\tau
−τ-\tau
4750Bk−1\tfrac{47}{50}B^{k-1}
−4750Bk−1-\tfrac{47}{50}B^{k-1}
h/2h/2
−h/2-h/2
小さい行列式 ∣t∣≤τ|t|\le\tau
SCS_C が入る範囲
SCS_C が入る範囲
(h=Bkh=B^k で割った余りの並び。長さは実際の比ではない)
図5概念図補題12の結論を、hh で割った余りを −h/2-h/2 から h/2h/2 まで並べた数直線に描いた図です。赤い範囲が「小さい行列式」∣t∣≤τ|t|\le\tau、緑の範囲が、ラベルが異なるときに行列式の代表 SCS_C が入る範囲です。長さは実際の比ではなく、並び順だけが意味をもちます。

このあと列を行列式 11 の行列でまぜ、余りを保ったまま整数の列に持ち上げます。どちらの操作でも行列式の hh で割った余りは det⁡C\det C のままなので、ラベルが異なる3列では ∣det⁡A∣>τ|\det A|>\tau です[2, §3, §6]。小さい行列式が起こりうるのは、3列のうち2つのラベルが一致したときだけになります。

全体の流れ

証明全体は、次の6段階に分かれます[2, §3–§8]。

  1. 1
    パラメータを決める

    d=41d=41 と kk を固定し、大きくしていく素数 rr に合わせて、素数 BB、h=Bkh=B^k、別の素数 qq、箱の大きさ N=(hq)10N=(hq)^{10} を選ぶ(§3、§5、§6)。

  2. 2
    列をランダムに作る

    各列に、ラベルと桁から決まる hh で割った余り(共通のランダムな行列式 11 の行列でまぜる)と、qq で割った余り(一直線に3点が並ばない集合から選ぶ)を与え、箱の中の整数の列に一様に持ち上げる(§4〜§6)。

  3. 3
    悪い組は少ないと数える

    同じ点を表す2列と、行列式の絶対値が τ\tau 以下の3列を「悪い組」と呼ぶ。悪い3列はラベルの一致を必要とし、その確率は 3r−d3r^{-d} 以下である。行列式がちょうど 00 の場合は qq での条件で別に数える(§6、§7)。

  4. 4
    悪い組を消す

    2nr2n_r 個の列を作ると、悪い組の個数 ZZ の期待値は nrn_r よりずっと小さい。Z<nrZ\lt n_r となる結果を選び、悪い組ごとに1列ずつ消して nrn_r 列を残す(§8。命題5と同じ形)。

  5. 5
    面積に戻す

    残った列はどの3つも ∣det⁡A∣>τ|\det A|>\tau なので、命題10から最小面積は τ/(16N3)\tau/(16N^3) 以上である。点の個数と合わせると nr2Δ(Pr)≥r2/64n_r^2\Delta(P_r)\ge r^2/64 となる(§8)。

  6. 6
    指数を読む

    点の個数は nr≈rαn_r\approx r^{\alpha}(α=45435k+16\alpha=45435k+16)の大きさなので、r2≈nr2/αr^2\approx n_r^{2/\alpha} となり、η=2/α\eta=2/\alpha が得られる。素数のあいだの nn にはベルトランの仮説の形の補題で対応する(§8、付録A)。

改良の源は、段階3の r−dr^{-d} です。点の個数を nr=⌊rN3/τ⌋n_r=\lfloor r\sqrt{N^3/\tau}\rfloor と選ぶと、悪い組の個数の期待値を nrn_r で割った EZ/nr\mathbb{E}Z/n_r の評価(式 (8.2))で、第2項が (log⁡2N)2(h/τ) r2−d(\log 2N)^2(h/\tau)\,r^{2-d} 以下になります[2, §8]。h/τh/\tau は r30r^{30} の定数倍以下で、d=41d=41 なので、この項は r→∞r\to\infty で 00 に近づきます。第1項(同じ点を表す2列の分)も 00 に近づくので、EZ\mathbb{E}Z は nrn_r よりずっと小さくなります。

論文はこれを、改良が「サンプルした点の算術的な分布」から来ると説明しています[2, §1.1]。最後の削除は命題5と同じ初等的な議論で、コムローシュ・ピンツ・セメレディが使ったハイパーグラフの独立集合の補題は使わない、とも述べています[2, §8]。

専門的な話体のノルムを行列式の和に分け、桁に書き込む

K=FrdK=\mathbb{F}_{r^d} の Fr\mathbb{F}_r 上の基底 β1,…,βd\beta_1,\dots,\beta_d をとり、3つの変数の組 X(1),X(2),X(3)X^{(1)},X^{(2)},X^{(3)}(それぞれ 3d3d 個の変数)から K3K^3 の列 Z(X(j))Z(X^{(j)}) を作ります。その行列式を D\mathcal{D} とし、係数にフロベニウス写像を施したものの積を N\mathcal{N} とします[2, §3.1]。

N=∏j=0d−1σj(D)\mathcal{N}=\prod_{j=0}^{d-1}\sigma^j(\mathcal{D})N=∏j=0d−1σj(D)\begin{aligned}\mathcal{N} &= \prod_{j=0}^{d-1}\sigma^j(\mathcal{D})\end{aligned}

N\mathcal{N} は係数が Fr\mathbb{F}_r に入る多項式で、Fr\mathbb{F}_r の値を代入すると行列式の体のノルムになります。dd が奇数なので、2つの変数の組を入れ替えると符号が変わります(交代式)。

交代式であることから、N\mathcal{N} は T=(M3)T=\binom M3 個の 3×33\times3 の「単項式の行列式」の和に書けます(Lemma 3.1。MM は 3d3d 変数の dd 次単項式の個数)。ラベルの列 (1,ξ,ξ2)(1,\xi,\xi^2) を代入すると、ヴァンデルモンドの行列式から次が成り立ちます。

N(X(ξ1),X(ξ2),X(ξ3))=Nm⁡K/Fr(∏i<j(ξj−ξi))≠0\begin{aligned} &\mathcal{N}\bigl(X(\xi_1),X(\xi_2),X(\xi_3)\bigr)\\ &\quad=\operatorname{Nm}_{K/\mathbb{F}_r}\Bigl(\prod_{i\lt j}(\xi_j-\xi_i)\Bigr)\ne0\end{aligned}

次に、和の aa 番目の項を、3つの行の決まった桁 ia=a, ja=Ta, ℓa=T2−(T+1)ai_a=a,\ j_a=Ta,\ \ell_a=T^2-(T+1)a に置きます。この3つの桁の和が k−1=T2k-1=T^2 になるのは、3つの添字がすべて同じ aa のときだけです(式 (3.7))。そのため、行列式を BB 進法で展開したときの Bk−1B^{k-1} の係数が、rr で割った余りで N\mathcal{N} と一致します。

桁 civc_{iv} は 11 以上 L=r10L=r^{10} 以下の整数から、必要な余りをもつものを一様に選びます。B>100k2L3B>100k^2L^3 なので、低い桁からの繰り上がりは Bk−1B^{k-1} に届きません。サレムとスペンサー、ベーレントの、3項等差数列を含まない集合の構成でも、大きな基数で繰り上がりを抑える同じ種類の工夫が使われている、と論文は述べています[2, §1.1]。

行列式がちょうど 00 の3列(3点が一直線)は、この補題では扱えません。そこで別の素数 qq で、楕円2次曲面のアフィン部分(一直線上に3点が並ばない集合、キャップ)から余りを選び、短い整数の一次関係を除きます(§5、Lemma 5.1、Lemma 5.2、Proposition 5.3。一次関係を除くのは主に Lemma 5.2)。行列式が 00 の3列は、それをみたす原始的な整数ベクトルごとに格子の点を数えて評価します(§7)。

どこまで確かなのか

  • 論文は2026年9月25日付のプレプリントで、査読はまだです[2]。
  • 2026年10月7日付の更新では、この論文は修正も取り下げもされていません[3]。
  • リポジトリには照合(Comparator)の設定 ComparatorChallenges/HeilbronnTriangle.json があり、2つの定理 heilbronn_power_lower_bound と almost_n_minus_two_refuted が検査の対象になっています[13]。許可された公理は、Leanで標準的な propext、Quot.sound、Classical.choice の3つです。
  • 形式化の範囲の説明によると、形式化されているのは「限りなく大きくなる点の個数の列 njn_j と、それぞれの点の配置」の存在です。論文の「十分大きいすべての nn」という主張のほうが広い、と明記されています[14]。
  • Leanの主張の指数は 1/(100000k)1/(100000k)(kk は上と同じ数)で、論文の η=2/(45435k+16)\eta=2/(45435k+16) より小さい、つまり弱い指数です。それでも 00 より大きいので、ほぼ n−2n^{-2} の上界の否定には十分です。
  • 形式化された論文の一覧 lean/formalization.yaml の主な結果の欄には、この論文は載っていません(10月9日に取得した版)。同じ欄には、照合の設定がある π\pi の無理数度の論文も載っていないので、この欄だけでは形式化の有無を判断できません。筆者はLeanの検査を手元で再実行していません。

Leanで照合される主張の前半は次のとおりです[13]。triangleAreasAtLeast (P j) a は「PjP_j の異なる3点が作る三角形の面積がすべて aa 以上」という意味です。

theorem heilbronn_power_lower_bound :
    0 < heilbronnExponent ∧ ∃ n : ℕ → ℕ, ∃ P : ℕ → Finset Point,
      Filter.Tendsto n Filter.atTop Filter.atTop ∧ ∀ j : ℕ,
        3 ≤ n j ∧ (P j).card = n j ∧ pointsInUnitSquare (P j) ∧
        triangleAreasAtLeast (P j)
          (Real.rpow (n j : ℝ) (-2 + heilbronnExponent))

後半の almost_n_minus_two_refuted は、ε=\varepsilon= heilbronnExponent / 2 で予想7の条件が成り立たないことを述べています。系9の証明と同じ形です。

なお論文の序論は、今回より強いべきの下界を主張する先行のプレプリントが2つあることにも触れています[2, §1]。論文によれば、エルマンの第12版[15]は、証明されていない局所一様性の仮定があることを認め、議論を発見的なものと述べています。

アガマの第13版[16]については、中心と境界の点を一緒に数えると一直線に並ぶ3点ができる、と論文は指摘しています。中心を除いても、その議論が評価しているのは中心と隣り合う境界の点が作る三角形なので、すべての3点についての評価は証明されないまま残る、とも書いています。

論文は、これらの指摘は引用した議論に対するもので、主張された評価が不可能だという意味ではないと断っています。

残された問題

  • 本当の指数:分かっているのは、指数が −2+η-2+\eta と −7/6-7/6 のあいだにあることだけです(図4)。Δ(n)\Delta(n) がどちらの端に近いのかは分かっていません。
  • 指数 η\eta の改良:論文はパラメータの最適化をしていないと書いています[2, §8]。3×10−2363\times10^{-236} よりずっと大きな η\eta が得られるかどうかは、これからの問題です。
  • 上からの評価:ザハロフは、指数 7/67/6 は今の手法の自然な限界で、それを超える改良はとても興味深いと書いています[5, §6]。

OpenAIの数学論文の解説には、ほかにπ\pi の無理数度は2と平面は5色では塗れないがあります。

参考文献

  1. OpenAI, openai/math: mathematical manuscripts and supporting proof artifacts, GitHub, 2026.
  2. OpenAI, A power improvement in the Heilbronn triangle lower bound, OpenAI Math Release preprint, September 25, 2026.
  3. OpenAI, openai/math: History, GitHub, 2026.
  4. K. F. Roth, On a problem of Heilbronn, Journal of the London Mathematical Society 26 (1951), 198–204.
  5. D. Zakharov, Small triangles, Journal of the London Mathematical Society 113 (2026), e70447.
  6. J. Komlós, J. Pintz, E. Szemerédi, A lower bound for Heilbronn's problem, Journal of the London Mathematical Society (2) 25 (1982), 13–24.
  7. W. M. Schmidt, On a problem of Heilbronn, Journal of the London Mathematical Society (2) 4 (1972), 545–550.
  8. K. F. Roth, On a problem of Heilbronn, II, Proceedings of the London Mathematical Society (3) 25 (1972), 193–212.
  9. K. F. Roth, On a problem of Heilbronn, III, Proceedings of the London Mathematical Society (3) 25 (1972), 543–549.
  10. J. Komlós, J. Pintz, E. Szemerédi, On Heilbronn's triangle problem, Journal of the London Mathematical Society (2) 24 (1981), 385–396.
  11. A. Cohen, C. Pohoata, D. Zakharov, A new upper bound for the Heilbronn triangle problem, arXiv:2305.18253, 2023.
  12. A. Cohen, C. Pohoata, D. Zakharov, Lower bounds for incidences, Inventiones mathematicae 240 (2025), 1045–1118.
  13. OpenAI, openai/math: lean/ComparatorChallenges/HeilbronnTriangle.lean, GitHub, 2026.
  14. OpenAI, openai/math: lean/docs/191.md(形式化の範囲), GitHub, 2026.
  15. G. Ellmann, A lower bound for Heilbronn's triangle-problem, arXiv:1703.03297v12, 2025.
  16. T. Agama, New bounds for the Heilbronn triangle problem, arXiv:2006.05269v13, 2026.