MathLabs

算術と数論

円周法

加法方程式の解の個数を数えるために単位円上のフーリエ積分を用いる解析的手法。

直観直感:円を一周しながら数える

整数 nn を、整数の集合 AA(たとえば素数や kk 乗数)から取った kk 個の要素の和として書く方法の個数を数えたいとする。AA を単位円上をちょうど一周する「波」として符号化する:[0,1)[0,1) の各実数 α\alpha に対して指数和 FA(α)=∑a∈Ae(aα)F_A(\alpha) = \sum_{a \in A} e(a\alpha) を作る。ここで e(x):=e2πixe(x) := e^{2\pi i x}。α\alpha が 00 から 11 まで動くと FA(α)kF_A(\alpha)^k は振動する。分母 qq が小さい有理数 a/qa/q という少数の「共鳴」点の近くでは和の各項が揃って足し合わさり(メジャーアーク)、それ以外のほとんどの場所では位相がほぼランダムな方向を向いて打ち消し合う(マイナーアーク)。円周法はこの幾何学的描像を厳密な公式に変え、弧ごとに評価する。

円周法の指数和を構成するために使われる、角度シータの点を示すインタラクティブな単位円。
θ\theta をドラッグすると点 e2πiθ/360e^{2\pi i\theta/360} が単位円上を動く。円周法は θ/360=α\theta/360=\alpha が [0,1)[0,1) を走るときの FA(α)kF_A(\alpha)^k に対してこのような点を積分する。

大学定義:円上の母関数

定義: 指数和と表現数

有限集合 A⊂Z≥0A\subset\mathbb Z_{\ge0} に対して FA(α)=∑a∈Ae(aα)F_A(\alpha) = \sum_{a \in A} e(a\alpha)(ただし e(x):=e2πixe(x) := e^{2\pi i x})とおく。n≥0n\ge0、k≥1k\ge1 について、rk(n)=#{(a1,…,ak)∈Ak:a1+⋯+ak=n}r_k(n)=\#\{(a_1,\dots,a_k)\in A^k : a_1+\cdots+a_k=n\} を、nn を AA の kk 個の要素の和として表す順序付き表現の個数とする。

FA(α)=∑a∈Ae(aα),e(x):=e2πixF_A(\alpha) = \sum_{a \in A} e(a\alpha), \qquad e(x) := e^{2\pi i x}

FAF_A を kk 乗して展開すると、e(nα)e(n\alpha) の係数は「本来」ちょうど rk(n)r_k(n) になるはずである。以下の抽出恒等式(定理1)がこれを厳密にする:rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha。ここで FA(α)kF_A(\alpha)^k は上の図の波そのものであり、積分は我々が関心を持つ一つの振動数 nn をちょうど取り出す——これがこの方法のすべてである:数え上げ問題を積分評価問題に置き換えるのだ。

rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha
メジャーアークとマイナーアーク
特徴メジャーアーク(小さい qq の a/qa/q 付近)マイナーアーク(それ以外)
位置q≤Qq\le Q である各有理数 a/qa/q の周りの短い区間すべてのメジャーアークを除いた後の [0,1)[0,1) の残り
FA(α)F_A(\alpha) の大きさ自明な最大値 ∣A∣|A| に近い:項が強め合う∣A∣|A| よりはるかに小さいと期待される:項が打ち消し合う
評価における役割主要項(「特異級数」)を生む上から評価して誤差項に吸収させる必要がある

大学二つの基礎定理

任意の有限集合 A⊂Z≥0A\subset\mathbb Z_{\ge0}、整数 k≥1k\ge1、n≥0n\ge0 に対して rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha が成り立つ。ここで rk(n)r_k(n) は和が nn になる AA からの順序付き kk 組の個数である。

なぜ正しいのか?

これにより、純粋に組合せ論的な数え上げ問題を、ある解析的積分の大きさに関する問いに置き換えることができる:その積分が正であることを示せれば、表現が存在することになる。

証明

まず直交関係を用いる:∫01e(mα) dα\int_0^1 e(m\alpha)\, d\alpha=1=1(m=0m=0 のとき)、=0=0(m≠0m\neq 0 のとき)。e(mα)=cos⁡(2πmα)+isin⁡(2πmα)e(m\alpha)=\cos(2\pi m\alpha)+i\sin(2\pi m\alpha) と書けば、m≠0m\neq0 のときこれは [0,1)[0,1) 上でサイン・コサイン波のちょうど整数周期分であり積分は 00 になる。m=0m=0 のとき被積分関数は定数 11 である。

次に母関数の kk 乗を展開する:FA(α)k=(∑a∈Ae(aα))k=∑(a1,…,ak)∈Ake((a1+⋯+ak)α)F_A(\alpha)^k=\Big(\sum_{a\in A}e(a\alpha)\Big)^k=\sum_{(a_1,\dots,a_k)\in A^k} e\big((a_1+\cdots+a_k)\alpha\big)、これは指数 m=a1+⋯+akm=a_1+\cdots+a_k ごとにまとめた、AA からの順序付き kk 組全体にわたる有限和である。

両辺に e(−nα)e(-n\alpha) を掛け、[0,1)[0,1) 上で項別に積分する(有限和なので積分と和の交換は正当である):∫01FA(α)ke(−nα) dα=∑(a1,…,ak)∈Ak∫01e((a1+⋯+ak−n)α) dα\int_0^1 F_A(\alpha)^k e(-n\alpha)\,d\alpha=\sum_{(a_1,\dots,a_k)\in A^k}\int_0^1 e\big((a_1+\cdots+a_k-n)\alpha\big)\,d\alpha。

直交関係により、右辺の各項は a1+⋯+ak=na_1+\cdots+a_k=n のときちょうど 11、それ以外は 00 である。したがって和全体は a1+⋯+ak=na_1+\cdots+a_k=n を満たす組の個数、すなわち rk(n)r_k(n) にちょうど一致する。これにより恒等式が近似なしに厳密に証明される:これは漸近式ではなく恒等式である。

任意の実数 α\alpha と任意の整数 Q≥1Q\ge1 に対して、1≤q≤Q1\le q\le Q、gcd⁡(a,q)=1\gcd(a,q)=1 を満たす有理数 a/qa/q が存在して ∣α−aq∣≤1q(Q+1)\left|\alpha - \dfrac{a}{q}\right| \le \dfrac{1}{q(Q+1)} が成り立つ。

なぜ正しいのか?

これにより、円上のすべての点が分母の小さいある有理数の近くにあることが保証される。これがまさに [0,1)[0,1) を、FAF_A が大きくなる小さい qq の良い近似 a/qa/q の周りの短い区間であるメジャーアークと、残りのマイナーアークに分割することを可能にする。

証明

0,{α},{2α},…,{Qα},10,\{\alpha\},\{2\alpha\},\dots,\{Q\alpha\},1 という Q+2Q+2 個の数を考える。ここで {x}\{x\} は xx の小数部分であり、すべて [0,1][0,1] に属する。[0,1][0,1] を長さ 1/(Q+1)1/(Q+1) の Q+1Q+1 個の等しい小区間に分割する:[0,1Q+1),[1Q+1,2Q+1),…[0,\tfrac1{Q+1}), [\tfrac1{Q+1},\tfrac2{Q+1}),\dots。

Q+2Q+2 個の数に対して小区間は Q+1Q+1 個しかないので、鳩の巣原理により、あるふたつの数 {jα}\{j\alpha\} と {iα}\{i\alpha\}(0≤i<j≤Q0\le i<j\le Q、i=0i=0 すなわち {iα}=0\{i\alpha\}=0 も許す)が同じ小区間に入り、1/(Q+1)1/(Q+1) 未満の差になる:∣{jα}−{iα}∣<1Q+1|\{j\alpha\}-\{i\alpha\}|<\tfrac{1}{Q+1}。

q=j−iq=j-i とおくと 1≤q≤Q1\le q\le Q。{jα}−{iα}=(jα−iα)−(⌊jα⌋−⌊iα⌋)=qα−a\{j\alpha\}-\{i\alpha\} = (j\alpha - i\alpha) - (\lfloor j\alpha\rfloor - \lfloor i\alpha\rfloor) = q\alpha - a(a=⌊jα⌋−⌊iα⌋a=\lfloor j\alpha\rfloor-\lfloor i\alpha\rfloor は整数)なので、∣qα−a∣<1Q+1|q\alpha - a| < \tfrac{1}{Q+1}、すなわち ∣α−aq∣<1q(Q+1)\left|\alpha-\dfrac aq\right| < \dfrac{1}{q(Q+1)} が得られる。

最後に gcd⁡(a,q)=d>1\gcd(a,q)=d>1 ならば aa と qq を dd で割る:得られる分数は分母がさらに小さく、α\alpha との距離も同じか小さいので、一般性を失うことなく gcd⁡(a,q)=1\gcd(a,q)=1 とできる。以上で証明が完了する。これは有限で構成的な鳩の巣論法であり、未証明の評価には一切頼っていない。

大学実世界での応用と具体例

数論を超えて、「波を足し合わせて共鳴を探す」という同じ発想は、信号処理や電気工学で広く使われる離散フーリエ変換の動作原理でもある。また、ハーディとラマヌジャンが円周法の初期版(分割数関数 p(n)p(n) に適用)で導いた漸近式 p(n)∼14n3exp⁡ ⁣(π2n3)p(n) \sim \dfrac{1}{4n\sqrt3}\exp\!\left(\pi\sqrt{\dfrac{2n}{3}}\right) は、統計力学において、固定した全エネルギー nn を持つ区別できないボース励起系の微視状態数――したがってエントロピー――を見積もるのに使われる。

例: 円周法の小規模計算

A={1,2,…,8}A=\{1,2,\dots,8\} とする。抽出恒等式を用いて、a+b=9a+b=9 を満たす順序対 (a,b)∈A2(a,b)\in A^2 の個数 r2(9)r_2(9) を計算せよ。

解答

定理1により、FA(α)=∑a=18e(aα)F_A(\alpha)=\sum_{a=1}^{8}e(a\alpha) として r2(9)=∫01FA(α)2e(−9α) dαr_2(9)=\int_0^1 F_A(\alpha)^2 e(-9\alpha)\,d\alpha である。実際に積分を解析的に計算する必要はない――恒等式によりそれは直接数えた個数に等しいことが保証されているので、組合せ論的に数を数えてその恒等式を信頼すればよい。

a,b∈{1,…,8}a,b\in\{1,\dots,8\}、a+b=9a+b=9 を満たす順序対 (a,b)(a,b) をすべて挙げる:(1,8),(2,7),(3,6),(4,5),(5,4),(6,3),(7,2),(8,1)(1,8),(2,7),(3,6),(4,5),(5,4),(6,3),(7,2),(8,1)。11 から 88 までの各項は b=9−ab=9-a を一意に定め、bb は常に {1,…,8}\{1,\dots,8\} に戻る(1≤a≤8⇒1≤9−a≤81\le a\le8 \Rightarrow 1\le 9-a\le8 より)ので、aa の 88 通りすべてが有効である。

したがって r2(9)r_2(9) は 88 に等しい。この小さな例は定理1の恒等式がまさに働いている場面である:解析的な積分と組合せ的な個数は、構成上、まったく同一の数である――この方法の本当の内容は、AA が(素数のように)無限または増大する集合になり、列挙ではなく評価が必要になって初めて現れる。

例: ディリクレの定理による最良有理近似

α=2\alpha=\sqrt2、Q=5Q=5 としてディリクレの近似定理を適用し、定理2の評価を満たす分数 a/qa/q(1≤q≤51\le q\le5)を一つ示し、数値的に確かめよ。

解答

1≤q≤51\le q\le5 かつ ∣2−a/q∣<1/(q⋅6)|\sqrt2-a/q|<1/(q\cdot6)(定理の評価で Q+1=6Q+1=6 とする)を満たす a/qa/q を探す。q=5q=5 を試す:52≈7.07115\sqrt2\approx7.0711 に最も近い整数は a=7a=7 で、7/5=1.47/5=1.4 を与える。

評価を確認する:∣2−7/5∣=∣1.41421…−1.4∣≈0.01421|\sqrt2-7/5|=|1.41421\ldots-1.4|\approx0.01421 であり、定理が保証するのは 1/(5⋅6)=1/30≈0.03331/(5\cdot6)=1/30\approx0.0333。実際 0.01421<0.03330.01421<0.0333 なので、評価は余裕を持って成り立ち、保証どおりである。

この 7/57/5 は実は 2=[1;2,2,2,… ]\sqrt2=[1;2,2,2,\dots] の連分数の二段階目の近似分数であり、これほど良く近似する理由でもある――ディリクレの鳩の巣論法自体は連分数を必要としないが、常に同程度に良い近似を与える。したがって 75\dfrac{7}{5} は有効な例である。

抽出恒等式 rk(n)=∫01FA(α)k e(−nα) dαr_k(n) = \int_0^1 F_A(\alpha)^k\, e(-n\alpha)\, d\alpha において、なぜこの積分はちょうど rk(n)r_k(n) に等しいのか?

A={1,2,…,8}A=\{1,2,\dots,8\} のとき、a+b=9a+b=9 を満たす順序対 (a,b)∈A2(a,b)\in A^2 の個数 r2(9)r_2(9) はいくつか?

ディリクレの近似定理によれば、任意の実数 α\alpha と整数 Q≥1Q\ge1 に対して、何の存在が保証されるか?

統計力学において、分割数関数のハーディ・ラマヌジャン漸近式 p(n)∼14n3exp⁡ ⁣(π2n3)p(n) \sim \dfrac{1}{4n\sqrt3}\exp\!\left(\pi\sqrt{\dfrac{2n}{3}}\right) は何を見積もるのに役立つか?

参考文献

  1. G. H. Hardy, S. Ramanujan (1918). Asymptotic formulae in combinatory analysis · DOI:10.1112/plms/s2-17.1.75
  2. J. Bourgain, C. Demeter, L. Guth (2016). Proof of the main conjecture in Vinogradov's Mean Value Theorem for degrees higher than three · DOI:10.4007/annals.2016.184.2.7 · arXiv:1512.01565
  3. B. Green, T. Tao (2008). The primes contain arbitrarily long arithmetic progressions · DOI:10.4007/annals.2008.167.481 · arXiv:math/0404188