Matematické Fórum

↑ Arabela:↑ Arabela:
Děkuji mockrát za odpověď. Budu se tedy muset spokojit s tím, co mám :)

↑ kaja.marik:
Na výpis využívám příkazu makelist(b,b,1,N). Nevím jaký jiný příkaz použít, protože pak se využívá příkazu delete pro vymazání složených čísel.

Zdravím,
podle Wikipedie jsem sestavil tento kód:

Eratos(n):=block([i, j, L, m],
    m: ceiling(sqrt(n)),
    L: makelist(if oddp(i) then i else 0, i, 1, n),
    L[1]: 0,
    L[2]: 2,
    for i: 3 thru m step 2 do (
        if is(L[i] # 0) then (
            L[2*i]: 0,
            for j: 3*i thru n step 2*i do (
                L[j]: 0
            )
        )
    ),
    delete(0, L)
)$
compile(Eratos);

Doba výpočtu je velmi závislá na velikosti .
Pro trvá výpočet 3 sekundy, pro již 35 sekund.

Pro větší jsem to raději nezkoumal.

Pro porovnání:

Prime_List(n):=block([i, L],
    L: makelist(if primep(i) then i else 0, i, 1, n),
    delete(0, L)
)$
compile(Prime_List);

Tento kód pro zobrazení prvočísel do 100000 trvá pouze 0.3 sekundy, ale pracuje na jiném principu.

Pokud chceš řešit proč tvůj kód nefunguje jak má, měl by si jej sem umístit, jinak budeme jen hádat...