login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A077854 Expansion of 1/((1-x)*(1-2*x)*(1+x^2)). 12
1, 3, 6, 12, 25, 51, 102, 204, 409, 819, 1638, 3276, 6553, 13107, 26214, 52428, 104857, 209715, 419430, 838860, 1677721, 3355443, 6710886, 13421772, 26843545, 53687091, 107374182, 214748364, 429496729, 858993459, 1717986918, 3435973836, 6871947673 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,2
COMMENTS
Partial sums of A007910. - Mircea Merca, Dec 27 2010
This is the decimal representation of the middle column of "Rule 54" elementary cellular automaton. - Karl V. Keller, Jr., Sep 26 2021
This same sequence (except that the offset is changed to 4) is 2^n with the final digit chopped off. - J. Lowell, May 11 2022
LINKS
FORMULA
a(n) = 3*a(n-1) - 3*a(n-2) + 3*a(n-3) - 2*a(n-4), with initial values a(0) = 1, a(1) = 3, a(2) = 6, a(3) = 12.
a(n) = (1/10)*(2^(n+4) + (-1)^floor(n/2) - 2*(-1)^floor((n+1)/2) - 5).
Row sums of A130306. - Gary W. Adamson, May 20 2007
a(n) = floor(2^(n+3)/5). - Gary Detlefs, Sep 06 2010
a(n) = round((2^(n+4)-5)/10) = floor((2^(n+3)-1)/5) = ceiling((2^(n+3)-4)/5) = round((2^(n+3)-2)/5); a(n) = a(n-4) + 3*2^(n-1), n > 3. - Mircea Merca, Dec 27 2010
a(n) = 2^(n+1) - 1 - a(n-2); a(n) = a(n-1)/2 for n == 2, 3 (mod 4); a(n) = (a(n-1)-1)/2 for n == 0, 1 (mod 4). - Arie Bos, Apr 06 2013
a(n) = floor(A000975(n+2)*3/5). - Armands Strazds, Oct 18 2014
a(n) = Sum_{k=1..n+3} floor(1 + sin(k*Pi/2 + 3*Pi/4))*2^(n-k+3). - Andres Cicuttin, Mar 28 2016
a(n) = (-15 + 3*2^(3+n) + 2^(1 + n - 4*floor((1+n)/4)) + 2^(2 + n - 4*floor((2+n)/4)))/15. - Andres Cicuttin, Mar 28 2016
a(n) = (16*2^n+(-1)^((2*n-1+(-1)^n)/4)-2*(-1)^((2*n+1-(-1)^n)/4)-5)/10. - Wesley Ivan Hurt, Apr 01 2016
EXAMPLE
The sequence in hexadecimal shows the pattern
1, 3, 6, c,
19, 33, 66, cc,
199, 333, 666, ccc,
1999, 3333, 6666, cccc,
19999, 33333, 66666, ccccc,
199999, 333333, 666666, cccccc,
1999999, 3333333, 6666666, ccccccc,
19999999, 33333333, 66666666, cccccccc,
... - Armands Strazds, Oct 09 2014
MAPLE
a := proc(n) option remember; if n=0 then RETURN(1); fi; if n=1 then RETURN(3); fi; if n=2 then RETURN(6); fi; if n=3 then RETURN(12); fi; 3*a(n-1)-3*a(n-2)+3*a(n-3)-2*a(n-4); end;
seq(iquo(2^n, 5), n=3..35); # Zerinvary Lajos, Apr 20 2008
MATHEMATICA
CoefficientList[Series[1/((1 - x) (1 - 2 x) (1 + x^2)), {x, 0, 32}], x] (* Michael De Vlieger, Mar 29 2016 *)
LinearRecurrence[{3, -3, 3, -2}, {1, 3, 6, 12}, 40] (* Harvey P. Dale, Feb 06 2019 *)
PROG
(Magma) [Round((2^(n+4)-5)/10): n in [0..40]]; // Vincenzo Librandi, Jun 25 2011
(PARI) a(n)=(16<<n)\10 \\ Charles R Greathouse IV, Sep 23 2012
(PARI) Vec(1/(1-3*x+3*x^2-3*x^3+2*x^4)+O(x^99)) \\ Derek Orr, Oct 26 2014
(Haskell)
import Data.Bits (xor)
a077854 n = a077854_list !! n
a077854_list = scanl1 xor $ tail a000975_list :: [Integer]
-- Reinhard Zumkeller, Jan 04 2013
(Python) print([2**(n+3)//5 for n in range(50)]) # Karl V. Keller, Jr., Sep 26 2021
CROSSREFS
Equals A007909(n+3) - [n congruent 2, 3 mod 4].
Cf. A130306, A043291 (subsequence); A000975, A007910, A133872, A259661 (binary).
Sequence in context: A068425 A329355 A136444 * A265700 A293313 A267539
KEYWORD
nonn,easy
AUTHOR
N. J. A. Sloane, Nov 17 2002
STATUS
approved

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

License Agreements, Terms of Use, Privacy Policy. .

Last modified March 28 18:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)