%I #10 Jun 25 2013 12:33:27
%S 1,1,1,3,1,1,3,1,3,5,1,46,1,5,66,1,1,3,3,3,10,3,1,1,172,1,1,398,5,132,
%T 255,1,336,41,1,3,615,1,462,918,283,1,1,1,453,1,11,713,1503,311,1,112,
%U 1,308,2272,1,788,481,2301,1,2686,1358,3,3,2651,3346,1
%N Numbers A226770((prime(n))^2-1).
%C a(n) = 1 iff (prime(n))^2 + prime(n) - 1 is prime (A053184).
%H Peter J. C. Moses, <a href="/A226856/b226856.txt">Table of n, a(n) for n = 1..400</a>
%t Table[(div=Most[Divisors[n+1]]; Length[FixedPoint[Union[Flatten[AppendTo[div, Map[Most[Rest[Divisors[n+#]]]&, #]]]]&, div]]-1), {n, Prime[Range[40]]^2-1}] (* _Peter J. C. Moses_, Jun 20 2013 *)
%Y Cf. A226770, A053184.
%K nonn
%O 1,4
%A _Vladimir Shevelev_, Jun 19 2013
%E More terms from _Peter J. C. Moses_, Jun 20 2013
|