login
Numbers that are not the sum of the nontrivial factors (excluding 1 and n) of some natural number.
8

%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_