dc (程序)

原作者 | Robert Morris (於AT&T貝爾實驗室) Lorinda Cherry |
---|---|
開發者 | 各種開源和商業開發者 |
程式語言 | B |
作業系統 | Unix, 類Unix, Plan 9 |
平台 | 跨平台 |
類型 | 命令 |
dc(desk calculator:桌面計數機)是採用逆波蘭表示法的跨平台計數機,它支援任意精度算術[1]。它是Robert Morris在貝爾實驗室工作期間書寫的[2],作為最老的Unix實用工具,先於C語言的發明。像那個年代的其他實用工具一樣,它有着一組強力的特徵和簡潔的語法[3][4]。傳統上,採用中綴表示法的bc計數機程式是在dc之上實現的。
歷史
[編輯]dc是倖存的最老的Unix語言[2]。在貝爾實驗室收到第一台PDP-11的時候,用B語言寫成的dc是在這個新機器上執行的第一個語言,甚至在組譯器之前[5]。
基本運算
[編輯]在dc中要做4和5的乘法:
$ dc
4 5 *
p
20
q
這可轉譯為「把4和5壓入棧頂,通過乘法算符,從棧中彈出兩個元素,將二者相乘並把結果壓回棧頂」。接着使用p
命令列印棧頂的元素。使用q
命令退出此次呼叫的dc實例。注意數值相互間必須以空白分隔,但某些算符可以不必如此。
還可以用如下命令得到這個結果:
$ dc -e '4 5 * p'
20
$ echo "4 5 * p" | dc
20
$ dc -
4 5*pq
20
$ cat <<EOF > cal.txt
4 5 *
p
EOF
$ dc cal.txt
20
使用命令k
來變更算術精度,它設置算術運算的小數碼數。因為預設精度是0
,例如:
$ dc -e '10946 6765 / p'
1
通過使用命令k
調整精度,可以產生任意數目的小數數碼,例如:
$ dc -e '7k 10946 6765 / p'
1.6180339
dc有科學計數機的基本運算功能,比如求黃金分割率的值:
$ dc -e '7k 5 v 1 + 2 / p'
1.6180339
dc將輸入數值的前導_
辨識為負號,命令^
計算冪,命令v
計算平方根。
d
命令用於複製棧頂元素。r
命令用於對棧頂和僅次棧頂的兩個元素進行對換。z
命令用於壓入當前棧深度,即執行z
命令前棧中元素的數目。
輸入/輸出
[編輯]使用?
命令,從stdin讀取一行並執行它。這允許從宏中向用戶要求輸入,故而此輸入必須是語法上正確的,並且這有潛在的安全問題,因為dc的!
命令可以執行任意系統命令。
前面提及過,p
命令列印棧頂元素,帶有隨後的一個換行。n
命令彈出棧頂元素並輸出它,沒有尾隨換行。f
命令列印整個棧,從棧頂到棧底並且一項一行。
dc還支援控制輸入和輸出的底數。o
命令設置輸出底數,輸出底數必須大於等於2。i
命令彈出棧頂元素並將它用作輸入底數,輸入底數必須在2和16之間,十六進制數字必須大寫以避免和dc命令衝突。要記住輸入底數將影響對後面的所有數值的分析,所以通常建議在設置輸入底數之前先設置輸出底數。例如將二進制轉換成十六進制:
$ echo '16o2i 1001101010111100110111101111 p' | dc
9ABCDEF
要讀取設置的這些數值,K
、I
和O
命令將壓入當前精度、輸入基數和輸出基數到棧頂。
語言特徵
[編輯]除了上述的基本算術和棧操作,dc包括了對宏、條件和儲存結果用於以後檢索的支援。
暫存器
[編輯]暫存器在dc中是有着單一字元名字的存貯位置,它可以通過命令來儲存和檢索,它是宏和條件的底層機制:sc
彈出棧頂元素並將它儲存入暫存器c
,而lc
將暫存器c
的值壓入棧頂。例如:
3 sc 4 lc * p
暫存器還被當作次要棧,可以使用S
和L
命令在它們和主要棧之間壓入和彈出數值。儲存棧頂元素到暫存器中並把這個元素留在棧頂,需要聯合使用ds
命令。
字串
[編輯]字串是包圍在[
和]
之中的字元,可以被壓入棧頂和存入暫存器。使用x
命令從棧頂彈出字串並執行它,使用P
命令從棧頂彈出並列印字串,無尾隨換行。a
命令可以把數值的低位位元組轉換成ASCII字元,或者在棧頂是字串時把它替換為這個字串的第一個字元。此外沒有方法去建造字串或進行字串操縱。
#
字元開始一個註釋直到此行結束。
宏
[編輯]通過允許暫存器和棧專案像數值一樣儲存字串,從而實現了宏。一個字串可以被列印,也可以被執行,就是說作為dc命令的序列而傳遞。例如可以把一個宏「加1
並接着乘以2
」儲存到一個暫存器m
中:
[1 + 2 *]sm
通過使用x
命令彈出棧頂的字串並執行之,如下這樣使用儲存的宏:
3 lmx p
Q
命令從棧頂彈出一個值作為退出宏的層數,比如2Q
命令退出2層宏,它永不導致退出dc。q
命令退出2層宏,如果宏少於2層則退出dc。
條件
[編輯]最後提供了有條件執行宏的機制。命令=r
將從棧頂彈出兩個值,如果二者相等,則執行儲存在暫存器r
中的宏。如下命令序列將在原棧頂元素等於5的條件下列印字串equal
。
[[equal]p]sm d5=m
這裏使用了d
命令保留原棧頂元素。其他條件有>
、!>
、<
、!<
、!=
,如果棧頂元素分別大於、不大於(小於等於)、小於、不小於(大於等於)、不等於僅次於棧頂的元素,則執行指定的宏。注意不同於Forth、PostScript和Factor,在不等式比較中的運算元的次序同在算術中的次序相反,5 3 -
等價於中綴表示法的5 - 3
,然而5 3 <t
在3 < 5
時執行暫存器t
的內容。下面是其執行範例:
$ echo 5 | dc -e '? [[equal]p]sm d5=m'
equal
$ echo 5 | dc -e '? [[[equal]pq]sm d5=m [not equal]p]x'
equal
$ echo 4 | dc -e '? [[[equal]pq]sm d5=m [not equal]p]x'
not equal
迭代範例
[編輯]階乘
[編輯]# F(x):
# if x > 1
# x * F(x-1)
# else
# x
它在x > 0
時有效,它可採用互遞歸寫為:
# F(x):
# x > 1 ? G(x) : x
# G(x):
# x * F(x-1)
[d 1- lFx *]sG [d1<G]dsFx
它可以進一步採用尾遞歸實現為:
[d 1- lFx]sG [d1<G]dsFx 1 [* z1<H]dsHx
這裏的運算由兩個步驟串接而成,首先將遞減的整數壓入堆疊形成整數集上偏序關係的區間,然後在這個數列上進行應用乘法運算的右歸約得出一個最終的結果值。這裏的命令z
,將當前的堆疊深度即堆疊中元素數目壓入棧頂。還可以進一步增加針對x ≤ 0
情況的預處理。下面是其執行範例:
$ echo 0 | dc -e '? [d-1+]sad0!<a [d 1- lFx]sG [d1<G]dsFx 1 [* z1<H]dsHx p'
1
$ echo 9 | dc -e '? [d-1+]sad0!<a [d 1- lFx]sG [d1<G]dsFx 1 [* z1<H]dsHx p'
362880
# n := x
# i := 0
# F(x):
# i := i + 1
# x := x * i
# print(x)
# i < n ? F(x) : x
# F(1)
這裏迭代的是需要初始化的堆疊中的自動變數x
,它是,其初始值1
是;i
既是迭代的遞增計數器,也是迭代二元運算的運算元,其遞增形成整數集上偏序關係的區間,同步地在這個數列上進行輸出中間值的應用乘法運算的左歸約得出結果數列,這裏的x
所起到的作用也被稱為累加器,其含義為「累計」或「累積」而不必然採用加法或乘法。它可實現為:
sn 0si 1 [li1+dsi *p liln>F]dsFx
下面是其執行範例:
$ echo 6 | dc -e '? sn 0si 1 [li1+dsi *p liln>F]dsFx'
1
2
6
24
120
720
可以將它改為計算單個的階乘:
$ echo 0 | dc -e '? sn 0si 1 [li1+dsi * liln>F]dsFx p'
1
$ echo 9 | dc -e '? sn 0si 1 [li1+dsi * liln>F]dsFx p'
362880
斐波那契數
[編輯]# F(x):
# if x > 1
# F(x-1) + F(x-2)
# else
# x
它在x ≥ 0
時有效,它可採用互遞歸寫為:
# F(x):
# x > 1 ? G(x) : x
# G(x):
# F(x-1) + F(x-2)
它可實現為:
[d 2- lFx r 1- lFx +]sG [d1<G]dsFx
這裏的r
命令反轉(交換)棧頂元素和僅次於棧頂的元素二者的次序。它可以進一步採用記憶化而寫為:
# F(x):
# x > 1 ? M(x) : x
# M(x):
# if x in a
# a[x]
# else
# temp := G(x)
# a[x] := temp
# temp
# G(x):
# F(x-1) + F(x-2)
基於陣列儲存命令:
和訪問命令;
,它可實現為:
0Sa 1si [d 2- lFx r 1- lFx +]sG [;aq]sN [dli!<N lGx dli1+dsi:a]sM [d1<M]dsFx Lasa
或者更一般化的寫為:
# a[0] := 0
# a[1] := 1
# F(x):
# if x in a
# a[x]
# else
# temp := G(x)
# a[x] := temp
# temp
# G(x):
# F(x-1) + F(x-2)
它可實現為:
0Sa 0 0:a 1 1:a 1si [d 2- lFx r 1- lFx +]sG [;aq]sN [dli!<N lGx dli1+dsi:a]dsFx Lasa
下面是其執行範例,展示了通過給它加列印斷點來觀察其遞歸呼叫的規模和詳細步驟:
$ echo 0 | dc -e '? [d 2- lFx r 1- lFx +]sG [d1<G]dsFx p'
0
$ echo 20 | dc -e '? [d 2- lFx r 1- lFx +]sG [d1<G]dsFx p'
6765
$ echo 21 | dc -e '? [d 2- lFx r 1- lFx +]sG [d1<G]dsFx p'
10946
$ echo 20 | dc -e '? [p d 2- lFx r 1- lFx +]sG [d1<G]dsFx' | wc -l
10945
$ echo 20 | dc -e '? 0Sa 0 0:a 1 1:a 1si [d 2- lFx r 1- lFx +]sG [;aq]sN [p dli!<N lGx dli1+dsi:a]dsFx Lasa' | wc -l
39
$ echo 6 | dc -e '? 0Sa 1si [d 2- lFx r 1- lFx +]sG [;aq]sN [zn[ F]np dli!<N lGx dli1+dsi:a]sM [d1<M]dsFx Lasa'
1 F6
2 F4
3 F2
3 F3
4 F2
2 F5
3 F3
3 F4
$ echo 6 | dc -e '? 0Sa 1si [d 1- lFx r 2- lFx +]sG [;aq]sN [zn[ F]np dli!<N lGx dli1+dsi:a]sM [d1<M]dsFx Lasa'
1 F6
2 F5
3 F4
4 F3
5 F2
4 F2
3 F3
2 F4
在遞歸定義的遞歸呼叫中,隨着的索引n
的遞減,而遞歸下降至被計算出來並儲存的第一項,然後沿着呼叫鏈逐級返回;在遞歸的每個步驟中都要計算並儲存,它所需要的和這兩項之中,如果首先進行遞歸呼叫並從其返回,則也需要計算並儲存,為此需要訪問已儲存的和;如果首先進行遞歸呼叫並從其返回,則需要訪問已儲存的;這種遞歸返回形成了遞推關係,每個的值被用到一次或兩次,首次是從遞歸呼叫直接返回,再次是對其已儲存的值進行訪問,每個步驟要麼訪問兩項並儲存兩項,要麼訪問一項並儲存一項。故而它可以進一步寫為:
1si [d 2- lFx r 1- lFx +]sG [2%;aq]sN [dli!<N lGx dli1+dsi2%:a]sM [d1<M]dsFx
在遞歸定義的上述兩種求值次序中,首先進行遞歸呼叫的這種形式,其記憶化實現從一開始就持續進行遞歸呼叫直到首次遞歸返回,然後就逐級遞歸返回而不再進行遞歸呼叫,故而它可以進一步寫為:
# a := 0
# F(x):
# x > 1 ? G(x) : x
# G(x):
# curt := F(x-1)
# prev := a
# a := curt
# curt + prev
這裏的a
被初始化為的值0
,x
在步入遞歸呼叫之後,除了在被遞減至的值1
之時作為呼叫結果來返回之外,沒有參與到後續的加法運算之中,而主要起到了遞減計數器的作用。它可實現為:
0sa [1- lFx la r dsa +]sG [d1<G]dsFx
下面的例子採用迭代法列印出斐波那契數列的不含第0
項的前n
項:
# i := x
# a := 1
# F(x):
# i > 0 ? G(x) : x
# G(x):
# prev := a
# a := x
# x := x + prev
# print(x)
# i := i - 1
# F(x)
# F(0)
針對遞推關係,它總共進行n
次迭代運算得出到。這裏的i
是進行迭代的遞減計數器,所迭代的是需要初始化的兩個變數x
和a
:在堆疊中的自動變數x
,先是被迭代的,再是迭代出來的,其初始值0
是,它所起到的作用也被稱為累加器;作為全域變數的a
,以所得的1
作為初始值,在首次迭代之後它是始於的。它可實現為:
si 1sa 0 [lardsa +p li1-si lFx]sG [li0<G]dsFx
下面是其執行範例:
$ echo 6 | dc -e '? si 1sa 0 [lardsa +p li1-si lFx]sG [li0<G]dsFx'
1
1
2
3
5
8
可以將它改為計算單個的斐波那契數:
$ echo 0 | dc -e '? si 1sa 0 [lardsa + li1-si lFx]sG [li0<G]dsFx p'
0
$ echo 9 | dc -e '? si 1sa 0 [lardsa + li1-si lFx]sG [li0<G]dsFx p'
34
參見
[編輯]參照
[編輯]- ^ Linux用戶命令(User Commands)手冊頁 : an arbitrary precision calculator –
- ^ 2.0 2.1 Brian Kernighan and Ken Thompson. A nerdy delight for any Vintage Computer Fest 2019 attendee: Kernighan interviewing Thompson about Unix. YouTube. 事件發生在 29m45s. [September 3, 2019]. (原始內容存檔於2022-02-01).
- ^ The sources for the manual page for 7th Edition Unix dc. [2020-09-25]. (原始內容存檔於2019-09-24).
- ^ Ritchie, Dennis M. The Evolution of the Unix Timesharing System. Sep 1979 [2019-05-31]. (原始內容存檔於2010-05-06).
- ^ McIlroy, M. D. A Research Unix reader: annotated excerpts from the Programmer's Manual, 1971–1986 (PDF) (技術報告). CSTR. Bell Labs. 1987 [2019-05-31]. 139. (原始內容存檔 (PDF)於2019-11-30).