[go: up one dir, main page]

login
A138836
Non-Mersenne numbers A001348.
4
1, 2, 4, 5, 6, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70
OFFSET
1,2
COMMENTS
Numbers that are not in A001348.
a(1) to a(2042) equals A133398, then a(2043)=2048 <> A133398(2043)=2047.
LINKS
Chai Wah Wu, Algorithms for complementary sequences, arXiv:2409.05844 [math.NT], 2024.
PROG
(Python)
from sympy import primepi, prime
def A138836(n): return n+(k:=int(primepi((n).bit_length())-1))+int(n+k+1>=1<<prime(k+1)) if n>1 else 1 # Chai Wah Wu, Sep 10 2024
KEYWORD
nonn
AUTHOR
Omar E. Pol, Apr 05 2008
STATUS
approved