Board logo

标题: [数值计算] 求Collatz猜想代码 [打印本页]
作者: email10t    时间: 2015-6-26 11:57     标题: 求Collatz猜想代码

求Collatz猜想代码
作者: pcl_test    时间: 2015-6-26 12:16

  1. @echo off

  2. set n=106

  3. :loop

  4. set /a t=%n%%%2

  5. if %t% equ 0 (echo,%n%&set /a n/=2) else (echo,%n%&set /a n=%n%*3+1)

  6. if %n% neq 1 (goto loop) else (echo,%n%&pause&exit)
复制代码
作者: email10t    时间: 2015-6-26 13:10

回复 2# pcl_test
作者: aa77dd@163.com    时间: 2015-6-26 15:24

本帖最后由 aa77dd@163.com 于 2015-6-28 00:57 编辑 
  1. @echo off

  2. for /l %%# in (1 1 99999) do call :Collatz %%#

  3. exit

  4. :Collatz

  5. set "X=%1" & set "out="

  6. :loop

  7. (

  8.     set "out=%out%%X%,"

  9.     if %X% equ 1 echo,%out%1 & echo, & exit /b

  10.     set /a "t = (X & 1) << 31 >> 31, X = t & (X << 1) + X + 1 | ~t & X >> 1"

  11.     goto :loop

  12. )

  13. exit /b
复制代码
作者: email10t    时间: 2015-7-1 13:27

  1. 1

  2. 

  3. 2,1

  4. 

  5. 3,10,5,16,8,4,2,1

  6. 

  7. 4,2,1

  8. 

  9. 5,16,8,4,2,1

  10. 

  11. 6,3,10,5,16,8,4,2,1

  12. 

  13. 7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  14. 

  15. 8,4,2,1

  16. 

  17. 9,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  18. 

  19. 10,5,16,8,4,2,1

  20. 

  21. 11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  22. 

  23. 12,6,3,10,5,16,8,4,2,1

  24. 

  25. 13,40,20,10,5,16,8,4,2,1

  26. 

  27. 14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  28. 

  29. 15,46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1

  30. 

  31. 16,8,4,2,1

  32. 

  33. 17,52,26,13,40,20,10,5,16,8,4,2,1

  34. 

  35. 18,9,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  36. 

  37. 19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  38. 

  39. 20,10,5,16,8,4,2,1

  40. 

  41. 21,64,32,16,8,4,2,1

  42. 

  43. 22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  44. 

  45. 23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1

  46. 

  47. 24,12,6,3,10,5,16,8,4,2,1

  48. 

  49. 25,76,38,19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  50. 

  51. 26,13,40,20,10,5,16,8,4,2,1

  52. 

  53. 27,82,41,124,62,31,94,47,142,71,214,107,322,161,484,242,121,364,182,91,274,137,4

  54. 12,206,103,310,155,466,233,700,350,175,526,263,790,395,1186,593,1780,890,445,133

  55. 6,668,334,167,502,251,754,377,1132,566,283,850,425,1276,638,319,958,479,1438,719

  56. ,2158,1079,3238,1619,4858,2429,7288,3644,1822,911,2734,1367,4102,2051,6154,3077,

  57. 9232,4616,2308,1154,577,1732,866,433,1300,650,325,976,488,244,122,61,184,92,46,2

  58. 3,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1

  59. 

  60. 28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  61. 

  62. 29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  63. 

  64. 30,15,46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1

  65. 

  66. 31,94,47,142,71,214,107,322,161,484,242,121,364,182,91,274,137,412,206,103,310,1

  67. 55,466,233,700,350,175,526,263,790,395,1186,593,1780,890,445,1336,668,334,167,50

  68. 2,251,754,377,1132,566,283,850,425,1276,638,319,958,479,1438,719,2158,1079,3238,

  69. 1619,4858,2429,7288,3644,1822,911,2734,1367,4102,2051,6154,3077,9232,4616,2308,1

  70. 154,577,1732,866,433,1300,650,325,976,488,244,122,61,184,92,46,23,70,35,106,53,1

  71. 60,80,40,20,10,5,16,8,4,2,1

  72. 

  73. 32,16,8,4,2,1

  74. 

  75. 33,100,50,25,76,38,19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  76. 

  77. 34,17,52,26,13,40,20,10,5,16,8,4,2,1

  78. 

  79. 35,106,53,160,80,40,20,10,5,16,8,4,2,1

  80. 

  81. 36,18,9,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  82. 

  83. 37,112,56,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  84. 

  85. 38,19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  86. 

  87. 39,118,59,178,89,268,134,67,202,101,304,152,76,38,19,58,29,88,44,22,11,34,17,52,

  88. 26,13,40,20,10,5,16,8,4,2,1

  89. 

  90. 40,20,10,5,16,8,4,2,1

  91. 

  92. 41,124,62,31,94,47,142,71,214,107,322,161,484,242,121,364,182,91,274,137,412,206

  93. ,103,310,155,466,233,700,350,175,526,263,790,395,1186,593,1780,890,445,1336,668,

  94. 334,167,502,251,754,377,1132,566,283,850,425,1276,638,319,958,479,1438,719,2158,

  95. 1079,3238,1619,4858,2429,7288,3644,1822,911,2734,1367,4102,2051,6154,3077,9232,4

  96. 616,2308,1154,577,1732,866,433,1300,650,325,976,488,244,122,61,184,92,46,23,70,3

  97. 5,106,53,160,80,40,20,10,5,16,8,4,2,1

  98. 

  99. 42,21,64,32,16,8,4,2,1

  100. 

  101. 43,130,65,196,98,49,148,74,37,112,56,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,

  102. 8,4,2,1

  103. 

  104. 44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  105. 

  106. 45,136,68,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  107. 

  108. 46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1

  109. 

  110. 47,142,71,214,107,322,161,484,242,121,364,182,91,274,137,412,206,103,310,155,466

  111. ,233,700,350,175,526,263,790,395,1186,593,1780,890,445,1336,668,334,167,502,251,

  112. 754,377,1132,566,283,850,425,1276,638,319,958,479,1438,719,2158,1079,3238,1619,4

  113. 858,2429,7288,3644,1822,911,2734,1367,4102,2051,6154,3077,9232,4616,2308,1154,57

  114. 7,1732,866,433,1300,650,325,976,488,244,122,61,184,92,46,23,70,35,106,53,160,80,

  115. 40,20,10,5,16,8,4,2,1

  116. 

  117. 48,24,12,6,3,10,5,16,8,4,2,1

  118. 

  119. 49,148,74,37,112,56,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1

  120. 

  121. ...
复制代码
作者: 老刘1号    时间: 2019-4-21 12:30

水，汇编
  1. ;Collatz猜想

  2. ;Code By 老刘

  3. ;环境：MASM32 SDK

  4. ;编译指令：ml /coff Collatz猜想.asm /link /subsystem:console

  5. ;参数：需要验证猜想的数。

  6. ;没写高精度，所以输入或计算中途的超过2^32 - 1的数会有溢出。

  7. ;参考：MHL批处理标准教程之3n+1猜想。

  8. 

  9. 

  10. Include masm32rt.inc

  11. Main Proto

  12. 

  13. .Data?

  14. 	bBuffer DB 11 Dup (?)

  15. 	dwInput DD ?

  16. 	dwCrLfZero DD ?

  17. 	dwJumpPlusOne DD ?

  18. ;End Data?

  19. 

  20. .Code

  21. 	Main Proc

  22. 		;获得参数

  23. 		Invoke ArgClC,1,Offset bBuffer

  24. 		Sub Esp,4

  25. 		Invoke atodw_ex,Offset bBuffer

  26. 		Add Esp,4

  27. 		Mov dwInput,Eax

  28. 		

  29. 		Push Eax

  30. 		Invoke StdOut,cfm$("{")

  31. 		Pop Eax

  32. 		;开始标记

  33. 		Xor Ecx,Ecx	;计数器

  34. 		.While Eax > 1

  35. 			PushAD

  36. 			Invoke dw2a,Eax,Offset bBuffer

  37. 			Invoke StdOut,Offset bBuffer

  38. 			Invoke StdOut,cfm$(",")

  39. 			PopAD

  40. 			Inc Ecx

  41. 			Mov Ebx,Eax

  42. 			Shr Eax,1

  43. 			Mov Edx,Eax	;除以2并取整数部分的结果

  44. 			Shl Eax,1

  45. 			.If Eax == Ebx

  46. 				Mov Eax,Edx

  47. 			.Else

  48. 				Lea Eax,[2*Ebx+Ebx+1]

  49. 			.EndIf

  50. 		.EndW

  51. 		

  52. 		Push Ecx

  53. 		Invoke StdOut,cfm$("1}\n")

  54. 		Invoke StdOut,cfm$("表长：")

  55. 		Pop Ecx

  56. 		Inc Ecx

  57. 		Invoke dw2a,Ecx,Offset bBuffer

  58. 		Invoke StdOut,Offset bBuffer

  59. 		Invoke StdOut,cfm$("\n")

  60. 		

  61. 		Ret

  62. 	Main EndP

  63. 	

  64. 	Start:

  65. 		Call Main

  66. 		Invoke ExitProcess,NULL

  67. 	End Start

  68. End
复制代码
  1. {1}

  2. 表长：1

  3. {2,1}

  4. 表长：2

  5. {3,10,5,16,8,4,2,1}

  6. 表长：8

  7. {4,2,1}

  8. 表长：3

  9. {5,16,8,4,2,1}

  10. 表长：6

  11. {6,3,10,5,16,8,4,2,1}

  12. 表长：9

  13. {7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  14. 表长：17

  15. {8,4,2,1}

  16. 表长：4

  17. {9,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  18. 表长：20

  19. {10,5,16,8,4,2,1}

  20. 表长：7

  21. {11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  22. 表长：15

  23. {12,6,3,10,5,16,8,4,2,1}

  24. 表长：10

  25. {13,40,20,10,5,16,8,4,2,1}

  26. 表长：10

  27. {14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  28. 表长：18

  29. {15,46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  30. 表长：18

  31. {16,8,4,2,1}

  32. 表长：5

  33. {17,52,26,13,40,20,10,5,16,8,4,2,1}

  34. 表长：13

  35. {18,9,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  36. 表长：21

  37. {19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  38. 表长：21

  39. {20,10,5,16,8,4,2,1}

  40. 表长：8

  41. {21,64,32,16,8,4,2,1}

  42. 表长：8

  43. {22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  44. 表长：16

  45. {23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  46. 表长：16

  47. {24,12,6,3,10,5,16,8,4,2,1}

  48. 表长：11

  49. {25,76,38,19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  50. 表长：24

  51. {26,13,40,20,10,5,16,8,4,2,1}

  52. 表长：11

  53. {27,82,41,124,62,31,94,47,142,71,214,107,322,161,484,242,121,364,182,91,274,137,412,206,103,310,155,466,233,700,350,175,526,263,790,395,1186,593,1780,890,445,1336,668,334,167,502,251,754,377,1132,566,283,850,425,1276,638,319,958,479,1438,719,2158,1079,3238,1619,4858,2429,7288,3644,1822,911,2734,1367,4102,2051,6154,3077,9232,4616,2308,1154,577,1732,866,433,1300,650,325,976,488,244,122,61,184,92,46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  54. 表长：112

  55. {28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  56. 表长：19

  57. {29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  58. 表长：19

  59. {30,15,46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  60. 表长：19

  61. {31,94,47,142,71,214,107,322,161,484,242,121,364,182,91,274,137,412,206,103,310,155,466,233,700,350,175,526,263,790,395,1186,593,1780,890,445,1336,668,334,167,502,251,754,377,1132,566,283,850,425,1276,638,319,958,479,1438,719,2158,1079,3238,1619,4858,2429,7288,3644,1822,911,2734,1367,4102,2051,6154,3077,9232,4616,2308,1154,577,1732,866,433,1300,650,325,976,488,244,122,61,184,92,46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  62. 表长：107

  63. {32,16,8,4,2,1}

  64. 表长：6

  65. {33,100,50,25,76,38,19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  66. 表长：27

  67. {34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  68. 表长：14

  69. {35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  70. 表长：14

  71. {36,18,9,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  72. 表长：22

  73. {37,112,56,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  74. 表长：22

  75. {38,19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  76. 表长：22

  77. {39,118,59,178,89,268,134,67,202,101,304,152,76,38,19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  78. 表长：35

  79. {40,20,10,5,16,8,4,2,1}

  80. 表长：9

  81. {41,124,62,31,94,47,142,71,214,107,322,161,484,242,121,364,182,91,274,137,412,206,103,310,155,466,233,700,350,175,526,263,790,395,1186,593,1780,890,445,1336,668,334,167,502,251,754,377,1132,566,283,850,425,1276,638,319,958,479,1438,719,2158,1079,3238,1619,4858,2429,7288,3644,1822,911,2734,1367,4102,2051,6154,3077,9232,4616,2308,1154,577,1732,866,433,1300,650,325,976,488,244,122,61,184,92,46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  82. 表长：110

  83. {42,21,64,32,16,8,4,2,1}

  84. 表长：9

  85. {43,130,65,196,98,49,148,74,37,112,56,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  86. 表长：30

  87. {44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  88. 表长：17

  89. {45,136,68,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  90. 表长：17

  91. {46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  92. 表长：17

  93. {47,142,71,214,107,322,161,484,242,121,364,182,91,274,137,412,206,103,310,155,466,233,700,350,175,526,263,790,395,1186,593,1780,890,445,1336,668,334,167,502,251,754,377,1132,566,283,850,425,1276,638,319,958,479,1438,719,2158,1079,3238,1619,4858,2429,7288,3644,1822,911,2734,1367,4102,2051,6154,3077,9232,4616,2308,1154,577,1732,866,433,1300,650,325,976,488,244,122,61,184,92,46,23,70,35,106,53,160,80,40,20,10,5,16,8,4,2,1}

  94. 表长：105

  95. {48,24,12,6,3,10,5,16,8,4,2,1}

  96. 表长：12

  97. {49,148,74,37,112,56,28,14,7,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  98. 表长：25

  99. {50,25,76,38,19,58,29,88,44,22,11,34,17,52,26,13,40,20,10,5,16,8,4,2,1}

  100. 表长：25
复制代码
作者: codegay    时间: 2019-4-21 16:11

https://zh.wikipedia.org/wiki/%E ... 9%E7%8C%9C%E6%83%B3

http://www.rosettacode.org/mw/in ... nce&redirect=no



欢迎光临 批处理之家 (http://bathome.net./) Powered by Discuz! 7.2