来写几段天书式的dc程序

星条旗

用Ramanujan算法快速逼近圆周率真实值
程序接受输入：圆周率精度   迭代次数

1. dc -e '?s$k 1 [ s~ l~ 4* [q]s.[d1=.d1-l_x*]ds_x 1103 26390 l~ * + * l~ [q]s![d1=!d1-l|x*]ds|x 4^ 396 4 l~ * ^ * / l~1+ dl$><]s<l<xs# z [s@ + l@ 1- d1<<]s<l<xs# 1103+ 8v * 9801 / 1 r / f'
   
2. 
   
3. 1000 150
   
4. 
   
5. 3.1415926535897932384626433832795028841971693993751058209749445923078\
   
6. 164062862089986280348253421170679821480865132823066470938446095505822\
   
7. 317253594081284811174502841027019385211055596446229489549303819644288\
   
8. 109756659334461284756482337867831652712019091456485669234603486104543\
   
9. 266482133936072602491412737245870066063155881748815209209628292540917\
   
10. 153643678925903600113305305488204665213841469519415116094330572703657\
    
11. 595919530921861173819326117931051185480744623799627495673518857527248\
    
12. 912279381830119491298336733624406566430860213949463952247371907021798\
    
13. 609437027705392171762931767523846748184676694051320005681271452635608\
    
14. 277857713427577896091736371787214684409012249534301465495853710507922\
    
15. 796892589235420199561121290219608640344181598136297747713099605187072\
    
16. 113499999983729780499510597317328160963185950244594553469083026425223\
    
17. 082533446850352619311881710100031378387528865875332083814206171776691\
    
18. 473035982534904287554687311595628638823537875937519577818577805321712\
    
19. 268066130019278766111959092164201997

复制代码

与bc运算结果之差小于10^-998.

1. echo "scale=1000;4*a(1)"|bc -l

复制代码

blackold

NB，不得不学习!!

星条旗

回复 22# blackold

过奖了！

下面贴个老外的dc作品，是计算汉诺塔的。
运行方法：dc -e n hanoi.dc

1. #
   
2. # The Towers Of Hanoi
   
3. # dc (Unix Desk Calculator) version
   
4. # Copyright (C) 2003 Amit Singh. All Rights Reserved.
   
5. # http://hanoi.kernelthread.com
   
6. #
   
7. # Tested under GNU dc 1.3
   
8. #
   
9. 
   
10. # Usage: dc -e n hanoi.dc
    
11. #        where n is the number of disks
    
12. 
    
13. [ # move(from, to)
    
14.     n           # print from
    
15.     [ --> ]n    # print " --> "
    
16.     p           # print to\n
    
17.     sw          # p doesn't pop, so get rid of the value
    
18. ]sm
    
19. 
    
20. [ # init(n)
    
21.     sw          # tuck n away temporarily
    
22.     9           # sentinel as bottom of stack
    
23.     lw          # bring n back
    
24.     1           # "from" tower's label
    
25.     3           # "to" tower's label
    
26.     0           # processed marker
    
27. ]si
    
28. 
    
29. [ # Move()
    
30.     lt          # push to
    
31.     lf          # push from
    
32.     lmx         # call move(from, to)
    
33. ]sM
    
34. 
    
35. [ # code block <d>
    
36.     ln          # push n
    
37.     lf          # push from
    
38.     lt          # push to
    
39.     1           # push processed marker 1
    
40.     ln          # push n
    
41.     1           # push 1
    
42.     -           # n - 1
    
43.     lf          # push from
    
44.     ll          # push left
    
45.     0           # push processed marker 0
    
46. ]sd
    
47. 
    
48. [ # code block <e>
    
49.     ln          # push n
    
50.     1           # push 1
    
51.     -           # n - 1
    
52.     ll          # push left
    
53.     lt          # push to
    
54.     0           # push processed marker 0
    
55. ]se
    
56. 
    
57. [ # code block <x>
    
58.     ln 1 =M
    
59.     ln 1 !=d
    
60. ]sx
    
61. 
    
62. [ # code block <y>
    
63.     lMx
    
64.     lex
    
65. ]sy
    
66. 
    
67. [ # quit()
    
68.     q           # exit the program
    
69. ]sq
    
70. 
    
71. [ # run()
    
72.     d 9 =q      # if stack empty, quit()
    
73.     sp          # processed
    
74.     st          # to
    
75.     sf          # from
    
76.     sn          # n
    
77.     6           #
    
78.     lf          #
    
79.     -           #
    
80.     lt          #
    
81.     -           # 6 - from - to
    
82.     sl          #
    
83.     lp 0 =x     #
    
84.     lp 0 !=y    #
    
85.     lrx         # loop
    
86. ]sr
    
87. 
    
88. lix # init(n)
    
89. lrx # run()
    

复制代码

星条旗

dc or bc ?

那么到底用dc还是bc呢：）

先考察一下各自的性能

计算大数的乘方

1. 
   
2. 
   
3. time dc -e '2 88888 ^ p'
   
4. 
   
5. real 0m0.241s
   
6. user 0m0.216s
   
7. sys 0m0.008s
   
8. 
   
9. 
   
10. 
    
11. time echo "2^88888" |bc
    
12. 
    
13. real 0m0.240s
    
14. user 0m0.216s
    
15. sys 0m0.004s
    

复制代码

计算小数点后一千位圆周率的值

1. 
   
2. time dc -e '1000k 150 s$ 1 [ s~ l~ 4* [q]s.[d1=.d1-l_x*]ds_x 1103 26390 l~ * + * l~ [q]s![d1=!d1-l|x*]ds|x 4^ 396 4 l~ * ^ * / l~1+ dl$><]s<l<xs# z [s@ + l@ 1- d1<<]s<l<xs# 1103+ 8v * 9801 / 1 r / f'
   
3. 
   
4. real 0m0.710s
   
5. user 0m0.704s
   
6. sys 0m0.000s
   
7. 
   
8. 
   
9. 
   
10. time echo "scale=1000;4*a(1)"|bc -l
    
11. 
    
12. real 0m0.600s
    
13. user 0m0.596s
    
14. sys 0m0.000s
    

复制代码

计算10000的乘阶

1. 
   
2. time echo "
   
3. define f (x) {
   
4. if (x <= 1) return (1);
   
5. return (f(x-1) * x);
   
6. }
   
7. f(10000)
   
8. " | bc
   
9. 
   
10. real 0m2.270s
    
11. user 0m2.236s
    
12. sys 0m0.008s
    
13. 
    
14. time dc -e '10000[d1-d1<r*]dsrxp'
    
15. 
    
16. real 0m2.272s
    
17. user 0m2.232s
    
18. sys 0m0.004s
    
19. 
    

复制代码

就上面三个测试结果来看，dc和bc是平分秋色的，额就没有进行其他的测试了。

但是在程序的易写性和易读性来说bc等工具远胜于dc ，这也就是为什么shell

论坛中基本上没有人喜欢讨论这个古老的unix命令的原因。

看来user friendly终究才是趋势，晦涩难懂的语法结构是dc和汇编等低级语言谈出

人们视线的根本原因。

全文完^_*

rdcwayx

我考虑给这帖子加精。 非常不错， 谢谢共享。

waker

星条旗 发表于 2012-12-03 17:49
dc or bc ?

那么到底用dc还是bc呢：）

曾经:bc只是个dc编译器
现代:很多dc使用bc的库来实现

所以dc~=bc?

星条旗

waker 发表于 2012-12-04 14:00
曾经:bc只是个dc编译器
现代:很多dc使用bc的库来实现

我进行过很多测试，应该说从运行速度上来看dc和bc的性能是完全一样的。
所以我觉得虽然从代码层面dc和bc风格迥异但是，在编译过后所生成的
低级机械语言层面他们的性能是完全一样的。只不过dc的语言显得更加接近
低级语言而已。
其实仔细想想的话dc对语言这种原始的处理方式在uinx命令家族里面并不孤单。
sed对脚本的处理也有类似的地方，都是在模仿cpu处理数据的方式。
dc脚本是对主栈和寄存器进行操作。
sed是对pattern space和hold space进行操作。
cpu是对存储器和寄存器进行操作。
我觉得cpu对数据处理的模型结构，很大程度上影响了以前那些老unix程序员
的思维，这可能也就是为什么他们把dc弄的这么复杂的原因吧^_^

hbmhalley

waker 发表于 2012-12-04 14:00
曾经:bc只是个dc编译器

{:3_182:} 真的假的 ....

hbmhalley

星条旗 发表于 2012-12-03 17:45
下面贴个老外的dc作品，是计算汉诺塔的。
运行方法：dc -e n hanoi.dc

求出处

星条旗

回复 29# hbmhalley

http://www.kernelthread.com/projects/hanoi/
程序来自于一个美国geek，他用各式各样的语言和方法编写了hanoi 包括dc ,sendmail ，tex ，sed ，postscript等
甚至还有世嘉机Dreamcast版本。
你在网页下面的列表里面就可以找到dc 。