OFFSET
0,2
COMMENTS
First differs from A025709 at a(84).
Positions of zeros in A025658. - R. J. Mathar, Jul 06 2025
LINKS
Robert Israel, Table of n, a(n) for n = 0..10000
FORMULA
a(n) = n + 1 + Sum_{j = 1 .. n} floor(j*log_6(4)). - Robert Israel, May 10 2015
From Amiram Eldar, Nov 13 2025: (Start)
A025618(a(n)) = 4^n.
a(n) ~ c * n^2 / 2, where c = log_6(4) (A153102). (End)
MAPLE
N:= 200: # to get a(0) to a(N)
[$1..N+1] + ListTools:-PartialSums([seq(floor(j*log[6](4)), j=0..N)]); # Robert Israel, May 10 2015
MATHEMATICA
Accumulate[Table[Ceiling[n * Log[6, 4]], {n, 0, 60}]] + 1 (* Amiram Eldar, Nov 13 2025 *)
PROG
(PARI) a(n)=my(N=1); n+1+sum(i=1, n, logint(N*=4, 6)); \\ Charles R Greathouse IV, Jan 11 2018
(PARI) first(n)=my(s, N=1/4); vector(n+1, i, s+=logint(N*=4, 6)+1) \\ Charles R Greathouse IV, Jan 11 2018
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
EXTENSIONS
Edited and offset corrected by Robert Israel, May 10 2015
STATUS
approved
