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 A048257 #27 May 10 2018 02:47:33 %S A048257 1,93,127,11811,112890,120054,124338,127330,132770,133998,134090, %T A048257 137058,138754,139962,146710,148665,148810,149534,153986,155510, %U A048257 160215,161194,164985,167134,170986,173098,183687,184682,187143,191913,198485,206823,206965,207687 %N A048257 Integers whose sum of divisors is a 7th power. %C A048257 If m and n are coprime members of the sequence, then m*n is also a member. - _Robert Israel_, May 10 2018 %H A048257 Donovan Johnson, <a href="/A048257/b048257.txt">Table of n, a(n) for n = 1..1000</a> %H A048257 Frits Beukers, Florian Luca and Frans Oort, <a href="http://www.jstor.org/stable/10.4169/amer.math.monthly.119.05.373">Power Values of Divisor Sums</a>, The American Mathematical Monthly, Vol. 119, No. 5 (May 2012), pp. 373-380. %F A048257 sigma(a(n)) = x^7, where the initial values of x are 1, 2, 4, 6 (48 times), ... %e A048257 Divisors(11811) = {1,3,31,93,127,381,3937,11811} and sigma(11811) = 16384 = 4^7. %p A048257 filter:= n -> type(map(t -> t[2]/7, ifactors(numtheory:-sigma(n))[2]),list(integer)): %p A048257 select(filter, [$1..21*10^4]); # _Robert Israel_, May 09 2018 %t A048257 Select[Range[210000],IntegerQ[Surd[DivisorSigma[1,#],7]]&] (* _Harvey P. Dale_, Jun 09 2017 *) %o A048257 (PARI) isok(n) = ispower(sigma(n), 7); \\ _Michel Marcus_, Dec 20 2013 %Y A048257 Cf. A000203, A006532, A020477, A019422, A019423, A019424, A048251-A048256, A048258. %K A048257 nonn %O A048257 1,2 %A A048257 _Labos Elemer_