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 A213401 #7 Aug 27 2013 10:47:43 %S A213401 3,7,6,22,19,14,12,21,22,27,43,33,198,28,24,66,67,30,98,57,44,197,367, %T A213401 42,343,63,85,91,463,54,48,93,86,202,76,66,511,111,99,120,163,60,1285, %U A213401 129,88,274,751,105,364,199,134,198,211,102,763,84,147,346,1765 %N A213401 Smallest number k such that the sum of divisors of k equals n times a square > 1. %C A213401 Smallest k such that sigma(k) = n*q^2. %H A213401 Michel Lagneau, <a href="/A213401/b213401.txt">Table of n, a(n) for n = 1..1500</a> %e A213401 a(8) = 21 because the sum of the divisors of 21 is 1 + 3 + 7 +21 = 32 = 8*4 where 4 is a square. %p A213401 a:= proc(n) local k, p; %p A213401 for k from 1 while irem(sigma(k), n, 'p')>0 or %p A213401 sqrt(p)<>floor(sqrt(p)) or p=1 do od; k %p A213401 end: %p A213401 seq (a(n), n=1..100); %t A213401 snk[n_]:=Module[{k=2,c},c=Sqrt[DivisorSigma[1,k]/n];While[!IntegerQ[c] || c==1,k++;c=Sqrt[DivisorSigma[1,k]/n]];k]; Array[snk,60] (* _Harvey P. Dale_, Aug 27 2013 *) %Y A213401 Cf. A000203, A090688, A213386. %K A213401 nonn %O A213401 1,1 %A A213401 _Michel Lagneau_, Jun 10 2012