Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
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!)

Revision History for A127666

(Underlined text is an addition; strikethrough text is a deletion.)

Showing entries 1-10 | older changes
A127666 Odd infinitary abundant numbers.
(history; published version)
#23 by Peter Luschny at Fri Sep 09 04:19:39 EDT 2022
STATUS

reviewed

approved

#22 by Joerg Arndt at Fri Sep 09 04:03:35 EDT 2022
STATUS

proposed

reviewed

#21 by Amiram Eldar at Fri Sep 09 03:29:43 EDT 2022
STATUS

editing

proposed

#20 by Amiram Eldar at Fri Sep 09 02:51:53 EDT 2022
COMMENTS

The numbers of terms not exceeding 10^k, for k = 4, 5, ..., are 1, 77, 473, 5703, 53569, 561610, 5525461, 54979537, ... . Apparently, the asymptotic density of this sequence exists and equals 0.0005... . - Amiram Eldar, Sep 09 2022

STATUS

approved

editing

#19 by Sean A. Irvine at Sun Jun 09 18:43:50 EDT 2019
STATUS

reviewed

approved

#18 by Michel Marcus at Sun Jun 09 12:57:48 EDT 2019
STATUS

proposed

reviewed

#17 by Joerg Arndt at Sun Jun 09 11:56:12 EDT 2019
STATUS

editing

proposed

#16 by Joerg Arndt at Sun Jun 09 11:56:09 EDT 2019
MATHEMATICA

ExponentList[n_Integer, factors_List]:={#, IntegerExponent[n, # ]}&/@factors; InfinitaryDivisors[1]:={1}; InfinitaryDivisors[n_Integer?Positive]:=Module[ { factors=First/@FactorInteger[n], d=Divisors[n] }, d[[Flatten[Position[ Transpose[ Thread[Function[{f, g}, BitOr[f, g]==g][ #, Last[ # ]]]&/@ Transpose[Last/@ExponentList[ #, factors]&/@d]], _?(And@@#&), {1}]] ]] ] Null; properinfinitarydivisorsum[k_]:=Plus@@InfinitaryDivisors[k]-k; Select[Range[1, 50000, 2], properinfinitarydivisorsum[ # ]># &][ # ]># &] (* end of program *)

STATUS

proposed

editing

#15 by Jinyuan Wang at Sun Jun 09 07:52:20 EDT 2019
STATUS

editing

proposed

#14 by Jinyuan Wang at Sun Jun 09 07:51:13 EDT 2019
PROG

(PARI) A049417(n) = {my(b, f=factorint(n)); prod(k=1, #f[, 2], b = binary(f[k, 2]); prod(j=1, #b, if(b[j], 1+f[k, 1]^(2^(#b-j)), 1)))}

isok(k) = A049417(k)>2*k&&k%2==1; \\ Jinyuan Wang, Jun 09 2019

STATUS

reviewed

editing

Discussion
Sun Jun 09 07:52
Jinyuan Wang: could anyone finish some of my draft edits?

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 18 23:05 EDT 2024. Contains 375284 sequences. (Running on oeis4.)