テトレーション

テトレーション（英: tetration）とは、冪乗の次の 4 番目のハイパー演算である。つまり、自らの冪乗を指定された回数反復する二項演算である。超冪（ちょうべき）ともいうが、この語は n ≥ 4 番目の一般のハイパー演算を総称するために用いられることもある。 第1から第4のハイパー演算は次のとおり。 加算 (hyper1) a + b = a + 1 + 1 + ⋯ + 1 ⏟ b {\displaystyle a+b=a+\underbrace {1+1+\cdots +1} _{b}} 乗算 (hyper2) a × b = a + a + ⋯ + a ⏟ b {\displaystyle a\times b=\underbrace {a+a+\cdots +a} _{b}} 冪乗 (hyper3) a b = a ↑ b = a × a × ⋯ × a ⏟ b {\displaystyle a^{b}=a\uparrow b=\underbrace {a\times a\times \cdots \times a} _{b}} テトレーション (hyper4) b a = a ↑↑ b = a ↑ a ↑ ⋯ ↑ a ⏟ b = a a ⋅ ⋅ a ⏟ b {\displaystyle {}^{b}a=a\uparrow \uparrow b=\underbrace {a\uparrow a\uparrow \cdots \uparrow a} _{b}=\underbrace {a^{a^{\cdot ^{\cdot ^{a}}}}} _{b}} 最初の 3 つ（加算、乗算、冪乗）は正則な複素関数へと一般化できるが、テトレーションは正則な複素関数への一般化が確立できず、また初等関数とも見なせない。 テトレーションは b を固定するごとに初等帰納的関数であるが、 b を変数と見なすと初等的ではない。 a を底（てい）、b を高さという。 任意の正の実数 a > 0 および非負整数 n ≥ 0 に対し、次のようにテトレーション na を再帰的に定める。 n a := { 1 if  n = 0 a [ ( n − 1 ) a ] if  n > 0 {\displaystyle {}^{n}a:={\begin{cases}1&{\text{if }}n=0\\a^{\left[{}^{(n-1)}a\right]}&{\text{if }}n>0\end{cases}}} この定義から明らかなように、冪の計算は右側（指数タワー表示における右上）から行われる。これはプログラミングにおける右結合に相当する。 4 2 = 2 2 2 2 = 2 [ 2 ( 2 2 ) ] = 2 ( 2 4 ) = 2 16 = 65536 {\displaystyle {}^{4}2=2^{2^{2^{2}}}=2^{\left[2^{\left(2^{2}\right)}\right]}=2^{\left(2^{4}\right)}=2^{16}=65536} 冪は結合法則を満たさないため、計算の順序を変えると値が変わってしまう。 4 2 = 2 2 2 2 ≠ [ ( 2 2 ) 2 ] 2 = 2 2 ⋅ 2 ⋅ 2 = 2 8 = 256 {\displaystyle {}^{4}2=2^{2^{2^{2}}}\neq \left[{\left(2^{2}\right)}^{2}\right]^{2}=2^{2\cdot 2\cdot 2}=2^{8}=256} a と 10 が互いに素であるとき、na の最後の d 桁がオイラーの定理から求められる。 テトレーションには多数の表記が存在する。テトレーションに使われる表記の中にはペンテーションやヘキセーションなど、より高次のハイパー演算の表記にも使用できるものもいくつかある。 上の表にある反復指数関数は、一般に次のように定義されるものである。 exp a n ⁡ ( x ) = a a ⋅ ⋅ a x {\displaystyle \exp _{a}^{n}(x)=a^{a^{\cdot ^{\cdot ^{a^{x}}}}}} （n 個の a の上に x が乗っている） 反復指数関数の表記のうち幾つかを以下に列挙する。 以下の表では、大部分の値が指数表記による表記すら困難なほど巨大であるため、それらの表記には底を 10 とした反復指数関数を用いた。なお小数部を持つ値はすべて近似値である。 任意の正の整数 n に対し、 nx の微分は次のようになる。 d d x n x = 1 x ∑ k = 1 n ( ln ⁡ x ) k − 1 ∏ j = n − k n j x {\displaystyle {\frac {\mathrm {d} }{{\mathrm {d} }x}}{}^{n}x={\frac {1}{x}}\sum _{k=1}^{n}(\ln x)^{k-1}\prod _{j=n-k}^{n}{}^{j}x} 1/2x、 2x の 0 から 1 までの定積分は二年生の夢と呼ばれる。 ∫ 0 1 1 2 x d x = ∑ n = 1 ∞ 1 2 n ( = 1.29128599706266354040728259059560054149861936827 … ) ∫ 0 1 2 x d x = − ∑ n = 1 ∞ 2 ( − n ) ( = 0.78343051071213440705926438652697546940768199014 … ) {\displaystyle {\begin{aligned}\int _{0}^{1}{\frac {1}{{}^{2}x}}\mathrm {d} x&=\quad \sum _{n=1}^{\infty }{\frac {1}{{}^{2}n}}&{\scriptstyle (=1.29128599706266354040728259059560054149861936827\dots )}\\\int _{0}^{1}{{}^{2}x}\,\mathrm {d} x&=-\sum _{n=1}^{\infty }{{}^{2}(-n)}&{\scriptstyle (=0.78343051071213440705926438652697546940768199014\dots )}\end{aligned}}} 任意の正の整数 n に対し、nx の不定積分は次のようになる。 ∫ n x d x = ∑ j = 0 n ( − 1 ) j ( j + 1 ) j − 1 j ! b j ( x ) + ∑ j = n + 1 ∞ ( − 1 ) j a n , j b j ( x ) + C {\displaystyle \int {}^{n}x{\mathrm {d} }x=\sum _{j=0}^{n}{\frac {(-1)^{j}(j+1)^{j-1}}{j!}}b_{j}(x)+\sum _{j=n+1}^{\infty }(-1)^{j}a_{n,j}b_{j}(x)+C} ここで aj,k は a j , k = { 1 if  k = 0 1 k ! if  j = 0 1 k ∑ l = 1 k l a j , k − l a j − 1 , l − 1 otherwise {\displaystyle a_{j,k}={\begin{cases}1&{\text{if }}k=0\\{\frac {1}{k!}}&{\text{if }}j=0\\{\frac {1}{k}}\sum _{l=1}^{k}la_{j,k-l}a_{j-1,l-1}&{\text{otherwise}}\end{cases}}} で与えられる有理数であり、 bj (x) は第2種不完全ガンマ関数を用いて b j ( x ) = Γ ( j + 1 , ln ⁡ x ) {\displaystyle b_{j}(x)=\Gamma \left(j+1,\ln x\right)} で与えられる。 テトレーションは、高さが正の整数以外の場合に拡張できる。 0の0乗が単純には定義できないため、 k0 は直接定義できないが、極限が lim n → 0 + k n = { 1 k ∈ 2 Z 0 k ∈ 2 Z + 1 {\displaystyle \lim _{n\to 0+}{}^{k}n={\begin{cases}1&k\in 2\mathbb {Z} \\0&k\in 2\mathbb {Z} +1\end{cases}}} と収束するため、 k 0 = { 1 k ∈ 2 Z 0 k ∈ 2 Z + 1 {\displaystyle {}^{k}0={\begin{cases}1&k\in 2\mathbb {Z} \\0&k\in 2\mathbb {Z} +1\end{cases}}} と定義する。（ここで、 2Z 、 2Z+1 はそれぞれ偶数、奇数の集合を表す。） なお、ここで 00 が一意に決まらないにもかかわらず 20 (= 00) が定義できるのは、 ab の a と b が等しいという条件下で極限を取ったからである。 複素数の累乗が可能なことから、テトレーションは複素数の底に対しても定義できる。 例えばテトレーション ni は対数関数の主枝を用いて定められる。このときオイラーの公式から次の式が得られる。 i a + b i = e 1 2 π i ( a + b i ) = e − 1 2 π b ( cos ⁡ π a 2 + i sin ⁡ π a 2 ) {\displaystyle i^{a+bi}=e^{{\frac {1}{2}}\pi i(a+bi)}=e^{-{\frac {1}{2}}\pi b}\left(\cos {\frac {\pi a}{2}}+i\sin {\frac {\pi a}{2}}\right)} 従って任意の ni = a + bi に対して n + 1i = a′ + b′i が次のように再帰的に定義できる。 a ′ = e − 1 2 π b cos ⁡ π a 2 b ′ = e − 1 2 π b sin ⁡ π a 2 {\displaystyle {\begin{aligned}a'&=e^{-{\frac {1}{2}}{\pi b}}\cos {\frac {\pi a}{2}}\\b'&=e^{-{\frac {1}{2}}{\pi b}}\sin {\frac {\pi a}{2}}\end{aligned}}} ここから以下の近似値が導かれる。 同様に値を逆向きに求めていくことで 0i = 1、 −1i = 0 が得られる。 ni の値を複素平面上にプロットすると、点列は渦巻状に極限値 0.4383 + 0.3606i へと近づく。この値は n→∞ のときと解釈できる。 このようなテトレーションの列に関する研究はオイラーの時代から続けられてきているものの、そのカオス的な振る舞いのために不明な所が多い。これまでに発表された研究のほとんどは、無限反復指数関数の収束について焦点を当てたものである。現在の研究は高性能のコンピュータを用いたフラクタルと数式処理システムの出現に大きな恩恵を受けている。テトレーションについて分かっていることの多くは、複素力学系の一般的な知識と、指数写像の専門的な研究によるものである。 ある範囲の底 a に対して limn→∞ na は有限の値に収束するので、この範囲においてテトレーションは高さ無限大の場合へ拡張できる。例えば limn→∞ n(√2) は収束して、その値は 2 であるから、 ∞(√2) = 2 であると言える。 一般に無限反復指数関数 x x ⋅ ⋅ ⋅ {\displaystyle x^{x^{\cdot ^{\cdot ^{\cdot }}}}} は、 n が無限大に向かうときの nx の極限として定義される。これが e−e ≤ x ≤ e1/e （およそ 0.066 から 1.44）の範囲で収束することはオイラーによって示された。 極限値 ∞x が存在するとき、これは方程式 ∞ x = x ∞ x {\displaystyle {}^{\infty }x=x^{{}^{\infty }x}} (1) を満たす正の実数に等しい。 式(1)より x = (∞x)1/∞x であり、このとき右辺の最大値が e1/e であることから、 x > e1/e については極限値が存在しないことがわかる。また x ↦ ∞ x {\displaystyle x\mapsto {}^{\infty }x} は y ↦ y 1 / y {\displaystyle y\mapsto y^{1/y}} の逆関数（の分岐の下部）であることがわかる。 式(1)から、極限値 ∞z をランベルトのW関数 W を用いて次のように定義することで複素数の底 z に対しても拡張される。 ∞ z = z z ⋅ ⋅ ⋅ = W ( − ln ⁡ z ) − ln ⁡ z {\displaystyle {}^{\infty }z=z^{z^{\cdot ^{\cdot ^{\cdot }}}}={\frac {\mathrm {W} (-\ln {z})}{-\ln {z}}}} (2) 定義より k n = log n ⁡ { ( k + 1 ) n } {\displaystyle {}^{k}n=\log _{n}\left\{{}^{(k+1)}n\right\}} が成り立つので、この関係を k ≤ 0 に対しても帰納的に拡張し 1 n = log n ⁡ ( 2 n ) = log n ⁡ ( n n ) = n log n ⁡ n = n 0 n = log n ⁡ ( 1 n ) = log n ⁡ n = 1 − 1 n = log n ⁡ ( 0 n ) = log n ⁡ 1 = 0 {\displaystyle {\begin{array}{rclclcccc}{}^{1}n&=&\log _{n}\left({}^{2}n\right)&=&\log _{n}\left(n^{n}\right)&=&n\log _{n}n&=&n\\{}^{0}n&=&\log _{n}\left({}^{1}n\right)&=&\log _{n}n&&&=&1\\{}^{-1}n&=&\log _{n}\left({}^{0}n\right)&=&\log _{n}1&&&=&0\end{array}}} と定義する。 ただし、定義できるのは k = −1 までで、log 0 が存在しないため k = −2 に対しては定義できず、したがって k ≤ −2 に対して拡張できない。 テトレーションを高さ実数または複素数へ拡張する、という一般的な問題への広く受け入れられた解答は今のところ存在しない。いくつかのアプローチについて以下で述べる。 一般にこの問題は任意の実数 a > 0 に対し、実数 x > −2 で定義され、次の条件を満たす超指数関数 f (x) = xa を求めるものである。 f (0) = 1、 f (−1) = 0 任意の実数 x > −1 に対し f (x) = a↑↑(f (x−1)) 三つ目の条件は通常次の中のどれかである。 連続性（通常は x > 0 における a、 x についての連続性） 微分可能性（x についての 1、 2、 k 回または無限回微分可能性） 正則性（x についての 2 階微分可能性を含む） 任意の実数 x > 0 に対して d 2 d x 2 f ( x ) > 0 {\displaystyle \displaystyle {{\frac {\mathrm {d} ^{2}}{\mathrm {d} x^{2}}}f(x)>0}} が成立する。 三つ目の条件は著者およびアプローチによって異なる。実数高さへの拡張には二つの主要なアプローチが存在し、一つは正則性、もう一つは微分可能性に基づいたものである。これらの二つのアプローチは、相反する結果を導くことから互いに大きく異なるとされ、調和は難しいと考えられている。 長さ 1 の区間で xa が定義されれば、任意の x >−2 に対し容易に拡張される。 一次近似（連続性のもと微分可能性を近似）は次のように与えられる。 x a ≈ { log a ⁡ ( x + 1 a ) x ≤ − 1 1 + x − 1 < x ≤ 0 a ( x − 1 a ) 0 < x {\displaystyle {}^{x}a\approx {\begin{cases}\log _{a}({}^{x+1}a)&x\leq -1\\1+x&-1<x\leq 0\\a^{\left({}^{x-1}a\right)}&0<x\end{cases}}} ゆえに と以下続く。但しこの微分可能性はあくまで区分的なものである。整数 x を境に微分係数が ln a 倍されるため、（a = e の場合を除いて）整数 x において微分不可能となる。 以下は値の計算例である。 π / 2 e ≈ 5.868... − 4.3 0.5 ≈ 4.03335... − 5.264 0.6 ≈ − 5.35997... 3.1 0.7 ≈ 0.7580... {\displaystyle {\begin{aligned}{}^{\pi /2}e&\approx 5.868...\\{}^{-4.3}{0.5}&\approx 4.03335...\\{}^{-5.264}{0.6}&\approx -5.35997...\\{}^{3.1}{0.7}&\approx 0.7580...\end{aligned}}} フーシュマンドは ultra exponential function という関数を導入した。 これはテトレーションの一次近似を表し、 uxpa (x) と表記される。 uxpa は、次の定理によって一意に定められる。 定理 ―  0 < a ≠ 1 {\displaystyle 0<a\neq 1} に対し、連続な関数 f : ( − 2 , + ∞ ) → R {\displaystyle f\colon (-2,+\infty )\to {\mathbb {R} }} が条件 任意の実数 x > − 1 {\displaystyle x>-1} に対し f ( x ) = a f ( x − 1 ) {\displaystyle f(x)=a^{f(x-1)}} f ( 0 ) = 1 {\displaystyle f(0)=1} f {\displaystyle f} は ( − 1 , 0 ) {\displaystyle (-1,0)} で微分可能 f ′ {\displaystyle f'} は ( − 1 , 0 ) {\displaystyle (-1,0)} で広義単調増加か広義単調減少 lim x → 0 + f ′ ( x ) = ( ln ⁡ a ) lim x → 0 − f ′ ( x ) {\displaystyle \lim _{x\to 0+}f'(x)=(\ln a)\lim _{x\to 0-}f'(x)} を満たすとき、 f {\displaystyle f} は次の方程式によって一意に定まる。 f ( x ) = exp a [ x ] ⁡ ( a { x } ) = exp a [ x + 1 ] ⁡ ( { x } ) {\displaystyle f(x)=\exp _{a}^{[x]}\left(a^{\lbrace {x}\rbrace }\right)=\exp _{a}^{[x+1]}\left(\lbrace {x}\rbrace \right)} 但し [ x ] {\displaystyle [x]} はガウス記号であり、また { x } = x − [ x ] {\displaystyle \lbrace {x}\rbrace =x-[x]} である。 証明は、三番目から五番目の条件より f が [−1, 0] で線型となることから従う。 フーシュマンドはさらに次のような一意性定理を導いた。 定理 ―  連続な関数 f : ( − 2 , + ∞ ) → R {\displaystyle f\colon (-2,+\infty )\to {\mathbb {R} }} が条件 任意の x > − 1 {\displaystyle x>-1} に対し f ( x ) = e f ( x − 1 ) {\displaystyle f(x)=e^{f(x-1)}} f ( 0 ) = 1 {\displaystyle f(0)=1} f {\displaystyle f} は ( − 1 , 0 ) {\displaystyle (-1,0)} で下に凸 lim x → 0 − f ′ ( x ) ≤ lim x → 0 + f ′ ( x ) {\displaystyle \lim _{x\to 0-}f'(x)\leq \lim _{x\to 0+}f'(x)} を満たすとき、 f {\displaystyle f} は u x p ( = u x p e ) {\displaystyle \mathrm {uxp} (=\mathrm {uxp} _{e})} に等しい。 証明は先とほぼ同様である。漸化式より lim x → − 1 + 0 f ′ ( x ) = lim x → + 0 f ′ ( x ) {\displaystyle \lim _{x\to -1+0}f'(x)=\lim _{x\to +0}f'(x)} となること、三番目の条件より f が (−1, 0) で線型となることから従う。 定理より、 x >−1 に対し f (x) = exp (f (x−1)) および f (0) = 1 であって、かつ (−1, +∞) で下に凸であるような関数 f は唯一 uxp のみである。f が十分微分可能であるためには (−1, 0) で極値を持つ必要がある。 （微分可能性についての）二次近似は次のように与えられる。 x a ≈ { log a ⁡ ( x + 1 a ) x ≤ − 1 1 + 2 ln ⁡ ( a ) 1 + ln ⁡ ( a ) x − 1 − ln ⁡ ( a ) 1 + ln ⁡ ( a ) x 2 − 1 < x ≤ 0 a ( x − 1 a ) 0 < x {\displaystyle {}^{x}a\approx {\begin{cases}\log _{a}({}^{x+1}a)&x\leq -1\\1+{\frac {2\ln(a)}{1+\ln(a)}}x-{\frac {1-\ln(a)}{1+\ln(a)}}x^{2}&-1<x\leq 0\\a^{\left({}^{x-1}a\right)}&0<x\end{cases}}} これは任意の x > 0 について微分可能であるが二階微分可能でない。a = e のときこれは一次近似に等しくなる。 三次近似および高次への一般化は次のように与えられる。 次の条件を満たす関数 F が一意に定まる事が証明されている 。 F (z + 1) = exp (F (x)) F (0) = 1 z→±i∞ のとき F (z) が対数関数の不動点（およそ 0.318 ± 1.337i ）に近づく 実数 z < −2 を除く複素平面全域で正則 この関数 F を右図に示す。また、底が e ではない場合についても、底が e 1 / e {\displaystyle e^{1/e}} よりも大きい場合については同様に証明されている。倍精度浮動小数点数近似はオンラインで公開されている。 テトレーションを一意に定めるためには正則性の条件が重要となる。いま、関数 F に対し関数 S を次のように構成する。 A ( z ) = ∑ n = 1 ∞ sin ⁡ ( 2 π n z ) α n B ( z ) = ∑ n = 1 ∞ ( 1 − cos ⁡ ( 2 π n z ) ) β n S ( z ) = F ( z + A ( z ) + B ( z ) ) {\displaystyle {\begin{aligned}A(z)&=\sum _{n=1}^{\infty }\sin(2\pi nz)\alpha _{n}\\B(z)&=\sum _{n=1}^{\infty }{\bigl (}1-\cos(2\pi nz){\bigr )}\beta _{n}\\S(z)&=F\left(z+A(z)+B(z)\right)\end{aligned}}} ここで αn、 βn は十分速く減衰する実数列であり、少なくとも実軸の近くで A (z)、 B (z) を収束させるとする。 この関数 S は F と同様に最初の二つの条件 S (z + 1) = exp (S (z))、 S (0) = 1 を満たす。また αn、 βn が十分速く 0 に近づくとき、S は正の実軸近傍で解析的となる。しかし αn、 βn が全て 0 でない場合、S は新たに大量の特異点と不連続線を複素平面上に持つことになる。これは sin (z)、 cos (z) が虚軸に沿って指数関数的に増大するためである。これらの特異点は αn、 βn が小さければ小さいほど実軸から離れていくため、 S が正則であるためには全ての αn、 βn が 0となる、即ち S = F であればよい。 実解析上のテトレーションは一意的に定まらないので、複素平面への拡張は一意性に必要である。 nπ, ne が整数になるような正の整数 n は存在するか。特に 4π は整数か。 与えられた正の整数 n と正の（整数でない）有理数 q に対し、nq は整数か。 特に 4x = 2 の正の解 x は有理数か。 冪は冪根と対数の二つの逆関係を持つ。これに倣って、以下テトレーションの逆関係をそれぞれ超冪根（英: super-root）と超対数（英: super-logarithm）と呼ぶ。 超冪根はテトレーションの底に関する逆関係である。 超平方根（英: super square root）は 2x の逆であり、二つの等価な表記 ssrt(x), √xs を持つ。 この関数は次のようなランベルトのW関数による表示を持つ。 s s r t ( x ) = e W ( l n ( x ) ) = l n ( x ) W ( l n ( x ) ) {\displaystyle \mathrm {ssrt} (x)=e^{W(\mathrm {ln} (x))}={\frac {\mathrm {ln} (x)}{W(\mathrm {ln} (x))}}} またこの関数により、冪根と対数の間の鏡映的な関係が表れる。次の方程式は y = ssrt x のときに真となる。 x y = log y ⁡ x {\displaystyle {\sqrt[{y}]{x}}=\log _{y}x} 平方根と同様に超平方根は一つとは限らない。ただし、平方根と異なり超平方根の個数を決定するのは容易とは言えない。一般に、e−1/e < x < 1 のとき x は二つの正の超平方根を 0 と 1 の間に持ち、1 < x のとき x は 1 より大きい一つの正の超平方根を持ち、0 < x < e−1/e のとき x は超平方根を実数の範囲で持たない。しかし上の式より、任意の x (≠ 1) は可算無限個の超平方根を複素数の範囲で持つことが従う。 超平方根はネットワークのクラスタサイズを決定するのに使用される。 任意の整数 n > 2 に対して nx は定義され、x ≥ 1 のとき増加となり n1 = 1 を満たす。従って x ≥ 1 のとき n√xs は存在する。しかし上述した一次近似を用いた場合、−1 < y ≤ 0 のとき yx は x によらず y + 1 となり、従ってこの場合 x = y√y + 1s は存在しない。 超平方根のほか n 次の超冪根も同様の記号を用いて n√xs と表すことができる。 超冪根は高さが無限大の場合へと拡張することができ、これは 1/e ≤ x ≤ e の場合に限り問題なく定義される（#高さが無限大を参照）。∞x が存在するとき ∞x = x∞x が成り立つことから、無限次の超冪根は初等関数によって ∞√xs = x1/x と表すことができる。例えば ∞√2s = 21/2 = √2 となる。 n を任意の正の整数とするとゲルフォント＝シュナイダーの定理より、超平方根 √ns は整数または超越数となり、超立方根 3√ns は整数または無理数（これが超越数かどうかは知られていない）となる。 超対数はテトレーションの高さに関する逆関係である。 テトレーション xa を x に関して連続的に増加するものとして定義すると、任意の実数 x に対し超対数 sloga x (a > 1) が定義される。 この関数 sloga x は以下の式を満たす。 slog a ⁡ x a = x slog a ⁡ a x = 1 + slog a ⁡ x slog a ⁡ x = 1 + slog a ⁡ log a ⁡ x slog a ⁡ x > − 2 {\displaystyle {\begin{array}{lcl}\operatorname {slog} _{a}{^{x}a}&=&x\\\operatorname {slog} _{a}a^{x}&=&1+\operatorname {slog} _{a}x\\\operatorname {slog} _{a}x&=&1+\operatorname {slog} _{a}\log _{a}x\\\operatorname {slog} _{a}x&>&-2\end{array}}} アッカーマン関数 二重指数関数（英語版） ハイパー演算 ペンテーション 反復対数（英語版） en:Symmetric level-index arithmetic ^ Exploding Array Function ^ 配列表記（巨大数研究 Wiki） ^ Ioannis Galidakis. On Extending hyper4 and Knuth’s Up-arrow Notation to the Reals ^ “Power Verb”. J Vocabulary. J Software. 2011年10月28日閲覧。 ^ Edward Drake Roe, Jr. (1898年). “68”. The American Mathematical Monthly 5: 110. doi:10.2307/2971013. http://www.jstor.org/stable/2971013.  ^ I. N. Galidakis, (2004年). “On an Application of Lambert’s W Function to Infinite Exponentials”. Complex Variables Th. Appl. 49: 759–780. doi:10.1080/02781070412331298796. ISSN 0278-1077. http://www.math.usm.edu/lee/InfiniteExponentials.pdf.  ^ Euler, L., "De serie Lambertina Plurimisque eius insignibus proprietatibus." Acta Acad. Scient. Petropol. 2, 29–51, 1783. Reprinted in Euler, L. Opera Omnia, Series Prima, Vol. 6: Commentationes Algebraicae. Leipzig, Germany: Teubner, pp. 350–369, 1921. (facsimile) ^ M. H. Hooshmand, (2006年). “Ultra power and ultra exponential functions”. Integral Transforms and Special Functions（英語版） 17 (8): 549–558. doi:10.1080/10652460500422247.  ^ Andrew Robbins. Solving for the Analytic Piecewise Extension of Tetration and the Super-logarithm ^ W. Paulsen and S. Cowgill (2017年3月). “Solving F ( z + 1 ) = b F ( z ) {\displaystyle F(z+1)=b^{F(z)}} in the complex plane”. Advances in Computational Mathematics: 1-22. doi:10.1007/s10444-017-9524-1. http://link.springer.com/article/10.1007/s10444-017-9524-1.  ^ テトレーションおよびその導関数を計算・描画するMathematicaコード ^ Marshall, Ash J., and Tan, Yiren, "A rational number of the form aa with a irrational", Mathematical Gazette 96, March 2012, pp. 106-109. ^ a b Corless, R. M.; Gonnet, G. H.; Hare, D. E. G.; Jeffrey, D. J.; Knuth, D. E. (1996年). “On the Lambert W function” (PostScript). Advances in Computational Mathematics 5: 333. doi:10.1007/BF02124750. http://www.apmaths.uwo.ca/~rcorless/frames/PAPERS/LambertW/LambertW.ps.  ^ BOSTON UNIVERSITY COLLEGE OF ENGINEERING – EFFICIENT SELF-ORGANIZATION OF LARGE WIRELESS SENSOR NETWORKS Daniel Geisler's site on tetration Tetration Forum Tetration - TORI - Mizugadro, the research site by Dmitrii Kouznetsov Gottfried Helms' site on tetration