OFFSET
1,2
COMMENTS
It is clear that n+1 must be a power of 2. Hence n=2^k-1 for some k. Found k=1, 2, 4, 6, 8, 16, 32. No other k<150. - T. D. Noe, Nov 04 2005
FORMULA
Note that this is different from the sequence A019434 - 2.
MATHEMATICA
Select[Range[10^6], 2^DivisorSigma[0, # ] == # + 1 &]
2^Select[Range[150], DivisorSigma[0, 2^#-1]==#&] - 1 (Noe)
PROG
(Python)
from sympy import divisor_count as tau
def afind(klimit, kstart=1):
for k in range(kstart, klimit+1):
m = 2**k - 1
if 2**tau(m) == m + 1: print(m, end=", ")
afind(klimit=100) # Michael S. Branicky, Dec 16 2021
CROSSREFS
KEYWORD
nonn,hard,more
AUTHOR
Joseph L. Pe, Nov 02 2005
EXTENSIONS
One more term from T. D. Noe, Nov 04 2005
STATUS
approved