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.

A213401 Smallest number k such that the sum of divisors of k equals n times a square > 1.

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