Batcher

【转帖1】批处理实现素数堆垒及哥德巴赫猜想局部验证

:: Prime.bat - Generate a serial prime number
:: Dirk van Deun - Will Sort Modified 2004/11/18
::
:: 改进: 将isprime和divided函数并入主函数以及其他一些风格上的改进
:: 效果: 函数,变量和代码均减少, 速度继续提升; 测试运行时间约 1.6秒
@echo off
if [%1]==[$] goto %2
if [%1]==[] %comspec% /e:5000 /c %0 $ init
del ~prime.bat
goto end
:: 初始化: 产生素数2, 将它存为第一个素数, 设置循环起始值
:init
echo I I
set prime-num=I
set %prime-num%=I I
set prime-in=I
:: 对3~n的奇数 %prime-in% 与已产生的所有素数由小到大循环相除
:: 若全部未整除则显示此整数, 否则递增 %prime-in% 后继续循环 
:runloop
set prime-in=I I %prime-in%
set divisor-no=I
    :divideloop
    echo set divisor=%%%divisor-no%%%>~prime.bat
    call ~prime.bat
    call %0 $ loopminus %prime-in%
    if "%min-out%"=="" goto runloop
    if "%divisor-no%"=="%prime-num%" goto isprime
    set divisor-no=I%divisor-no%
    goto divideloop
:isprime
echo %prime-in%
set prime-num=I%prime-num%
if "%prime-num%"=="IIIIIIIIII" goto end
set %prime-num%=%prime-in%
goto runloop
:: 对传入的 %dividend%(被除数) %divisor%(除数) 循环相减
:: 若不足相减 (%2!=I) 则返回下溢错, 否则直接返回空
:loopminus
for %%n in (%divisor%) do shift
if not [%3]==[] goto loopminus
set min-out=
if [%2]==[] set min-out=underflow
goto end
:end
复制代码

原文地址：http://www.cn-dos.net/forum/viewthread.php?tid=14580

Batcher

【转帖2】批处理筛选质数

@echo off
setlocal enabledelayedexpansion
::::::::::::::::::::::::::::Find Prime Numbers::::::::::::::::::::::::::::
::::::::::::::::::::::::::::{s11ss  2007-9-20}::::::::::::::::::::::::::::
:r
echo Please input the upper limit number:
set /p n=
if not !n! geq 2 (echo 2 at least. & goto :r)
echo.
echo Calculating...
set /a i=2
for /l %%a in (2,1,!n!) do (
        set m%%a=0
)
:ci
set /a j=!i!
:cj
set /a m=!i!*!j!
if !m! leq !n! (
        set /a j+=1
        set m!m!=1
        goto :cj
) else (
        set /a i+=1
        set /a ii=!i!*!i!
        if !ii! leq !n! (goto :ci) else (goto :e)
)
:e
set /a counter=0
echo.
echo In [2,!n!],prime numbers are:
for /l %%a in (2,1,!n!) do (
        if !m%%a! equ 0 (
                echo %%a
                set /a counter+=1
                )
)
echo.
echo In total:!counter!
echo.
echo Press Any Key To Exit.
pause>nul
复制代码

原文地址：http://www.cn-dos.net/forum/viewthread.php?tid=33724

Batcher

【转帖3】批处理寻找大素数-32位正整数的素性判定

以下代码完全照搬了以下链接中的判定算法——米勒拉宾检验+二次检验
http://blog.csdn.net/bsrw/archive/2006/11/28/1419145.aspx

遗憾的是其中的二次检测会导致大素数判定时数据溢出
测试结果从46400开始大量出现因为数据溢出而导致的漏检

1000以内测试没有漏检
1000~46400 没有测试

@echo off
setlocal EnableDelayedExpansion

:loop_test
set testnum=
set /p testnum=请输入一个整数(按i使用内置测试集，直接按回车退出): 
if "%testnum%"=="" goto :eof
if /i not "%testnum%"=="i" (
    call :JudgePrime %testnum%
    if errorlevel 2 (echo 无效输入：%testnum%
    ) else if errorlevel 1 (echo %testnum% 是素数
    ) else (echo %testnum% 是合数)
    goto :loop_test
)

set time0=%time%
for /l %%i in (46001,2,48000) do (
    set /a testnum=%%i
    call :JudgePrime !testnum!
    if !errorlevel! equ 1 set /p=!testnum!  <nul & set /a iprime+=1
)
echo.
echo.
echo found: %iprime%
echo begin: %time0%
echo finish: %time%
pause
goto :eof

:JudgePrime
if [%1]==[] exit /b 2
set /a tmp1=%1
if not %tmp1%==%1 exit /b 2
if %1 lss 2 exit /b 0
if %1 equ 2 exit /b 1

set i=0
for %%i in (2,3,5,7,11) do (
   set prime_!i!=%%i
   set /a prime6p_!i!=%%i*%%i*%%i
   set /a prime6p_!i!*=prime6p_!i!
   set /a i+=1
)

set i=0
:loop1_JP
call set prime_i=%%prime_%i%%%
call set prime6p_i=%%prime6p_%i%%%
if %1 geq %prime6p_3% (
    if !prime_i! equ 3 (
        set /a i+=1
        goto :loop1_JP
    ) 
) else (
    if !prime_i! neq 2 if %1 lss !prime6p_i! exit /b 1    
)
call :LikePrime %1 %prime_i%
if not errorlevel 1 exit /b 0
set /a i+=1
if %i% lss 5 goto :loop1_JP
exit /b 1
goto :eof

:LikePrime
set /a x=result=1, tmp1=%1-1, bits=0

:loop1_LP
set /a bits+=1
set /a "tmp1>>=1" 
if %tmp1% gtr 0 goto :loop1_LP

set /a tmp1=%1-1
:loop2_LP
set /a bits-=1
set /a result=(x*x) %% %1  %=此句代码判断导致大素数时数据溢出=%
if %result% equ 1 if %x% neq 1 if %x% neq %tmp1% exit /b 0
set /a "tmp2=%tmp1% & (1 << %bits%)"
if %tmp2% neq 0 set /a result=(result*%2) %% %1
set x=%result%
if %bits% gtr 0 goto :loop2_LP
if %result% equ 1 exit /b 1
goto :eof
复制代码

原文地址：http://www.cn-dos.net/forum/viewthread.php?tid=27198

Batcher

999983应该怎样分解质因数？