|
%I #21 Dec 26 2020 03:49:32
%S 6,8,13,14,15,16,17,19,20,22,23,27,29,32,42,44,46,50,54,62,69,90,92,
%T 100,104,108,110,114,130,136,148,150,152,156,166,170,176,182,184,186,
%U 198,200,202,214,230,232,234,236,240,242,244,254,258,266,272,280,286
%N Positive integers k with exactly 2 aliquot sequence predecessors. In other words, there are exactly two solutions x for which s(x) = n. The function s(x) here is the sum of all proper divisors of x (A001065).
%H Amiram Eldar, <a href="/A057710/b057710.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from T. D. Noe)
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/RestrictedDivisorFunction.html">Restricted Divisor Function</a>.
%H Eric Weisstein's World of Mathematics, <a href="http://mathworld.wolfram.com/AliquotSequence.html">Aliquot sequence</a>.
%e 14 is a member of the sequence because s(22) = 14 and s(169) = 14 (and because no other integer x satisfies s(x) = 14).
%t len = max = 57; f[_List] := (s = Select[ Split[ Sort[ Table[ DivisorSigma[1, n] - n, {n, 1, max *= 2}]]], Length[#] == 2 & ][[All, 1]]; s [[1 ;; Min[len, Length[s]]]]); FixedPoint[f, {}] (* _Jean-François Alcover_, Oct 07 2011 *)
%Y Cf. A001065, A005114, A057709.
%K nonn
%O 1,1
%A _Jack Brennen_, Oct 24 2000
|
