cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-2 of 2 results.

A171861 Expansion of x*(1+x+x^2) / ( (x-1)*(x^3+x^2-1) ).

Original entry on oeis.org

1, 2, 4, 6, 9, 13, 18, 25, 34, 46, 62, 83, 111, 148, 197, 262, 348, 462, 613, 813, 1078, 1429, 1894, 2510, 3326, 4407, 5839, 7736, 10249, 13578, 17988, 23830, 31569, 41821, 55402, 73393, 97226, 128798, 170622, 226027, 299423, 396652, 525453, 696078, 922108
Offset: 1

Views

Author

Ed Pegg Jr, Oct 16 2010

Keywords

Comments

Number of wins in Penney's game if the two players start HHT and TTT and HHT beats TTT.
HHT beats TTT 70% of the time. - Geoffrey Critzer, Mar 01 2014

Examples

			a(n) enumerates length n+2 sequences on {H,T} that end in HHT but do not contain the contiguous subsequence TTT.
a(3)=4 because we have: TTHHT, THHHT, HTHHT, HHHHT.
a(4)=6 because we have: TTHHHT, THTHHT, THHHHT, HTTHHT, HTHHHT, HHHHHT. - _Geoffrey Critzer_, Mar 01 2014

Links

Crossrefs

Related sequences are A000045 (HHH beats HHT, HTT beats TTH), A006498 (HHH beats HTH), A023434 (HHH beats HTT), A000930 (HHH beats THT, HTH beats HHT), A000931 (HHH beats TTH), A077868 (HHT beats HTH), A002620 (HHT beats HTT), A000012 (HHT beats THH), A004277 (HHT beats THT), A070550 (HTH beats HHH), A000027 (HTH beats HTT), A097333 (HTH beats THH), A040000 (HTH beats TTH), A068921 (HTH beats TTT), A054405 (HTT beats HHH), A008619 (HTT beats HHT), A038718 (HTT beats THT), A128588 (HTT beats TTT).
Cf. A164315 (essentially the same sequence).

Programs

  • Maple
    A171861 := proc(n) option remember; if n <=4 then op(n,[1,2,4,6]); else procname(n-1)+procname(n-2)-procname(n-4) ; end if; end proc:
  • Mathematica
    nn=44;CoefficientList[Series[x(1+x+x^2)/(1-x-x^2+x^4),{x,0,nn}],x] (* Geoffrey Critzer, Mar 01 2014 *)
  • PARI
    a(n)=([0,1,0,0; 0,0,1,0; 0,0,0,1; -1,0,1,1]^(n-1)*[1;2;4;6])[1,1] \\ Charles R Greathouse IV, Oct 03 2016

Formula

a(n) = a(n-1) +a(n-2) -a(n-4) = A000931(n+10)-3 = A134816(n+6)-3 = A078027(n+12)-3.
a(n) = A164315(n-1). - Alois P. Heinz, Oct 12 2017

A023434 Dying rabbits: a(n) = a(n-1) + a(n-2) - a(n-4).

Original entry on oeis.org

0, 1, 1, 2, 3, 4, 6, 8, 11, 15, 20, 27, 36, 48, 64, 85, 113, 150, 199, 264, 350, 464, 615, 815, 1080, 1431, 1896, 2512, 3328, 4409, 5841, 7738, 10251, 13580, 17990, 23832, 31571, 41823, 55404, 73395, 97228, 128800, 170624, 226029, 299425, 396654, 525455
Offset: 0

Views

Author

N. J. A. Sloane

Keywords

Comments

Limit_{n->infinity} a(n)/a(n-1) = positive root of 1+x-x^3 (smallest Pisot-Vijayaraghavan number, A060006). - Gerald McGarvey, Sep 19 2004
a(n) is the number of distinct even run-types taken over nonempty subsets of [n+1]. The run-type of a set of positive integers is the sequence of lengths when the set is decomposed into maximal runs of consecutive integers and it is even if all its entries are even. For example, the set {2,3,5,6,9,10,11} has run-type (2,2,3) and a(6)=6 counts (2),(4),(6),(2,2),(2,4),(4,2). - David Callan, Jul 14 2006
Partial sums of the sequence obtained by deleting the first 2 terms of A000931. Example: 0+1+0+1+1 = 3 = a(4). - David Callan, Jul 14 2006
One less than the sequence obtained by deleting the first 7 terms of A000931. - Ira M. Gessel, May 02 2007
This sequence counts ordered partitions of (n-1) into parts less than or equal to 3, in which the order of 1's are unimportant. Alternately, the order of 2's and 3's are important (see example). - David Neil McGrath, Apr 26 2015
Interleaving of A289692 and A077855. - Bruce J. Nicholson, Apr 09 2018

Examples

			G.f. = x + x^2 + 2*x^3 + 3*x^4 + 4*x^5 + 6*x^6 + 8*x^7 + 11*x^8 + ...
a(7)=8, with (n-1)=6. The partially ordered partitions of 6 are (33),(321,312,132=one),(231,213,123=one),(3111,1311,1131,1113=one),(222),(2211,1122,1221,2112,1212,2121=one),(21111,12111,11211,11121,11112=one),(111111). - _David Neil McGrath_, Apr 26 2015

Links

Crossrefs

Cf. A000931, A054405, A060006, A077905.
Cf. A289692, A077855.

Programs

  • Magma
    [0,1] cat [ n le 4 select (n) else Self(n-1)+Self(n-2)-Self(n-4): n in [1..45] ]; // Vincenzo Librandi, Apr 27 2015
  • Maple
    f:= gfun:-rectoproc({a(n)=a(n-1)+a(n-2)-a(n-4),seq(a(i)=[0,1,1,2][i+1],i=0..3)},a(n),remember):
seq(f(i),i=0..100); # Robert Israel, May 04 2015
  • Mathematica
    a[ n_] := If[ n < 0, SeriesCoefficient[ -x^3 / (1 - x^2 - x^3 + x^4), {x, 0, -n}], SeriesCoefficient[ x / (1 - x - x^2 + x^4), {x, 0, n}]]; (* Michael Somos, Nov 29 2013 *)
LinearRecurrence[{1, 1, 0, -1}, {0, 1, 1, 2}, 50] (* Vincenzo Librandi, Apr 27 2015 *)
  • PARI
    {a(n) = polcoeff( if( n<0, -x^3 / (1 - x^2 - x^3 + x^4), x / (1 - x - x^2 + x^4)) + x * O(x^abs(n)), abs(n))}; /* Michael Somos, Nov 29 2013 */
  • PARI
    x='x+O('x^99); concat(0, Vec(x/((1-x)*(1-x^2-x^3)))) \\ Altug Alkan, Apr 09 2018

Formula

a(n) = A000931(n+7)-1.
a(0)=0, a(1)=1, a(2)=1 then for n>2 a(n)=ceiling(r*a(n-1)) where r is the positive root of x^3-x-1=0. - Benoit Cloitre, Jun 19 2004
G.f.: x/((1-x)*(1-x^2-x^3)). - Jon Perry, Jul 04 2004
For n>2 a(n) = floor(sqrt(a(n-3)*a(n-2) + a(n-2)*a(n-1) + a(n-1)*a(n-3))) + 1. - Gerald McGarvey, Sep 19 2004
a(n) = Sum_{k=1..floor((n+2)/3)} binomial(floor((n+2-k)/2),k). This formula counts even run-types by length. - David Callan, Jul 14 2006
a(n) = a(n-2) + a(n-3) + 1. - Mark Dols, Feb 01 2010
a(n) + a(n+1) = A054405(n). Partial sums is A054405. - Michael Somos, Dec 01 2013
a(-3-n) = -A077905(n) for all n in Z. - Michael Somos, Sep 25 2014
Showing 1-2 of 2 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.