|
%I #24 Nov 05 2019 06:51:11
%S 1,4,51,87,95,119,123,145,161,187,205,209,215,237,245,247,261,267,275,
%T 287,289,291,303,305,321,323,325,335,341,371,405,407,425,429,447,471,
%U 473,497,515,517,519,529,539,551,555,561,575,583,611,623,625,627,657
%N Numbers that are not the sum of the nontrivial factors (excluding 1 and n) of some natural number.
%H Amiram Eldar, <a href="/A048995/b048995.txt">Table of n, a(n) for n = 1..10000</a> (terms 1..1000 from Donovan Johnson)
%t a048995[n_Integer] := Module[{t = Table[i, {i, n}], a048050, k},
%t a048050[m_] := Plus @@ Take[Divisors[m], {2, -2}];
%t Do[
%t If[a048050[k] == 0 || a048050[k] > n, Null, t[[a048050[k]]] = 0],
%t {k, 2, n^2}];
%t Drop[DeleteDuplicates[t], {2}]
%t ]; a048995[660] (* _Michael De Vlieger_, Nov 30 2014 *)
%o (PARI) mx=8479; v=vector(mx); for(i=2, mx^2, ch=sigma(i)-i-1; if(ch<=mx, if(ch>0, v[ch]=1))); c=0; for(i=1, mx, if(v[i]==0, c++; write("b048995.txt", c " " i))) /* _Donovan Johnson_, Feb 11 2013 */
%Y Complement of A048050.
%Y Cf. A005114 (the same property with the sum of proper divisors) and A007369 (the same property with the sum of all divisors).
%K nonn,nice
%O 1,2
%A Bill Taylor (W.Taylor(AT)math.canterbury.ac.nz)
%E Corrected and extended by _Jud McCranie_
|
