A001834 -id:A001834 - cp's OEIS frontend

A098495 Array T(r,c) read by antidiagonals: values of triangle A098493 interpreted as polynomials, r >= 0.

1, 1, 0, 1, -1, -1, 1, -2, -1, -1, 1, -3, 1, 1, 0, 1, -4, 5, 1, 1, 1, 1, -5, 11, -7, -2, -1, 1, 1, -6, 19, -29, 9, 1, -1, 0, 1, -7, 29, -71, 76, -11, 1, 1, -1, 1, -8, 41, -139, 265, -199, 13, -2, 1, -1, 1, -9, 55, -239, 666, -989, 521, -15, 1, -1, 0, 1, -10, 71, -377, 1393, -3191, 3691, -1364, 17, 1, -1, 1, 1, -11, 89, -559

Offset: 0

Ralf Stephan, Sep 12 2004

			Array begins
  1,  0, -1,  -1,   0,    1,    1,   0, -1, ...
  1, -1, -1,   1,   1,   -1,   -1,   1,  1, ...
  1, -2,  1,   1,  -2,    1,    1,  -2,  1, ...
  1, -3,  5,  -7,   9,  -11,   13, -15, ...
  1, -4, 11, -29,  76, -199,  521, ...
  1, -5, 19, -71, 265, -989, 3691, ...
  ...

Alex Fink, Richard K. Guy, and Mark Krusemeyer, Partitions with parts occurring at most thrice, Contributions to Discrete Mathematics, Vol 3, No 2 (2008), pp. 76-114.

See A094954 (with negative k) for negative r and more formulas and programs.

Rows include (-1)^c times A005408, A002878, A001834, A030221, A002315. Columns include A028387. Antidiagonal sums are in A098496.

T[r_, 1] := 1; T[r_, 2] := -r - 1; T[r_, c_] := -r*T[r, c - 1] - T[r, c - 2]; Flatten[ Table[ T[n - i, i], {n, 0, 11}, {i, n + 1}]] (* Robert G. Wilson v, May 10 2005 *)

{ t(r,c)=if(c>r||c<0||r<0,0,if(c>=r-1,(-1)^r*if(c==r,1,-c),if(r==1,0,if(c==0,t(r-1,0)-t(r-2,0),t(r-1,c)-t(r-2,c)-t(r-1,c-1))))) }
T(r,c)=sum(i=0,c,t(c,i)*r^i);
matrix(5,5,n,k,T(n-1,k-1))

Recurrence: T(r, 1) = 1, T(r, 2) = -r-1, T(r, c) = -rT(r, c-1) - T(r, c-2). (Corrected Oct 19 2004)

More terms from Robert G. Wilson v, May 10 2005

A102206 a(0) = 3, a(1) = 8, a(n+2) = 4*a(n+1) - a(n) - 2.

Original entry on oeis.org

3, 8, 27, 98, 363, 1352, 5043, 18818, 70227, 262088, 978123, 3650402, 13623483, 50843528, 189750627, 708158978, 2642885283, 9863382152, 36810643323, 137379191138, 512706121227, 1913445293768, 7141075053843, 26650854921602, 99462344632563, 371198523608648

Offset: 0

Creighton Dement, Dec 30 2004

Colin Barker, Table of n, a(n) for n = 0..1000
Index entries for linear recurrences with constant coefficients, signature (5,-5,1).

Cf. A092184, A001834, A001353, A102207.

a[0] = 3; a[1] = 8; a[n_] := a[n] = 4a[n - 1] - a[n - 2] - 2; Table[a[n], {n, 0, 23}] (* Or *)
CoefficientList[ Series[(2x - 1)(x - 3)/((1 - x)(x^2 - 4x + 1)), {x, 0, 22}], x] (* Robert G. Wilson v, Jan 12 2005 *)
LinearRecurrence[{5,-5,1},{3,8,27},30] (* Harvey P. Dale, Jul 25 2012 *)

Vec((2*x-1)*(x-3)/((1-x)*(x^2-4*x+1)) + O(x^30)) \\ Colin Barker, Nov 03 2016

G.f.: (2x-1)(x-3)/((1-x)(x^2-4x+1)).
a(n) = A092184(n+1) + 2; a(n+1) - a(n) = A001834(n+1) (see comment).
a(0)=3, a(1)=8, a(2)=27, a(n) = 5*a(n-1) - 5*a(n-2) + a(n-3). - Harvey P. Dale, Jul 25 2012
a(n) = (2+(2-sqrt(3))^(1+n)+(2+sqrt(3))^(1+n))/2. - Colin Barker, Nov 03 2016

More terms from Robert G. Wilson v, Jan 12 2005

Recurrence in the definition corrected by R. J. Mathar, Aug 07 2008

A103772 Larger of two sides in a (k,k,k-1)-integer-sided triangle with integer area.

Original entry on oeis.org

1, 17, 241, 3361, 46817, 652081, 9082321, 126500417, 1761923521, 24540428881, 341804080817, 4760716702561, 66308229755041, 923554499868017, 12863454768397201, 179164812257692801, 2495443916839302017, 34757050023492535441, 484103256412056194161

Offset: 1

Zak Seidov, Feb 23 2005

Corresponding areas are 0, 120, 25080, 4890480, 949077360, 184120982760, ...
Values of (x^2 + y^2)/2, where the pair (x, y) satisfies x^2 - 3*y^2 = -2, i.e., a(n) = {(A001834(n))^2 + (A001835(n))^2}/2 = {(A001834(n))^2 + A046184(n)}/2. - Lekraj Beedassy, Jul 13 2006
The heights of these triangles are given in A028230. (A028230(n), A045899(n), A103772(n)) forms a primitive Pythagorean triple.
Shortest side of (k,k+2,k+3) triangle such that median to longest side is integral. Sequence of such medians is A028230. - James R. Buddenhagen, Nov 22 2013
Numbers n such that (n+1)*(3n-1) is a square. - James R. Buddenhagen, Nov 22 2013

Colin Barker, Table of n, a(n) for n = 1..850
J. B. Cosgrave, The Gauss-Factorial Motzkin connection (Maple worksheet, change suffix to .mw)
J. B. Cosgrave and K. Dilcher, An Introduction to Gauss Factorials, The American Mathematical Monthly, 118 (Nov. 2011), 812-829.
Project Euler, Problem 94: Almost Equilateral Triangles.
Index entries for linear recurrences with constant coefficients, signature (15,-15,1).

Cf. A102341, A103974, A016064, A011945, A028230.

I:=[1,17]; [n le 2 select I[n] else 14*Self(n-1)-Self(n-2)+4: n in [1..20]]; // Vincenzo Librandi, Mar 05 2016

a[1] = 1; a[2] = 17; a[3] = 241; a[n_] := a[n] = 15a[n - 1] - 15a[n - 2] + a[n - 3]; Table[ a[n] - 1, {n, 17}] (* Robert G. Wilson v, Mar 24 2005 *)
LinearRecurrence[{15,-15,1},{1,17,241},20] (* Harvey P. Dale, Jan 02 2016 *)
RecurrenceTable[{a[1] == 1, a[2] == 17, a[n] == 14 a[n-1] - a[n-2] + 4}, a, {n, 20}] (* Vincenzo Librandi, Mar 05 2016 *)

Vec(x*(1+x)^2/((1-x)*(1-14*x+x^2)) + O(x^25)) \\ Colin Barker, Mar 05 2016

a(n) = (4*A001570(n+1) - 1)/3, n > 0. - Ralf Stephan, May 20 2007
a(n) = A052530(n-1)*A052530(n) + 1. - Johannes Boot, May 21 2011
G.f.: x*(1+x)^2/((1-x)*(1-14*x+x^2)). - Colin Barker, Apr 09 2012
a(n) = 15*a(n-1) - 15*a(n-2) + a(n-3); a(1)=1, a(2)=17, a(3)=241. - Harvey P. Dale, Jan 02 2016
a(n) = (-1+(7-4*sqrt(3))^n*(2+sqrt(3))-(-2+sqrt(3))*(7+4*sqrt(3))^n)/3. - Colin Barker, Mar 05 2016
a(n) = 14*a(n-1) - a(n-2) + 4. - Vincenzo Librandi, Mar 05 2016
a(n) = A001353(n)^2 + A001353(n-1)^2. - Antonio Alberto Olivares, Apr 06 2020

More terms from Robert G. Wilson v, Mar 24 2005

A107903 Generalized NSW numbers.

Original entry on oeis.org

1, 10, 76, 568, 4240, 31648, 236224, 1763200, 13160704, 98232832, 733219840, 5472827392, 40849739776, 304906608640, 2275853910016, 16987204845568, 126794223124480, 946404965613568, 7064062832410624, 52726882796830720

Offset: 0

Paul Barry, May 27 2005

Counts total area under elevated Schroeder paths of length 2n+2, where horizontal steps can choose from three colors.
Case r=3 for family (1+(r-1)x)/(1-2(1+r)x+(1-r)^2*x^2). Case r=2 gives NSW numbers A002315 and case r=4 gives NSW numbers A096053.
Fifth binomial transform of (1+8x)/(1-16x^2), A107906.
If p is an odd prime, a((p-1)/2) == 1 mod p. - Altug Alkan, Mar 17 2016

Index entries for linear recurrences with constant coefficients, signature (8,-4).

Cf. A001834, A099156, A102591.

Table[Sum[Binomial[2 n + 1, 2 k] 3^k, {k, 0, n}], {n, 0, 20}] (* or *) CoefficientList[Series[(1 + 2 x)/(1 - 8 x + 4 x^2), {x, 0, 20}], x] (* Michael De Vlieger, Mar 17 2016 *)

Vec((1+2*x)/(1-8*x+4*x^2) + O(x^40)) \\ Michel Marcus, Mar 17 2016

G.f.: (1+2*x)/(1-8*x+4*x^2). [corrected by Ralf Stephan, Nov 30 2010]
a(n) = Sum_{k=0..n} binomial(2*n+1, 2*k)*3^k.
a(n) = ((1+sqrt(3))*(4+2*sqrt(3))^n+(1-sqrt(3))*(4-2*sqrt(3))^n)/2 = A099156(n+1)+2*A099156(n).
a(n) = 8*a(n-1) - 4*a(n-2); a(0) = 1, a(1) = 10. - Lekraj Beedassy, Apr 19 2020
a(n) = 2^n*A001834(n). - Philippe Deléham, Mar 18 2023

Typo corrected and link added by Johannes W. Meijer, Aug 07 2010

A108412 Expansion of (1 + x + x^2)/(1 - 4x^2 + x^4).

Original entry on oeis.org

1, 1, 5, 4, 19, 15, 71, 56, 265, 209, 989, 780, 3691, 2911, 13775, 10864, 51409, 40545, 191861, 151316, 716035, 564719, 2672279, 2107560, 9973081, 7865521, 37220045, 29354524, 138907099, 109552575, 518408351, 408855776, 1934726305

Offset: 0

Ralf Stephan, Jun 05 2005

This is the sequence of Lehmer numbers u_n(sqrt(R),Q) with the parameters R = 6 and Q = 1. It is a strong divisibility sequence, that is, gcd(a(n),a(m)) = a(gcd(n,m)) for all natural numbers n and m. The sequence satisfies a linear recurrence of order four. - Peter Bala, Apr 18 2014
The sequence of convergents of the 2-periodic continued fraction [0; 1, -6, 1, -6, ...] = 1/(1 - 1/(6 - 1/(1 - 1/(6 - ...)))) = 3 - sqrt(3) begins [0/1, 1/1, 6/5, 5/4, 24/19, 19/15, 90/71,...]. The present sequence is the sequence of denominators; the sequence of numerators of the continued fraction convergents [1, 6, 5, 24, 19, 90,...] is also a strong divisibility sequence. Cf. A005013 and A203976. - Peter Bala, May 19 2014
From Peter Bala, Mar 25 2018: (Start)
The following remarks assume an offset of 1.
Define a binary operation o on the real numbers by x o y = x*sqrt(1 + (1/2)*y^2) + y*sqrt(1 + (1/2)*x^2). The operation o is commutative and associative with identity 0. We have a(2*n + 1) = 1 o 1 o ... o 1 (2*n + 1 terms) and sqrt(6)*a(2*n) = (1 o 1 o ... o 1) (2*n terms). Cf. A005013 and A084068. For example, 1 o 1 = sqrt(6) and 1 o 1 o 1 = sqrt(6) o 1 = 5 = a(3).
From the obvious identity ( 1 o 1 o ... o 1 (2*n terms) ) o ( 1 o 1 o ... o 1 (2*m terms) ) = 1 o 1 o ... o 1 (2*n + 2*m terms) we find the relation a(2*n+2*m) = a(2*n)*sqrt(1 + 3*a(2*m)^2) + a(2*m)*sqrt(1 + 3*a(2*n)^2).
Similarly, from a(2*n+1) o a(2*m+1) = sqrt(6)*a(2*n+2*m+2) we find sqrt(6)*a(2*n+2*m+2) = a(2*n+1)*sqrt(1 + (1/2)*a(2*m+1)^2) + a(2*m+1)*sqrt(1 + (1/2)*a(2*n+1)^2). (End)

			G.f. = 1 + x + 5*x^2 + 4*x^3 + 19*x^4 + 15*x^5 + 71*x^6 + 56*x^7 + ...

Harvey P. Dale, Table of n, a(n) for n = 0..1000
P. Bala, Notes on 2-periodic continued fractions and Lehmer sequences
Seong Ju Kim, R. Stees, L. Taalman, Sequences of Spiral Knot Determinants, Journal of Integer Sequences, Vol. 19 (2016), #16.1.4.
Ryan Stees, Sequences of Spiral Knot Determinants, Senior Honors Projects, Paper 84, James Madison Univ., May 2016.
E. W. Weisstein, MathWorld: Lehmer Number
Index to divisibility sequences
Index entries for linear recurrences with constant coefficients, signature (0,4,0,-1).

Cf. A001353, A001834.

Cf. A026741, A005013, A084068.

a := proc (n) if `mod`(n, 2) = 1 then 1/sqrt(2)*( ((sqrt(6) + sqrt(2))/2 )^n - ( (sqrt(6) - sqrt(2))/2 )^n) else 1/sqrt(12)*( ((sqrt(6) + sqrt(2))/2 )^n - ( (sqrt(6) - sqrt(2))/2 )^n) end if;
end proc:
seq(simplify(a(n)), n = 1..30); # Peter Bala, Mar 25 2018

CoefficientList[Series[(1+x+x^2)/(1-4x^2+x^4),{x,0,40}],x] (* or *) LinearRecurrence[{0,4,0,-1},{1,1,5,4},40] (* Harvey P. Dale, Nov 15 2012 *)

{a(n) = my( w = quadgen(24)); simplify( polchebyshev( n, 2, w/2) / if( n%2, w, 1))}; /* Michael Somos, Feb 10 2015 */

a(0)=a(1)=1, a(2)=5, a(n)a(n+3) - a(n+1)a(n+2) = -1.
a(0)=1, a(1)=1, a(2)=5, a(3)=4, a(n) = 4*a(n-2)-a(n-4). - Harvey P. Dale, Nov 15 2012
a(n) = (alpha^n - beta^n)/(alpha - beta) for n odd, and a(n) = (alpha^n - beta^n)/(alpha^2 - beta^2) for n even, where alpha = (1/2)*(sqrt(6) + sqrt(2)) (A188887) and beta = (1/2)*(sqrt(6) - sqrt(2)) (A101263). Equivalently, a(n) = U(n-1,sqrt(6)/2) for n odd and a(n) = (1/sqrt(6))*U(n-1,sqrt(6)/2) for n even, where U(n,x) is the Chebyshev polynomial of the second kind. - Peter Bala, Apr 18 2014
a(2*n) = A001834(n). a(2*n + 1) = A001353(n+1). - Michael Somos, Feb 10 2015
a(n) = -a(-2-n) for all n in Z. - Michael Somos, Feb 10 2015

A227418 Array A(n,k) with all numbers m such that 3*m^2 +- 3^k is a square and their corresponding square roots, read by downward antidiagonals.

Original entry on oeis.org

0, 1, 1, 0, 2, 4, 3, 3, 7, 15, 0, 6, 12, 26, 56, 9, 9, 21, 45, 97, 209, 0, 18, 36, 78, 168, 362, 780, 27, 27, 63, 135, 291, 627, 1351, 2911, 0, 54, 108, 234, 504, 1086, 2340, 5042, 10864, 81, 81, 189, 405, 873, 1881, 4053, 8733, 18817, 40545

Offset: 0

Richard R. Forberg, Sep 02 2013

Array is analogous to A228405 in goal and structure, with key differences.
Left column is A001353.  Top row (not in OEIS) interleaves 0 with the powers of 3, as: 0, 1, 0, 3, 0, 9, 0, 27, 0, 81.
Either or both may be used as initializing values. See Formula section.
The left column is the second binomial transform of the top row. The intermediate transform sequence is A002605, not present in this array.
The columns of the array hold all values, in sequential order, of numbers m such that 3*m^2 + 3^k or 3*m^2 - 3^k are squares, and their corresponding square roots in the next column, which then form the "next round" of m values for column k+1.
For example: A(n,0) are numbers such that 3*m^2 + 1 are squares, the integer square roots of each are in A(n,1), which are then numbers m such that 3*m^2 - 3 are squares, with those square roots in A(n,2), etc.  The sign alternates for each increment of k, etc. No integer square roots exist for the opposite sign in a given column, regardless of n.
Also, A(n,1) are values of m such that floor(m^2/3) is square, with the corresponding square roots given by A(n,0).
A(n,k)/A(n,k-2) = 3; A(n,k)/A(n,k-1) converges to sqrt(3) for large n.
A(n,k)/A(n-1,k) converges to 2 + sqrt(3) for large n.
Several ways of combining the first few columns give OEIS sequences:
  A(n,0) + A(n,1) = A001835; A(n,1) + A(n,2)= A001834; A(n,2) + A(n,3) = A082841;
  A(n,0)*A(n,1)/2 = A007655(n); A(n+2,0)*A(n+1,1) = A001922(n);
  A(n,0)*A(n+1,1) = A001921(n); A(n,0)^2 + A(n,1)^2 = A103974(n);
  A(n,1)^2 - A(n,0)^2 = A011922(n); (A(n+2,0)^2 + A(n+1,1)^2)/2 = A122770(n) = 2*A011916(n).
The main diagonal (without initial 0) = 2*A090018.  The first subdiagonal =  abs(A099842). First superdiagonal = A141041.
A001353 (in left column) are the only initializing set of numbers where the recursive square root equation (see below) produces exclusively integer values, for all iterations of k.  For any other initial values only even iterations (at k = 2, 4, ...) produce integers.

			The array, A(n, k), begins as:
    0,    1,    0,    3,    0,     9,     0,    27, ... see A000244;
    1,    2,    3,    6,    9,    18,    27,    54, ... A038754;
    4,    7,   12,   21,   36,    63,   108,   189, ... A228879;
   15,   26,   45,   78,  135,   234,   405,   702, ...
   56,   97,  168,  291,  504,   873,  1512,  2619, ...
  209,  362,  627, 1086, 1881,  3258,  5643,  9774, ...
  780, 1351, 2340, 4053, 7020, 12159, 21060, 36477, ...
Antidiagonal triangle, T(n, k), begins as:
   0;
   1,  1;
   0,  2,   4;
   3,  3,   7,  15;
   0,  6,  12,  26,  56;
   9,  9,  21,  45,  97,  209;
   0, 18,  36,  78, 168,  362,  780;
  27, 27,  63, 135, 291,  627, 1351, 2911;
   0, 54, 108, 234, 504, 1086, 2340, 5042, 10864;
  81, 81, 189, 405, 873, 1881, 4053, 8733, 18817, 40545;

G. C. Greubel, Antidiagonals n = 0..50, flattened

Cf. A001353, A001075, A005320, A038754, A090018, A099842, A141041, A151961.

Cf. A000244, A084156, A228879.

function A(n,k)
  if k lt 0 then return 0;
  elif n eq 0 then return Round((1/2)*(1-(-1)^k)*3^((k-1)/2));
  elif k eq 0 then return Evaluate(ChebyshevSecond(n), 2);
  else return 2*A(n, k-1) - A(n-1, k-1);
  end if; return A;
end function;
A227418:= func< n,k | A(k, n-k) >;
[A227418(n,k): k in [0..n], n in [0..15]]; // G. C. Greubel, Oct 09 2022

A[n_, k_]:= If[k<0, 0, If[k==0, ChebyshevU[n-1, 2], 2*A[n, k-1] - A[n-1, k-1]]];
T[n_, k_]:= A[k, n-k];
Table[T[n, k], {n,0,15}, {k,0,n}]//Flatten (* G. C. Greubel, Oct 09 2022 *)

def A(n,k):
    if (k<0): return 0
    elif (k==0): return chebyshev_U(n-1,2)
    else: return 2*A(n, k-1) - A(n-1, k-1)
def A227418(n, k): return A(k, n-k)
flatten([[A227418(n,k) for k in range(n+1)] for n in range(15)]) # G. C. Greubel, Oct 09 2022

If using the left column and top row to initialize, then: A(n,k) = 2*A(n, k-1) - A(n-1, k-1).
If using only the top row to initialize, then: A(n,k) = 4*A(n-1,k) - A(n-2,k).
If using the left column to initialize, then: A(n,k) = sqrt(3*A(n,k-1) + (-3)^(k-1)), for all n, k > 0.
Other internal relationships that apply are: A(2*n-1, 2*k) = A(n,k)^2 - A(n-1,k)^2;
A(n+1,k) * A(n,k+1) - A(n+1, k+1) * A(n,k) = (-3)^k, for all n, k > 0.
A(n, 0) = A001353(n).
A(n, 1) = A001075(n).
A(n, 2) = A005320(n).
A(n, 3) = A151961(n).
A(1, k) = A038754(k).
A(n, n) = 2*A090018(n), for n > 0 (main diagonal).
A(n, n+1) = A141041(n-1) (superdiagonal).
A(n+1, n) = abs(A099842(n)) (subdiagonal).
From G. C. Greubel, Oct 09 2022: (Start)
T(n, 0) = (1/2)*(1-(-1)^n)*3^((n-1)/2).
T(n, 1) = A038754(n-1).
T(n, 2) = A228879(n-2).
T(2*n-1, n-1) = A141041(n-1).
T(2*n, n) = 2*A090018(n-1), n > 0.
T(n, n-4) = 3*A005320(n-4).
T(n, n-3) = 3*A001075(n-3).
T(n, n-2) = 3*A001353(n-2).
T(n, n-1) = A001075(n-1).
T(n, n) = A001353(n).
Sum_{k=0..n-1} T(n, k) = A084156(n).
Sum_{k=0..n} T(n, k) = A084156(n) + A001353(n). (End)

Offset corrected by G. C. Greubel, Oct 09 2022

A126866 a(n) = 13*a(n-1) - a(n-2).

Original entry on oeis.org

1, 14, 181, 2339, 30226, 390599, 5047561, 65227694, 842912461, 10892634299, 140761333426, 1819004700239, 23506299769681, 303762892305614, 3925411300203301, 50726584010337299, 655520180834181586, 8471035766834023319, 109467944788008121561

Offset: 0

Diego A. Penta (diego(AT)alum.mit.edu), Mar 15 2007

Nonnegative x values in solutions to the Diophantine equation 11*x^2 - 15*y^2 = -4. The corresponding y values are in A085260. Note that a(n+1)^2 - a(n)*a(n+2) = 15. - Klaus Purath, Mar 21 2025

Harvey P. Dale, Table of n, a(n) for n = 0..899
K. Andersen, L. Carbone, and D. Penta, Kac-Moody Fibonacci sequences, hyperbolic golden ratios, and real quadratic fields, Journal of Number Theory and Combinatorics, Vol 2, No. 3 pp 245-278, 2011. See Section 9.
Alex Fink, Richard K. Guy, and Mark Krusemeyer, Partitions with parts occurring at most thrice, Contributions to Discrete Mathematics, Vol 3, No 2 (2008), pp. 76-114. See Section 13.
Index entries for linear recurrences with constant coefficients, signature (13,-1).

Cf. A002878, A001834, A030221, A002315.

LinearRecurrence[{13,-1},{1,14},30] (* Harvey P. Dale, Mar 28 2013 *)

[(lucas_number2(n,13,1)-lucas_number2(n-1,13,1))/11 for n in range(1, 16)] # Zerinvary Lajos, Nov 10 2009

a(n) = 13*a(n-1) - a(n-2); a(0)=1, a(1)=14.

G.f.: (x+1)/(x^2-13*x+1). - Harvey P. Dale, Mar 28 2013

Corrected and extended by Harvey P. Dale, Mar 28 2013

A131887 Number of Khalimsky-continuous functions with a three-point codomain.

Original entry on oeis.org

3, 5, 11, 19, 41, 71, 153, 265, 571, 989, 2131, 3691, 7953, 13775, 29681, 51409, 110771, 191861, 413403, 716035, 1542841, 2672279, 5757961, 9973081, 21489003, 37220045, 80198051, 138907099, 299303201, 518408351, 1117014753, 1934726305, 4168755811, 7220496869, 15558008491, 26947261171

Offset: 1

Shiva Samieinia (shiva(AT)math.su.se), Oct 05 2007, Oct 09 2007

Neo Scott, Table of n, a(n) for n = 1..1500
Shiva Samieinia, Digital straight line segments and curves. Licentiate Thesis. Stockholm University, Department of Mathematics, Report 2007:6.
Index entries for linear recurrences with constant coefficients, signature (0,4,0,-1).

Cf. A001045, A000213, A131935, A001834 (bisection), A001835 (bisection)

LinearRecurrence[{0,4,0,-1},{3,5,11,19},40] (* Harvey P. Dale, Jan 01 2017 *)

a(2k) = a(2k-1) + a(2k-2) + a(2k-3) and a(2k-1) = a(2k-2) + 2a(2k-3).
The asymptotic behavior is a(2k) = t(2k) sqrt(3)(2 + sqrt(3))^k, a(2k-1) = t(2k-1)(2 + sqrt(3))^k where t(n) tends to 1/2 + sqrt(3)/6.
G.f.: -x*(-3-5*x+x^2+x^3) / ( 1-4*x^2+x^4 ). - R. J. Mathar, Nov 08 2013

A321118 T(n,k) = A321119(n) - (-1)^kA321119(n-2k)/2 for 0 < k < n, with T(0,0) = 0 and T(n,0) = T(n,n) = A002530(n+1) for n > 0, triangle read by rows; unreduced numerator of the weights of Holladay-Sard's quadrature formula.

Original entry on oeis.org

0, 1, 1, 3, 10, 3, 4, 11, 11, 4, 11, 32, 26, 32, 11, 15, 43, 37, 37, 43, 15, 41, 118, 100, 106, 100, 118, 41, 56, 161, 137, 143, 143, 137, 161, 56, 153, 440, 374, 392, 386, 392, 374, 440, 153, 209, 601, 511, 535, 529, 529, 535, 511, 601, 209

Offset: 0

Franck Maminirina Ramaharo, Nov 01 2018

The n-th row common denominator is factorized out and is given by A321119(n).

Given a continuous function f over the interval [0,n], the best quadrature formula in the sense of Holladay-Sard is given by Integral_{x=0..n} f(x) dx = Sum_{k=0..n} T(n,k)*f(k)/A321119(n). The formula is exact if f belongs to the class of natural cubic splines.

			Triangle begins (denominator is factored out):
    0;                                                 1/4
    1,   1;                                            1/2
    3,  10,   3;                                       1/8
    4,  11,  11,   4;                                  1/10
   11,  32,  26,  32,  11;                             1/28
   15,  43,  37,  37,  43,  15;                        1/38
   41, 118, 100, 106, 100, 118,  41;                   1/104
   56, 161, 137, 143, 143, 137, 161,  56;              1/142
  153, 440, 374, 392, 386, 392, 374, 440, 153;         1/388
  209, 601, 511, 535, 529, 529, 535, 511, 601, 209;    1/530
  ...
If f is a continuous function over the interval [0,3], then the quadrature formula yields Integral_{x=0..3} f(x) d(x) = (1/10)*(4*f(0) + 11*f(1) + 11*f(2) + 4*f(3)).

Harold J. Ahlberg, Edwin N. Nilson and Joseph L. Walsh, The Theory of Splines and Their Applications, Academic Press, 1967. See p. 47, Table 2.5.2.

Franck Maminirina Ramaharo, Rows n = 0..150 of triangle, flattened
Harold J. Ahlberg, Edwin N. Nilson and Joseph L. Walsh, Chapter II The Cubic Spline, Mathematics in Science and Engineering Volume 38 (1967), p. 9-74.
John C. Holladay, A smoothest curve approximation, Math. Comp. Vol. 11 (1957), 233-243.
Peter Köhler, On the weights of Sard's quadrature formulas, CALCOLO Vol. 25 (1988), 169-186.
Leroy F. Meyers and Arthur Sard, Best approximate integration formulas, J. Math. Phys. Vol. 29 (1950), 118-123.
Arthur Sard, Best approximate integration formulas; best approximation formulas, American Journal of Mathematics Vol. 71 (1949), 80-91.
Isaac J. Schoenberg, Spline interpolation and best quadrature formulae, Bull. Amer. Math. Soc. Vol. 70 (1964), 143-148.
Frans Schurer, On natural cubic splines, with an application to numerical integration formulae, EUT report. WSK, Dept. of Mathematics and Computing Science Vol. 70-WSK-04 (1970), 1-32.

Cf. A093735, A093736, A100641, A002176, A100640, A321119, A321120, A321121, A321122.

alpha = (Sqrt[2] + Sqrt[6])/2; T[0,0] = 0;
T[n_, k_] := If[n > 0 && k == 0 || k == n, (alpha^(n + 1) - (-alpha)^(-(n + 1)))/(2*Sqrt[6]*(alpha^n + (-alpha)^(-n))), 1 - (-1)^k*(alpha^(n - 2*k) + (-alpha)^(2*k - n))/(2*(alpha^n + (-alpha)^(-n)))];
a321119[n_] := 2^(-Floor[(n - 1)/2])*((1 - Sqrt[3])^n + (1 + Sqrt[3])^n);
Table[FullSimplify[a321119[n]*T[n, k]],{n, 0, 10}, {k, 0, n}] // Flatten

(b[0] : 0, b[1] : 1, b[2] : 1, b[3] : 3, b[n] := 4*b[n-2] - b[n-4])$ /* A002530 */
d(n) := 2^(-floor((n - 1)/2))*((1 - sqrt(3))^n + (1 + sqrt(3))^n) $ /* A321119 */
T(n, k) := if n = 0 and k = 0 then 0 else if n > 0 and k = 0 or k = n then b[n + 1] else d(n) - (-1)^k*d(n - 2*k)/2$
create_list(ratsimp(T(n, k)), n, 0, 10, k, 0, n);

T(n,k)/A321119(n) = (alpha^(n + 1) - (-alpha)^(-(n + 1)))/(2*sqrt(6)*(alpha^n + (-alpha)^(-n))) if k = 0 or k = n, and 1 - (-1)^k*(alpha^(n - 2*k) + (-alpha)^(2*k - n))/(2*(alpha^n + (-alpha)^(-n))) if 0 < k < n, where alpha = (sqrt(2) + sqrt(6))/2.
T(n,k) = T(n,n-k).
T(n,k) = 4*T(n-2,k) - T(n-4,k), n >= k + 4.
T(2*n+2,k)*A001834(n+1) = A001834(n)*T(2*n,k) + 2*A003500(n)*T(2*n+1,k) for k < 2*n.
T(2*n+3,k)*A003500(n+1) = A003500(n)*T(2*n+1,k) + 2*A001834(n+1)*T(2*n+2,k) for k < 2*n + 1.
Sum_{k=0..n} T(n,k)/A321119(n) = n.

A321119 a(n) = ((1 - sqrt(3))^n + (1 + sqrt(3))^n)/2^floor((n - 1)/2); n-th row common denominator of A321118.

Original entry on oeis.org

4, 2, 8, 10, 28, 38, 104, 142, 388, 530, 1448, 1978, 5404, 7382, 20168, 27550, 75268, 102818, 280904, 383722, 1048348, 1432070, 3912488, 5344558, 14601604, 19946162, 54493928, 74440090, 203374108, 277814198, 759002504, 1036816702, 2832635908, 3869452610

Offset: 0

Franck Maminirina Ramaharo, Nov 01 2018

			a(0) = ((1 - sqrt(3))^0 + (1 + sqrt(3))^0)/2^floor((0 - 1)/2) = 2*(1 + 1) = 4.

Harold J. Ahlberg, Edwin N. Nilson and Joseph L. Walsh, The Theory of Splines and Their Applications, Academic Press, 1967. See p. 47, Table 2.5.2.

Encyclopedia of Mathematics, Quadrature formula
John C. Holladay, A smoothest curve approximation, Math. Comp. Vol. 11 (1957), 233-243.
Peter Köhler, On the weights of Sard's quadrature formulas, CALCOLO Vol. 25 (1988), 169-186.
Leroy F. Meyers and Arthur Sard, Best approximate integration formulas, J. Math. Phys. Vol. 29 (1950), 118-123.
Arthur Sard, Best approximate integration formulas; best approximation formulas, American Journal of Mathematics Vol. 71 (1949), 80-91.
Isaac J. Schoenberg, Spline interpolation and best quadrature formulae, Bull. Amer. Math. Soc. Vol. 70 (1964), 143-148.
Frans Schurer, On natural cubic splines, with an application to numerical integration formulae, EUT report. WSK, Dept. of Mathematics and Computing Science Vol. 70-WSK-04 (1970), 1-32.
Index entries for linear recurrences with constant coefficients, signature (0,4,0,-1).

Cf. A002176 (common denominators of Cotesian numbers).
Cf. A321118, A321120.
Cf. A001834, A002530, A002531, A003500, A005246, A048788, A083336, A125905.

LinearRecurrence[{0, 4, 0, -1}, {4, 2, 8, 10}, 50]

a(n) := ((1 - sqrt(3))^n + (1 + sqrt(3))^n)/2^floor((n - 1)/2)$
makelist(ratsimp(a(n)), n, 0, 50);

a(n) = (((sqrt(2) - sqrt(6))/2)^n + ((sqrt(6) + sqrt(2))/2)^n)*((2 - sqrt(2))*(-1)^n + 2 + sqrt(2))/2.
a(-n) = (-1)^n*a(n).
a(n) = 2*A000034(n+1)*A002531(n).
a(2*n) = 2*A001834(n).
a(2*n+1) = 2*A003500(n).
a(n) = 4*a(n-2) - a(n-4) with a(0) = 4, a(1) = 2, a(2) = 8, a(3) = 10.
a(2*n+3) = a(2*n+1) + a(2*n+2).
a(2*n+2) = a(2*n) + 2*a(2*n+1).
G.f.: 2*(1 - x)*(2 + 3*x - x^2)/(1 - 4*x^2 + x^4).
E.g.f.: (1 + exp(-sqrt(6)*x))*((2 - sqrt(2))*exp(sqrt(2 - sqrt(3))*x) + (2 + sqrt(2))*exp(sqrt(2 + sqrt(3))*x))/2.
Lim_{n->infinity} a(2*n+1)/a(2*n) = (1 + sqrt(3))/2.

cp's OEIS Frontend

A098495 Array T(r,c) read by antidiagonals: values of triangle A098493 interpreted as polynomials, r >= 0.

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A102206 a(0) = 3, a(1) = 8, a(n+2) = 4*a(n+1) - a(n) - 2.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A103772 Larger of two sides in a (k,k,k-1)-integer-sided triangle with integer area.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

Extensions

A107903 Generalized NSW numbers.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

Formula

Extensions

A108412 Expansion of (1 + x + x^2)/(1 - 4x^2 + x^4).

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Formula

A227418 Array A(n,k) with all numbers m such that 3*m^2 +- 3^k is a square and their corresponding square roots, read by downward antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Magma

Mathematica

SageMath

Formula

Extensions

A126866 a(n) = 13*a(n-1) - a(n-2).

Views

Author

Keywords

A321118 T(n,k) = A321119(n) - (-1)^kA321119(n-2k)/2 for 0 < k < n, with T(0,0) = 0 and T(n,0) = T(n,n) = A002530(n+1) for n > 0, triangle read by rows; unreduced numerator of the weights of Holladay-Sard's quadrature formula.