OFFSET
1,1
COMMENTS
Primes p such that A006577(p) is prime. - Charles R Greathouse IV, May 12 2015
LINKS
Charles R Greathouse IV, Table of n, a(n) for n = 1..10000
MATHEMATICA
pnsQ[n_]:=PrimeQ[Length[NestWhileList[If[EvenQ[#], #/2, 3#+1]&, n, #>1&]]- 1]; Select[Prime[Range[500]], pnsQ] (* Harvey P. Dale, Nov 28 2016 *)
PROG
(PARI) A006577(n)=my(s); while(n>1, n=if(n%2, 3*n+1, n/2); s++); s
select(p->isprime(A006577(p)), primes(100)) \\ Charles R Greathouse IV, May 12 2015
CROSSREFS
KEYWORD
nonn
AUTHOR
Kent Edward Hays (hayske(AT)gmail.com), Apr 08 2010
STATUS
approved