OFFSET
1,1
COMMENTS
A006666: Number of halving steps to reach 1 in '3x+1' problem.
A006667: Number of tripling steps to reach 1 in '3x+1' problem.
We observe an interesting property: the subsequence {b(i)} of perfect squares is {49, 144, 196, 576, 3844, 12544, 15376, 51529, 61504, 246016, ...} with the property that b(3) = 4*b(1), b(4) = 4*b(2), b(7) = 4*b(5), b(9) = 4*b(7), b(10) = 4*b(9), ...
The primes of the sequence are 2, 2953, 3739, 9931, 38303, 44641, ...
EXAMPLE
MAPLE
nn:=10^6:U:=array(1..nn):V:=array(1..nn):
for i from 1 to nn do:
m:=i:it0:=0:it1:=0:
for j from 1 to nn while(m<>1) do:
if irem(m, 2)=0
then
m:=m/2:it0:=it0+1:
else
m:=3*m+1:it1:=it1+1:
fi:
od:
U[i]:=it0:V[i]:=it1:
od:
for n from 1 to 100 do:
ii:=0:
for k from 1 to nn while(ii=0) do:
if U[k]=U[n]+U[n+1] and V[k]=V[n]+V[n+1]
then
ii:=1:printf(`%d, `, k):
else
fi:
od:
od:
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Lagneau, Oct 11 2016
EXTENSIONS
Name edited by Michel Marcus, Sep 13 2017
STATUS
approved