This site is supported by donations to The OEIS Foundation.

Thanks to everyone who made a donation during our annual appeal!
To see the list of donors, or make a donation, see the OEIS Foundation home page.

 Hints (Greetings from The On-Line Encyclopedia of Integer Sequences!)
 A004128 a(n) = Sum_{k=1..n} floor(3n/3^k). 19
 0, 1, 2, 4, 5, 6, 8, 9, 10, 13, 14, 15, 17, 18, 19, 21, 22, 23, 26, 27, 28, 30, 31, 32, 34, 35, 36, 40, 41, 42, 44, 45, 46, 48, 49, 50, 53, 54, 55, 57, 58, 59, 61, 62, 63, 66, 67, 68, 70, 71, 72, 74, 75, 76, 80, 81, 82, 84, 85, 86, 88, 89, 90, 93, 94, 95, 97, 98, 99, 101, 102 (list; graph; refs; listen; history; text; internal format)
 OFFSET 0,3 COMMENTS 3-adic valuation of (3n)! - cf. A054861. Denominators of expansion of (1-x)^{-1/3} are 3^a(n). Numerators are in |A067622|. REFERENCES Gary W. Adamson, in "Beyond Measure, A Guided Tour Through Nature, Myth and Number", by Jay Kappraff, World Scientific, 2002, p. 356. LINKS T. D. Noe, Table of n, a(n) for n = 0..1000 FORMULA a(n) = n+[n/3]+[n/9]+[n/27]+... = n+a([n/3]) = n+A054861(n) = A054861(3n) = (3n-A053735(n))/2. - Henry Bottomley, May 01 2001 a(n) = sum{k>=0, floor(n/3^k)}. a(n)=sum{0<=k<=floor(log_3(n)), floor(n/3^k)}, n>=1. - Hieronymus Fischer, Aug 14 2007 Recurrence: a(n) = n+a(floor(n/3)); a(3n) = 3n+a(n); a(n*3^m) = 3*n*(3^m-1)/2 + a(n). - Hieronymus Fischer, Aug 14 2007 a(k*3^m) = k*(3^(m+1)-1)/2, 0<=k<3, m>=0. - Hieronymus Fischer, Aug 14 2007 Asymptotic behavior: a(n) = 3/2*n+O(log(n)), a(n+1)-a(n) = O(log(n)); this follows from the inequalities below. - Hieronymus Fischer, Aug 14 2007 a(n) <= (3n-1)/2; equality holds for powers of 3. - Hieronymus Fischer, Aug 14 2007 a(n) >= (3n-2)/2-floor(log_3(n)); equality holds for n=3^m-1, m>0. - Hieronymus Fischer, Aug 14 2007 lim inf (3n/2-a(n)) = 1/2, for n-->oo. - Hieronymus Fischer, Aug 14 2007 lim sup (3n/2-log_3(n)-a(n)) = 0, for n-->oo. - Hieronymus Fischer, Aug 14 2007 lim sup (a(n+1)-a(n)-log_3(n)) = 1, for n-->oo. - Hieronymus Fischer, Aug 14 2007 G.f.: sum{k>=0, x^(3^k)/(1-x^(3^k))}/(1-x). - Hieronymus Fischer, Aug 14 2007 a(n) = Sum_{k>=0} A030341(n,k)*A003462(k+1). - Philippe Deléham, Oct 21 2011 MAPLE A004128 := proc(n)     A054861(3*n) ; end proc: seq(A004128(n), n=0..100) ; # R. J. Mathar, Nov 04 2017 MATHEMATICA Table[Total[NestWhileList[Floor[#/3] &, n, # > 0 &]], {n, 0, 70}] (* Birkas Gyorgy, May 20 2012 *) A004128 = Log[3, CoefficientList[ Series[1/(1+x)^(1/3), {x, 0, 100}], x] // Denominator] (* Jean-François Alcover, Feb 19 2015 *) PROG (PARI) {a(n) = my(s, t=1); while(t<=n, s += n\t; t*=3); s}; /*Michael Somos, Feb 26 2004 */ (Sage) A004128 = lambda n: A004128(n//3) + n if n > 0 else 0 [A004128(n) for n in (0..70)]  # Peter Luschny, Nov 16 2012 (Haskell) a004128 n = a004128_list !! (n-1) a004128_list = scanl (+) 0 a051064_list -- Reinhard Zumkeller, May 23 2013 CROSSREFS Cf. A004117, A001511, A051064, A055457. A051064(n) = a(n+1) - a(n). - Alford Arnold, Jul 19 2000 Cf. A054861, A067080, A098844, A132027, A005187, A054899. Sequence in context: A267137 A095775 A035063 * A023717 A171599 A288174 Adjacent sequences:  A004125 A004126 A004127 * A004129 A004130 A004131 KEYWORD nonn AUTHOR EXTENSIONS Current definition suggested by Jason Earls, Jul 04 2001 STATUS approved

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

Last modified January 18 03:13 EST 2019. Contains 319260 sequences. (Running on oeis4.)