Logo

Greetings from The On-Line Encyclopedia of Integer Sequences!

Hints

Search: id:A118932
Displaying 1-1 of 1 results found. page 1
     Format: long | short | internal | text      Sort: relevance | references | number      Highlight: on | off
A118932 E.g.f.: A(x) = exp( Sum_{n>=0} x^(3^n)/3^((3^n-1)/2) ). +0
5
1, 1, 1, 3, 9, 21, 81, 351, 1233, 10249, 75841, 388411, 3733401, 33702813, 215375889, 1984583511, 19181083041, 141963117201, 1797976123393, 22534941675379, 202605151063081, 2992764505338021, 43182110678814801, 445326641624332623 (list; graph; listen)
OFFSET

0,4

COMMENT

Equals invariant column vector V that satisfies matrix product A118931*V = V, where A118931(n,k) = n!/[k!(n-3k)!*3^k] for n>=3*k>=0; thus a(n) = Sum_{k=0..[n/3]} A118931(n,k)*a(k), with a(0)=1.

FORMULA

a(n) = Sum_{k=0..[n/3]} n!/[k!*(n-3*k)!*3^k] * a(k), with a(0)=1.

EXAMPLE

E.g.f. A(x) = exp( x + x^3/3 + x^9/3^4 + x^27/3^13 + x^81/3^40 +...)

= 1 + 1*x + 1*x^2/2! + 3*x^3/3! + 9*x^4/4! + 21*x^5/5!+ 81*x^6/6!+...

PROGRAM

(PARI) {a(n)=if(n==0, 1, sum(k=0, n\3, n!/(k!*(n-3*k)!*3^k)*a(k)))} (PARI) /* Defined by E.G.F.: */ {a(n)=n!*polcoeff( exp(sum(k=0, ceil(log(n+1)/log(3)), x^(3^k)/3^((3^k-1)/2))+x*O(x^n)), n, x)}

CROSSREFS

Cf. A118931; variants: A118930, A118935.

Adjacent sequences: A118929 A118930 A118931 this_sequence A118933 A118934 A118935

Sequence in context: A073947 A062811 A001470 this_sequence A053499 A074000 A067645

KEYWORD

nonn

AUTHOR

Paul D. Hanna (pauldhanna(AT)juno.com), May 06 2006

page 1

Search completed in 0.002 seconds

Lookup | Welcome | Find friends | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Transforms | Puzzles | Hot | Classics
More pages | Superseeker | Maintained by N. J. A. Sloane (njas@research.att.com)

Last modified October 12 15:26 EDT 2008. Contains 144830 sequences.


AT&T Labs Research