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!)
A094423 a(n) = A045873(n)^2. 3
1, 4, 1, 144, 361, 484, 19321, 28224, 128881, 2427364, 1745041, 26501904, 285643801, 64995844, 4675961161, 31354493184, 149793121, 741117817924, 3178942795681, 545370434064, 107989070784841, 292105630845604 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
The g.f. is an example of a rational function with nonnegative integer coefficients that is not N-rational.
LINKS
Ira M. Gessel, Rational Functions With Nonnegative Integer Coefficients, slides, 50th Séminaire Lotharingien de Combinatoire, 2003.
FORMULA
G.f.: x*(1+5*x)/(1+x-5*x^2-125*x^3).
a(n) = A250102(n)/16.
a(n) = (1/16)*( 2*5^n - (1+2*i)^(2*n) - (1-2*i)^(2*n) ) = (1/4)*( Im(1+2*i)^n )^2 = (1/4) * 5^n * sin(n*arctan(2))^2.
From G. C. Greubel, Jan 11 2024: (Start)
a(n) = (1/8)*5^n*(1 - ChebyshevU(n, -3/5) - (3/5)*ChebyshevU(n-1, -3/5)).
a(n) = (1/8)*( 5^n - (-1)^n*A066771(n) ).
E.g.f.: (1/8)*exp(-3*x)*(exp(8*x) - cos(4*x)). (End)
MATHEMATICA
LinearRecurrence[{-1, 5, 125}, {1, 4, 1}, 40] (* G. C. Greubel, Jan 11 2024 *)
PROG
(PARI) Vec((x+5*x^2)/(1+x-5*x^2-125*x^3) + O(x^30)) \\ Michel Marcus, Aug 28 2015
(Magma) I:=[1, 4, 1]; [n le 3 select I[n] else -Self(n-1) +5*Self(n-2) +125*Self(n-3): n in [1..41]]; // G. C. Greubel, Jan 11 2024
(SageMath)
@CachedFunction
def a(n): # a = A094423
if (n<4): return (0, 1, 4, 1)[n]
else: return -a(n-1) + 5*a(n-2) + 125*a(n-3)
[a(n) for n in range(1, 41)] # G. C. Greubel, Jan 11 2024
CROSSREFS
Sequence in context: A262405 A152841 A125083 * A262404 A299522 A300140
KEYWORD
nonn,easy
AUTHOR
Ralf Stephan, May 04 2004
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 29 00:26 EDT 2024. Contains 371264 sequences. (Running on oeis4.)