The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
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!)
A272553 Numbers n whose sum of divisors equals the sum of divisors of 2n+1. 4
20, 464, 650, 2744, 3980, 5504, 5736, 5922, 7032, 8130, 10472, 18618, 24312, 27654, 38874, 39500, 43032, 45492, 56870, 64410, 71058, 79068, 85158, 89178, 92130, 97014, 109928, 117114, 118902, 127688, 130304, 175554, 180438, 187304, 188292, 208452, 224058, 244674, 249788, 269192, 294380, 305624, 347964 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Most terms are even; the first three odd ones are 1264545, 8770125, and 9346995, and these are the only odd terms among the first 10^7 numbers that include 135 terms.
For some n, 2n+1 is prime; for example, this is so for the first three terms, but this happens rarely with only 4 cases among the first 10^7 numbers.
All terms are abundant numbers (A005101): since sigma(x)>x for x>1, sigma(2n+1)>2n+1>2n for n>0, the defining formula, sigma(n)=sigma(2n+1), implies sigma(n)>2n, which proves that n is an abundant number.
Up to 6*10^9 there are 1151 terms, 46 of which are odd. All these odd terms are multiple of 3 and all are multiple of 5, except 1501989489 and 4242679749. The values n for which 2n+1 is a prime number are a subset of A088831, thus it is easy to verify that up to 10^13 there are only 4 such values (20, 464, 650, and 130304). - Giovanni Resta, May 03 2016
LINKS
FORMULA
A000203(n) = A000203(2n+1).
EXAMPLE
20 is a term as its sum of divisors, 42=1+2+4+5+10+20, is the same as the sum of divisors of 41=2*20+1; 41 has only two divisors 1 and 41.
MAPLE
select(t -> numtheory:-sigma(t) = numtheory:-sigma(2*t+1), [$1..10^6]); # Robert Israel, May 03 2016
MATHEMATICA
Select[Range@500000, DivisorSigma[1, #]==DivisorSigma[1, 2*#+1]&]
PROG
(PARI) for (n=1, 500000, (sigma(n)==sigma(2*n+1)) && print1(n ", "))
CROSSREFS
Cf. A000203 (sum of divisors), A074821 (similar sequence for the number of divisors), A005101 (abundant numbers, supersequence), A088831,
Sequence in context: A111158 A049382 A288034 * A268884 A324069 A065412
KEYWORD
nonn
AUTHOR
Waldemar Puszkarz, May 02 2016
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 May 14 01:40 EDT 2024. Contains 372528 sequences. (Running on oeis4.)