cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Previous Showing 11-17 of 17 results.

A354653 G.f. A(x) satisfies: -3 = Sum_{n=-oo..oo} (-1)^n * x^(n*(n+1)/2) * A(x)^(n*(n-1)/2).

Original entry on oeis.org

4, 63, 3024, 188688, 13492350, 1044853344, 85281392688, 7224776707896, 629288553814092, 56002675660109424, 5070000855941708292, 465454828626459320736, 43230859988456631732954, 4054827527508982869148392, 383529048423080768494135488, 36541031890621600233033859488
Offset: 0

Views

Author

Paul D. Hanna, Jun 02 2022

Keywords

Examples

			G.f.: A(x) = 4 + 63*x + 3024*x^2 + 188688*x^3 + 13492350*x^4 + 1044853344*x^5 + 85281392688*x^6 + 7224776707896*x^7 + 629288553814092*x^8 + ...
such that A = A(x) satisfies:
(1) -3 = ... + x^36*A^28 - x^28*A^21 + x^21*A^15 - x^15*A^10 + x^10*A^6 - x^6*A^3 + x^3*A - x + 1 - A + x*A^3 - x^3*A^6 + x^6*A^10 - x^10*A^15 + x^15*A^21 - x^21*A^28 + x^28*A^36 + ...
(2) -3 = (1-x) - (1-x^3)*A + x*(1-x^5)*A^3 - x^3*(1-x^7)*A^6 + x^6*(1-x^9)*A^10 - x^10*(1-x^11)*A^15 + x^15*(1-x^13)*A^21 - x^21*(1-x^15)*A^28 + ...
(3) -3 = (1-A) - (1-A^3)*x + A*(1-A^5)*x^3 - A^3*(1-A^7)*x^6 + A^6*(1-A^9)*x^10 - A^10*(1-A^11)*x^15 + A^15*(1-A^13)*x^21 - A^21*(1-A^15)*x^28 + ...
(4) -3 = (1 - x*A)*(1 - A)*(1-x) * (1 - x^2*A^2)*(1 - x*A^2)*(1 - x^2*A) * (1 - x^3*A^3)*(1 - x^2*A^3)*(1 - x^3*A^2) * (1 - x^4*A^4)*(1 - x^3*A^4)*(1 - x^4*A^3) * (1 - x^5*A^5)*(1 - x^4*A^5)*(1 - x^5*A^4) * ...
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[4]); for(i=1,n, A = concat(A,0);
    A[#A] = polcoeff(3 + sum(m=0,sqrtint(2*#A+9), (-1)^m * x^(m*(m-1)/2) * (1 - x^(2*m+1)) * Ser(A)^(m*(m+1)/2) ),#A-1) );A[n+1]}
    for(n=0,30,print1(a(n),", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) -3 = Sum_{n=-oo..oo} (-1)^n * x^(n*(n-1)/2) * A(x)^(n*(n+1)/2).
(2) -3 = Sum_{n>=0} (-1)^n * x^(n*(n-1)/2) * (1 - x^(2*n+1)) * A(x)^(n*(n+1)/2).
(3) -3 = Sum_{n>=0} (-1)^n * A(x)^(n*(n-1)/2) * (1 - A(x)^(2*n+1)) * x^(n*(n+1)/2).
(4) -3 = Product_{n>=1} (1 - x^n*A(x)^n) * (1 - x^(n-1)*A(x)^n) * (1 - x^n*A(x)^(n-1)), by the Jacobi triple product identity.
a(n) = (-1)^(n+1) * Sum_{k=0..2*n+1} A354649(n,k)*(-3)^k, for n >= 0.
a(n) = Sum_{k=0..2*n+1} A354650(n,k)*3^k, for n >= 0.

A354654 G.f. A(x) satisfies: -4 = Sum_{n=-oo..oo} (-1)^n * x^(n*(n+1)/2) * A(x)^(n*(n-1)/2).

Original entry on oeis.org

5, 124, 9300, 912520, 102616748, 12498655200, 1604505393140, 213790010204692, 29287693334340840, 4099332312599011100, 583685111605968443456, 84277588096627459702860, 12310921909740521584887824, 1816058097888803062860159620, 270156262107594683175523302780
Offset: 0

Views

Author

Paul D. Hanna, Jun 02 2022

Keywords

Examples

			G.f.: A(x) = 5 + 124*x + 9300*x^2 + 912520*x^3 + 102616748*x^4 + 12498655200*x^5 + 1604505393140*x^6 + 213790010204692*x^7 + 29287693334340840*x^8 + ...
such that A = A(x) satisfies:
(1) -4 = ... + x^36*A^28 - x^28*A^21 + x^21*A^15 - x^15*A^10 + x^10*A^6 - x^6*A^3 + x^3*A - x + 1 - A + x*A^3 - x^3*A^6 + x^6*A^10 - x^10*A^15 + x^15*A^21 - x^21*A^28 + x^28*A^36 + ...
(2) -4 = (1-x) - (1-x^3)*A + x*(1-x^5)*A^3 - x^3*(1-x^7)*A^6 + x^6*(1-x^9)*A^10 - x^10*(1-x^11)*A^15 + x^15*(1-x^13)*A^21 - x^21*(1-x^15)*A^28 + ...
(3) -4 = (1-A) - (1-A^3)*x + A*(1-A^5)*x^3 - A^3*(1-A^7)*x^6 + A^6*(1-A^9)*x^10 - A^10*(1-A^11)*x^15 + A^15*(1-A^13)*x^21 - A^21*(1-A^15)*x^28 + ...
(4) -4 = (1 - x*A)*(1 - A)*(1-x) * (1 - x^2*A^2)*(1 - x*A^2)*(1 - x^2*A) * (1 - x^3*A^3)*(1 - x^2*A^3)*(1 - x^3*A^2) * (1 - x^4*A^4)*(1 - x^3*A^4)*(1 - x^4*A^3) * (1 - x^5*A^5)*(1 - x^4*A^5)*(1 - x^5*A^4) * ...
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[5]); for(i=1,n, A = concat(A,0);
    A[#A] = polcoeff(4 + sum(m=0,sqrtint(2*#A+9), (-1)^m * x^(m*(m-1)/2) * (1 - x^(2*m+1)) * Ser(A)^(m*(m+1)/2) ),#A-1) );A[n+1]}
    for(n=0,30,print1(a(n),", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) -4 = Sum_{n=-oo..oo} (-1)^n * x^(n*(n-1)/2) * A(x)^(n*(n+1)/2).
(2) -4 = Sum_{n>=0} (-1)^n * x^(n*(n-1)/2) * (1 - x^(2*n+1)) * A(x)^(n*(n+1)/2).
(3) -4 = Sum_{n>=0} (-1)^n * A(x)^(n*(n-1)/2) * (1 - A(x)^(2*n+1)) * x^(n*(n+1)/2).
(4) -4 = Product_{n>=1} (1 - x^n*A(x)^n) * (1 - x^(n-1)*A(x)^n) * (1 - x^n*A(x)^(n-1)), by the Jacobi Triple Product identity.
a(n) = (-1)^(n+1) * Sum_{k=0..2*n+1} A354649(n,k)*(-4)^k, for n >= 0.
a(n) = Sum_{k=0..2*n+1} A354650(n,k)*4^k, for n >= 0.

A354661 G.f. A(x) satisfies: 1 = Sum_{n=-oo..oo} (-x)^(n*(n+1)/2) * A(x)^(n*(n-1)/2), with A(0) = 0.

Original entry on oeis.org

1, 0, 0, 2, 0, 0, 8, 0, 0, 44, 0, 6, 280, 0, 96, 1934, 0, 1124, 14088, 18, 11792, 106536, 648, 117626, 828360, 13416, 1142288, 6580780, 216000, 10921088, 53184864, 3019614, 103408416, 435930008, 38629656, 973041448, 3615741192, 465419760, 9118011128, 30298375236
Offset: 1

Views

Author

Paul D. Hanna, Jun 02 2022

Keywords

Examples

			G.f.: A(x) = x + 2*x^4 + 8*x^7 + 44*x^10 + 6*x^12 + 280*x^13 + 96*x^15 + 1934*x^16 + 1124*x^18 + 14088*x^19 + 18*x^20 + 11792*x^21 + ...
such that A = A(x) satisfies:
(1) 1 = ... + x^36*A^28 + x^28*A^21 - x^21*A^15 - x^15*A^10 + x^10*A^6 + x^6*A^3 - x^3*A - x + 1 + A - x*A^3 - x^3*A^6 + x^6*A^10 + x^10*A^15 - x^15*A^21 - x^21*A^28 + x^28*A^36 +--+ ...
(2) 1 = (1-x) + (1-x^3)*A - x*(1-x^5)*A^3 - x^3*(1-x^7)*A^6 + x^6*(1-x^9)*A^10 + x^10*(1-x^11)*A^15 - x^15*(1-x^13)*A^21 - x^21*(1-x^15)*A^28 + ...
(3) 1 = (1+A) - (1+A^3)*x - A*(1+A^5)*x^3 + A^3*(1+A^7)*x^6 + A^6*(1+A^9)*x^10 - A^10*(1+A^11)*x^15 - A^15*(1+A^13)*x^21 + A^21*(1+A^15)*x^28 + ...
(4) 1 = (1 + x*A)*(1 + A)*(1-x) * (1 - x^2*A^2)*(1 - x*A^2)*(1 + x^2*A) * (1 + x^3*A^3)*(1 + x^2*A^3)*(1 - x^3*A^2) * (1 - x^4*A^4)*(1 - x^3*A^4)*(1 + x^4*A^3) * (1 + x^5*A^5)*(1 + x^4*A^5)*(1 - x^5*A^4) * ...
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[0]); for(i=0,n, A = concat(A,0);
    A[#A] = -polcoeff(-1 + sum(m=0,sqrtint(2*#A+9), (-x)^(m*(m-1)/2) * (1 - x^(2*m+1)) * Ser(A)^(m*(m+1)/2) ),#A-1) );H=A;A[n+1]}
    for(n=1,50,print1(a(n),", "))

Formula

G.f. A(x) = Sum_{n>=1} a(n)*x^n satisfies:
(1) 1 = Sum_{n=-oo..oo} (-x)^(n*(n-1)/2) * A(x)^(n*(n+1)/2).
(2) 1 = Sum_{n>=0} (-x)^(n*(n-1)/2) * (1 - x^(2*n+1)) * A(x)^(n*(n+1)/2).
(3) 1 = Sum_{n>=0} (-1)^(n*(n+1)/2) * A(x)^(n*(n-1)/2) * (1 + A(x)^(2*n+1)) * x^(n*(n+1)/2).
(4) 1 = Product_{n>=1} (1 - (-x)^n*A(x)^n) * (1 + (-x)^(n-1)*A(x)^n) * (1 + (-x)^n*A(x)^(n-1)), by the Jacobi triple product identity.
(5) A(-A(-x)) = x.
a(n) = (-1)^n * Sum_{k=0..2*n+1} A354649(n,k), for n >= 0.
a(n) = -Sum_{k=0..2*n+1} A354650(n,k)*(-1)^k, for n >= 0.

A354663 G.f. A(x) satisfies: 3 = Sum_{n=-oo..oo} (-x)^(n*(n+1)/2) * A(x)^(n*(n-1)/2).

Original entry on oeis.org

2, 9, 108, 1848, 36306, 771768, 17280096, 401451192, 9587095686, 233892105912, 5804193409056, 146051807458320, 3717875447707254, 95571022734750600, 2477365983601721280, 64684289495622383472, 1699638032224106092368, 44909438746576707103608
Offset: 0

Views

Author

Paul D. Hanna, Jun 02 2022

Keywords

Examples

			G.f.: A(x) = 2 + 9*x + 108*x^2 + 1848*x^3 + 36306*x^4 + 771768*x^5 + 17280096*x^6 + 401451192*x^7 + 9587095686*x^8 + 233892105912*x^9 + ...
such that A = A(x) satisfies:
(1) 3 = ... + x^36*A^28 + x^28*A^21 - x^21*A^15 - x^15*A^10 + x^10*A^6 + x^6*A^3 - x^3*A - x + 1 + A - x*A^3 - x^3*A^6 + x^6*A^10 + x^10*A^15 - x^15*A^21 - x^21*A^28 + x^28*A^36 +--+ ...
(2) 3 = (1-x) + (1-x^3)*A - x*(1-x^5)*A^3 - x^3*(1-x^7)*A^6 + x^6*(1-x^9)*A^10 + x^10*(1-x^11)*A^15 - x^15*(1-x^13)*A^21 - x^21*(1-x^15)*A^28 + ...
(3) 3 = (1+A) - (1+A^3)*x - A*(1+A^5)*x^3 + A^3*(1+A^7)*x^6 + A^6*(1+A^9)*x^10 - A^10*(1+A^11)*x^15 - A^15*(1+A^13)*x^21 + A^21*(1+A^15)*x^28 + ...
(4) 3 = (1 + x*A)*(1 + A)*(1-x) * (1 - x^2*A^2)*(1 - x*A^2)*(1 + x^2*A) * (1 + x^3*A^3)*(1 + x^2*A^3)*(1 - x^3*A^2) * (1 - x^4*A^4)*(1 - x^3*A^4)*(1 + x^4*A^3) * (1 + x^5*A^5)*(1 + x^4*A^5)*(1 - x^5*A^4) * ...
		

Crossrefs

Programs

  • PARI
    {a(n) = my(A=[2]); for(i=1,n, A = concat(A,0);
    A[#A] = -polcoeff(-3 + sum(m=0,sqrtint(2*#A+9), (-x)^(m*(m-1)/2) * (1 - x^(2*m+1)) * Ser(A)^(m*(m+1)/2) ),#A-1) );H=A;A[n+1]}
    for(n=0,30,print1(a(n),", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) 3 = Sum_{n=-oo..oo} (-x)^(n*(n-1)/2) * A(x)^(n*(n+1)/2).
(2) 3 = Sum_{n>=0} (-x)^(n*(n-1)/2) * (1 - x^(2*n+1)) * A(x)^(n*(n+1)/2).
(3) 3 = Sum_{n>=0} (-1)^(n*(n+1)/2) * A(x)^(n*(n-1)/2) * (1 + A(x)^(2*n+1)) * x^(n*(n+1)/2).
(4) 3 = Product_{n>=1} (1 - (-x)^n*A(x)^n) * (1 + (-x)^(n-1)*A(x)^n) * (1 + (-x)^n*A(x)^(n-1)), by the Jacobi triple product identity.
a(n) = (-1)^n * Sum_{k=0..2*n+1} A354649(n,k)*3^k, for n >= 0.
a(n) = -Sum_{k=0..2*n+1} A354650(n,k)*(-3)^k, for n >= 0.

A354664 G.f. A(x) satisfies: 4 = Sum_{n=-oo..oo} (-x)^(n*(n+1)/2) * A(x)^(n*(n-1)/2).

Original entry on oeis.org

3, 28, 756, 28200, 1205228, 55731456, 2714642292, 137199520340, 7127794098792, 378292284479388, 20421818573265728, 1117886561607128940, 61904487399635790288, 3461693986652051482948, 195203095905903229325340, 11087371481682320212435332, 633751222047605882649272600
Offset: 0

Views

Author

Paul D. Hanna, Jun 02 2022

Keywords

Examples

			G.f.: A(x) = 3 + 28*x + 756*x^2 + 28200*x^3 + 1205228*x^4 + 55731456*x^5 + 2714642292*x^6 + 137199520340*x^7 + 7127794098792*x^8 + ...
such that A = A(x) satisfies:
(1) 4 = ... + x^36*A^28 + x^28*A^21 - x^21*A^15 - x^15*A^10 + x^10*A^6 + x^6*A^3 - x^3*A - x + 1 + A - x*A^3 - x^3*A^6 + x^6*A^10 + x^10*A^15 - x^15*A^21 - x^21*A^28 + x^28*A^36 +--+ ...
(2) 4 = (1-x) + (1-x^3)*A - x*(1-x^5)*A^3 - x^3*(1-x^7)*A^6 + x^6*(1-x^9)*A^10 + x^10*(1-x^11)*A^15 - x^15*(1-x^13)*A^21 - x^21*(1-x^15)*A^28 + ...
(3) 4 = (1+A) - (1+A^3)*x - A*(1+A^5)*x^3 + A^3*(1+A^7)*x^6 + A^6*(1+A^9)*x^10 - A^10*(1+A^11)*x^15 - A^15*(1+A^13)*x^21 + A^21*(1+A^15)*x^28 + ...
(4) 4 = (1 + x*A)*(1 + A)*(1-x) * (1 - x^2*A^2)*(1 - x*A^2)*(1 + x^2*A) * (1 + x^3*A^3)*(1 + x^2*A^3)*(1 - x^3*A^2) * (1 - x^4*A^4)*(1 - x^3*A^4)*(1 + x^4*A^3) * (1 + x^5*A^5)*(1 + x^4*A^5)*(1 - x^5*A^4) * ...
		

Crossrefs

Programs

  • Mathematica
    (* Calculation of constant d: *) 1/r /. FindRoot[{r*s * QPochhammer[1/r, -r*s] * QPochhammer[-1/s, -r*s] * QPochhammer[-r*s]/((-1 + r)*(1 + s)) == 4, -4*(Log[-r*s] - (1 + s)*QPolyGamma[0, 1, -r*s] + (1 + s) * QPolyGamma[0, -Log[-s]/Log[-r*s], -r*s]) / (s*Log[-r*s]) + 4*r*(1 + s) * Derivative[0, 1][QPochhammer][1/r, -r*s] / QPochhammer[1/r, -r*s] + r^2*s*QPochhammer[1/r, -r*s]*QPochhammer[-r*s] * Derivative[0, 1][QPochhammer][-1/s, -r*s]/(-1 + r) + 4*r*(1 + s)*Derivative[0, 1][QPochhammer][-r*s, -r*s] / QPochhammer[-r*s] == 0}, {r, 1/50}, {s, 2}, WorkingPrecision -> 70] (* Vaclav Kotesovec, Jan 19 2024 *)
  • PARI
    {a(n) = my(A=[3]); for(i=1,n, A = concat(A,0);
    A[#A] = -polcoeff(-4 + sum(m=0,sqrtint(2*#A+9), (-x)^(m*(m-1)/2) * (1 - x^(2*m+1)) * Ser(A)^(m*(m+1)/2) ),#A-1) );H=A;A[n+1]}
    for(n=0,30,print1(a(n),", "))

Formula

G.f. A(x) = Sum_{n>=0} a(n)*x^n satisfies:
(1) 4 = Sum_{n=-oo..oo} (-x)^(n*(n-1)/2) * A(x)^(n*(n+1)/2).
(2) 4 = Sum_{n>=0} (-x)^(n*(n-1)/2) * (1 - x^(2*n+1)) * A(x)^(n*(n+1)/2).
(3) 4 = Sum_{n>=0} (-1)^(n*(n+1)/2) * A(x)^(n*(n-1)/2) * (1 + A(x)^(2*n+1)) * x^(n*(n+1)/2).
(4) 4 = Product_{n>=1} (1 - (-x)^n*A(x)^n) * (1 + (-x)^(n-1)*A(x)^n) * (1 + (-x)^n*A(x)^(n-1)), by the Jacobi triple product identity.
a(n) = (-1)^n * Sum_{k=0..2*n+1} A354649(n,k)*4^k, for n >= 0.
a(n) = -Sum_{k=0..2*n+1} A354650(n,k)*(-4)^k, for n >= 0.
a(n) ~ c * d^n / n^(3/2), where d = 62.81220628370975097276726417958831026998790927499386157136003... and c = 0.71771306470564419436314253512374835316192083855385416486... - Vaclav Kotesovec, Jun 08 2022
Formula (4) can be rewritten as the functional equation QPochhammer(-x*y) * QPochhammer(1/x, -x*y)/(1 - 1/x) * QPochhammer(-1/y, -x*y)/(1 + 1/y) = 4. - Vaclav Kotesovec, Jan 19 2024

A356500 Coefficients T(n,k) of x^n*y^k in A(x,y) for n >= 0, k = 0..3*n+1, where A(x,y) satisfies: y = Sum_{n=-oo..+oo} (-x)^(n^2) * A(x,y)^((n-1)^2), as an irregular triangle read by rows.

Original entry on oeis.org

0, 1, 1, 0, 0, 0, 1, 0, 0, 0, 4, 0, 0, 0, 4, 0, 0, 6, 0, 0, 0, 28, 0, 0, 0, 22, 0, 3, 0, 0, 0, 84, 0, 0, 0, 219, 0, 0, 0, 140, 0, 0, 0, 0, 135, 0, 0, 0, 981, 0, 0, 0, 1807, 0, 0, 0, 969, 0, 0, 0, 120, 0, 0, 0, 2568, 0, 0, 0, 10764, 0, 0, 0, 15368, 0, 0, 0, 7084, 0, 0, 54, 0, 0, 0, 4284, 0, 0, 0, 38896, 0, 0, 0, 114240, 0, 0, 0, 133266, 0, 0, 0, 53820, 0, 9, 0, 0, 0, 4662, 0, 0, 0, 94390, 0, 0, 0, 525980, 0, 0, 0, 1187433, 0, 0, 0, 1171390, 0, 0, 0, 420732
Offset: 0

Views

Author

Paul D. Hanna, Aug 09 2022

Keywords

Comments

T(n, 3*n+1) = [x^n*y^(3*n+1)] A(x,y) = binomial(4*n, n)/(3*n + 1) = A002293(n) for n >= 0, where g.f. G(x) of A002293 satisfies: G(x) = 1 + x*G(x)^4.
T(4*n, 1) = A000716(n) for n >= 0 (nonzero terms in column 1).
T(4*n+3, 2) = [x^(4*n+3)*y^2] A(x,y) = 2 * A354655(n+1) for n >= 0, where A354655 equals column 2 of triangle A354650.
T(4*n+2, 3) = [x^(4*n+2)*y^3] A(x,y) = 4 * A354656(n+1) for n >= 0, where A354656 equals column 3 of triangle A354650.
T(2*n, 2*n+1) = A356504(n), for n >= 0.
T(2*n+1, 2*n) = A356505(n) for n >= 0.
T(3*n, n+1) = A356506(n) for n >= 0.
T(3*n+1, n) = A355365(n) where A355365 is the central terms of A355360 = A355360(2*n,n).
Sum_{k=0..3*n+1} T(n, k) = A354248(n) for n >= 0 (row sums).
Sum_{k=0..3*n+1} T(n, k) * 2^k = A356502(n) for n >= 0.
Sum_{k=0..3*n+1} T(n, k) * 3^k = A356503(n) for n >= 0.
Sum_{k=0..3*n+1} T(4*n+1-k, k) = A355872(n+1) for n >= 0 (nonzero antidiagonal sums).
SPECIFIC VALUES.
(V.1) 1 = A(x,y) at x = -exp(-Pi) and y = Pi^(1/4)/gamma(3/4).
(V.2) 1 = A(x,y) at x = -exp(-2*Pi) and y = Pi^(1/4)/gamma(3/4) * (6 + 4*sqrt(2))^(1/4)/2.
(V.3) 1 = A(x,y) at x = -exp(-3*Pi) and y = Pi^(1/4)/gamma(3/4) * (27 + 18*sqrt(3))^(1/4)/3.
(V.4) 1 = A(x,y) at x = -exp(-4*Pi) and y = Pi^(1/4)/gamma(3/4) * (8^(1/4) + 2)/4.
(V.5) 1 = A(x,y) at x = -exp(-sqrt(3)*Pi) and y = gamma(4/3)^(3/2)*3^(13/8)/(Pi*2^(2/3)).

Examples

			G.f.: A(x,y) = y + x*(1 + y^4) + x^2*(4*y^3 + 4*y^7) + x^3*(6*y^2 + 28*y^6 + 22*y^10) + x^4*(3*y + 84*y^5 + 219*y^9 + 140*y^13) + x^5*(135*y^4 + 981*y^8 + 1807*y^12 + 969*y^16) + x^6*(120*y^3 + 2568*y^7 + 10764*y^11 + 15368*y^15 + 7084*y^19) + x^7*(54*y^2 + 4284*y^6 + 38896*y^10 + 114240*y^14 + 133266*y^18 + 53820*y^22) + x^8*(9*y + 4662*y^5 + 94390*y^9 + 525980*y^13 + 1187433*y^17 + 1171390*y^21 + 420732*y^25) + x^9*(3250*y^4 + 160965*y^8 + 1670942*y^12 + 6640711*y^16 + 12167001*y^20 + 10399545*y^24 + 3362260*y^28) + ...
such that A = A(x,y) satisfies
y = ... + x^16*A^25 - x^9*A^16 + x^4*A^9 - x*A^4 + A - x + x^4*A - x^9*A^4 + x^16*A^9 - x^25*A^16 +- ... + (-x)^(n^2) * A(x,y)^((n-1)^2) + ...
This irregular triangle of coefficients T(n,k) of x^n*y^k in A(x,y) for n >= 0, k = 0..3*n+1, begins:
  n = 0: [0, 1];
  n = 1: [1, 0, 0, 0, 1];
  n = 2: [0, 0, 0, 4, 0, 0, 0, 4];
  n = 3: [0, 0, 6, 0, 0, 0, 28, 0, 0, 0, 22];
  n = 4: [0, 3, 0, 0, 0, 84, 0, 0, 0, 219, 0, 0, 0, 140];
  n = 5: [0, 0, 0, 0, 135, 0, 0, 0, 981, 0, 0, 0, 1807, 0, 0, 0, 969];
  n = 6: [0, 0, 0, 120, 0, 0, 0, 2568, 0, 0, 0, 10764, 0, 0, 0, 15368, 0, 0, 0, 7084];
  n = 7: [0, 0, 54, 0, 0, 0, 4284, 0, 0, 0, 38896, 0, 0, 0, 114240, 0, 0, 0, 133266, 0, 0, 0, 53820];
  n = 8: [0, 9, 0, 0, 0, 4662, 0, 0, 0, 94390, 0, 0, 0, 525980, 0, 0, 0, 1187433, 0, 0, 0, 1171390, 0, 0, 0, 420732];
  n = 9: [0, 0, 0, 0, 3250, 0, 0, 0, 160965, 0, 0, 0, 1670942, 0, 0, 0, 6640711, 0, 0, 0, 12167001, 0, 0, 0, 10399545, 0, 0, 0, 3362260];
  ...
Reading this triangle by nonzero antidiagonals [x^(4*n+1-k)*y^k] A(x,y) for n >= 0, k = 0..3*n+1, yields triangle A356501:
  [1, 1];
  [0, 3, 6, 4, 1];
  [0, 9, 54, 120, 135, 84, 28, 4];
  [0, 22, 294, 1360, 3250, 4662, 4284, 2568, 981, 219, 22];
  [0, 51, 1260, 10120, 41405, 103020, 170324, 196172, 160965, 94390, 38896, 10764, 1807, 140];
  [0, 108, 4590, 58380, 368145, 1404102, 3587696, 6515712, 8715465, 8763645, 6684744, 3863496, 1670942, 525980, 114240, 15368, 969];
  ...
		

Crossrefs

Programs

  • PARI
    {T(n,k) = my(A=[y],M); for(i=1,n, A = concat(A,0); M = ceil(sqrt(n+1));
    A[#A] = -polcoeff( sum(m=-M,M, (-x)^(m^2)*Ser(A)^((m-1)^2)), #A-1)); polcoeff(A[n+1],k,y) }
    for(n=0,9, for(k=0,3*n+1, print1(T(n,k),", "));print(""))

Formula

G.f. A(x,y) = Sum_{n>=0} Sum_{k=0..3*n+1} T(n,k) * x^n * y^k satisfies:
(1) y = Sum_{n=-oo..+oo} (-x)^(n^2) * A(x,y)^((n+1)^2).
(2) y = A(x,y) * Product_{n>=1} (1 - x^(2*n)*A(x,y)^(2*n)) * (1 - x^(2*n-1)*A(x,y)^(2*n+1)) * (1 - x^(2*n-1)*A(x,y)^(2*n-3)), by the Jacobi triple product identity.
(3) y = (-x) * Product_{n>=1} (1 - x^(2*n)*A(x,y)^(2*n)) * (1 - x^(2*n+1)*A(x,y)^(2*n-1)) * (1 - x^(2*n-3)*A(x,y)^(2*n-1)), by the Jacobi triple product identity.
(4) y = A(x, F(x,y)) where F(x,y) = Sum_{n=-oo..+oo} (-x)^(n^2) * y^((n-1)^2).
(5) 1 = A(x, theta_4(x)) where theta_4(x) = 1 + 2*Sum_{n>=1} (-1)^n * x^(n^2) is a Jacobi theta function.

A354657 a(n) = A354655(n)/3, for n >= 1.

Original entry on oeis.org

1, 9, 49, 210, 765, 2492, 7434, 20700, 54420, 136360, 327789, 760102, 1707342, 3728025, 7935525, 16507152, 33624045, 67186077, 131891825, 254710260, 484474753, 908538081, 1681364124, 3073166600, 5551851375, 9919925145, 17541289017, 30714092066, 53279031420
Offset: 1

Views

Author

Paul D. Hanna, Jun 02 2022

Keywords

Crossrefs

Programs

  • PARI
    {A354650(n,k) = my(A=[1+y]); for(i=1,n, A = concat(A,0);
    A[#A] = polcoeff(y + sum(m=0,sqrtint(2*#A+9), (-1)^m * x^(m*(m-1)/2) * (1 - x^(2*m+1)) * Ser(A)^(m*(m+1)/2) ),#A-1) );
    polcoeff(A[n+1],k,y)}
    for(n=1,30,print1(A354650(n,2)/3,", "))

Formula

a(n) = A354655(n)/3, for n >= 1.
a(n) = A354650(n,2)/3, for n >= 1.
a(n) = (-1)^(n+1) * A354649(n,2)/3, for n >= 1.
Previous Showing 11-17 of 17 results.