ie>Вчера вечером о выводах как-то недумалось — спать хотелось. Сегодня утром сделал вывод, что надо было плюнуть на все и идти таки спать

ie>Вообщем вывод который я для себя сделал, пора завязывать с пенесометрией! Вот только alpha-blend на J сделаю, уже пообещал Хропову

Lisp  - от 2 до 8,5 часов
C/C++ - от 3 до 25 часов
Java  - от 4 до 63 часов

RK>>Ну будет не так кратко, но что быстрее в 4-8 раз — почти не сомневаюсь

ie>    def eratosthenes = (p, i) => 
ie>        foreach (j in [i,i+p..arr.Length-1])
ie>            arr[j] = true;
ie>    def addCulculated = (s, e) => 
ie>        foreach (j in [s..e-1])
ie>            when (!arr[j])
ie>                lst.Add(num(j));
ie>

genMap 2 3 = [5,7] -- special case, even p0^2
genMap p0 p1 = arrayFilter $ accum (\_ e -> e)
    (listArray (1, bound-1) $ repeat True :: UArray Int Bool)
    [ (j, False)
    | i <- takeWhile (<p1) $ tail primes
    , j <- takeWhile (<bound) $
        iterate (+i) $ (correct (lo-rem lo i+i) i - lo) `quot` 2
    ]
    where
    lo = p0^2; hi = p1^2
    bound = (hi-lo) `quot` 2
    correct n i
        | even n = n+i
        | otherwise = n
    arrayFilter a = intern (lo+2) 1
        where intern n i
            | i==bound = []
            | a!i = n : intern (n+2) (i+1)
            | otherwise = intern (n+2) (i+1)

import Data.Array.IO
import System.IO.Unsafe

genMap 2 3 = [5,7]
genMap p0 p1 = unsafePerformIO $ do
    a <- newArray (1, bound) True :: IO (IOUArray Int Bool)
    let erase i = intern $ (correct (lo-rem lo i+i) i - lo) `quot` 2
        where
        intern j
            | j>bound = return ()
            | otherwise = do writeArray a j False; intern $! j+i
        correct n i
            | even n = n+i
            | otherwise = n
    mapM_ erase $ takeWhile (<p1) $ tail primes
    let extract n i acc
        | i==0 = return acc
        | otherwise = do
            f <- readArray a i
            if f
                then ((extract $! n-2) $! i-1) $ n:acc
                else ((extract $! n-2) $! i-1) acc
    extract (hi-2) bound [] >>= return
    where
    lo = p0^2; hi = p1^2
    bound = (hi-lo) `quot` 2 - 1