A085260 -id:A085260 - cp's OEIS frontend

A108299 Triangle read by rows, 0 <= k <= n: T(n,k) = binomial(n-[(k+1)/2],[k/2])*(-1)^[(k+1)/2].

1, 1, -1, 1, -1, -1, 1, -1, -2, 1, 1, -1, -3, 2, 1, 1, -1, -4, 3, 3, -1, 1, -1, -5, 4, 6, -3, -1, 1, -1, -6, 5, 10, -6, -4, 1, 1, -1, -7, 6, 15, -10, -10, 4, 1, 1, -1, -8, 7, 21, -15, -20, 10, 5, -1, 1, -1, -9, 8, 28, -21, -35, 20, 15, -5, -1, 1, -1, -10, 9, 36, -28, -56, 35, 35, -15, -6, 1, 1, -1, -11, 10, 45, -36, -84, 56, 70

Offset: 0

Reinhard Zumkeller, Jun 01 2005

Matrix inverse of A124645.
Let L(n,x) = Sum_{k=0..n} T(n,k)*x^(n-k) and Pi=3.14...:
L(n,x) = Product_{k=1..n} (x - 2*cos((2*k-1)*Pi/(2*n+1)));
Sum_{k=0..n} T(n,k) = L(n,1) = A010892(n+1);
Sum_{k=0..n} abs(T(n,k)) = A000045(n+2);
abs(T(n,k)) = A065941(n,k), T(n,k) = A065941(n,k)*A087960(k);
T(2*n,k) + T(2*n+1,k+1) = 0 for 0 <= k <= 2*n;
T(n,0) = A000012(n) = 1; T(n,1) = -1 for n > 0;
T(n,2) = -(n-1) for n > 1; T(n,3) = A000027(n)=n for n > 2;
T(n,4) = A000217(n-3) for n > 3; T(n,5) = -A000217(n-4) for n > 4;
T(n,6) = -A000292(n-5) for n > 5; T(n,7) = A000292(n-6) for n > 6;
T(n,n-3) = A058187(n-3)*(-1)^floor(n/2) for n > 2;
T(n,n-2) = A008805(n-2)*(-1)^floor((n+1)/2) for n > 1;
T(n,n-1) = A008619(n-1)*(-1)^floor(n/2) for n > 0;
T(n,n) = L(n,0) = (-1)^floor((n+1)/2);
L(n,1) = A010892(n+1); L(n,-1) = A061347(n+2);
L(n,2) = 1; L(n,-2) = A005408(n)*(-1)^n;
L(n,3) = A001519(n); L(n,-3) = A002878(n)*(-1)^n;
L(n,4) = A001835(n+1); L(n,-4) = A001834(n)*(-1)^n;
L(n,5) = A004253(n); L(n,-5) = A030221(n)*(-1)^n;
L(n,6) = A001653(n); L(n,-6) = A002315(n)*(-1)^n;
L(n,7) = A049685(n); L(n,-7) = A033890(n)*(-1)^n;
L(n,8) = A070997(n); L(n,-8) = A057080(n)*(-1)^n;
L(n,9) = A070998(n); L(n,-9) = A057081(n)*(-1)^n;
L(n,10) = A072256(n+1); L(n,-10) = A054320(n)*(-1)^n;
L(n,11) = A078922(n+1); L(n,-11) = A097783(n)*(-1)^n;
L(n,12) = A077417(n); L(n,-12) = A077416(n)*(-1)^n;
L(n,13) = A085260(n);
L(n,14) = A001570(n); L(n,-14) = A028230(n)*(-1)^n;
L(n,n) = A108366(n); L(n,-n) = A108367(n).
Row n of the matrix inverse (A124645) has g.f.: x^floor(n/2)*(1-x)^(n-floor(n/2)). - Paul D. Hanna, Jun 12 2005
From L. Edson Jeffery, Mar 12 2011: (Start)
Conjecture: Let N=2*n+1, with n > 2. Then T(n,k) (0 <= k <= n) gives the k-th coefficient in the characteristic function p_N(x)=0, of degree n in x, for the n X n tridiagonal unit-primitive matrix G_N (see [Jeffery]) of the form
  G_N=A_{N,1}=
  (0 1 0 ... 0)
  (1 0 1 0 ... 0)
  (0 1 0 1 0 ... 0)
  ...
  (0 ... 0 1 0 1)
  (0 ... 0 1 1),
with solutions phi_j = 2*cos((2*j-1)*Pi/N), j=1,2,...,n. For example, for n=3,
  G_7=A_{7,1}=
  (0 1 0)
  (1 0 1)
  (0 1 1).
We have {T(3,k)}=(1,-1,-2,1), while the characteristic function of G_7 is p(x) = x^3-x^2-2*x+1 = 0, with solutions phi_j = 2*cos((2*j-1)*Pi/7), j=1,2,3. (End)
The triangle sums, see A180662 for their definitions, link A108299 with several sequences, see the crossrefs. - Johannes W. Meijer, Aug 08 2011
The roots to the polynomials are chaotic using iterates of the operation (x^2 - 2), with cycle lengths L and initial seeds returning to the same term or (-1)* the seed. Periodic cycle lengths L are shown in A003558 such that for the polynomial represented by row r, the cycle length L is A003558(r-1). The matrices corresponding to the rows as characteristic polynomials are likewise chaotic [cf. Kappraff et al., 2005] with the same cycle lengths but substituting 2*I for the "2" in (x^2 - 2), where I = the Identity matrix. For example, the roots to x^3 - x^2 - 2x + 1 = 0 are 1.801937..., -1.246979..., and 0.445041... With 1.801937... as the initial seed and using (x^2 - 2), we obtain the 3-period trajectory of 8.801937... -> 1.246979... -> -0.445041... (returning to -1.801937...). We note that A003558(2) = 3. The corresponding matrix M is: [0,1,0; 1,0,1; 0,1,1,]. Using seed M with (x^2 - 2*I), we obtain the 3-period with the cycle completed at (-1)*M. - Gary W. Adamson, Feb 07 2012

			Triangle begins:
  1;
  1,  -1;
  1,  -1,  -1;
  1,  -1,  -2,   1;
  1,  -1,  -3,   2,   1;
  1,  -1,  -4,   3,   3,  -1;
  1,  -1,  -5,   4,   6,  -3,  -1;
  1,  -1,  -6,   5,  10,  -6,  -4,   1;
  1,  -1,  -7,   6,  15, -10, -10,   4,   1;
  1,  -1,  -8,   7,  21, -15, -20,  10,   5,  -1;
  1,  -1,  -9,   8,  28, -21, -35,  20,  15,  -5,  -1;
  1,  -1, -10,   9,  36, -28, -56,  35,  35, -15,  -6,   1;
  ...

Friedrich L. Bauer, 'De Moivre und Lagrange: Cosinus eines rationalen Vielfachen von Pi', Informatik Spektrum 28 (Springer, 2005).
Jay Kappraff, S. Jablan, G. Adamson, & R. Sazdonovich: "Golden Fields, Generalized Fibonacci Sequences, & Chaotic Matrices"; FORMA, Vol 19, No 4, (2005).

Reinhard Zumkeller, Rows n = 0..150 of triangle, flattened
Henry W. Gould, A Variant of Pascal's Triangle, Corrections, The Fibonacci Quarterly, Vol. 3, Nr. 4, Dec. 1965, p. 257-271.
L. Edson Jeffery, Unit-primitive matrices.
Ju, Hyeong-Kwan On the sequence generated by a certain type of matrices. Honam Math. J. 39, No. 4, 665-675 (2017).
Michelle Rudolph-Lilith, On the Product Representation of Number Sequences, with Application to the Fibonacci Family, arXiv preprint arXiv:1508.07894 [math.NT], 2015.
Frank Ruskey and Carla Savage, Gray codes for set partitions and restricted growth tails, Australasian Journal of Combinatorics, Volume 10(1994), pp. 85-96. See Table 1 p. 95.

Cf. A049310, A039961, A124645 (matrix inverse).
Triangle sums (see the comments): A193884 (Kn11), A154955 (Kn21), A087960 (Kn22), A000007 (Kn3), A010892 (Fi1), A134668 (Fi2), A078031 (Ca2), A193669 (Gi1), A001519 (Gi3), A193885 (Ze1), A050935 (Ze3). - Johannes W. Meijer, Aug 08 2011
Cf. A003558.
Cf. A033999, A059841.

a108299 n k = a108299_tabl !! n !! k
a108299_row n = a108299_tabl !! n
a108299_tabl = [1] : iterate (\row ->
   zipWith (+) (zipWith (*) ([0] ++ row) a033999_list)
               (zipWith (*) (row ++ [0]) a059841_list)) [1,-1]
-- Reinhard Zumkeller, May 06 2012

A108299 := proc(n,k): binomial(n-floor((k+1)/2), floor(k/2))*(-1)^floor((k+1)/2) end: seq(seq(A108299 (n,k), k=0..n), n=0..11); # Johannes W. Meijer, Aug 08 2011

t[n_, k_?EvenQ] := I^k*Binomial[n-k/2, k/2]; t[n_, k_?OddQ] := -I^(k-1)*Binomial[n+(1-k)/2-1, (k-1)/2]; Table[t[n, k], {n, 0, 12}, {k, 0, n}] // Flatten (* Jean-François Alcover, May 16 2013 *)

{T(n,k)=polcoeff(polcoeff((1-x*y)/(1-x+x^2*y^2+x^2*O(x^n)),n,x)+y*O(y^k),k,y)} (Hanna)

T(n,k) = binomial(n-floor((k+1)/2),floor(k/2))*(-1)^floor((k+1)/2).
T(n+1, k) = if sign(T(n, k-1))=sign(T(n, k)) then T(n, k-1)+T(n, k) else -T(n, k-1) for 0 < k < n, T(n, 0) = 1, T(n, n) = (-1)^floor((n+1)/2).
G.f.: A(x, y) = (1 - x*y)/(1 - x + x^2*y^2). - Paul D. Hanna, Jun 12 2005
The generating polynomial (in z) of row n >= 0 is (u^(2*n+1) + v^(2*n+1))/(u + v), where u and v are defined by u^2 + v^2 = 1 and u*v = z. - Emeric Deutsch, Jun 16 2011
From Johannes W. Meijer, Aug 08 2011: (Start)
abs(T(n,k)) = A065941(n,k) = abs(A187660(n,n-k));
T(n,n-k) = A130777(n,k); abs(T(n,n-k)) = A046854(n,k) = abs(A066170(n,k)). (End)

Corrected and edited by Philippe Deléham, Oct 20 2008

A094954 Array T(k,n) read by antidiagonals. G.f.: x(1-x)/(1-kx+x^2), k>1.

Original entry on oeis.org

1, 1, 1, 1, 2, 1, 1, 3, 5, 1, 1, 4, 11, 13, 1, 1, 5, 19, 41, 34, 1, 1, 6, 29, 91, 153, 89, 1, 1, 7, 41, 169, 436, 571, 233, 1, 1, 8, 55, 281, 985, 2089, 2131, 610, 1, 1, 9, 71, 433, 1926, 5741, 10009, 7953, 1597, 1, 1, 10, 89, 631, 3409, 13201, 33461, 47956, 29681

Offset: 1

Ralf Stephan, May 31 2004

Also, values of polynomials with coefficients in A098493 (see Fink et al.). See A098495 for negative k.

Number of dimer tilings of the graph S_{k-1} X P_{2n-2}.

			1,1,1,1,1,1,1,1,1,1,1,1,1,1, ...
1,2,5,13,34,89,233,610,1597, ...
1,3,11,41,153,571,2131,7953, ...
1,4,19,91,436,2089,10009,47956, ...
1,5,29,169,985,5741,33461,195025, ...
1,6,41,281,1926,13201,90481,620166, ...

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.
Elizabeth Wilmer, A note on Stephan's conjecture 87
Elizabeth Wilmer, A note on Stephan's conjecture 87 [cached copy]

Rows are first differences of rows in array A073134.
Rows 2-14 are A000012, A001519, A079935/A001835, A004253, A001653, A049685, A070997, A070998, A072256, A078922, A077417, A085260, A001570. Other rows: A007805 (k=18), A075839 (k=20), A077420 (k=34), A078988 (k=66).
Columns include A028387. Diagonals include A094955, A094956. Antidiagonal sums are A094957.

max = 14; row[k_] := Rest[ CoefficientList[ Series[ x*(1-x)/(1-k*x+x^2), {x, 0, max}], x]]; t = Table[ row[k], {k, 2, max+1}]; Flatten[ Table[ t[[k-n+1, n]], {k, 1, max}, {n, 1, k}]] (* Jean-François Alcover, Dec 27 2011 *)

PARI
```
T(k,n)=polcoeff(x*(1-x)/(1-k*x+x*x),n)
```

Recurrence: T(k, 1) = 1, T(k, 2) = k-1, T(k, n) = kT(k, n-1) - T(k, n-2).
For n>3, T(k, n) = [k(k-2) + T(k, n-1)T(k, n-2)] / T(k, n-3).
T(k, n+1) = S(n, k) - S(n-1, k) = U(n, k/2) - U(n-1, k/2), with S, U = Chebyshev polynomials of second kind.
T(k+2, n+1) = Sum[i=0..n, k^(n-i) * C(2n-i, i)] (from comments by Benoit Cloitre).

A238379 Expansion of (1 - x)/(1 - 36*x + x^2).

Original entry on oeis.org

1, 35, 1259, 45289, 1629145, 58603931, 2108112371, 75833441425, 2727895778929, 98128414600019, 3529895029821755, 126978092658983161, 4567681440693572041, 164309553772309610315, 5910576254362452399299, 212616435603275976764449

Offset: 0

Bruno Berselli, Feb 25 2014

First bisection of A041611.

Bruno Berselli, Table of n, a(n) for n = 0..100
Tanya Khovanova, Recursive Sequences
Index entries for linear recurrences with constant coefficients, signature (36,-1).

Cf. similar sequences with g.f. (1-x)/(1-k*x+x^2): A122367 (k=3), A079935 (k=4), A004253 (k=5), A001653 (k=6), A049685 (k=7), A070997 (k=8), A070998 (k=9), A138288 (k=10), A078922 (k=11), A077417 (k=12), A085260 (k=13), A001570 (k=14), A160682 (k=15), A157456 (k=16), A161595 (k=17). From 18 to 38, even k only, except k=27 and k=31: A007805 (k=18), A075839 (k=20), A157014 (k=22), A159664 (k=24), A153111 (k=26), A097835 (k=27), A159668 (k=28), A157877 (k=30), A111216 (k=31), A159674 (k=32), A077420 (k=34), this sequence (k=36), A097315 (k=38).

[n le 2 select 35^(n-1) else 36*Self(n-1)-Self(n-2): n in [1..20]];

R:=PowerSeriesRing(Integers(), 20); Coefficients(R!( (1 - x)/(1 - 36*x + x^2))); // Marius A. Burtea, Jan 14 2020

CoefficientList[Series[(1 - x)/(1 - 36 x + x^2), {x, 0, 20}], x] (* or *) LinearRecurrence[{36, -1}, {1, 35}, 20]

a(n)=([0,1; -1,36]^n*[1;35])[1,1] \\ Charles R Greathouse IV, May 10 2016

m = 20; L. = PowerSeriesRing(ZZ, m); f = (1-x)/(1-36*x+x^2)
print(f.coefficients())

G.f.: (1 - x)/(1 - 36*x + x^2).
a(n) = a(-n-1) = 36*a(n-1) - a(n-2).
a(n) = ((19-sqrt(323))/38)*(1+(18+sqrt(323))^(2*n+1))/(18+sqrt(323))^n.
a(n+1) - a(n) = 34*A144128(n+1).
323*a(n+1)^2 - ((a(n+2)-a(n))/2)^2 = 34.
Sum_{n>0} 1/(a(n) - 1/a(n)) = 1/34.
See also Tanya Khovanova in Links field:
a(n) = 35*a(n-1) + 34*Sum_{i=0..n-2} a(i).
a(n+2)*a(n) - a(n+1)^2 = 36-2 = 34 = 34*1,
a(n+3)*a(n) - a(n+1)*a(n+2) = 36*(36-2) = 1224 = 34*36.
Generalizing:
a(n+4)*a(n) - a(n+1)*a(n+3) = 44030 = 34*1295,
a(n+5)*a(n) - a(n+1)*a(n+4) = 1583856 = 34*46584,
a(n+6)*a(n) - a(n+1)*a(n+5) = 56974786 = 34*1675729, etc.,
where 1, 36, 1295, 46584, 1675729, ... is the sequence A144128, which is the second bisection of A041611.
a(n)^2 - 36*a(n)*a(n+1) + a(n+1)^2 + 34 = 0 (see comments by Colin Barker in similar sequences).

A165253 Triangle T(n,k), read by rows given by [1,0,1,0,0,0,0,0,0,...] DELTA [0,1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.

Original entry on oeis.org

1, 1, 0, 1, 1, 0, 1, 3, 1, 0, 1, 6, 5, 1, 0, 1, 10, 15, 7, 1, 0, 1, 15, 35, 28, 9, 1, 0, 1, 21, 70, 84, 45, 11, 1, 0, 1, 28, 126, 210, 165, 66, 13, 1, 0, 1, 36, 210, 462, 495, 286, 91, 15, 1, 0, 1, 45, 330, 924, 1287, 1001, 455, 120, 17, 1, 0, 1, 55, 495, 1716, 3003, 3003, 1820, 680

Offset: 0

Philippe Deléham, Sep 10 2009

Mirror image of triangle in A121314.

			Triangle begins:
  1;
  1,    0;
  1,    1,    0;
  1,    3,    1,    0;
  1,    6,    5,    1,    0;
  1,   10,   15,    7,    1,    0;
  1,   15,   35,   28,    9,    1,    0;
  1,   21,   70,   84,   45,   11,    1,    0;
  1,   28,  126,  210,  165,   66,   13,    1,    0;
  1,   36,  210,  462,  495,  286,   91,   15,    1,    0,
  1,   45,  330,  924, 1287, 1001,  455,  120,   17,    1,    0;

Indranil Ghosh, Rows 0..125 of triangle, flattened

Cf. A054142, A085478, A121314.

m = 13;
(* DELTA is defined in A084938 *)
DELTA[Join[{1, 0, 1}, Table[0, {m}]], Join[{0, 1}, Table[0, {m}]], m] // Flatten (* Jean-François Alcover, Feb 19 2020 *)

T(0,0)=1, T(n,k) = binomial(n-1+k,2k) for n >= 1.
Sum {k=0..n} T(n,k)*x^k = A000012(n), A001519(n), A001835(n), A004253(n), A001653(n), A049685(n-1), A070997(n-1), A070998(n-1), A072256(n), A078922(n), A077417(n-1), A085260(n), A001570(n) for x = 0,1,2,3,4,5,6,7,8,9,10,11,12 respectively.
Sum_{k=0..n} T(n,k)*x^(n-k) = A000007(n), A001519(n), A047849(n), A165310(n), A165311(n), A165312(n), A165314(n), A165322(n), A165323(n), A165324(n) for x= 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively. - Philippe Deléham, Sep 26 2009
T(n,k) = 2*T(n-1,k) + T(n-1,k-1) - T(n-2,k), T(0,0)=T(1,0)=1, T(1,1)=0. - Philippe Deléham, Feb 18 2012
G.f.: (1-x-y*x)/((1-x)^2-y*x). - Philippe Deléham, Feb 19 2012

A160682 The list of the A values in the common solutions to 13k+1 = A^2 and 17k+1 = B^2.

Original entry on oeis.org

1, 14, 209, 3121, 46606, 695969, 10392929, 155197966, 2317576561, 34608450449, 516809180174, 7717529252161, 115246129602241, 1720974414781454, 25699370092119569, 383769576967012081, 5730844284413061646, 85578894689228912609, 1277952576054020627489

Offset: 1

Paul Weisenhorn, May 23 2009

This summarizes the case C=13 of common solutions to C*k+1=A^2, (C+4)*k+1=B^2.
The 2 equations are equivalent to the Pell equation x^2-C*(C+4)*y^2=1,
with x=(C*(C+4)*k+C+2)/2; y=A*B/2 and with smallest values x(1) = (C+2)/2, y(1)=1/2.
Generic recurrences are:
A(j+2)=(C+2)*A(j+1)-A(j) with A(1)=1; A(2)=C+1.
B(j+2)=(C+2)*B(j+1)-B(j) with B(1)=1; B(2)=C+3.
k(j+3)=(C+1)*(C+3)*( k(j+2)-k(j+1) )+k(j) with k(1)=0; k(2)=C+2; k(3)=(C+1)*(C+2)*(C+3).
x(j+2)=(C^2+4*C+2)*x(j+1)-x(j) with x(1)=(C+2)/2; x(2)=(C^2+4*C+1)*(C+2)/2;
Binet-type of solutions of these 2nd order recurrences are:
R=C^2+4*C; S=C*sqrt(R); T=(C+2); U=sqrt(R); V=(C+4)*sqrt(R);
A(j)=((R+S)*(T+U)^(j-1)+(R-S)*(T-U)^(j-1))/(R*2^j);
B(j)=((R+V)*(T+U)^(j-1)+(R-V)*(T-U)^(j-1))/(R*2^j);
x(j)+sqrt(R)*y(j)=((T+U)*(C^2*4*C+2+(C+2)*sqrt(R))^(j-1))/2^j;
k(j)=(((T+U)*(R+2+T*U)^(j-1)+(T-U)*(R+2-T*U)^(j-1))/2^j-T)/R. [Paul Weisenhorn, May 24 2009]
.C -A----- -B----- -k-----
01 A001519 A002878 A058038
02 A001653 A002315 A045899/2
03 A004253 A030221 A160695
04 A001653 A002315 A078522/4
05 A049685 A033890 A161582
06 A070997 A057080 A159683/2
07 A070998 A057081 A161585
08 A072256 A054320 A045502/4
09 A078922 A097783 A161586
10 A077417 A077416 A159681/2
11 A085260 A126816 A161588
12 A001570 A028230 A059989/4
13 A160682 A161591 A161584
14 A157456 A159678 A159679/2
15 A161595 A161599 A161583
16 A007805 A049629 A157459/4
For n>=2, a(n) equals the permanent of the (2n-2)X(2n-2) tridiagonal matrix with sqrt(13)'s along the main diagonal, and 1's along the superdiagonal and the subdiagonal. [John M. Campbell, Jul 08 2011]
Positive values of x (or y) satisfying x^2 - 15xy + y^2 + 13 = 0. - Colin Barker, Feb 11 2014

Vincenzo Librandi, Table of n, a(n) for n = 1..200
J.-C. Novelli, J.-Y. Thibon, Hopf Algebras of m-permutations,(m+1)-ary trees, and m-parking functions, arXiv preprint arXiv:1403.5962 [math.CO], 2014.
Index entries for the Pell equation
Index entries for linear recurrences with constant coefficients, signature (15,-1).

Cf. similar sequences listed in A238379.

I:=[1,14]; [n le 2 select I[n] else 15*Self(n-1)-Self(n-2): n in [1..30]]; // Vincenzo Librandi, Feb 12 2014

LinearRecurrence[{15,-1},{1,14},20] (* Harvey P. Dale, Oct 08 2012 *)
CoefficientList[Series[(1 - x)/(1 - 15 x + x^2), {x, 0, 40}], x] (* Vincenzo Librandi, Feb 12 2014 *)

a(n) = round((2^(-1-n)*((15-sqrt(221))^n*(13+sqrt(221))+(-13+sqrt(221))*(15+sqrt(221))^n))/sqrt(221)) \\ Colin Barker, Jul 25 2016

a(n) = 15*a(n-1)-a(n-2).
G.f.: (1-x)*x/(1-15*x+x^2).
a(n) = (2^(-1-n)*((15-sqrt(221))^n*(13+sqrt(221))+(-13+sqrt(221))*(15+sqrt(221))^n))/sqrt(221). - Colin Barker, Jul 25 2016

Edited, extended by R. J. Mathar, Sep 02 2009

First formula corrected by Harvey P. Dale, Oct 08 2012

A133607 Triangle read by rows: T(n, k) = qStirling2(n, k, q) for q = -1, with 0 <= k <= n.

Original entry on oeis.org

1, 0, 1, 0, 1, -1, 0, 1, -1, -1, 0, 1, -1, -2, 1, 0, 1, -1, -3, 2, 1, 0, 1, -1, -4, 3, 3, -1, 0, 1, -1, -5, 4, 6, -3, -1, 0, 1, -1, -6, 5, 10, -6, -4, 1, 0, 1, -1, -7, 6, 15, -10, -10, 4, 1, 0, 1, -1, -8, 7, 21, -15, -20, 10, 5, -1, 0, 1, -1, -9, 8, 28, -21, -35, 20, 15, -5, -1

Offset: 0

Philippe Deléham, Dec 27 2007

Previous name: Triangle T(n,k), 0<=k<=n, read by rows given by [0, 1, 0, 0, 0, 0, 0, 0, 0, ...] DELTA [1, -2, 1, 0, 0, 0, 0, 0, 0, 0, ...] where DELTA is the operator defined in A084938.

			Triangle begins:
  1;
  0, 1;
  0, 1, -1;
  0, 1, -1, -1;
  0, 1, -1, -2, 1;
  0, 1, -1, -3, 2, 1;
  0, 1, -1, -4, 3, 3, -1;
  0, 1, -1, -5, 4, 6, -3, -1;
  0, 1, -1, -6, 5, 10, -6, -4, 1;
  0, 1, -1, -7, 6, 15, -10, -10, 4, 1;
  0, 1, -1, -8, 7, 21, -15, -20, 10, 5, -1;
  0, 1, -1, -9, 8, 28, -21, -35, 20, 15, -5, -1;
  0, 1, -1, -10, 9, 36, -28, -56, 35, 35, -15, -6, 1;
  ...
Triangle A103631 begins:
  1;
  0, 1;
  0, 1, 1;
  0, 1, 1, 1;
  0, 1, 1, 2, 1;
  0, 1, 1, 3, 2, 1;
  0, 1, 1, 4, 3, 3, 1;
  0, 1, 1, 5, 4, 6, 3, 1;
  0, 1, 1, 6, 5, 10, 6, 4, 1;
  0, 1, 1, 7, 6, 15, 10, 10, 4, 1;
  0, 1, 1, 8, 7, 21, 15, 20, 10, 5, 1;
  0, 1, 1, 9, 8, 28, 21, 35, 20, 15, 5, 1;
  0, 1, 1, 10, 9, 36, 28, 56, 35, 35, 15, 6, 1;
  ...
Triangle A108299 begins:
  1;
  1, -1;
  1, -1, -1;
  1, -1, -2, 1;
  1, -1, -3, 2, 1;
  1, -1, -4, 3, 3, -1;
  1, -1, -5, 4, 6, -3, -1;
  1, -1, -6, 5, 10, -6, -4, 1;
  1, -1, -7, 6, 15, -10, -10, 4, 1;
  1, -1, -8, 7, 21, -15, -20, 10, 5, -1;
  1, -1, -9, 8, 28, -21, -35, 20, 15, -5, -1;
  1, -1, -10, 9, 36, -28, -56, 35, 35, -15, -6, 1;
  ...

Another version is A108299.

Unsigned version is A103631 (T(n,k) = A103631(n,k)*A057077(k)).

m = 13
(* DELTA is defined in A084938 *)
DELTA[Join[{0, 1}, Table[0, {m}]], Join[{1, -2, 1}, Table[0, {m}]], m] // Flatten (* Jean-François Alcover, Feb 19 2020 *)
qStirling2[n_, k_, q_] /; 1 <= k <= n := q^(k-1) qStirling2[n-1, k-1, q] + Sum[q^j, {j, 0, k-1}] qStirling2[n-1, k, q];
qStirling2[n_, 0, _] := KroneckerDelta[n, 0];
qStirling2[0, k_, _] := KroneckerDelta[0, k];
qStirling2[, , _] = 0;
Table[qStirling2[n, k, -1], {n, 0, 11}, {k, 0, n}] // Flatten (* Jean-François Alcover, Mar 10 2020 *)

from sage.combinat.q_analogues import q_stirling_number2
for n in (0..9):
    print([q_stirling_number2(n,k).substitute(q=-1) for k in [0..n]])
# Peter Luschny, Mar 09 2020

Sum_{k, 0<=k<=n}T(n,k)*x^(n-k)= A057077(n), A010892(n), A000012(n), A001519(n), A001835(n), A004253(n), A001653(n), A049685(n-1), A070997(n-1), A070998(n-1), A072256(n), A078922(n), A077417(n-1), A085260(n), A001570(n-1) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14 respectively .
Sum_{k, 0<=k<=n}T(n,k)*x^k = A000007(n), A010892(n), A133631(n), A133665(n), A133666(n), A133667(n), A133668(n), A133669(n), A133671(n), A133672(n) for x = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9 respectively .
G.f.: (1-x+y*x)/(1-x+y^2*x^2). - Philippe Deléham, Mar 14 2012
T(n,k) = T(n-1,k) - T(n-2,k-2), T(0,0) = T(1,1) = T(2,1) = 1, T(1,0) = T(2,0) = 0, T(2,2) = -1 and T(n,k) = 0 if k<0 or if k>n. - Philippe Deléham, Mar 14 2012

New name from Peter Luschny, Mar 09 2020

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

A123971 Triangle T(n,k), read by rows, defined by T(n,k)=3*T(n-1,k)-T(n-1,k-1)-T(n-2,k), T(0,0)=1, T(1,0)=2, T(1,1)=-1, T(n,k)=0 if k<0 or if k>n.

Original entry on oeis.org

1, 2, -1, 5, -5, 1, 13, -19, 8, -1, 34, -65, 42, -11, 1, 89, -210, 183, -74, 14, -1, 233, -654, 717, -394, 115, -17, 1, 610, -1985, 2622, -1825, 725, -165, 20, -1, 1597, -5911, 9134, -7703, 3885, -1203, 224, -23, 1, 4181, -17345, 30691, -30418, 18633, -7329

Offset: 0

Gary W. Adamson and Roger L. Bagula, Oct 30 2006

This entry is the result of merging two sequences, this one and a later submission by Philippe Deléham, Nov 29 2013 (with edits from Ralf Stephan, Dec 12 2013). Most of the present version is the work of Philippe Deléham, the only things remaining from the original entry are the sequence data and the Mathematica program. - N. J. A. Sloane, May 31 2014
Subtriangle of the triangle given by (0, 2, 1/2, 1/2, 0, 0, 0, 0, 0, 0, 0, ...) DELTA (1, -2, 0, 0, 0, 0, 0, 0, 0, ...) where DELTA is the operator defined in A084938.
Apart from signs, equals A126124.
Row sums = 1.
Sum_{k=0..n} T(n,k)*(-x)^k = A001519(n+1), A079935(n+1), A004253(n+1), A001653(n+1), A049685(n), A070997(n), A070998(n), A072256(n+1), A078922(n+1), A077417(n), A085260(n+1), A001570(n+1) for x=0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11 respectively.

			Triangle begins:
  1
  2, -1
  5, -5, 1
  13, -19, 8, -1
  34, -65, 42, -11, 1
  89, -210, 183, -74, 14, -1
  233, -654, 717, -394, 115, -17, 1
Triangle (0, 2, 1/2, 1/2, 0, 0, ...) DELTA (1, -2, 0, 0, ...) begins:
  1
  0, 1
  0, 2, -1
  0, 5, -5, 1
  0, 13, -19, 8, -1
  0, 34, -65, 42, -11, 1
  0, 89, -210, 183, -74, 14, -1
  0, 233, -654, 717, -394, 115, -17, 1

Cf. A094954, A098495, A123971, A126124, A152063, A001519, A079935, A004253, A001653, A049685, A070997, A070998, A072256, A078922, A077417, A085260, A001570, A001870, A126124.

Mathematica ( general k th center) Clear[M, T, d, a, x, k] k = 3 T[n_, m_, d_] := If[ n == m && n < d && m < d, k, If[n == m - 1 || n == m + 1, -1, If[n == m == d, k - 1, 0]]] M[d_] := Table[T[n, m, d], {n, 1, d}, {m, 1, d}] Table[M[d], {d, 1, 10}] Table[Det[M[d]], {d, 1, 10}] Table[Det[M[d] - x*IdentityMatrix[d]], {d, 1, 10}] a = Join[{M[1]}, Table[CoefficientList[ Det[M[d] - x*IdentityMatrix[d]], x], {d, 1, 10}]] Flatten[a] MatrixForm[a] Table[NSolve[Det[M[d] - x*IdentityMatrix[d]] == 0, x], {d, 1, 10}] Table[x /. NSolve[Det[M[d] - x*IdentityMatrix[d]] == 0, x][[d]], {d, 1, 10}]

T(n,k)=polcoeff(polcoeff(Ser((1-x)/(1+(y-3)*x+x^2)),n,x),n-k,y) \\ Ralf Stephan, Dec 12 2013

@CachedFunction
def A123971(n,k): # With T(0,0) = 1!
    if n< 0: return 0
    if n==0: return 1 if k == 0 else 0
    h = 2*A123971(n-1,k) if n==1 else 3*A123971(n-1,k)
    return A123971(n-1,k-1) - A123971(n-2,k) - h
for n in (0..9): [A123971(n,k) for k in (0..n)] # Peter Luschny, Nov 20 2012

T(n,k) = (-1)^n*A126124(n+1,k+1).
T(n,k) = (-1)^k*Sum_{m=k..n} binomial(m,k)*binomial(m+n,2*m). - Wadim Zudilin, Jan 11 2012
G.f.: (1-x)/(1+(y-3)*x+x^2).
T(n,0) = A001519(n+1) = A000045(2*n+1).
T(n+1,1) = -A001870(n).

Edited by N. J. A. Sloane, May 31 2014

cp's OEIS Frontend

A108299 Triangle read by rows, 0 <= k <= n: T(n,k) = binomial(n-[(k+1)/2],[k/2])*(-1)^[(k+1)/2].

Views

Author

Keywords

Comments

Examples

References

Links

Crossrefs

Programs

Haskell

Maple

Mathematica

PARI

Formula

Extensions

A094954 Array T(k,n) read by antidiagonals. G.f.: x(1-x)/(1-kx+x^2), k>1.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Formula

A238379 Expansion of (1 - x)/(1 - 36*x + x^2).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Magma

Mathematica

PARI

Sage

Formula

A165253 Triangle T(n,k), read by rows given by [1,0,1,0,0,0,0,0,0,...] DELTA [0,1,0,0,0,0,0,0,0,...] where DELTA is the operator defined in A084938.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

A160682 The list of the A values in the common solutions to 13*k+1 = A^2 and 17*k+1 = B^2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Magma

Mathematica

PARI

Formula

Extensions

A133607 Triangle read by rows: T(n, k) = qStirling2(n, k, q) for q = -1, with 0 <= k <= n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

Sage

Formula

Extensions

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

A160682 The list of the A values in the common solutions to 13k+1 = A^2 and 17k+1 = B^2.

A161588 The list of the k values in the common solutions to the 2 equations 11k+1=A^2, 15k+1=B^2.