在數學裡面，**迭代冪次** （亦作**超-4運算**或**四級運算**），或可理解為迭代乘方、冪塔運算和超冪運算等等，是專指冪的下一個超運算級別，用以表示極大的數字。以下列舉了首四個超運算級別，其中**迭代冪次**為第四級，（后继函数，例如a' = a + 1即將a加上一，可理解為第零級運算，相關解釋參見皮亞諾公理）。範例如下：

1. 加法
a + n = a\!\underbrace{''{}^{\cdots}{}'}\_n
*a* 的連續*n* 次後繼
2. 乘法
a \times n = \underbrace{a + a + \cdots + a}\_n
*a* 連續加上自己*n* 次
3. 冪
a^n = \underbrace{a \times a \times \cdots \times a}\_n
*a* 連續乘以自己*n* 次
4. **迭代冪次**
{^{n}a} = \underbrace{a^{a^{\cdot^{\cdot^{a}}}}}\_n
*a* 連續取冪於自己*n* 次

以上每一個運算級別皆被定義為對上一運算級別的迭代（**迭代冪次**的下一個運算級別為五級運算（超-5運算））。**迭代冪次**跟首三個超運算級別的一大不同之處在於首三個超運算級別中的*n* 可以是任意複數，而*n* 為任意複數的**迭代冪次**目前則未有一個概括的定義。另外，**迭代冪次**不屬於初等函數。

加法(*a* + *n*)是最基本的運算級別；乘法(*a* × *n*)亦是其中一種初等函數，在自然數的域當中，它可被視為*a* 的*n* 次鏈式加法；冪(*a^n*)則可被視為*a* 的*n* 次鏈式乘法。如此類推，**迭代冪次**(*{^{n}a}*)可被視為*a* 的*n* 次鏈式冪。當中，變量*a* 將會在下文被稱為底數，而變量*n* 則是此函數的*高度值*，在下文有時會被稱為**上標數**（早段提及的**上標數**皆為整數，而後則會擴展到分數、實數以及複數，如下所示）。

定義
--

對於任何正實數*a* ( a > 0 )及非負整數*n* ( n \ge 0 )，\,\! {^{n}a} 被定義為：

{^{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^{\left[2^{\left(2^2\right)}\right]} = 2^{\left(2^4\right)} = 2^{16} = 65,\!536
要注意，冪是不遵從結合律的，因此以其他順序來計算上述表達式將會出現不一樣的答案，例如：

\,\! 2^{2^{2^2}} \ne \left[{\left(2^2\right)}^2\right]^2 = 2^{2 \cdot 2 \cdot2} = 256
因此，**冪塔**一定要從上而下（或從右至左）來運算。在電腦程式中，此制式稱為**右結合律**。

當*a*與*n*為互質時，我們可以透過歐拉定理來計算\scriptstyle\underbrace{a^{a^{\cdot^{\cdot^{a}}}}}\_n的最後*m*個小數位值。

專門用語
----

**迭代冪次**在英文裡面稱作*tetration*，有時亦會被稱為*superexponentiation*及*hyperpower*（中文意譯**超冪**）等，這些詞語也可被用來表示這種運算模式。

**迭代冪次**有時會跟一些相關的函數及表達式混淆，這是因為在這些函數及表達式當中的大部分專門用語均適用於**迭代冪次**。以下列舉了一些相關用語：

| 形式 | 用語 |
| --- | --- |
| a^{a^{\cdot^{\cdot^{a^a}}}} | 迭代冪次 |
| a^{a^{\cdot^{\cdot^{a^x}}}} | 迭代指數 |
| a\_1^{a\_2^{\cdot^{\cdot^{a\_n}}}} | 指數群 （亦作指數塔） |
| a\_1^{a\_2^{a\_3^{\cdot^{\cdot^\cdot}}}} | 無窮指數 （亦作無窮指數塔） |

在首兩種表達式當中的*a* 是*底數*，而*a* 出現的數目則是*高度值*（*x* 的出現使*高度值* 加1）。在第三種表達式當中，*n* 是*高度值*，但每一個底皆不相同。

要注意的是，**迭代指數**的形式有時也會被稱為**迭代冪次**。這是模稜兩可的，因為這可以指迭代乘方或迭代指數。

符號標示法
-----

可以用來表示**迭代冪次**的符號有很多，當中有一些符號可用來表示更高級的迭代運算（hyper-5、hyper-6 等等）。

| 名稱 | 形式 | 描述 |
| --- | --- | --- |
| 標準符號記法 | {\displaystyle \,{}^{n}a,a[4]n} | Maurer及Goodstein分別於1901年及1947年使用此記法，其後由美國數學家鲁迪·拉克於其著作"Infinity and the Mind"中將這個記法普及化。 |
| 高德納箭號表示法 | a {\uparrow\uparrow} n | 允許加上更多箭號，乃至在箭號的右上方標上正整數*n*（↑*n*），以作推廣。 |
| 康威鏈式箭號表示法 | a \rightarrow n \rightarrow 2 | 允許把*2* 改成更大的整數以作推廣（同上述推廣法），亦可透過延長鏈式來作出推廣。 |
| 阿克曼函數 | {}^{n}2 = \operatorname{A}(4, n - 3) + 3 | 底數為2（a=2）的情況下可寫成阿克曼函數式。 |
| 迭代指數表示法 | {}^{n}a = \exp\_a^n(1) | 迭代指數可以為1以外的數字。 |
| Hooshmand符號記法 | \operatorname{uxp}\_a n, \, a^{\frac{n}{}}  |  |
| 超運算符號 | a^{(4)}n, \, \operatorname{hyper}\_4(a,n) | 允許把4改為更大的整數，以表示更高級數的超運算 |
| ASCII符號 | `a^^n` | 由於上箭號的用法跟脫字符（`^`）一樣，**迭代冪次**運算符可寫成（`^^`）。同样允许更多^连续表示； |

上述的迭代指數表示法中使用的迭代指數記號，一般被定義成：

\exp\_a^n(x) = a^{a^{\cdot^{\cdot^{a^x}}}}，當中包含*n* 個*a*。
以下是一些用以表示迭代指數的符號：

| 名稱 | 形式 | 描述 |
| --- | --- | --- |
| 標準符號記法 | \exp\_a^n(x) | 歐拉創造了符號\exp\_a(x)來表示a^x，設f(x)=\exp\_a(x)，則\exp\_a^n(x)可表示成迭代函數f^n(x)。 |
| 高德納箭號表示法 | (a{\uparrow})^n(x) | 允許增加箭號的數目用以表示**超冪**（**迭代冪次**）和超指數（迭代指數），在大數一條目中常被用到。 |
| Ioannis Galidakis符號記法 | \,{}^{n}(a, x) | 底數部分可以改成較長的表達式。 |
| ASCII（輔助） | `a^^n@x` | 以迭代指數作為迭代冪次的輔助函數。 |
| ASCII（標準） | `exp_a^n(x)` | 基於標準符號記法。 |
| J符號記法 | `x^^:(n-1)x` | 重複冪次。詳見J語言 |

例子
--

在下表，大部分數值大得連科學記數法也難以表示，因此使用了迭代指數記號，設底數為10來表示。包含小數點的數值是近似值。

| x | {}^{2}x | {}^{3}x | {}^{4}x |
| --- | --- | --- | --- |
| 1 | 1 | 1 | 1 |
| 2 | 4 | 16 | 65,536 |
| 3 | 27 | 7,625,597,484,987 | \exp\_{10}^3(1.09902) |
| 4 | 256 | \exp\_{10}^2(2.18788) | \exp\_{10}^3(2.18726) |
| 5 | 3,125 | \exp\_{10}^2(3.33931) | \exp\_{10}^3(3.33928) |
| 6 | 46,656 | \exp\_{10}^2(4.55997) | \exp\_{10}^3(4.55997) |
| 7 | 823,543 | \exp\_{10}^2(5.84259) | \exp\_{10}^3(5.84259) |
| 8 | 16,777,216 | \exp\_{10}^2(7.18045) | \exp\_{10}^3(7.18045) |
| 9 | 387,420,489 | \exp\_{10}^2(8.56784) | \exp\_{10}^3(8.56784) |
| 10 | 10,000,000,000 | \exp\_{10}^3(1) | \exp\_{10}^4(1) |

以較原始的函數來作逼近法
------------

### 多項式*逼近法*

#### 線性逼近法

以下是對迭代冪次函數的*線性逼近法*（以滿足連續函數的需要性，逼近法基於此函數的可微性質）的定義：

{}^{x}a \approx \begin{cases}
\log\_a(^{x+1}a) & x \le -1 \\
1 + x & -1 < x \le 0 \\
a^{\left(^{x-1}a\right)} & 0 < x
\end{cases}
由此可得：

| 逼近法 | 定義域 |
| --- | --- |
| \,{}^{x}a \approx x+1 | for -1<x<0 |
| \,{}^{x}a \approx a^x | for 0<x<1 |
| \,{}^{x}a \approx a^{a^{(x-1)}} | for 1<x<2 |

及其他逼近值。不過，這個函數只是分段可微的；在x為整數的時候，函數的導數要乘以\ln{a}。

##### 例子

\begin{align}
 {}^{\frac{1}{2}\pi}e &\approx 5.868...\\
 {}^{-4.3}0.5 &\approx 4.03335...
\end{align}

Hooshmand的手稿中有一個重要定理：設 0 <a \neq 1。若f:(-2,+\infty)\rightarrow \mathbb{R} 是連續的並滿足以下條件：

* f(x)=a^{f(x-1)} \; \; \mbox{for all} \; \; x>-1, \; f(0)=1,
* f 於  (-1, 0)  之上可微；
* f^{\prime } 在  (-1,0), 之上是一個單調函數；
* f^\prime (0^+) = (\ln a) f^\prime (0^-) \mbox{ or } f^\prime (-1^+) = f^\prime (0^-).

由此，f 可於以下方程式中獨特地定義出來：

f(x)=\exp^{[x]}\_a (a^{x})=\exp^{[x+1]}\_a(x) \quad \mbox{for all} \; \; x > -2，
當中， (x)=x-[x]  標示x的分數部分，以及 \exp^{[x]}\_a  是函數 \exp\_a  的 [x] -迭代函數。

以上四個條件中的第二個條件僅當f 在[-1, 0]之上是線性函數，由此可作為上述對f 的定義的證明。

對自然**迭代冪次**函數{}^xe 的線性逼近法是連續可微的，但其二階導數的輻角並不是整數值。Hooshmand為此導出了以下的這一個獨特定理：

若 f: (-2, +\infty)\rightarrow \mathbb{R}
是一個連續函數並滿足以下條件:

* f(x)=e^{f(x-1)} \; \; \mbox{for all} \; \; x>-1, \; f(0)=1；
* f 於(-1,0) 之上是凸函數；
* f^\prime (0^-)\leq f^\prime (0^+).

那麼f=\mbox{uxp}。（這裡的f=\mbox{uxp} 是Hooshmand給予自然迭代冪次函數的線性逼近法的名稱。）

這個定理的證明跟之前提到的證明法十分相似；遞迴方程式保證了f^\prime (-1^+) = f^\prime (0^+), ，而f 的凸函數的性質僅當f 在(-1, 0)之上是線性的。

所以，自然迭代函數的線性逼近法是凸於(-1,+\infty) 之上的方程式 f(x)=e^{f(x-1)} \; \; (x>-1) 的唯一解。所有其他充足可微的解在區間(-1, 0)之上必定存在一個拐點。

#### 更高次的逼近法

以下是對a ≠ e的迭代冪次函數的*二次逼近法*（逼近法基於此函數的可微性質）的定義：

{}^{x}a \approx \begin{cases}
\log\_a({}^{x+1}a) & x \le -1 \\
\frac{-\log(a) \;+\; \sqrt{x [1 \;-\; \log(a)^2] \;+\; 1}}{1 \;-\; \log(a)} & -1 < x \le 0 \\
a^{\left({}^{x-1}a\right)} & 0 < x
\end{cases}
這對所有x>0 可微，但並不二次可微。若a = e，應採用線性逼近法。

有關於三次逼近法，以及一個能歸納出*n* 次逼近法的方法，詳見：。

推廣
--

**迭代冪次**能被推廣至定義 {^n 0}  乃至其他定義域。

### 對底數定義域的推廣

#### 推廣至底數為0

指數0^{0} 是不連續定義的。所以，迭代冪次\,{^n0} 於早期提出的公式中亦並不被清晰定義。不過，\lim\_{x\rightarrow0}{^n}x 是定義良好的，並存在：

\lim\_{x\rightarrow0} {}^{n}x = \begin{cases} 1, & n \mbox{ even} \\ 0, & n \mbox{ odd} \end{cases} 
所以我們能連續地定義{}^{n}0 = \lim\_{x\rightarrow0} {}^{n}x。這等價於定義 0^0 = 1 。

在這推廣之下，{}^{0}0 = 1，所以最初定義出來的法則 {^{0}a} = 1  依然成立。

#### 推廣至複數底

由於複數可以作為指數，迭代冪次的底數可以為\scriptstyle z \;=\; a + bi 的形式，當中的\scriptstyle i 是−1的平方根。舉例來說，設\scriptstyle z \;=\; i，對於\scriptstyle {}^{n}z，其迭代冪次可由自然對數中的主枝（英语：）來達成，並用歐拉公式得出以下關係：

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)
這表明了\scriptstyle {}^{(n+1)}i \;=\; a'+b'i 在任何\scriptstyle {}^{n}i \;=\; a+bi 的情況下的遞迴定義為：

  

\begin{align}
 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{align}
從而導出以下的逼近值：

| {}^{n}i | 逼近值 |
| --- | --- |
| {}^{1}i = i | *i* |
| {}^{2}i = i^{\left({}^{1}i\right)} | 0.2079 |
| {}^{3}i = i^{\left({}^{2}i\right)} | 0.9472 + 0.3208i |
| {}^{4}i = i^{\left({}^{3}i\right)} | 0.0501 + 0.6021i |
| {}^{5}i = i^{\left({}^{4}i\right)} | 0.3872 + 0.0305i |
| {}^{6}i = i^{\left({}^{5}i\right)} | 0.7823 + 0.5446i |
| {}^{7}i = i^{\left({}^{6}i\right)} | 0.1426 + 0.4005i |
| {}^{8}i = i^{\left({}^{7}i\right)} | 0.5198 + 0.1184i |
| {}^{9}i = i^{\left({}^{8}i\right)} | 0.5686 + 0.6051i |

根據上一部分對於迭代冪次的逆向關係的定義，可得\scriptstyle \,{}^{0}i \;=\; 1 及\scriptstyle \,{}^{(-1)}i \;=\; 0，當中負值的*n* 在虛數軸上會得出無窮的結果。在複平面當中，整個序列成螺旋形地趨向於極限0.4383 + 0.3606i，這個極限可理解為*n* 為無窮時，函數相對應的值。

這樣的迭代冪次序列由歐拉時期已開始被研究，但是由於序列的雜亂性而難以被理解。歷史上大部分有正式發表的研究皆集中於冪塔函數的收歛性。高運算效率電腦連同計算機代數系統和分形幾何系統的出現大大地促進了近代對於*迭代冪次*的研究。現時對迭代冪次的研究均建基於複動力學的普遍知識及對指數映射的專門研究。

### 對高度值定義域的推廣

#### 推廣至無窮高

迭代冪次可被推廣至無窮高（{}^{n}a 當中的*n*）。這是因為當底數在一個特定的區間之內而高度值趨向於無窮時，迭代冪次會收歛於一個有限的數值。舉例來說，\sqrt{2}^{\sqrt{2}^{\sqrt{2}^{\cdot^{\cdot^{\cdot}}}}} 收歛於2，因此可以說是等於2。對2的趨向性可從對以下小型有限冪塔的計算而看出來：

\begin{align}
\sqrt{2}^{\sqrt{2}^{\sqrt{2}^{\sqrt{2}^{\sqrt{2}^{1.414}}}}} &= \sqrt{2}^{\sqrt{2}^{\sqrt{2}^{\sqrt{2}^{1.63}}}} \\
&= \sqrt{2}^{\sqrt{2}^{\sqrt{2}^{1.76}}} \\
&= \sqrt{2}^{\sqrt{2}^{1.84}} \\
&= \sqrt{2}^{1.89} \\
&= 1.93
\end{align}
一般來說，有限冪塔x^{x^{\cdot^{\cdot^{\cdot}}}} （定義為當*n* 趨向於無窮時{}^{n}x 的極限）收歛於*e*−*e* ≤ *x* ≤ *e*1/*e*，大約是位於0.066和1.44之間的區間，這是由萊昂哈德·歐拉所證明的。如果存在一個極限，這會是一個對於方程式*y* = *x**y* 的正實數解。所以，*x* = *y*1/*y*。根據這個極限的定義，當*x* > *e*1/*e* 時，*x* 的無窮迭代冪次不具收歛性，因為*y*1/*y* 的最大值為*e*1/*e*。

以上特性可以被推廣至複數底*z*，定義如下：

{}^{\infty}z = z^{z^{\cdot^{\cdot^{\cdot}}}} = \frac{\mathrm{W}(-\ln{z})}{-\ln{z}} 
當中的*W*（*z*）表示朗伯W函數。

由於極限*y* = ∞*x* （如果存在的話，即當*e*−*e* < *x* < *e*1/*e*）必定滿足*x**y* = *y*，因而得出*x* ↦ *y* = ∞*x* 是*y* ↦ *x* = *y*1/*y* 的反函數（較低枝）。

#### （有限地）推廣至負高

為了維持原有法則：

 {^{(k+1)}a} = a^{({^{k}a})} 
當k 為負值時，必須用到以下的遞迴關係：

 {^{k}a} = \log\_a \left( {^{(k+1)}a} \right) 
所以：

 {}^{(-1)}a = \log\_{a} \left( {}^{0}a \right) = \log\_{a} 1 = 0 
不過，當高度值為更小的負值時就不能以此方法良好地定義出來了，因為

 {}^{(-2)}a = \log\_{a} \left( {}^{-1}a \right) = \log\_{a} 0 
這是定義不良好的。

更要注意的是，當 a = 1  時，任何根據法則對於\,\! {^{(-1)}1}  的定義都是一致的，因為

 {^{0}1} = 1 = 1^n  對於任何 \,\! n = {^{(-1)}1} 。
#### 推廣至實高

當前還未有對於推廣迭代冪次至n 為實數或複數值的共識解。以下提到了兩種不同的逼近法。

一般來說，問題在於對任何實數*a* > 0，找出一個能滿足以下條件的*超冪函數* \,f(x) = {}^{x}a ，當中x > -2 並為實數：

* \,{}^{(-1)}a = 0
* \,{}^{0}a = 1
* \,{}^{x}a = a^{\left({}^{(x-1)}a\right)} 對所有實數*x* > -1。
* 第四個條件通常為以下其中一個：

* A的*連續性* 需要（通常是指在*a* 和*x* 皆可變時，{}^{x}a 為連續，當中x>0）。
* A的*可微性* 需要（可以是對於*x* 一次、二次、*k* 次，或是無窮可微）。
* A的*規律性* 需要（僅當對於*x* 二次可微），即是：

\left( \frac{d^2}{dx^2}f(x) > 0\right) 對於所有 x > 0 
對於第四個條件，不同的編者有不同的說法，而且亦視乎於採用何種逼近法。對於把迭代冪次推廣至實高，有兩種主要的逼近法，一種是建基於*規律性* 需要，另一種則建基於*可微性* 需要。這兩種逼近法似乎十分相異，皆因它們所得出的結果並不相符。

幸運地，任何在一段長度的區間內滿足到其中一種逼近法的解，皆能被推廣為一個對於所有正實數高度值的迭代冪次的通解。當\,{}^{x}a 在一段長度的區間內被定義，對任何x > -2，整個函數的後續將能被輕易地定義出來。

其中一个简单的推广方式为：

{\displaystyle a[n]b=c}

当且仅当{\displaystyle c[n]({\frac {1}{b}})=a,n>1} 。可以判断，当{\displaystyle n=2}、{\displaystyle n=3}时，乘法和幂次成立：{\displaystyle a\times b=c} 当且仅当{\displaystyle c\times ({\frac {1}{b}})=a};{\displaystyle a^{b}=c} 当且仅当{\displaystyle c^{\frac {1}{b}}=a}

举例，当{\displaystyle n=4} 时, 计算 {\displaystyle {}^{\frac {1}{2}}256} 和 {\displaystyle {}^{\frac {3}{4}}256} ：
{\displaystyle {}^{2}4=4^{4}=256}, 则 {\displaystyle {}^{\frac {1}{2}}256=4}。
{\displaystyle {}^{3}2={2}^{{2}^{2}}=16},则{\displaystyle {}^{\frac {1}{3}}16=2}；
由 256= 16^2 = 16^(16^^(1/3)) = 16^^(1+1/3) = 16^^(4/3) 得 256^^(3/4) = 16。

#### 推廣至複高

以下為有關猜想：函數*F* 為方程式*F*(*z*+1)=exp(*F*(*z*)) 的解並滿足以下附加條件：當*z*逼近於±*i* ∞及*F* 在整個複數*z* 平面當中為全純函數，*F*(0)=1 及*F*(*z*)逼近於對數的不動點 （大約為 0.31813150520476413531 ± 1.33723570143068940890*i*）。

相關條目
----

* 阿克曼函數
* 超運算

參考資料
----

1. 1 2  M. H. Hooshmand,.. Integral Transforms and Special Functions. 2006, **17** (8): 549–558. doi:10.1080/10652460500422247.
2. ↑  Ioannis Galidakis. On Extending hyper4 and Knuth’s Up-arrow Notation to the Reals （页面存档备份，存于）
3. ↑ . J Vocabulary. J Software.  [28 October 2011]. （原始内容存档于2021-05-06）.
4. ↑  Andrew Robbins. Solving for the Analytic Piecewise Extension of Tetration and the Super-logarithm （页面存档备份，存于）
5. ↑  崔雷.. 科协论坛. 2011, (1): 96.
6. ↑  D. Kouznetsov.  (PDF). Mathematics of Computation. July 2009, **78** (267): 1647–1670. doi:10.1090/S0025-5718-09-02188-7.

* Daniel Geisler, *tetration.org （页面存档备份，存于）*
* Ioannis Galidakis, *On extending hyper4 to nonintegers* (undated, 2006 or earlier) *(A simpler, easier to read review of the next reference)*
* Ioannis Galidakis, *On Extending hyper4 and Knuth's Up-arrow Notation to the Reals* (undated, 2006 or earlier).
* Robert Munafo, *Extension of the hyper4 function to reals （页面存档备份，存于）* *(An informal discussion about extending tetration to the real numbers.)*
* Lode Vandevenne, *Tetration of the Square Root of Two （页面存档备份，存于）*, (2004). *(Attempt to extend tetration to real numbers.)*
* Ioannis Galidakis, *Mathematics*, *(Definitive list of references to tetration research. Lots of information on the Lambert W function, Riemann surfaces, and analytic continuation.)*
* Galidakis, Ioannis and Weisstein, Eric W. Power Tower （页面存档备份，存于）
* Joseph MacDonell, *Some Critical Points of the Hyperpower Function （页面存档备份，存于）*.
* Dave L. Renfro, *Web pages for infinitely iterated exponentials （页面存档备份，存于）* *(Compilation of entries from questions about tetration on sci.math.)*
* R. Knobel. "Exponentials Reiterated." *American Mathematical Monthly* **88**, (1981), p. 235-252.
* Hans Maurer. "Über die Funktion y=x^{[x^{[x(\cdots)]}]} für ganzzahliges Argument (Abundanzen)." *Mittheilungen der Mathematische Gesellschaft in Hamburg* **4**, (1901), p. 33-50. *(Reference to usage of \ {^{n} a} from Knobel's paper.)*
* Marco Ripà, "La strana coda della serie n^n^...^n", Trento (2011). ISBN 9788861787896

外部連結
----

* Andrew Robbins' site on tetration
* Daniel Geisler's site on tetration （页面存档备份，存于）
* Tetration Forum （页面存档备份，存于）
* tetration at citizendium （页面存档备份，存于）
* Gottfried Helms' site on tetration （页面存档备份，存于）