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!)
A346119 Expansion of the e.g.f. sqrt(2*x*exp(x) - 2*exp(x) + 3). 3
1, 0, 1, 2, 0, -16, -35, 342, 2779, -6424, -239382, -822460, 22393657, 278844084, -1553468891, -68399947042, -275025888900, 15302175612416, 243541868882077, -2463105309082902, -121649966081262521, -473088821582805820, 50905612811064360006, 945133249101683013812, -15321255878414345388335 (list; graph; refs; listen; history; text; internal format)
OFFSET
0,4
LINKS
FORMULA
E.g.f. y(x) satisfies y*y' = x*exp(x).
a(0)=1, a(n) = Sum_{k=1..floor(n/2)} (-1)^(k-1)*A006677(k)*A008306(n,k) for n > 0.
For all p prime, a(p) == -1 (mod p).
For n > 1, a(n) == 0 (mod (n-1)).
Conjecture: a(n) = 0 for only n = 1 and n = 4.
EXAMPLE
sqrt(2*x*exp(x)-2*exp(x)+3) = 1 + x^2/2! + 2*x^3/3! - 16*x^5/5! - 35*x^6/6! + 342*x^7/7! + 2779*x^8/8! - 6424*x^9/9! + ...
a(11) = Sum_{k=1..5} (-1)^(k-1)*A006677(k)*A008306(11,k) = -822460.
For k=1, (-1)^(1-1)*A006677(1)*A008306(11,1) == -1 (mod 11), because A006677(1) = 1 and A008306(11,1) = (11-1)!
For k>=2, (-1)^(k-1)*A006677(k)*A008306(11,k) == 0 (mod 11), because A008306(11,k) == 0 (mod 11), result a(11) == -1 (mod 11).
a(8) = Sum_{k=1..4} (-1)^(k-1)*A006677(k)*A008306(8,k) = 2779.
a(8) == 0 (mod (8-1)), because for k >= 1, A008306(8,k) == 0 (mod 7).
MAPLE
stirtr:= proc(p) proc(n) add(p(k)*Stirling2(n, k), k=0..n) end end: f:= n-> `if`(n=0, 1, (2*n-2)!/ (n-1)!/ 2^(n-1)): A006677:= stirtr(f): # Alois P. Heinz, 2008.
A008306 := proc(n, k): if k=1 then (n-1)! ; elif n<=2*k-1 then 0; else (n-1)*procname(n-1, k)+(n-1)*procname(n-2, k-1) ; end if; end proc:
a:= n-> add(((-1)^(k-1)*A006677(k)*A008306(n, k)), k=1..iquo(n, 2)):a(0):=1 ; seq(a(n), n=0..24);
# second program:
a := series(sqrt(2*x*exp(x)-2*exp(x)+3), x=0, 25):seq(n!*coeff(a, x, n), n=0..24);
MATHEMATICA
CoefficientList[Series[Sqrt(2*x*E^x-2*E^x+3), {x, 0, 24}], x] * Range[0, 24]!
PROG
(PARI) my(x='x+O('x^30)); Vec(serlaplace(sqrt(2*x*exp(x) - 2*exp(x) + 3))) \\ Michel Marcus, Jul 05 2021
CROSSREFS
Sequence in context: A009045 A060313 A015154 * A009266 A009246 A358629
KEYWORD
sign
AUTHOR
Mélika Tebni, Jul 05 2021
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 28 22:04 EDT 2024. Contains 371254 sequences. (Running on oeis4.)