来写几段天书式的dc程序

星条旗

喜欢用dc写程序的有两种人：

1.那些热衷搞机械指令集的老派的程序员。

2.喜欢让别人和自己都读不懂自己所写程序的偏执狂。

那段出名到被维基百科的dc programming页码中列出的，神秘难懂，用来求质数的dc程序，额估计就是后者写的。

http://bbs.chinaunix.net/thread-316197-2-1.html

但是由于dc 语言独特的语法模型结构，花一点时间深入一下还是有好处的。

对dc程序命令和逆波兰一窍不通的童鞋先去看入门介绍：

http://bbs.chinaunix.net/thread-2303294-1-1.html

http://en.wikipedia.org/wiki/Dc_(computer_program)

http://www.mkssoftware.com/docs/man1/dc.1.asp

1. ok，首先，先来做做热身：
   
2. 
   
3. dc为逆波兰式运算，k命令将主栈栈顶的值弹出，并设置为计算精度 ，p命令打印栈顶值到终端
   
4. 
   
5. dc -e '2k 1 2 3 * / p'
   
6. 
   
7. 
   
8. 
   
9. r命令 交换主栈栈顶头两个值的位置。
   
10. 
    
11. dc -e '2k 2 3 * 1 r / p'
    
12. 
    
13. 
    
14. 
    
15. ？ 表示用户可键入标准输入，sA和 lA为对寄存器A的操作，用来释放和储存临时变量，寄存器也是一个栈。
    
16. 
    
17. 此程序用来计算欧拉常数E，输入值为迭代次数。
    
18. 
    
19. dc -e '100k ? sAlA1r/1+lA^ p'
    
20. 
    
21. 
    
22. 
    
23. 下面展示了一个基本的dc程序循环流程，类似于c中的for或while。
    
24. 
    
25. z命令取到主栈长度并将该值压入主栈，d命令为复制栈顶值并将该值压入主栈。
    
26. 
    
27. <J为条件判别，当主栈顶的头两个值a，b符合b<a时弹出 a ,b 并跳转到寄存器J，执行J的顶的命令。
    
28. 
    
29. x为执行主栈顶部[ ... ]中的命令，并将[ ... ]弹出。
    
30. 
    
31. 此程序计算了主栈中所有数值的和。
    
32. 
    
33. dc -e '1 2 3 0 4 0.5 2 k z [sR + lR 1- d1<J]sJlJx sK p'
    
34. 
    
35. 
    
36. 
    
37. 注意 dc所接受的输入也必须是逆波兰形式。
    
38. 
    
39. f命令表示打印整个主栈。
    
40. 
    
41. dc -e '?  f'
    
42. 
    
43. 1 0 1 -
    
44. 
    
45. -1
    
46. 
    
47. 1
    
48. 
    
49. dc -e '?  f'
    
50. 
    
51. 1 2 3 + 1
    
52. 
    
53. 1
    
54. 
    
55. 5
    
56. 
    
57. 1
    
58. 
    
59. 
    
60. 
    
61. dc -e '1[d 1+ d3>R ]sRlRx f '
    
62. 
    
63. 下面分析一下这个程序，列个表来表述最清楚。。
    
64. 
    
65. 命令           1         [d 1+ d3>R ]sR       lR                   x            d 1+ d3       >R        。。。
    
66. 
    
67. 主栈            1             1               1  [d 1+ d3>R ]       1           1 2 2 3      1 2       。。。
    
68. 
    
69. 寄存器R       空          [d1+ d3>R ]      [d..不变]           [d..不变]       [d..不变]   [d..不变]         
    
70. 
    
71. 
    
72. 
    
73. 继续上面的表格，程序继续运行寄存器R的栈顶的 [d1+ d3>R]内容     
    
74. 
    
75. 命令           d 1+ d3          >R         循环结束      f
    
76. 
    
77. 主栈          1 2 3  3  3      1 2 3                   1 2 3
    
78. 
    
79. 寄存器R      [d..不变]       [d..不变]                [d..不变]
    
80. 
    
81. 所以最后到f命令打印主栈内容的结果到终端，最后是 3 2 1
    
82. 
    
83. 
    
84. 
    
85. 我的运行环境
    
86. 
    
87. simone@Infinite-Inspiron:~/shell$ uname -a
    
88. Linux Infinite-Inspiron 3.5.0-17-generic #28-Ubuntu SMP Tue Oct 9 19:31:23 UTC 2012 x86_64 x86_64 x86_64 GNU/Linux
    
89. 
    
90. simone@Infinite-Inspiron:~/shell$ bash --version
    
91. GNU bash, version 4.2.37(1)-release (x86_64-pc-linux-gnu)
    
92. Copyright (C) 2011 Free Software Foundation, Inc.
    
93. 
    
94. simone@Infinite-Inspiron:~/shell$ uname -a
    
95. Linux Infinite-Inspiron 3.5.0-17-generic #28-Ubuntu SMP Tue Oct 9 19:31:23 UTC 2012 x86_64 x86_64 x86_64 GNU/Linux
    
96. 
    
97. 如果大家还感兴趣的话，额会继续来写几段稍微复杂点的：）

复制代码

hbmhalley

回复 1# 星条旗


    yo ~~~~~~~~~
    wel - come ~

星条旗

回复 2# hbmhalley

拜读过您的帖子，很有启发！

星条旗

好，先拿projecteuler上的一道简单题目开刀：
The series, 1^1 + 2^2 + 3^3 + ... + 10^10 = 10405071317.

Find the last ten digits of the series, 1^1 + 2^2 + 3^3 + ... + 1000^1000.
（计算上面和式的后十位）

http://projecteuler.net/problem=48

简单的用dc加起来就行，没有什么算法：

1. dc -e ' 1[ s! l!d^  l!1+ l!1000>w]swlwxs# z [s@ + l@ 1- d1<<]s<l<xs# f'

复制代码

hbmhalley

回复 4# 星条旗


    你真的 ... 算了 1000^1000 ...
    其实有个 | 运算符，专门算模意义下的幂

1. dc -e '?10 10^sw[d1-d0<a rdlw|+lw%]dsaxp'

复制代码

星条旗

回复 5# hbmhalley


汗，学习了！
以后请多多指教！


   

星条旗

用级数展开公式，计算正弦函数值：

1. dc -e '100k ? s~ 1[ s! l~ 2 l! * 1 + ^ 2 l! * 1 + [q]s.[d1=.d1-l_x*]ds_x / 0 1 - l! ^ * l!1+ l!20>w]swlwxs= l~ z [s@ + l@ 1- d1<<]s<l<xs= p'

复制代码

输入 1 回车 即得到 sin（1）的值
由于 dc所能接受的输入必须是逆波兰式输入，所以要计算 sin（-1）的话需键入：
0 1 -

hbmhalley

回复 7# 星条旗


    0 1 -    ====     _1

waker

牛!学习了!

ecjtubaowp

不实用，感觉。