dc (UNIX)

dc におけるレジスタは、マクロや条件式の基本となる仕組みである。それぞれのレジスタは 1 文字の名前を持っている。レジスタを利用して値を保存したりスタックに取り出したりすることができる。sc コマンドはスタックからポップした要素をレジスタ c に保存する。lc コマンドは c レジスタの値をスタックにプッシュする。以下に示すのはレジスタを利用する例である:

3 sc 4 lc * p

レジスタをスタックとして利用することもできる。メインのスタックとレジスタのスタックとの間でプッシュ・ポップするには S・L コマンドを利用する。

dc では [ と ] で囲まれたものを文字列値として扱う。文字列も数値と同様にスタックにプッシュしたりレジスタに格納したりできる。a コマンドはスタックから要素をポップし、それが数値である場合は下位バイトを ASCII 文字に変換してプッシュし、文字列の場合は最初の 1 文字をプッシュする。文字列を x コマンドでマクロとして実行したり P コマンドで表示することはできるが、文字列を構築したり文字列操作を行う方法はない。

また、dc では # から行末までをコメントとして扱い、実行には影響を及ぼさない。

dc ではレジスタやスタックに数値だけでなく、文字列を格納できるようにすることでマクロ機能を実現している。文字列は表示するだけでなく、それ自身を dc コマンドとして実行することもできる。以下に示すのは 1 を加算して 2 を掛けるマクロをレジスタ m に格納する例である:

[1 + 2 *] sm

そして次のようにマクロを呼び出す (x コマンドはスタックからポップした値をマクロとして実行する):

3 lm x p

マクロの仕組みを利用して条件分岐を行うことができる。=r コマンドはスタックから要素を 2 つポップし、それらの値が等しいときにレジスタ r に格納されているマクロを実行する。以下に示すのはスタックの先頭に 5 があるときに Equal と表示する例である:

[[Equal]p]sm 5 =m

その他にも、> (もともとの先頭が大きい場合にマクロ実行)、!> (大きくない場合)、< (小さい場合)、!< (小さくない場合)、!= (等しくない場合) がある。

dc でループそのものを行う方法はないが、マクロを再帰的に実行することで再現できる。スタックの先頭の要素の階乗を求めるには次のようにする:

# F(x): return x!
# if x-1 > 1
#     return x * F(x-1)
# otherwise
#     return x
[d1-d1<F*]dsFxp

1Q コマンドはマクロから脱出する。これに対して q コマンドはそのマクロとそのマクロを呼び出したマクロから脱出する。マクロの呼び出し深度が 2 未満の場合は dc プログラムをも終了する。z コマンドは現在の (z コマンドの実行前の) スタックの深さをプッシュする。

これはレジスタ a に格納されたマクロを再帰的に呼び出してスタック全体の合計を計算する。マクロの動作は、まずスタックの先頭 2 を足し合わせ (+)、スタックの深さが 1 より大きい場合は再帰する (1z<a) ものである。マクロ文直前の 0d は一見必要ないように思えるがが、加算演算子がスタックの先頭 2 の要素を確実に取り出せることを保証するために必要である。

1 2 4 8 16 100  # 加算したい数をスタックに積む
0d[+z1<a]dsaxp  # すべて加算して表示

実行結果は 131 となる。

ファイルの各行に記述された dc コマンドを実行し、その合計値を計算する。単なる数値は dc コマンドとして有効であるので、行ごとに記述された数値の合計を得ることもできる。

これも上の例と同様にレジスタ a に格納されたマクロを再帰的に呼び出して合計を計算している。

$ cat file | dc -e '0d[?+z1<a]dsaxp'

? コマンドは標準入力から dc コマンドを読み込み、結果をスタックにプッシュする。ファイルの終端に達した場合はコマンドは読み込まれず、スタックには何もプッシュされない。

$ { echo "5"; echo "7"; } | dc -e '0d[?+z1<a]dsaxp'

この結果は 12 となる。 入力の dc コマンドは複雑でも良い。

$ { echo "3 5 *"; echo "4 3 *"; echo "5dd++"; } | dc -e '0d[?+z1<a]dsaxp'

この結果は ${\displaystyle (3\times 5)+(4\times 3)+(5+5+5)=42}$ より 42 となる。

この手法を利用することの長所は、dc の演算は無限精度であるので、 AWK などの簡明な手法を用いる場合に起こりうるオーバーフローや精度落ちが発生しないことである。

対して短所もあり、空行 (技術的には要素をプッシュしない行) に遭遇するとループが停止すること、負値の符号 ${\displaystyle -}$ を dc の負符号 _ に置換する必要があることが挙げられる。dc の ? には空行とファイル終端を見分けるきれいな方法が提供されていない。

以下に示すのは、ワンライナーで書かれた較的簡単な単位変換プログラムの例である:

$ dc -e '[[Enter a number (metres), or 0 to exit]psj]sh[q]sz[lhx?d0=z10k39.370079*.5+0k12~1/rn[ feet ]Pn[ inches]P10Pdx]dx'

この例はメートルで表された距離をフィートとインチに変換する。ユーザに入力を促し、適切な形式で変換結果を出力し、別の数値を変換するためにループしている。

以下に示すのは互除法を利用して最大公約数を求める例である:

$ dc -e '??[dSarLa%d0<a]dsax+p'                  # 最短実装
$ dc -e '[a=]P?[b=]P?[dSarLa%d0<a]dsax+[GCD:]Pp' # 読みやすくしたもの

入力された値 ${\displaystyle n}$ の階乗 ${\displaystyle n!=\prod _{k=1}^{n}{k}}$ を計算する:

$ dc -e '?[q]sQ[d1=Qd1-lFx*]dsFxp'

dc にも quine が存在する。

$ dc -e '[91Pn[dx]93Pn]dx'
$ dc -e '[91PP93P[dx]P]dx'

$ echo '2p3p[dl!d2+s!%0=@l!l^!<#]s#[s/0ds^]s@[p]s&[ddvs^3s!l#x0<&2+l.x]ds.x' | dc

このプログラムは Michel Charpentier によって書かれた、すべての素数を表示するプログラムである。このプログラムは次のように短縮でき、恐らくこれが最短であると思われる。

$ echo '2p3p[dl!d2+s!%0=@l!l^!<#]s#[0*ds^]s@[p]s&[ddvs^3s!l#x0<&2+l.x]ds.x' | dc

これもまた Michel Charpentier によって書かれた[6]:

$ dc -e '[n=]P?[p]s2[lip/dli%0=1dvsr]s12sid2%0=13sidvsr[dli%0=1lrli2+dsi!>.]ds.xd1<2'

短くすると:

$ dc -e "[n=]P?[lfp/dlf%0=Fdvsr]sF[dsf]sJdvsr2sf[dlf%0=Flfdd2%+1+sflr<Jd1<M]dsMx"

より高速なものでは (200 bit の数 ${\displaystyle 2^{200}-1}$ で試すなら 2 200^1- を入力する):

$ dc -e "[n=]P?[lfp/dlf% 0=Fdvsr]sFdvsr2sfd2%0=F3sfd3%0=F5sf[dlf%0=Flfd4+sflr>M]sN[dlf%0=Flfd2+sflr>N]dsMx[p]sMd1<M"

これを更に高速化するには、定数へのアクセスをレジスタへのアクセスに置き換える:

$ dc -e "[n=]P?[lfp/dlf%l0=Fdvsr]sF2s2dvsr2sf4s4d2%0=F3sfd3%0=F5sf[dlf%l0=Flfdl4+sflr>M]sN[dlf%l0=Flfdl2+sflr>N]dsMx[p]sMd1<M"

より複雑な例では、Perl スクリプトに組み込まれている Diffie–Hellman 鍵交換がある。これは ITAR 議論のときにサイファーパンクたちの間で人気であった。このスクリプトは Unix ライクな OS にはどこにでもある Perl と dc だけで実行できる[7]:

#!/usr/bin/perl -- -export-a-crypto-system-sig Diffie-Hellman-2-lines
($g, $e, $m) = @ARGV, $m || die "$0 gen exp mod\n";
print `echo "16dio1[d2%Sa2/d0<X+d*La1=z\U$m%0]SX$e"[$g*]\EszlXx+p | dc`

これのコメント付きのバージョンは少しわかりやすく、ループや条件分岐、マクロから脱出するための q コマンドなどの使い方が示されている。GNU の dc では、以下のコード中でレジスタ X に格納されるマクロを使う代わりに | コマンドで冪剰余を演算できる。

#!/usr/bin/perl

my ($g, $e, $m) = map { "\U$_" } @ARGV;
die "$0 gen exp mod\n" unless $m;

print `echo $g $e $m | dc -e '
# Hex input and output
16dio
# Read m, e and g from stdin on one line
?SmSeSg

# Function z: return g * top of stack
[lg*]sz

# Function Q: remove the top of the stack and return 1
[sb1q]sQ

# Function X(e): recursively compute g^e % m
# It is the same as Sm^Lm%, but handles arbitrarily large exponents.
# Stack at entry: e
# Stack at exit: g^e % m
# Since e may be very large, this uses the property that g^e % m == 
#     if( e == 0 )
#         return 1
#     x = (g^(e/2)) ^ 2
#     if( e % 2 == 1 )
#         x *= g
#     return x %
[
    d 0=Q   # return 1 if e==0 (otherwise, stack: e)
    d 2% Sa # Store e%2 in a (stack: e)
    2/      # compute e/2
    lXx     # call X(e/2)
    d*      # compute X(e/2)^2
    La1=z   # multiply by g if e%2==1
    lm %    # compute (g^e) % m
] SX

le          # Load e from the register
lXx         # compute g^e % m
p           # Print the result
'`;