Erlang: Running time = 8.65s
+-%echo
echo(Message,Sender)->Sender ! Message.
+-%prime_list
prime_list()->
case whereis(primeList) of
undefined->
Pid=spawn(euler,prime_list,[initialize]),
register(primeList, Pid),
Pid;
Exists->
Exists
end.
prime_list(initialize)->
put(primes,array:from_list([2,3,5,7,11,13,17,19])),
put(total,8),
prime_list(idle);
prime_list(idle)->
Size=get(total),
receive
{get,FromZero,Reply} when FromZero < Size ->
Reply ! {prime,FromZero,array:get(FromZero,get(primes))};
{get,FromZero,Reply} ->
echo({get,FromZero,Reply},self()),
prime_list(gen);
{all_below,What,Reply} ->
case array:get(Size-1,get(primes)) < What of
true->
prime_list(gen),
echo({all_below,What,Reply},self());
false->
Reply ! {primes_below,What,prime_list(below,What,0,[])}
end;
{is_prime,What,Reply}->
Skirt=math:sqrt(What),
Largest=array:get(Size-1,get(primes)),
case Largest < What of
true->
case Largest < Skirt of
true->
prime_list(gen),
echo({is_prime,What,Reply},self());
false->
Reply ! {is_prime,What,prime_list(check_factors,0,What,Skirt)}
end;
false->
Reply ! {is_prime,What,prime_list(binary_search,What,0,Size-1)}
end
end,
prime_list(idle);
prime_list(gen)->
Size=get(total),
Start=array:get(Size-1,get(primes))+2,
End=Start*10+1,
Sqrt=math:sqrt(End),
AsList=array:to_list(get(primes)),
[2|Odds]=AsList,
RelevantOdds=lists:takewhile(fun(N)->N =< Sqrt end,Odds),
Newbs=prime_list(gen,Size,Start,End,RelevantOdds,ordsets:from_list(lists:seq(Start,End,2))),
put(primes,array:from_list(lists:append(AsList,ordsets:to_list(Newbs)))),
put(total,get(total)+ordsets:size(Newbs)).
prime_list(check_factors,NextIndex,N,Sqrt)->
NextPrime=array:get(NextIndex,get(primes)),
case NextPrime > Sqrt of
true->true;
false->
case N rem NextPrime ==0 of
true->false;
false->
prime_list(check_factors,NextIndex+1,N,Sqrt)
end
end;
prime_list(binary_search,What,Low,High) when High-Low < 2 ->
H=array:get(High,get(primes)),
L=array:get(Low,get(primes)),
if
H == What; L == What ->
true;
true->false
end;
prime_list(binary_search,What,Low,High)->
Mid=Low+((High-Low) div 2),
Got=array:get(Mid,get(primes)),
if
What > Got->prime_list(binary_search,What,Mid+1,High);
What < Got->prime_list(binary_search,What,Low,Mid-1);
true->true
end;
prime_list(below,What,Index,SoFar)->
Next=array:get(Index,get(primes)),
if
Next >= What->
lists:reverse(SoFar);
true->
prime_list(below,What,Index+1,[Next|SoFar])
end.
prime_list(gen,_,_,_,[],Survivors)->Survivors;
prime_list(gen,Total,Lowest,Highest,[MyPrime|RestPrimes],Survivors)->
prime_list(gen,Total,Lowest,Highest,RestPrimes,
ordsets:subtract(Survivors,
ordsets:from_list(
lists:seq(Lowest-(Lowest rem MyPrime),Highest,MyPrime)))).
p50()->
prime_list() ! {all_below,1000000,self()},
put(primes,receive {primes_below,1000000,X}->X end),
L=p50(initial,get(primes),0,0),
put(prime_set,sets:from_list(get(primes))),
p50(L).
p50(Length)->
{First,Rest}=lists:split(Length,get(primes)),
case p50(search,Rest,queue:from_list(First),lists:sum(First)) of
{found,Ans}->
io:format("~w~n",[Ans]);
false->
p50(Length-1)
end.
p50(search,_,_,Total)when Total >= 1000000->false;
p50(search,[NextPrime|Others],Current,Total)->
case sets:is_element(Total,get(prime_set)) of
true->
{found,Total};
false->
{{value,Pop},Future1}=queue:out(Current),
Future2=queue:in(NextPrime,Future1),
p50(search,Others,Future2,Total-Pop+NextPrime)
end;
p50(initial,[NextPrime|_],Sum,L) when Sum+NextPrime >= 1000000 -> L;
p50(initial,[NextPrime|Others],Sum,L)->p50(initial,Others,Sum+NextPrime,L+1).