Board logo

标题: [数值计算] 批处理快速计算勾股素数 [打印本页]

作者: plp626    时间: 2009-3-7 20:18     标题: 批处理快速计算勾股素数

  1. :::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
  2. ::算法原理:http://zh.wikipedia.org/wiki/%E5%8B%BE%E8%82%A1%E6%95%B0
  3. ::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::::
  4. @echo off||Code By plp626
  5. ::勾股数中最小的一个数 < MAX^2 这里暂定MAX=30
  6. set MAX=30 & if not %1.==. goto%*
  7. echo    稍等。。。
  8. for /f "delims=" %%a in ('%~s0 :rtnum^|sort')do echo.%%a
  9. pause&exit/b
  10. :rtnum
  11. setlocal enabledelayedexpansion&set "s=        "
  12. for /l %%n in (1 1 %MAX%)do (set/a tp=%%n+1
  13. for /l %%m in (!tp! 2 %MAX%)do (call:ifprime %%m %%n prm
  14.     if !prm! equ 1 (
  15.       set/a mm=%%m*%%m,nn=%%n*%%n,a=mm-nn,b=2*%%m*%%n,c=mm+nn
  16.         set a=%s%!a!&set b=%s%!b!&set c=%s%!c!
  17.         if !a! lss !b! (echo !a:~-7! !b:~-7! !c:~-7!) else echo !b:~-7! !a:~-7! !c:~-7!
  18. )))
  19. exit/b
  20. :ifprime
  21. set/a x=%1,y=%2
  22. :loop
  23. set/a r=x%%y,x=y,y=r
  24. if %r% neq 0 (goto:loop) else set %3=%x%
复制代码





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