设为首页收藏本站
切换到宽版

 找回密码
 注册
搜索
批处理之家»论坛 › 批处理专区: BAT CMD DOS › BAT求助&讨论 › 求Collatz猜想代码
[新手上路]批处理新手入门导读[视频教程]批处理基础视频教程[视频教程]VBS基础视频教程[批处理精品]批处理版照片整理器
[批处理精品]纯批处理备份&还原驱动[批处理精品]CMD命令50条不能说的秘密[在线下载]第三方命令行工具[在线帮助]VBScript / JScript 在线参考
返回列表 发新帖
查看: 28682|回复: 6

[数值计算] 求Collatz猜想代码

[复制链接]
email10t
发表于 2015-6-26 11:57:06 | 显示全部楼层 |阅读模式
求Collatz猜想代码
回复

举报

pcl_test
发表于 2015-6-26 12:16:06 | 显示全部楼层
  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:47 | 显示全部楼层
回复 2# pcl_test
回复

举报

aa77dd@163.com
发表于 2015-6-26 15:24:38 | 显示全部楼层
本帖最后由 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:32 | 显示全部楼层
  1. 1

  3. 2,1

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

  7. 4,2,1

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

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

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

  15. 8,4,2,1

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

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

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

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

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

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

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

  31. 16,8,4,2,1

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

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

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

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

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

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

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

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

  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

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

  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

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

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

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

  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

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

  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

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

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

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

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

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

  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

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

  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

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

  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

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

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

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

  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

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

  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

  121. ...
复制代码
回复

举报

老刘1号
发表于 2019-4-21 12:30:55 | 显示全部楼层
水，汇编
  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猜想。


  10. Include masm32rt.inc
  11. Main Proto

  13. .Data?
  14.         bBuffer DB 11 Dup (?)
  15.         dwInput DD ?
  16.         dwCrLfZero DD ?
  17.         dwJumpPlusOne DD ?
  18. ;End Data?

  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:20 | 显示全部楼层
https://zh.wikipedia.org/wiki/%E ... 9%E7%8C%9C%E6%83%B3

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

举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|手机版|小黑屋|批处理之家 ( 渝ICP备10000708号 )

GMT+8, 2026-3-19 05:12 , Processed in 0.019466 second(s), 8 queries , File On.

Powered by Discuz! X3.5

© 2001-2026 Discuz! Team.

快速回复 返回顶部 返回列表