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 A264154 #12 Feb 14 2019 05:10:58 %S A264154 1,1,2,1,3,3,3,1,3,3,2,2,2,3,4,5,1,3,2,2,7,1,2,4,3,3,2,3,3,2,5,2,2,3, %T A264154 3,5,2,3,7,3,3,3,5,3,4,2,5,3,2,7,2,3,3,3,5,2,7,2,6,2,5,3,2,3,3,2,3,1, %U A264154 3,3,4,3,11,4,7,3,2,2,5,3,3,5,3,4,4,7,4 %N A264154 For numbers m such that rad(n) divides sigma(n), this sequence gives the minimum exponent k such that sigma(m)^k divides m. %H A264154 Michel Marcus, <a href="/A264154/b264154.txt">Table of n, a(n) for n = 1..10000</a> %e A264154 A175200(2) is 6, and for 6, sigma(6)^k/6 is already an integer with k=1, so a(2)=6. %e A264154 A175200(3) is 24, and for 24, sigma(24)/24 is not an integer while sigma(24)^2/24 is an integer, so a(3)=2. %o A264154 (PARI) fk(s, m) = {j = 1; while(denominator(s^j/m) != 1, j++); j;} %o A264154 rad(n) = factorback(factorint(n)[, 1]); %o A264154 lista(nn) = {for (n=1, nn, if (denominator(sigma(n)/rad(n)) == 1, k = fk(sigma(n), n); print1(k, ", ");););} %Y A264154 Cf. A000203 (sigma(n)), A007947 (rad(n)), A175200. %K A264154 nonn %O A264154 1,3 %A A264154 _Michel Marcus_, Nov 06 2015