Svoboda | Graniru | BBC Russia | Golosameriki | Facebook
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A181487 Numbers k such that Sum_{d|k, d<k, d not occurring before} d > k. 7

%I #13 Jan 06 2024 19:00:50

%S 12,18,20,30,42,48,56,66,70,72,78,80,84,88,90,102,104,108,114,120,138,

%T 150,162,174,180,186,192,196,200,210,220,222,246,252,258,260,264,270,

%U 272,280,282,288,294,300,304,308,312,318,320,330,336,340,354,364,366

%N Numbers k such that Sum_{d|k, d<k, d not occurring before} d > k.

%C This is the complement of the set S occurring in S-perfect numbers A118372.

%C From _Amiram Eldar_, Aug 11 2023: (Start)

%C Sometimes called S-abundant numbers, since they are analogous to abundant numbers (A005101) as S-perfect numbers (A118372) are analogous to perfect numbers (A000396).

%C De Koninck and Ivić conjectured that this sequence has an asymptotic density.

%C The numbers of terms not exceeding 10^k, for k = 2, 3, ..., are 15, 152, 1567, 15336, 154301, 1541445, 15392073, ... . Apparently, the asymptotic density of this sequence exists and equals 0.15... . (End)

%D Elena Deza, Perfect and Amicable Numbers, World Scientific, 2023, pp. 325-327.

%H Donovan Johnson, <a href="/A181487/b181487.txt">Table of n, a(n) for n = 1..10000</a>

%H Jean-Marie De Koninck and Aleksandar Ivić, <a href="https://www.emis.de/journals/PIMB/078/2.html">On a sum of divisors problem</a>, Publications de l'Institut Mathématique (Beograd), New Series, Vol. 64 (78) (1998), pp. 9-20.

%H Gérard Villemin, <a href="http://villemin.gerard.free.fr/aNombre/TYPDIVIS/ParfaitS.htm">Nombres S-PARFAITS ou Nombres de Granville</a>, NOMBRES - Curiosités, théorie et usages, 2019 (in French).

%H Wikipedia, <a href="https://en.wikipedia.org/wiki/Granville_number">Granville number</a>.

%t seq[kmax_] := Module[{s = {1}, a = {}, sum}, Do[sum = Total[Select[Divisors[k], MemberQ[s, #] &]]; If[sum <= k, AppendTo[s, k], AppendTo[a, k]], {k, 2, kmax}]; a]; seq[400] (* _Amiram Eldar_, Aug 11 2023 *)

%o (PARI) A181487(Nmax) = { my(C=0); for(n=2,Nmax, sumdiv(n,d,!bittest(C,d)*d)>2*n & !print1(n", ") & C+=1<<n )}

%Y Cf. A000396, A005101, A118372.

%K nonn

%O 1,1

%A _M. F. Hasler_, Oct 28 2010

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 22 00:18 EDT 2024. Contains 375353 sequences. (Running on oeis4.)