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

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A318263 Expansion of Product_{k>=1} (1 + Lucas(k)*x^k). 3
1, 1, 3, 7, 11, 30, 62, 129, 235, 541, 1034, 2101, 4140, 8129, 15984, 31903, 60398, 117646, 228808, 433768, 836552, 1601282, 3031299, 5736396, 10899112, 20466182, 38556342, 72522116, 135662847, 253047629, 473785878, 878655661, 1634304062, 3033385668, 5608183925 (list; graph; refs; listen; history; text; internal format)
OFFSET

0,3

LINKS

Vaclav Kotesovec, Table of n, a(n) for n = 0..4500

FORMULA

From Vaclav Kotesovec, Aug 24 2018: (Start)

a(n) ~ c * A000032(n) * A000009(n) ~ c * phi^n * exp(Pi*sqrt(n/3)) / (4 * 3^(1/4) * n^(3/4)), where c = Product_{k>=1} ((1 + Lucas(k)/phi^k)/2) = 0.8503149035690839100210269103058319341315494385103929947491... and phi = A001622 = (1+sqrt(5))/2 is the golden ratio.

Equivalently, c = Product_{k>=1} (1 + (-1)^k/(2*phi^(2*k))),

c = 2/3 * QPochhammer[-1/2, -1/GoldenRatio^2]. (End)

MATHEMATICA

nmax = 40; CoefficientList[Series[Product[1+LucasL[k]*x^k, {k, 1, nmax}], {x, 0, nmax}], x]

nmax = 40; poly = ConstantArray[0, nmax + 1]; poly[[1]] = 1; poly[[2]] = 1; Do[Do[poly[[j + 1]] += LucasL[k]*poly[[j - k + 1]], {j, nmax, k, -1}]; , {k, 2, nmax}]; poly

CROSSREFS

Cf. A000032, A022629, A318248, A318264.

Sequence in context: A233517 A125879 A238673 * A213740 A267357 A152084

Adjacent sequences:  A318260 A318261 A318262 * A318264 A318265 A318266

KEYWORD

nonn

AUTHOR

Vaclav Kotesovec, Aug 22 2018

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 April 18 14:42 EDT 2019. Contains 322209 sequences. (Running on oeis4.)