XProfan

Kompilieren Markieren Separieren

p.specht

Hab´s etwas weiterentwickelt zu einem Primzahlenfinder auf FLOAT-Basis...

WindowTitle "Zufallsprimzahlengenerator"
'(D)Demo 2018-05 by P.Specht/Wien, nach einer Idee von Jörg Sellmeyer

proc frac :parameters x!:var s!=(x!>0)-(x!<0):x!=abs(x!)

    x!=x!-round(x!,0):case x!<0:x!=1+x!:return s!*x!

endproc

proc intf :parameters x!

    var s!=(x!>0)-(x!<0)
    x!=abs(x!)
    x!=x!-frac(x!)
    return s!*x!

endproc

Proc IsIntf

    Parameters f!
    Return (f! - Intf(f!) = 0)

EndProc

Proc NaheMultiplikatoren

    Parameters Wert![],Zahl!
    Declare a!,f_0!,f_1!
    f_0! = Sqrt(Zahl!)
    Case Zahl!=0:Return 0

    If IsIntf(f_0!)

        Wert![0] = f_0!
        Wert![1] = f_0!

    Else

        WhileLoop 0,Intf(f_0!)

            a! = Intf(f_0!) - &Loop
            f_1! = Zahl! / a!

            If IsIntf(f_1!)

                Wert![0] = a!
                Wert![1] = Zahl! / a!
                Break

            EndIf

        Endwhile

    EndIf

    Return f_0!

EndProc

Cls
Declare Test![1],ZufallsZahl!
Randomize:set("decimals",1)

While 1

    ZufallsZahl! = intf(10000000.00000*Rnd()+rnd())
    NaheMultiplikatoren(Test![],ZufallsZahl!)

    If (Test![0]=1.0) and (Test![1]>1.0)

        cls
        Print "\n Nah-Multiplikatoren:",Test![0],"*",Test![1],"=",ZufallsZahl!
        '   If (Test![0]=1.0) and (Test![1]>1)
        font 2
        print "\n Diese Zufallszahl ist eine Primzahl!"
        sound 60,16
        font 0
        waitinput 1000

    Endif

Endwhile

XProfan 11
Computer: Gerät, daß es in Mikrosekunden erlaubt, 50.000 Fehler zu machen, zB 'daß' statt 'das'...

19.05.2018 ▲

Jörg
Sellmeyer

Da würde ich ja eher ~~Doubles~~ Quadint (meinte ich natürlich!) als Übergabeparameter nehmen.
Und wenn es nur um die Primzahlen geht, kann man das ja noch sehr beschleunigen, wenn man 0 zurückgibt, sobald der erste Multiplikant <> 1 gefunden ist.

Windows XP SP2 XProfan X4
... und hier mal was ganz anderes als Profan ...

p.specht

Stimmt, danke! Bin aber noch auf XProfan-11, warte auf die Vollversion von X4. Ausserdem läuft die Schleife dzt. nur mit Integers. Wenn´s bei Float zu lange dauert, müsste man auf Inline-Assembler ausweichen. Deine Idee ist aber super!

Computer: Gerät, daß es in Mikrosekunden erlaubt, 50.000 Fehler zu machen, zB 'daß' statt 'das'...

Jörg
Sellmeyer

Und man kann natürlich vornedran auch noch ein Sieb einbauen, in dem alle Zahlen hängenbleiben, die durch 2,3,5,7,11,13,17,23,... primzahl&& teilbar sind und einen Wert > 1 ergeben.

Windows XP SP2 XProfan X4
... und hier mal was ganz anderes als Profan ...

p.specht

..und mit Mersenne-Primzahlen arbeiten... da gibt es Formeln. Ich brauche aber echte Zufallsprimeln für RSA, da kann es ruhig etwas länger dauern

XProfan 11
Computer: Gerät, daß es in Mikrosekunden erlaubt, 50.000 Fehler zu machen, zB 'daß' statt 'das'...

Kompilieren Markieren Separieren

Jörg
Sellmeyer

Das Sieb ist ja nur, um den Vorgang zu beschleunigen, indem Zahlen schneller rausgeworfen werden, die keine Primzahlen sind.

Also:

If (Eingang&& \ Prim&&) > 1

    Print "kann keine Primzahl sein"
    return 0

ElseIf (Eingang&& \ NaechstePrim&&) > 1

    ...

EndIf

XProfan X3
Windows XP SP2 XProfan X4
... und hier mal was ganz anderes als Profan ...

p.specht

Alles klar.
So, bin dann wieder auf Paules PC-Forum / Xprofan zu finden...

XProfan 11
Computer: Gerät, daß es in Mikrosekunden erlaubt, 50.000 Fehler zu machen, zB 'daß' statt 'das'...