cp's OEIS Frontend

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.

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.

This page as a plain text file.
%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