The OEIS mourns the passing of Jim Simons and is grateful to the Simons Foundation for its support of research in many branches of science, including the OEIS.
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!)
A322226 a(n) = A322227(n) / (n*(n+1)/2), where A322227(n) = [x^(n-1)] Product_{k=1..n} (k + x - k*x^2), for n >= 1. 5

%I #14 Jan 05 2019 02:28:20

%S 1,1,-2,-14,36,736,-1819,-88227,163185,19160769,-15294993,-6611719113,

%T -4120456053,3331391605595,7909978546731,-2311176691053557,

%U -10230565685787877,2114443392338548619,14290732101459857168,-2468208492961535459212,-23089123066195736850322,3581650102772343613724618,43704098963536055443651815,-6326399008631824968253597825,-96742680150222152446019734205

%N a(n) = A322227(n) / (n*(n+1)/2), where A322227(n) = [x^(n-1)] Product_{k=1..n} (k + x - k*x^2), for n >= 1.

%H Paul D. Hanna, <a href="/A322226/b322226.txt">Table of n, a(n) for n = 1..301</a>

%e The irregular triangle A322225 formed from coefficients of x^k in Product_{m=1..n} (m + x - m*x^2), for n >= 0, k = 0..2*n, begins

%e 1;

%e 1, 1, -1;

%e 2, 3, -3, -3, 2;

%e 6, 11, -12, -21, 12, 11, -6;

%e 24, 50, -61, -140, 75, 140, -61, -50, 24;

%e 120, 274, -375, -1011, 540, 1475, -540, -1011, 375, 274, -120;

%e 720, 1764, -2696, -8085, 4479, 15456, -5005, -15456, 4479, 8085, -2696, -1764, 720;

%e 5040, 13068, -22148, -71639, 42140, 169266, -50932, -221389, 50932, 169266, -42140, -71639, 22148, 13068, -5040; ...

%e in which the central terms equal A322228.

%e RELATED SEQUENCES.

%e Note that the terms in the secondary diagonal A322227 in the above triangle

%e [1, 3, -12, -140, 540, 15456, -50932, -3176172, 7343325, 1053842295, ...]

%e may be divided by triangular numbers to obtain this sequence

%e [1, 1, -2, -14, 36, 736, -1819, -88227, 163185, 19160769, -15294993, ...].

%t a[n_] := SeriesCoefficient[Product[k + x - k x^2, {k, 1, n}], {x, 0, n - 1}]/(n (n + 1)/2);

%t Array[a, 25] (* _Jean-François Alcover_, Dec 29 2018 *)

%o (PARI) {A322225(n, k) = polcoeff( prod(m=1, n, m + x - m*x^2) +x*O(x^k), k)}

%o /* Print the irregular triangle */

%o for(n=0, 10, for(k=0, 2*n, print1( A322225(n, k), ", ")); print(""))

%o /* Print this sequence */

%o for(n=1, 30, print1( A322225(n, n-1)/(n*(n+1)/2), ", "))

%Y Cf. A322227, A322228.

%Y Cf. A322236 (variant), A322894 (variant).

%K sign

%O 1,3

%A _Paul D. Hanna_, Dec 15 2018

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 May 13 12:32 EDT 2024. Contains 372519 sequences. (Running on oeis4.)