A241602 Erroneous version of A000287.
0, 0, 0, 0, 0, 1, 0, 4, 6, 24, 214
Offset: 1
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.
A(x;t) = t^3*x^3 + (4*t^4 + 3*t^5)*x^4 + (3*t^4 + 24*t^5 + 33*t^6 + 13*t^7)*x^5 + ... Triangle starts: n\k [1] [2] [3] [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [1] 0; [2] 0, 0; [3] 0, 0, 1; [4] 0, 0, 0, 4, 3; [5] 0, 0, 0, 3, 24, 33, 13; [6] 0, 0, 0, 0, 33, 188, 338, 252, 68; [7] 0, 0, 0, 0, 13, 338, 1705, 3580, 3740, 1938, 399; [8] 0, 0, 0, 0, 0, 252, 3580, 16980, 39525, 51300, 38076, 15180, 2530; [9] ...
T(n,k) = { if (n < 3 || k < 3, return(0)); sum(i=0, k-1, sum(j=0, n-1, (-1)^((i+j+1)%2) * binomial(i+j, i)*(i+j+1)*(i+j+2)/2* (binomial(2*n, k-i-1) * binomial(2*k, n-j-1) - 4 * binomial(2*n-1, k-i-2) * binomial(2*k-1, n-j-2)))); }; N=10; concat(concat([0,0,0], apply(n->vector(2*n-3, k, T(n,k)), [3..N]))) \\ test 1: N=100; y=x*Ser(vector(N, n, sum(i=1+(n+2)\3, (2*n)\3-1, T(i,n-i)))); 0 == x*(x+1)^2*(x+2)*(4*x-1)*y' + 2*(x^2-11*x+1)*(x+1)^2*y + 10*x^6 /* \\ test 2: x='x; t='t; N=44; y=Ser(apply(n->Polrev(vector(2*n-3, k, T(n, k)), 't), [3..N+2]), 'x) * t*x^3; 0 == (t + 1)^3*(x + 1)^3*(t + x + t*x)^3*y^4 + t*(t + 1)^2*x*(x + 1)^2*((4*t^4 + 12*t^3 + 12*t^2 + 4*t)*x^4 + (12*t^4 + 16*t^3 - 4*t^2 - 8*t)*x^3 + (12*t^4 - 4*t^3 - 49*t^2 - 30*t + 3)*x^2 + (4*t^4 - 8*t^3 - 30*t^2 - 21*t)*x + 3*t^2)*y^3 + t^2*(t + 1)*x^2*(x + 1)*((6*t^5 + 18*t^4 + 18*t^3 + 6*t^2)*x^5 + (18*t^5 + 12*t^4 - 30*t^3 - 24*t^2)*x^4 + (18*t^5 - 30*t^4 - 123*t^3 - 58*t^2 + 17*t)*x^3 + (6*t^5 - 24*t^4 - 58*t^3 + 25*t^2 + 56*t)*x^2 + (17*t^3 + 56*t^2 + 48*t + 3)*x + 3*t)*y^2 + t^3*x^3*((4*t^6 + 12*t^5 + 12*t^4 + 4*t^3)*x^6 + (12*t^6 - 36*t^4 - 24*t^3)*x^5 + (12*t^6 - 36*t^5 - 99*t^4 - 26*t^3 + 25*t^2)*x^4 + (4*t^6 - 24*t^5 - 26*t^4 + 81*t^3 + 80*t^2)*x^3 + (25*t^4 + 80*t^3 + 44*t^2 - 14*t)*x^2 + (-14*t^2 - 17*t)*x + 1)*y + t^6*x^6*((t^4 + 2*t^3 + t^2)*x^4 + (2*t^4 - 7*t^3 - 9*t^2)*x^3 + (t^4 - 9*t^3 + 11*t)*x^2 + (11*t^2 + 13*t)*x - 1) */
Q(n,k) = { \\ c-nets with n-edges, k-vertices (see A290326) if (k < 2+(n+2)\3 || k > 2*n\3, return(0)); sum(i=2, k, sum(j=k, n, (-1)^((i+j+1-k)%2)*binomial(i+j-k,i)*i*(i-1)/2* (binomial(2*n-2*k+2,k-i)*binomial(2*k-2, n-j) - 4*binomial(2*n-2*k+1, k-i-1)*binomial(2*k-3, n-j-1)))); }; a(n)={sum(k=2+(n+2)\3, 2*n\3, k!*Q(n,k))/(4*n)} \\ Andrew Howroyd, May 05 2021
max = 24; Clear[a, eq, s]; gf = Sum[a[k]*x^k, {k, 0, max}]; a[0] = 0; a[1] = 1; a[2] = 2; coes = CoefficientList[(x^4 - 2*x^3 + x^2)*gf^5 + (8*x^4 - 14*x^3 + 8*x^2 - 2*x)*gf^4 + (25*x^4 - 16*x^3 - 14*x^2 + 8*x + 1)*gf^3 + (38*x^4 + 15*x^3 - 30*x^2 - x + 2)*gf^2 + (28*x^4 + 36*x^3 - 5*x^2 - 12*x + 1)*gf + 8*x^4 + 17*x^3 + 8*x^2 - x, x]; eq[n_] := eq[n] = If[n == 1, Thread[Drop[coes, 3] == 0], eq[n-1] /. s[n-1] // First]; s[n_] := s[n] = (Print["n = ", n]; Solve[eq[n][[n]], a[n+2]]); sol = Table[s[n], {n, 1, max-2}] // Flatten; Table[a[n], {n, 1, max}] /. sol (* Jean-François Alcover, Apr 15 2014 *)
R:=proc(n) option remember; if n=1 then -1 elif n=2 then 2 else (1/(2*n))*((7*n-22)*R(n-1)+2*(2*n-1)*R(n-2)); fi; end; [seq(R(n),n=1..40)];
a[n_] := a[n] = Which[n == 1, -1, n == 2, 2, True, (1/(2*n))*((7*n-22)*a[n-1]+2*(2*n-1)*a[n-2])]; Table[a[n], {n, 1, 40}] (* Jean-François Alcover, Mar 06 2014, after Maple *)
Comments