login
This site is supported by donations to The OEIS Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A180002 a(n) red and b(n) blue balls in an urn; draw 5 balls without replacement; Probability(5 red balls) = Probability(3 red and 2 blue balls). 3
4, 8, 24, 43, 179, 783, 1504, 6668, 29604, 56983, 253079, 1124043, 2163724, 9610208, 42683904, 82164403, 364934699, 1620864183, 3120083464, 13857908228, 61550154924, 118481007103, 526235577839, 2337285022803, 4499158186324, 19983094049528 (list; graph; refs; listen; history; text; internal format)
OFFSET

1,1

COMMENTS

This is equivalent to the Pell equation A(n)^2 - 10*B(n)^2 = -9 with a(n) = (A(n)+7)/2, b(n) = (B(n)+1)/2, and the 3 fundamental solutions (1,1), (9,3), (41,13), and the solution (19,6) for the unit form.

LINKS

G. C. Greubel, Table of n, a(n) for n = 1..1000

Index entries for linear recurrences with constant coefficients, signature (1,0,38,-38,0,-1,1).

FORMULA

G.f.: x*(4 +4*x +16*x^2 -133*x^3 -16*x^4 -4*x^5 +3*x^6)/((1-x)*(1 -38*x^3 +x^6)).

a(n+9) = 39*a(n+6) - 39*a(n+3) + a(n).

Let r = sqrt(10) then:

a(3*n+1) = (14 + (1+r)*(19+6*r)^n + (1-r)*(19-6*r)^n)/4.

a(3*n+2) = (14 + 3*(3+r)*(19+6*r)^n + 3*(3-r)*(19-6*r)^n)/4.

a(3*n+3) = (14 + (41+13*r)*(19+6*r)^n + (41-13*r)*(19-6*r)^n)/4.

a(n) = a(n-1) + 38*a(n-3) - 38*a(n-4) - a(n-6) + a(n-7). - G. C. Greubel, Mar 20 2019

EXAMPLE

For n=3: a(3)=24 and b(3)=7 since binomial(24,5) = binomial(24,3)*binomial(7,2) = 42504.

MATHEMATICA

Rest[CoefficientList[Series[x*(4+4*x+16*x^2-133*x^3-16*x^4-4*x^5 +3*x^6 )/((1-x)*(1-38*x^3+x^6)), {x, 0, 30}], x]] (* G. C. Greubel, Mar 20 2019 *)

PROG

(PARI) my(x='x+O('x^30)); Vec(x*(4+4*x+16*x^2-133*x^3-16*x^4-4*x^5+3*x^6) /((1-x)*(1-38*x^3+x^6))) \\ G. C. Greubel, Mar 20 2019

(MAGMA) R<x>:=PowerSeriesRing(Integers(), 30); Coefficients(R!( x*(4+4*x+ 16*x^2-133*x^3-16*x^4-4*x^5+3*x^6)/((1-x)*(1-38*x^3+x^6)) )); // G. C. Greubel, Mar 20 2019

(Sage) a=(x*(4+4*x+16*x^2-133*x^3-16*x^4-4*x^5+3*x^6)/((1-x)*(1-38*x^3 +x^6))).series(x, 30).coefficients(x, sparse=False); a[1:] # G. C. Greubel, Mar 20 2019

CROSSREFS

Cf. A180003 (b(n)).

Sequence in context: A006640 A212686 A316961 * A266821 A306484 A208901

Adjacent sequences:  A179999 A180000 A180001 * A180003 A180004 A180005

KEYWORD

nonn

AUTHOR

Paul Weisenhorn, Aug 05 2010

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.

License Agreements, Terms of Use, Privacy Policy. .

Last modified October 15 04:33 EDT 2019. Contains 328026 sequences. (Running on oeis4.)