OFFSET
1,2
COMMENTS
If 1 is excluded (of form [e^0]) then complement of A000149. - Michel Marcus, Jun 16 2013
REFERENCES
J. Roberts, Lure of the Integers, Math. Assoc. America, 1992, p. 11.
N. J. A. Sloane and Simon Plouffe, The Encyclopedia of Integer Sequences, Academic Press, 1995 (includes this sequence).
LINKS
Reinhard Zumkeller, Table of n, a(n) for n = 1..10000
J. Lambek and L. Moser, Inverse and complementary sequences of natural numbers, Amer. Math. Monthly, 61 (1954), 454-458.
FORMULA
a(n) = n + [ log (n + 1 + [ log (n + 1) ]) ].
a(n) = n+m if n+m>=exp(m) and a(n) = n+m-1 otherwise, where m = 1+floor(log(n)). - Chai Wah Wu, Nov 09 2025
MATHEMATICA
Table[n + Floor@ Log[n + 1 + Floor@ Log[n + 1]], {n, 50}] (* Michael De Vlieger, Oct 06 2017 *)
PROG
(PARI) a(n) = n + floor( log (n + 1 + floor( log (n + 1) )) ) \\ Michel Marcus, Jun 16 2013
(Haskell)
a003619 n = n + floor (log (x + fromIntegral (floor $ log x)))
where x = fromIntegral n + 1
-- Reinhard Zumkeller, Mar 17 2015
(Python)
from sympy import exp, log, floor
def A003619(n): return n+(m:=int(1+floor(log(n))))-bool(n+m<exp(m)) # Chai Wah Wu, Nov 09 2025
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved
