This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.
%I A072086 #22 Feb 09 2018 09:12:51 %S A072086 0,1,2,1,2,2,3,1,2,2,3,2,3,3,2,1,2,2,3,2,3,3,2,2,2,3,2,3,2,2,3,1,3,2, %T A072086 3,2,3,3,3,2,3,3,2,3,2,2,3,2,3,2,2,3,2,2,3,3,3,2,3,2,3,3,3,1,3,3,3,2, %U A072086 2,3,2,2,3,3,2,3,3,3,3,2,2,3,2,3,2,2,2,3,2,2,3,2,3,3,3,2 %N A072086 Number of steps to reach 1, starting with n and applying the A072084-map repeatedly. %H A072086 Reinhard Zumkeller, <a href="/A072086/b072086.txt">Table of n, a(n) for n = 1..10000</a> %e A072086 n=21: '11'*'111' -> '10'*'11' -> '11'*'1' -> '1'; i.e., 21=3*7 -> 6=2*3 -> 2*1 -> 1*1=1, therefore a(21)=3. %p A072086 A072086 := proc(n) local c,i,j; c :=0; i := n; %p A072086 while i > 1 do i:=A072084(i); c:=c+1 od; c end: %p A072086 # Note that this gives A072086(0)=0 if desired %p A072086 # without any additional case discrimination. %p A072086 # _Peter Luschny_, Jan 16 2010 %t A072086 b[1] = 1; b[n_] := Times @@ Power @@@ (FactorInteger[n] /. {p_Integer, e_} :> {DigitCount[p, 2, 1], e}); a[n_] := Length[FixedPointList[b, n]] - 2; Array[a, 100] (* _Jean-François Alcover_, Feb 09 2018 *) %o A072086 (Haskell) %o A072086 a072086 n = fst $ %o A072086 until ((== 1) . snd) (\(i, x) -> (i + 1, a072084 x)) (0, n) %o A072086 -- _Reinhard Zumkeller_, Feb 10 2013 %Y A072086 Cf. A072085. %K A072086 nonn %O A072086 1,3 %A A072086 _Reinhard Zumkeller_, Jun 14 2002