批处理快速计算勾股素数

plp626

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

复制代码