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 A066301 #27 Jan 05 2024 02:50:08 %S A066301 0,0,0,1,0,0,0,2,1,0,0,1,0,0,0,3,0,1,0,1,0,0,0,2,1,0,2,1,0,0,0,4,0,0, %T A066301 0,1,0,0,0,2,0,0,0,1,1,0,0,3,1,1,0,1,0,2,0,2,0,0,0,1,0,0,1,5,0,0,0,1, %U A066301 0,0,0,2,0,0,1,1,0,0,0,3,3,0,0,1,0,0,0,2,0,1 %N A066301 a(n) = 0 if n is squarefree, otherwise 1 + a(n/rad(n)) where rad = A007947 (squarefree kernel). %C A066301 This sequence is not the same as A046660. %H A066301 Reinhard Zumkeller, <a href="/A066301/b066301.txt">Table of n, a(n) for n = 1..10000</a> %F A066301 a(n) = A051903(n)-1 for n > 1, a(1) = 0. %F A066301 Asymptotic mean: Limit_{m->oo} (1/m) * Sum_{k=1..m} a(k) = A033150 - 1 = 0.705211... . - _Amiram Eldar_, Jan 05 2024 %e A066301 a(24) = 1 + a(24/rad(24)) = 1 + a(24/6) = 1 + a(4) = 1 + (1+a(4/rad(4))) = 1 + (1+a(4/2)) = 2 + a(2) = 2 + 0 = 2. %t A066301 a[n_] := Max[FactorInteger[n][[;;, 2]]] - 1; Array[a, 100] (* _Amiram Eldar_, Jan 05 2024 *) %o A066301 (PARI) a(n)=if(n>1, vecmax(factor(n)[,2])-1, 0) \\ _Charles R Greathouse IV_, Jul 15 2013 %o A066301 (Haskell) %o A066301 a066301 1 = 0 %o A066301 a066301 n = a051903 n - 1 -- _Reinhard Zumkeller_, Jul 23 2013 %Y A066301 Cf. A005117, A007947, A003557, A008966, A033150, A046660, A051903. %K A066301 nonn,easy %O A066301 1,8 %A A066301 _Reinhard Zumkeller_, Jan 01 2002