%I #7 Jun 08 2013 01:44:15
%S 10,2,1,1,9,1,1,4,2,2,5,28,3,8,1,1,1,17,2,1,27,1,1,34,7,2,4,12,4,3,2,
%T 16,2,2,1,1,1,1,12,4,9,1,33,1,6,12,1,26,2,16,11,5,21,4,2,2,6,8,15,2,3,
%U 6,1,11,3,27,2,4,1,15,2,1,1,3,12,2,2,1,8,2,7,3,6,3,16,11,4,2,25,8,4,10
%N Number of steps needed to reach a prime when the following map is repeatedly applied to n: if n is even then 2n + int(sqrt(n)) + 1, otherwise 2n - int(sqrt(n)) - 1; or -1 if no prime is ever reached.
%e For n=3, a(3)=2 because 3 -> 4 -> 11
%o (UBASIC)
%o 10 cls
%o 30 for I=2 to 100
%o 32 H=I
%o 40 if odd(H)=1 then goto 90 else goto 50
%o 50 A=2*H+int(sqrt(H))+1:K=K+1
%o 60 if prmdiv(A)=A then print I,K:goto 120
%o 65 if K>1000 then print I,0:goto 120
%o 70 H=A:goto 40
%o 90 A=2*H-int(sqrt(H))-1:K=K+1
%o 100 if prmdiv(A)=A then print I,K:goto 120
%o 105 if K>1000 then print I,0:goto 120
%o 110 H=A:goto 40
%o 120 K=0
%o 130 next
%K easy,nonn
%O 2,1
%A _Felice Russo_, Sep 02 2002
|