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.

Showing 1-5 of 5 results.

A131074 Triangular array T read by rows: T(j,1) = 1 for ((j-1) mod 8) < 4, else 0; T(j,k) = T(j-1,k-1) + T(j,k-1) for 2 <= k <= j.

Original entry on oeis.org

1, 1, 2, 1, 2, 4, 1, 2, 4, 8, 0, 1, 3, 7, 15, 0, 0, 1, 4, 11, 26, 0, 0, 0, 1, 5, 16, 42, 0, 0, 0, 0, 1, 6, 22, 64, 1, 1, 1, 1, 1, 2, 8, 30, 94, 1, 2, 3, 4, 5, 6, 8, 16, 46, 140, 1, 2, 4, 7, 11, 16, 22, 30, 46, 92, 232, 1, 2, 4, 8, 15, 26, 42, 64, 94, 140, 232, 464, 0, 1, 3, 7, 15, 30, 56, 98
Offset: 1

Views

Author

Klaus Brockhaus, following a suggestion of Paul Curtz, Jun 14 2007

Keywords

Comments

All columns are periodic with period length 8. The (4+8*i)-th row equals the first (4+8*i) terms of the main diagonal (i >= 0). Main diagonal and eighth subdiagonal agree; generally j-th subdiagonal equals (j+8)-th subdiagonal.

Examples

			First seven rows of T are
[ 1 ]
[ 1, 2 ]
[ 1, 2, 4 ]
[ 1, 2, 4, 8 ]
[ 0, 1, 3, 7, 15 ]
[ 0, 0, 1, 4, 11, 26 ]
[ 0, 0, 0, 1, 5, 16, 42 ].

Crossrefs

Cf. A131022, A129961 (main diagonal of T), A131075 (first subdiagonal of T), A131076 (row sums of T), A131077 (antidiagonal sums of T). First through sixth column of T are in A131078, A131079, A131080, A131081, A131082, A131083 resp.

Programs

  • Magma
    m:=13; M:=ZeroMatrix(IntegerRing(), m, m); for j:=1 to m do if (j-1) mod 8 lt 4 then M[j, 1]:=1; end if; end for; for k:=2 to m do for j:=k to m do M[j, k]:=M[j-1, k-1]+M[j, k-1]; end for; end for; &cat[ [ M[j, k]: k in [1..j] ]: j in [1..m] ];
  • Mathematica
    T[j_, 1] := If[Mod[j-1, 8]<4, 1, 0]; T[j_, k_] := T[j, k] = T[j-1, k-1]+T[j, k-1]; Table[T[j, k], {j, 1, 13}, {k, 1, j}] // Flatten (* Jean-François Alcover, Mar 06 2014 *)
  • PARI
    {m=13; M=matrix(m, m); for(j=1, m, M[j, 1]=if((j-1)%8<4, 1, 0)); for(k=2, m, for(j=k, m, M[j, k]=M[j-1, k-1]+M[j, k-1])); for(j=1, m, for(k=1, j, print1(M[j, k], ",")))}

A131075 First subdiagonal of triangular array T: T(j,1) = 1 for ((j-1) mod 8) < 4, else 0; T(j,k) = T(j-1,k-1) + T(j,k-1) for 2 <= k <= j.

Original entry on oeis.org

1, 2, 4, 7, 11, 16, 22, 30, 46, 92, 232, 628, 1652, 4096, 9544, 21000, 43912, 87824, 169120, 315952, 578096, 1048576, 1913440, 3567072, 6874336, 13748672, 28384384, 59797312, 126906176, 268435456, 561834112, 1158971520, 2353246336, 4706492672, 9292452352
Offset: 1

Views

Author

Klaus Brockhaus, following a suggestion of Paul Curtz, Jun 14 2007

Keywords

Comments

Also first differences of main diagonal A129961.

Examples

			For first seven rows of T see A131074 or A129961.

Links

Crossrefs

Cf. A131074 (T read by rows), A129961 (main diagonal of T), A131076 (row sums of T), A131077 (antidiagonal sums of T). First through sixth column of T are in A131078, A131079, A131080, A131081, A131082, A131083 resp.

Programs

  • Magma
    m:=34; M:=ZeroMatrix(IntegerRing(), m, m); for j:=1 to m do if (j-1) mod 8 lt 4 then M[j, 1]:=1; end if; end for; for k:=2 to m do for j:=k to m do M[j, k]:=M[j-1, k-1]+M[j, k-1]; end for; end for; [ M[n+1, n]: n in [1..m-1] ];
  • PARI
    lista(m) = my(v=concat([1, 2, 4, 7, 11], vector(m-5))); for(n=6, m, v[n]=6*v[n-1]-14*v[n-2]+16*v[n-3]-10*v[n-4]+4*v[n-5]); v

Formula

a(1) = 1, a(2) = 2, a(3) = 4, a(4) = 7, a(5) = 11; for n > 5, a(n) = 6*a(n-1)-14*a(n-2)+16*a(n-3)-10*a(n-4)+4*a(n-5).
G.f.: x*(1-4*x+6*x^2-5*x^3+3*x^4)/((1-2*x)*(1-4*x+6*x^2-4*x^3+2*x^4)).

A131076 Row sums of triangular array T: T(j,1) = 1 for ((j-1) mod 8) < 4, else 0; T(j,k) = T(j-1,k-1) + T(j,k-1) for 2 <= k <= j.

Original entry on oeis.org

1, 3, 7, 15, 26, 42, 64, 93, 139, 231, 463, 1092, 2744, 6840, 16384, 37383, 81295, 169119, 338239, 654192, 1232288, 2280864, 4194304, 7761375, 14635711, 28384383, 56768767, 116566080, 243472256, 511907712, 1073741824, 2232713343, 4585959679, 9292452351
Offset: 1

Views

Author

Klaus Brockhaus, following a suggestion of Paul Curtz, Jun 14 2007

Keywords

Comments

Sum of n-th row equals (n+1)-th term of main diagonal minus (n+1)-th term of first column: a(n) = A129961(n+1) - A131078(n+1).

Examples

			For first seven rows of T see A131074 or A129961.

Links

Crossrefs

Cf. A131074 (T read by rows), A129961 (main diagonal of T), A131075 (first subdiagonal of T), A131077 (antidiagonal sums of T). First through sixth column of T are in A131078, A131079, A131080, A131081, A131082, A131083 resp.

Programs

  • Magma
    m:=32; M:=ZeroMatrix(IntegerRing(), m, m); for j:=1 to m do if (j-1) mod 8 lt 4 then M[j, 1]:=1; end if; end for; for k:=2 to m do for j:=k to m do M[j, k]:=M[j-1, k-1]+M[j, k-1]; end for; end for; [ &+[ M[j, k]: k in [1..j] ]: j in [1..m] ];
  • Mathematica
    LinearRecurrence[{7,-20,30,-27,21,-24,30,-26,14,-4},{1,3,7,15,26,42,64,93,139,231},40] (* Harvey P. Dale, Jun 23 2025 *)
  • PARI
    lista(m) = my(M=matrix(m, m)); for(j=1, m, M[j, 1]=if((j-1)%8<4, 1, 0)); for(k=2, m, for(j=k, m, M[j, k]=M[j-1, k-1]+M[j, k-1])); for(j=1, m, print1(sum(k=1, j, M[j, k]), ", "))

Formula

G.f.: x*(1-4*x+6*x^2-4*x^3-2*x^4+10*x^5-10*x^6+5*x^7-x^8)/((1-x)*(1-2*x)*(1+x^4)*(1-4*x+6*x^2-4*x^3+2*x^4)).

A131077 Antidiagonal sums of triangular array T: T(j,1) = 1 for ((j-1) mod 8) < 4, else 0; T(j,k) = T(j-1,k-1) + T(j,k-1) for 2 <= k <= j.

Original entry on oeis.org

1, 1, 3, 3, 6, 5, 11, 8, 20, 14, 35, 24, 59, 41, 100, 77, 178, 162, 341, 364, 705, 837, 1542, 1915, 3458, 4282, 7741, 9280, 17021, 19461, 36482, 39559, 76042, 78218, 154261, 151184, 305445, 287509, 592954, 542223, 1135178, 1023210, 2158389, 1949312, 4107701
Offset: 1

Views

Author

Klaus Brockhaus, following a suggestion of Paul Curtz, Jun 14 2007

Keywords

Examples

			For first seven rows of T see A131074 or A129961.

Links

Crossrefs

Cf. A131074 (T read by rows), A129961 (main diagonal of T), A131075 (first subdiagonal of T), A131076 (row sums of T). First through sixth column of T are in A131078, A131079, A131080, A131081, A131082, A131083 resp.

Programs

  • Magma
    m:=44; M:=ZeroMatrix(IntegerRing(), m, m); for j:=1 to m do if (j-1) mod 8 lt 4 then M[j, 1]:=1; end if; end for; for k:=2 to m do for j:=k to m do M[j, k]:=M[j-1, k-1]+M[j, k-1]; end for; end for; [ &+[ M[j-k+1, k]: k in [1..(j+1) div 2] ]: j in [1..m] ];
  • PARI
    lista(m) = my(M=matrix(m, m)); for(j=1, m, M[j, 1]=if((j-1)%8<4, 1, 0)); for(k=2, m, for(j=k, m, M[j, k]=M[j-1, k-1]+M[j, k-1])); for(j=1, m, print1(sum(k=1, (j+1)\2, M[j-k+1, k]), ", "))

Formula

G.f.: x*(1-4*x^2+6*x^4-x^5-4*x^6+3*x^7+x^8-3*x^9+x^10+2*x^11-x^12) / ((1-x)*(1-2*x^2)*(1+x^4)*(1-4*x^2+6*x^4-4*x^6+2*x^8)).

A130668 Diagonal of A129819.

Original entry on oeis.org

0, 0, 1, -2, 5, -11, 23, -48, 102, -220, 476, -1024, 2184, -4624, 9744, -20480, 42976, -90048, 188352, -393216, 819328, -1704192, 3539200, -7340032, 15203840, -31456256, 65010688, -134217728, 276826112, -570429440, 1174409216
Offset: 0

Views

Author

Paul Curtz, Jun 27 2007

Keywords

Comments

This sequence is connected to A124072. To see this, change the sign of every negative term and consider the differences of every line. Hence for the second line, and following lines, the four terms form periodic sequences:
0 1 0 1 0
1 0 0 1 1
1 0 1 2 1
1 1 3 3 1
2 4 6 4 2
6 10 10 6 6
16 20 16 12 16
36 36 28 28 36
72 64 56 64 72
136 120 120 136 136
256 240 256 272 256.
The lines are connected as seen by the examples: (3rd line connected to 2nd, from right to left) 1+1=2, 1+0=1, 0+0=0, 0+1=1; (11th line connected to 10th) 136+136=272, 136+120=256, 120+120=240, 120+136=256.
The 4 columns are almost known (must the first line be suppressed?): A038503 (without the first 1), A000749 (without the first 0), A038505, A038504. Like the present sequence, every sequence of A124072 beginning with a negative number (-2, -11, ...) is a "twisted" sequence (see A129339 comments, A129961 and the present 4 columns). Periodic with period 2^n.
Inverse binomial transform of A129819. - R. J. Mathar, Feb 25 2009

Links

Programs

  • GAP
    a:=[-2,5,-11,23];; for n in [5..30] do a[n]:=-6*a[n-1]+-14*a[n-2] -16*a[n-3]-8*a[n-4]; od; Concatenation([0,0,1], a); # G. C. Greubel, Mar 24 2019
  • Magma
    I:=[-2,5,-11,23]; [0,0,1] cat [n le 4 select I[n] else -6*Self(n-1) - 14*Self(n-2)-16*Self(n-3)-8*Self(n-4): n in [1..30]]; // G. C. Greubel, Mar 24 2019
  • Mathematica
    gf = x^2*(1+x)*(1+3*x+4*x^2+3*x^3)/((1+2*x+2*x^2)*(1+2*x)^2); CoefficientList[Series[gf, {x, 0, 30}], x] (* Jean-François Alcover, Dec 16 2014, after R. J. Mathar *)
Join[{0, 0, 1}, LinearRecurrence[{-6,-14,-16,-8}, {-2,5,-11,23}, 30]] (* Jean-François Alcover, Feb 15 2016 *)
  • PARI
    my(x='x+O('x^30)); concat([0,0], Vec(x^2*(1+x)*(1+3*x+4*x^2+3*x^3 )/((1+2*x +2*x^2)*(1+2*x)^2))) \\ G. C. Greubel, Mar 24 2019
  • Sage
    (x^2*(1+x)*(1+3*x+4*x^2+3*x^3)/((1+2*x+2*x^2)*(1+2*x)^2 )).series(x, 30).coefficients(x, sparse=False) # G. C. Greubel, Mar 24 2019

Formula

From R. J. Mathar, Feb 25 2009: (Start)
G.f.: x^2*(1+x)*(1 + 3*x + 4*x^2 + 3*x^3)/((1 + 2*x + 2*x^2)*(1+2*x)^2).
a(n) = ((-1)^n*A001787(n+1) - 4*A108520(n) + 4*A122803(n))/32, n > 2. (End)
a(n) = -6*a(n-1) - 14*a(n-2) - 16*a(n-3) - 8*a(n-4) for n >= 7. - G. C. Greubel, Mar 24 2019

Extensions

Extended by R. J. Mathar, Feb 25 2009
Showing 1-5 of 5 results.

Started in 1964 by Neil J. A. Sloane | Maintained by The OEIS Foundation Inc.

Content is available under The OEIS End-User License Agreement.