OFFSET
1,1
COMMENTS
If x is even f(x) = x/2 else f(x) = 3x + 1.
LINKS
Abhiram R Devesh, Table of n, a(n) for n = 1..10000
EXAMPLE
a(1) = p = 173; q = 179
c(p) = c(q) = 32
MATHEMATICA
t = Table[Length@ NestWhileList[If[EvenQ@ #, #/2, 3 # + 1] &, n, # != 1 &] - 1, {n, 10^4}]; Prime@ Flatten@ Position[#, k_ /; Length@ k == 1] &@ Map[Union@ Part[t, #] &, #] &@ Partition[#, 2, 1] &@ Prime@ Range@ 410 (* Michael De Vlieger, Sep 01 2016 *)
PROG
(PARI) A008908(n)=my(c=1); while(n>1, n=if(n%2, 3*n+1, n/2); c++); c
t=A008908(p=2); forprime(q=3, 1e4, tt=A008908(q); if(t==tt, print1(p", ")); p=q; t=tt) \\ Charles R Greathouse IV, Sep 01 2016
(Python)
import sympy
def lcs(n):
....a=1
....while n>1:
........if n%2==0:
............n=n//2
........else:
............n=(3*n)+1
........a=a+1
....return(a)
m=2
while m>0:
....n=sympy.nextprime(m)
....if lcs(m)==lcs(n):
........print(m, )
....m=n
# Abhiram R Devesh, Sep 02 2016
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Abhiram R Devesh, Aug 15 2016
STATUS
approved