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!)
A254603 Nonprime numbers n such that sum of the divisors of n is a power of 2. 0
1, 21, 93, 217, 381, 651, 889, 2667, 3937, 11811, 24573, 27559, 57337, 82677, 172011, 253921, 393213, 761763, 917497, 1040257, 1572861, 1777447, 2752491, 3120771, 3670009, 4063201, 5332341, 7281799, 11010027, 12189603, 16252897, 16646017, 21845397, 28442407 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
a(1)=1; for n>=2, a(n) = composite numbers that are a product of distinct Mersenne primes (A046528).
Also nonprime numbers n such that A051027(n) = sigma(sigma(n)) = 2*sigma(n)-1 = 2^(k+1)-1 for some k. If n is composite number (product of distinct Mersenne primes) then k is the sum of Mersenne exponents (A000043) of these distinct Mersenne primes. Example: 651 = 3*7*31 = (2^2-1)*(2^3-1)*(2^5-1); k=2+3+5=10; A051027(651) = sigma(sigma(651)) = 2^(10+1)-1 = 2047.
Complement of A000668 (Mersenne primes) with respect to A046528.
LINKS
EXAMPLE
651 = 3*7*31 (product of three distinct Mersenne primes); sigma(651) = 1024 = 2^10.
PROG
(Magma) [n: n in [1..10^6] | not IsPrime(n) and SumOfDivisors(SumOfDivisors(n)) eq 2*SumOfDivisors(n) - 1]
(Magma)[n: n in[1..10000], k in [0..100] | not IsPrime(n) and SumOfDivisors(n) eq 2^k]
CROSSREFS
Sequence in context: A326164 A143843 A119109 * A144856 A065522 A218844
KEYWORD
nonn
AUTHOR
Jaroslav Krizek, Feb 02 2015
STATUS
approved

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 April 24 14:18 EDT 2024. Contains 371960 sequences. (Running on oeis4.)