OFFSET
1,1
COMMENTS
This is to A214408 as deficient numbers are to abundant numbers.
Differs from A097218, which does not contain 105, for example.
The deficient numbers which are *not* in the sequence are 2, 4, 8, 10, 16, 32, 44, 64, 128, 136, 152, 184, 256, 512, 752, 884, 1024, 2048, 2144, 2272, 2528, 4096, 8192, 8384, 12224, 16384, 17176, 18632, 18904, 32768, 32896, 33664, ... the union of powers of 2 and the terms of A060326. - M. F. Hasler, Jul 21 2012
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
EXAMPLE
7 is in the sequence because 7 is deficient, and its abundance is -6, and |-6| = 6 does not divide 7.
MAPLE
filter:= proc(n) local t;
t:= 2*n-numtheory:-sigma(n);
t > 0 and n mod t <> 0
end proc:
select(filter, [$1..200]); # Robert Israel, Nov 13 2019
MATHEMATICA
q[n_] := Module[{def = 2*n - DivisorSigma[1, n]}, def > 0 && !Divisible[n, def]]; Select[Range[120], q] (* Amiram Eldar, Apr 07 2024 *)
PROG
(PARI) is_A214547(n)={sigma(n)<2*n & n%(2*n-sigma(n))} \\ M. F. Hasler, Jul 21 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jonathan Vos Post, Jul 20 2012
EXTENSIONS
Given terms double-checked with the PARI script by M. F. Hasler, Jul 21 2012
STATUS
approved