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.

User: Gary W. Adamson

Gary W. Adamson's wiki page.

Gary W. Adamson has authored 2826 sequences. Here are the ten most recent ones:

A362455 Squarefree positive integers d such that the dimension of the space of cuspidal harmonic automorphic forms for SL(2, O_{-d}) is zero, where O_{-d} is the ring of integers in Q(sqrt(-d)).

Original entry on oeis.org

1, 2, 3, 5, 6, 7, 11, 15, 19, 23, 31, 39, 47, 71
Offset: 1

Author

Gary W. Adamson and N. J. A. Sloane, May 06 2023, following a suggestion from Roger L. Bagula

Keywords

Comments

Reid, Alan W., and Colin MacLachlan. "The Arithmetic of hyperbolic 3-manifolds." Graduate Texts in Mathematics 219, Springer (2003).

A358549 Triangle read by rows where row n is reversed partial sums of row n of the Sierpinski triangle (A047999).

Original entry on oeis.org

1, 2, 1, 2, 1, 1, 4, 3, 2, 1, 2, 1, 1, 1, 1, 4, 3, 2, 2, 2, 1, 4, 3, 3, 2, 2, 1, 1, 8, 7, 6, 5, 4, 3, 2, 1, 2, 1, 1, 1, 1, 1, 1, 1, 1, 4, 3, 2, 2, 2, 2, 2, 2, 2, 1, 4, 3, 3, 2, 2, 2, 2, 2, 2, 1, 1, 8, 7, 6, 5, 4, 4, 4, 4, 4, 3, 2, 1, 4, 3, 3, 3, 3, 2, 2, 2, 2, 1, 1, 1, 1
Offset: 0

Author

Gary W. Adamson, Nov 21 2022

Keywords

Comments

Row reversal of A261363 (which is the main entry).
These sums can be formed by taking A047999 as a lower triangular matrix times an all-1's lower triangular matrix.

Examples

			Triangle begins:
      k=0  1  2  3  4  5  6  7  8
  n=0:  1;
  n=1:  2, 1;
  n=2:  2, 1, 1;
  n=3:  4, 3, 2, 1;
  n=4:  2, 1, 1, 1, 1;
  n=5:  4, 3, 2, 2, 2, 1;
  n=6:  4, 3, 3, 2, 2, 1, 1;
  n=7:  8, 7, 6, 5, 4, 3, 2, 1;
  n=8:  2, 1, 1, 1, 1, 1, 1, 1, 1;
For n=5, row 5 here and row 5 of A047999 are:
  row      4, 3, 2, 2, 2, 1
  sums of  1, 1, 0, 0, 1, 1
		

Crossrefs

Cf. A047999, A261363 (rows reversed).
Cf. A001316 (column k=0), A000012 (main diagonal).

Programs

  • Mathematica
    row[n_] := Reverse[Accumulate[Array[Boole[0 == BitAnd[n-#, #]] &, n + 1, 0]]]; Array[row, 13, 0] // Flatten (* Amiram Eldar, May 13 2025 *)

Formula

T(n,k) = Sum_{i=k..n} A047999(n,i).

A347029 a(n) = ceiling(e^(n*(Pi/2))).

Original entry on oeis.org

1, 5, 24, 112, 536, 2576, 12392, 59610, 286752, 1379411, 6635624, 31920520, 153552936, 738662923, 3553321281, 17093171649, 82226315586, 395547831245, 1902773895293, 9153250784395, 44031505860633, 211812562992414, 1018919543279305, 4901489415968643, 23578503968558227
Offset: 0

Author

Gary W. Adamson, Aug 11 2021

Keywords

Comments

Alternative formula for e^(n*Pi/2) is i^(-n*i), where i = sqrt(-1). Substitute 2i for n in each identity, resulting in e^(Pi*i) = -1; Euler's formula.
A121905 is the bisection of the sequence, ceiling(e^(n*Pi)).

Examples

			a(5) = ceiling(e^(5*Pi/2)) = ceiling(i^(-5*i)) = 2576.
		

Crossrefs

Cf. A121905 (even bisection), A124507 (floor), A042972.

Programs

  • Mathematica
    a[n_]:=Ceiling[Exp[n Pi/2]]; Table[a[n],{n,0,24}] (* Stefano Spezia, Aug 12 2021 *)
  • PARI
    a(n) = ceil(exp(n*Pi/2)); \\ Michel Marcus, Aug 12 2021

Formula

a(n) = ceiling(e^(n*Pi/2)) = ceiling(i^(-n*i)).

A345226 Irregular triangle, row sums equal A000041; in the format of A233932.

Original entry on oeis.org

1, 1, 1, 2, 1, 2, 1, 2, 4, 1, 2, 4, 5, 2, 8, 5, 2, 8, 5, 2, 7, 16, 5, 2, 7, 16, 17, 2, 7, 30, 17, 2, 7, 30, 17, 23, 7, 54, 17, 23, 7, 54, 51, 23, 7, 95, 51, 23, 7, 95, 51, 23, 7, 55, 161, 51, 23, 7, 55, 161, 139, 23, 7, 55, 266, 139, 23, 7, 55, 266, 139, 160, 7, 55, 431, 139, 160, 7, 55
Offset: 1

Author

Gary W. Adamson, Jun 12 2021

Keywords

Comments

The format of A233932 has a Gray code property of one term change in the next row. Using the production matrix shown below, we can obtain an array with row sums of any target sequence.

Examples

			The first few rows of the array equal A000041 with offset 1 as to sums:
    1;
    1,   1;
    2,   1;
    2,   1,  2;
    4,   1,  2;
    4,   5,  2;
    8,   5,  2;
    8,   5,  2, 7;
   16,   5,  2, 7;
   16,  17,  2, 7;
   30,  17,  2, 7;
   30,  17, 23, 7;
   54,  17, 23, 7;
   54,  51, 23, 7;
   95,  51, 23, 7;
   95,  51, 23, 7, 55;
  161,  51, 23, 7, 55;
  161, 139, 23, 7, 55;
  ...
The leftmost column is (1, 1, 2, 2, 4, 4, 8, 8, ...), being the partial sums of the first column in matrix M:  (1, 0, 1, 0, 2, 0, 4, ...).
		

Formula

Let P equal the infinite lower triangular matrix with 1's in every row: [(1), (1, 1), (1, 1, 1), ...]. Begin with the following matrix format such that M[n, A001511(k)] = 1, otherwise 0:
1
0, 1
1
0, 0, 1
1
0, 1
1
0, 0, 0, 1
...
Replace the 1's with A187219 (the first difference row of A000041), getting M:
1
0, 1
1
0, 0, 2
2
0, 4
4
0, 0, 0, 7
...
Perform the operation P * M, getting A345226 as an irregular matrix. The operation P * M is equivalent to taking partial sums of column terms from top to bottom.

A343233 Triangle read by rows: Riordan triangle T = (1 - x*c(x), x), with the generating function c of A000108 (Catalan).

Original entry on oeis.org

1, -1, 1, -1, -1, 1, -2, -1, -1, 1, -5, -2, -1, -1, 1, -14, -5, -2, -1, -1, 1, -42, -14, -5, -2, -1, -1, 1, -132, -42, -14, -5, -2, -1, -1, 1, -429, -132, -42, -14, -5, -2, -1, -1, 1, -1430, -429, -132, -42, -14, -5, -2, -1, -1, 1
Offset: 0

Author

Gary W. Adamson and Wolfdieter Lang, Apr 12 2021

Keywords

Comments

As an unsigned sequence a(n) this is identical with the one of A155586(n+1), for n >= 0, but the triangle is not a simple signed version of A155586. See the formula.
This lower triangular Riordan matrix T of Toeplitz type is the inverse of the Riordan matrix (c(x), x) = |A106270|, also of Toeplitz type.

Examples

			The triangle matrix T begins:
  n/m     0    1    2   3   4   5   6   7   8   9 ...
  --------------------------------------------------
  0:      1
  1:     -1    1
  2:     -1   -1    1
  3:     -2   -1   -1   1
  4:     -5   -2   -1  -1   1
  5:    -14   -5   -2  -1  -1   1
  6:    -42  -14   -5  -2  -1  -1   1
  7:   -132  -42  -14  -5  -2  -1  -1   1
  8:   -429 -132  -42 -14  -5  -2  -1  -1   1
  9:  -1430 -429 -132 -42 -14  -5  -2  -1  -1   1
  ...
		

Crossrefs

Cf. A106270 (unsigned), A155586.

Formula

The lower triangular matrix T satisfies: T = I - L^{tr}*|A106270|, also for the finite N X N version, with the unit matrix I and the lower triangular matrix L^{tr}(i, j) = delta_{i, j-1} (Kronecker symbol delta) with first lower diagonal of 1s and 0 otherwise.
T(n, n) = 1, and for T(n, m) = -C_{n - 1 - m } = - |A106270(n-1, m)|, for 0 <= m <= n-1, with the Catalan numbers C(n) = A000108, and T(n, m) = 0 for n < m.
O.g.f. of column m: (1/c(x))*x^m = (1 - x*c(x))*x^m (Riordan matrix of Toeplitz type), with the o.g.f. c of A000108.
O.g.f. row polynomials R(n, x) = Sum_{m=0..n} T(n, m)*x^m, that is the o.g.f. of the triangle. G(z, x) = c(z)/(1 - x*z).

A339050 Triangle read by rows T(n, m) = F(2*m-1)*(n-m) + F(2*m), for 1 <= m <= n, where F = A000045 (Fibonacci).

Original entry on oeis.org

1, 2, 3, 3, 5, 8, 4, 7, 13, 21, 5, 9, 18, 34, 55, 6, 11, 23, 47, 89, 144, 7, 13, 28, 60, 123, 233, 377, 8, 15, 33, 73, 157, 322, 610, 987, 9, 17, 38, 86, 191, 411, 843, 1597, 2584, 10, 19, 43, 99, 225, 500, 1076, 2207, 4181, 6765
Offset: 1

Author

Gary W. Adamson and Wolfdieter Lang, Jan 15 2021

Keywords

Comments

This is the partial sum triangle of triangle A143929.
The main diagonal is the INVERT transform of the first column (offset 1 in both sequences).

Examples

			The triangle T(n, m) begins:
n\m   1  2  3  4   5   6    7    8    9   10 ...
1:    1
2:    2  3
3:    3  5  8
4:    4  7 13 21
5:    5  9 18 34  55
6:    6 11 23 47  89 144
7:    7 13 28 60 123 233  377
8:    8 15 33 73 157 322  610  987
9:    9 17 38 86 191 411  843 1597 2584
10:  10 19 43 99 225 500 1076 2207 4181 6765
...
		

Crossrefs

The first columns (without leading zeros) are A001477(n), A005408(n+1), A005408(n+1), for n >= 1.
The first (sub)diagonals are A001906(m), A001519(m+1), A005248(m), for m >= 1.

Formula

T(n, m) = Sum_{k=1..m} A143929(n, k), n >=1, m = 1, 2, ..., n, otherwise 0.
T(n, m) = A(m)*n + B(m), with A(m) = A(m-1) + F(2*(m-1)), for m >= 2 and A(1) = 1, and B(m) = B(m-1) + (m-1)*F(2*(m-1)), for m >= 2 and B(1) = 0, where F(2*m) =A001906(m) and F(2*m-1) = A001519(m).
T(n, 1) = n, for n >= 1; T(n, m) = F(2*(m-1))*(n-m+1), if m >= 2 and n >= m, and 0 otherwise.
G.f. of column m: G(m,x) = x^m*(x*F(2*m-1)/(1-x)^2 + F(2*m)/(1-x)), for m >= 1.
G.f. of row polynomials R(n, x) := Sum{m=1..n} T(n, m)*x^m, that is g.f. of the triangle: G(z,x) = (x*z)*(1 - x*z^2)/((1- 3*x*z + (x*z)^2)*(1 - z)^2).
G.f. of (sub)diagonal k: D(k,x) = x*((k-1)*(1-x) + 1)/(1 - 3*x + x^2), for k >= 1.

A337712 Irregular triangle read by rows: row n gives the complete system of cycles of the doubling sequences modulo N = 2*n+1, for n >= 0.

Original entry on oeis.org

1, 2, 1, 2, 4, 3, 1, 2, 4, 3, 6, 5, 1, 2, 4, 8, 7, 5, 1, 2, 4, 8, 5, 10, 9, 7, 3, 6, 1, 2, 4, 8, 3, 6, 12, 11, 9, 5, 10, 7, 1, 2, 4, 8, 7, 14, 13, 11, 1, 2, 4, 8, 16, 15, 13, 9, 3, 6, 12, 7, 14, 11, 5, 10, 1, 2, 4, 8, 16, 13, 7, 14, 9, 18, 17, 15, 11, 3, 6, 12, 5, 10
Offset: 0

Author

Gary W. Adamson and Wolfdieter Lang, Oct 14 2020

Keywords

Comments

The length of row n is A037225(n), for n >= 0.
The doubling sequence modulo N = 2*n+1, for n >= 0, has entries DS(N, s(N,i), j) = s(N,i)*2^j (mod N), with j >= 0, and certain positive odd integer seeds s(N, i), for i = 1, 2, ..., S(N) = A037226((N-1)/2), where gcd(s(N, i), N) = 1 (restricted seeds modulo N). These doubling sequences are periodic with period length P(N) = A002326((N-1)/2) (order of 2 modulo N). Only the periods (cycles) {DS(N, s(N, i), j)}_{j=0..P(N)-1}, for i = 1, 2, ..., S(N), are listed.
N = 1 (n=0) is special: one takes here the restricted residue system modulo N not as [0] but as [1]. The order of 2 modulo 1 is 1, because 2^1 == 1 (mod 1) (== 0 (mod 1)).
In order to obtain the complete system of doubling sequences one starts with seed s(N, 1) = 1, and if all numbers from the smallest positive reduced residue system modulo N (called RRS(N), given in row N of A038566) are obtained, i.e., if P(N) = #RRS(N) = phi(N) = A000010(N), then the system is complete. Otherwise the smallest missing number from RRS(N) is taken as new seed s(N, 2), etc. until the system is complete. This means that the number of seeds needed is S(N) = phi(N)/P(N) = A037226((N-1)/2)).
The irregular subtriangle where only seed s(N, 1) = 1 has been used is given in A201908. But there 0 (not 1) for N = 1 has been used.
From Gary W. Adamson and Wolfdieter Lang, Dec 15 2020: (Start)
The cycles in row n, for N = 2*n + 1, of period length P(N) = A002326((N-1)/2) give the periods of the iterated doubling function D(x) = frac(2*x) with seeds x = s(N, i)/N, for i = 1, 2, ..., S(N) = A037226((N-1)/2), after multiplication with N. This is the doubling function used in the Devaney reference, pp. 24-25, 27, 125. 132, 171,289.
Each cycle in row n can also be used to find from the base 2 version of its first entry (the seed s = s(N, i)) divided by N the other entries by repeated application of a cyclic left shift by one step (called sigma operation) to the period of the base 2 expression of s/N. E.g., n = 7, N = 15, P(N) = 4, s = 1: (1/15){10->2} = .repeat(0001), then (.repeat(0010)){2->10} = 2/10, (.repeat(0100)){2->10} = 4/10 and (.repeat(1000)){2->10} = 8/15. Similarly for s = 7: from (7/15)_{10->2} = .repeat(0111) one obtains by repeated sigma operations 14/15, 13/15 and 11/15. The proof uses the elementary formulas for the conversion from base 10 to base 2, and the reverse one, from base 2 to base 10. See also a comment on the period length P(N) given in A002326. (End)

Examples

			The irregular triangle T(n, k) begins (cycles are separated by a vertical bar)
n,  N \ k 1 2 3 4  5  6  7  8  9 10 11 12 13 14 15 16 17 18 19 20 21 22 ...
0,  1:    1
1,  3:    1 2
2,  5:    1 2 4 3
3,  7:    1 2 4|3  6  5
4,  9:    1 2 4 8  7  5
5,  11:   1 2 4 8  5 10  9  7  3  6
6,  13:   1 2 4 8  3  6 12 11  9  5 10  7
7,  15:   1 2 4 8| 7 14 13 11
8,  17:   1 2 4 8 16 15 13  9| 7 14 11  5 10  3  6 12
9,  19:   1 2 4 8 16 13  7 14  9 18 17 15 11  3  6 12 5 10
10, 21:   1 2 4 8 16 11| 5 10 20 19 17 13
11, 23:   1 2 4 8 16  9 18 13  3  6 12| 5 10 20 17 11 22 21 19 15  7 14
12, 25:   1 2 4 8 16  7 14  3  6 12 24 23 21 17  9 18 11 22 19 13
13, 27:   1 2 4 8 16  5 10 20 13 26 25 23 19 11 22 17  7 14
...
n = 14, N = 29:  1 2 4 8 16  3  6 12 24 19  9 18  7 14 28 27 25 21 13 26 23 17  5 10 20 11 22 15,
n = 15, N = 31: 1 2 4 8 16|3 6 12 24 17|5 10 20 9 18|7 14 28 25 19|11 22 13 26 21|15 30 29 27 23.
		

References

  • Robert L. Devaney, A First Course in Chaotic Dynamical Systems, Addison-Wesley., 1992. pp. 24-25, 27, 125, 132, 171, 289. Second edition 2020.

Crossrefs

Cf. A000010, A002326, A037225, A037226, A201908, A038566, A334430 (modified doubling), A337936 (tripling), A339046 (quadrupling).

Programs

  • Mathematica
    Array[Block[{a = {}, k = 2, n = 2 # + 1, m}, m = EulerPhi[n]; While[Length@ Flatten@ a < m, AppendTo[a, Most@ NestWhileList[Mod[2 #, n] &, If[Length@ a == 0, 1, k], UnsameQ, All]]; Set[k, SelectFirst[Complement[Range[n], Union@ Flatten@ a], GCD[#, n] == 1 &] ]]; a] &, 9] // Flatten (* Michael De Vlieger, Nov 06 2020 *)

Formula

T(n, k) gives the k-th entry in the complete doubling system modulo N = 2*n+1, for n >= 0, with the S(N) = A037226((N-1)/2) cycles of length A002326((N-1)/2) written in row n. See the comment above for DS(N,s(N,i)), i = 1, 2, ..., S(N).

A334430 Irregular triangle read by rows. Row n gives the cycles of positive integers of the complete modified doubling sequence MDS(b) for b = 2*n + 1, for n >= 1.

Original entry on oeis.org

1, 2, 1, 2, 3, 1, 2, 4, 1, 2, 4, 3, 5, 1, 2, 4, 5, 3, 6, 1, 2, 4, 7, 1, 2, 4, 8, 1, 6, 5, 7, 3, 2, 4, 8, 3, 6, 7, 5, 9, 1, 2, 4, 8, 5, 10, 1, 2, 4, 8, 7, 9, 5, 10, 3, 6, 11, 1, 2, 4, 8, 9, 7, 11, 3, 6, 12, 1, 2, 4, 8, 11, 5, 10, 7, 13, 1, 2, 4, 8, 13, 3, 6, 12, 5, 10, 9, 11, 7, 14, 1
Offset: 1

Author

Gary W. Adamson and Wolfdieter Lang, Aug 11 2020

Keywords

Comments

The length of row n is phi(2*n+1)/2 with phi = A000010, for n > = 1.
There are c(2*n+1) = A135303(n) cycles in row n, and the length of each cycle is k(2*n+1) = A003558(n). c(2*n+1)*k(2*n+1) = phi(n)/2 (in Hilton and Pedersen [HP], the coach theorem, p. 262).
In the construction of the cycles a modified modular equivalence relation, called mod* b, for b = 2*n + 1, with n >= 1, proposed in the Brändli and Beyne [BB] paper, is used. It is defined as mod*(a, b) = mod(a, b) (with values in I(b) = {0, 1, ..., b-1}) for integer numbers a, if mod(a, b) <= (b-1)/2 and mod*(a, b) = mod(-a, b) if mod(a, b) > (b-1)/2. The equivalence relation k ~ m if mod*(k, b) = mod*(m,b) is multiplicative (not additive).
The MDS(b) cycles are obtained from the doubling sequence {a(b,i)*2^j}A003558(n)%20for%20each%20odd%20a(b,%20i)%20with%20gcd(a(b,%20i),%20b)%20=%201%20and%20a(b,%20i)%20%3C=%20(b-1)/2.%20The%20first%20cycle%20is%20Cy*(b,1)%20obtained%20for%20a(b,%201)%20=%201.%20If%20not%20all%20odd%20numbers%20relatively%20prime%20to%20b%20and%20%3C=%20(b-1)/2%20are%20present,%20the%20next%20cycle%20Cy*(b,2)%20uses%20the%20smallest%20missing%20reduced%20odd%20number,%20etc.%20The%20number%20of%20such%20inputs%20a(b,%20i),%20hence%20cycles,%20is%20c*(b),%20coinciding%20with%20c(2*n+1)%20=%20A135303(n),%20and%20all%20odd%20and%20even%20positive%20numbers%20with%20gcd(m,%20b)%20=%201%20and%20m%20%3C=%20(b-1)/2%20appear%20just%20once%20in%20the%20complete%20cycle%20system%20for%20b,%20called%20Cy*(b)%20=%20%7BCy*(b,i)%7D">{j=1..P(b)} evaluated mod*b: Cy*(b,i)_j = mod*(a(b,i)*2^j, b), with certain inputs (seeds) a(b, i). P(b), the length of the period independent, coincides with k(2*n+1) = A003558(n) for each odd a(b, i) with gcd(a(b, i), b) = 1 and a(b, i) <= (b-1)/2. The first cycle is Cy*(b,1) obtained for a(b, 1) = 1. If not all odd numbers relatively prime to b and <= (b-1)/2 are present, the next cycle Cy*(b,2) uses the smallest missing reduced odd number, etc. The number of such inputs a(b, i), hence cycles, is c*(b), coinciding with c(2*n+1) = A135303(n), and all odd and even positive numbers with gcd(m, b) = 1 and m <= (b-1)/2 appear just once in the complete cycle system for b, called Cy*(b) = {Cy*(b,i)}{i=1..c(b)}.
Such modified doubling sequences have been considered in the Kappraff-Adamson paper using iterations of x^2 - 2, and also in comments and examples by Adamson like the Aug 25 2019 comment in A065941, where it is named "r-t table" (for roots trajectory).
The recurrence relation for the cycle Cy*(b,i) with elements {c_j(b,i)} is
The complete cycle system Cy*(b) is equivalent to the complete coach system of Hilton and Pedersen [HP], with the coach number c(b) = c*(b) and k(b) = P(b). The odd numbers of each cycle Cy*(b,i), read backwards (anticyclically), give the first row of the coach; and the number of steps to go from one odd number to the next odd number, ending with the starting number, gives the second row of the coach with sum k(b). See [HP] pp. 102-103 for the modified coach system and the extended list, eq. (7. 12), in the proof of their quasi-order theorem.
The complete cycle system of Schick (unsigned) and Brändli-Beyne SBB(b) is also equivalent to Cy*(b). Each cycle Cy*(b,i) gives the SBB(b,i) cycle (q(b, i)j) with elements obtained from q(b, i)_j = b - 2*c(b,i){P(b)-1+j}, for j = 0, 1, 2, ..., P(b)-1. Cyclicity is used: c(b,i){P(b)+k} = c(b,i){k}.
For more details and proofs see the arxiv paper by W. Lang.

Examples

			The irregular triangle a(n,m) begins (cycles are enclosed by brackets; blanks are used to fit the index m):
n,   b \ m 1 2 3  4  5   6  7  8  9 10  11 12 13 14 15 16 17 18 19 20 ...
1,   3:   (1)
2,   5:   (2 1)
3,   7:   (2 3 1)
4,   9:   (2 4 1)
5,  11:   (2 4 3  5  1)
6,  13:   (2 4 5  3  6   1)
7,  15:   (2 4 7  1)
8,  17:   (2 4 8  1)(6   5  7  3)
9,  19:   (2 4 8  3  6   7  5  9  1)
10, 21:   (2 4 8  5 10   1)
11, 23:   (2 4 8  7  9   5 10  3  6 11  1)
12, 25:   (2 4 8  9  7  11  3  6 12  1)
13, 27:   (2 4 8 11  5  10  7 13  1)
14, 29:   (2 4 8 13  3   6 12  5 10  9  11  7 14  1)
15, 31:   (2 4 8 15  1) (6 12  7 14  3)(10 11  9 13  5)
16, 33:   (2 4 8 16  1)(10 13  7 14  5)
17, 35:   (2 4 8 16  3   6 12 11 13  9  17  1)
18, 37:   (2 4 8 16  5  10 17  3  6 12  13 11 15  7 14  9 18  1)
19, 39:   (2 4 8 16  7  14 11 17  5 10  19  1)
20, 41:   (2 4 8 16  9  18  5 10 20  1) (6 12 17  7 14 13 15 11 19  3)
...
-------------------------------------------------------------------------
n = 8, b = 17: (2, 4, 8, 16) = (2, 4, 8, 16) (mod(17)) = (2, 4, 8, 1) (mod*(17), and the first odd (reduced) number missing is 3, with cycle (3*2, 3*4, 3*8, 3*16) = (6, 12, 7, 14) (mod(17)) = (6, 5, 7, 3) (mod*(17).
From this the complete coach system Sigma(17) becomes {[[1], [4]]; [[3 7 5], [1 1 2]]]}.
The complete SBB(17) cycle systen is (unsigned, offset 0) {(1 15 13 9); (3 11 5 7)}, from  17 - 2*8 = 1, 17 - 2*1 = 15, 17 - 2*2 = 13,  17 - 2*4 = 9, and from  17 - 2*7 = 3, 17 - 2*3 = 11, 17 - 2*6 = 5, 17 - 2*5 = 7.
		

References

  • Peter Hilton and Jean Pederson, A Mathematical Tapestry, 2010 (3rd printing 2012), Cambridge University Press, pp. 102-103, 260-264.
  • Carl Schick, Trigonometrie und unterhaltsame Zahlentheorie, Bokos Druck, Zürich, 2003 (ISBN 3-9522917-0-6). Tables 3.1 to 3.10, for odd p = 3 to 113 (with some gaps), pp. 158-166.

Formula

See the above given recurrence relation for the cycles Cy*(b, i), i = 1, 2, ..., c(b) = A135303((b-1)/2), for b = 2*n + 1 >= 3, and the procedure to choose the inputs a(b, i).

A332602 Tridiagonal matrix M read by antidiagonals: main diagonal is 1,2,2,2,2,..., two adjacent diagonals are 1,1,1,1,1,...

Original entry on oeis.org

1, 1, 1, 0, 2, 0, 0, 1, 1, 0, 0, 0, 2, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 2, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0, 0, 0, 0, 0
Offset: 0

Author

N. J. A. Sloane, Mar 06 2020, following a suggestion from Gary W. Adamson

Keywords

Comments

From Gary W. Adamson, Mar 11 2020: (Start)
The upper left entry of M^n gives the Catalan numbers A000108. Extracting 2 X 2, 3 X 3, and 4 X 4 submatrices from M; then generating sequences from the upper left entries of M^n, we obtain the following sequences:
1, 1, 2, 5, 13, ... = A001519 and the convergent is 2.61803... = 2 + 2*cos(2*Pi/5) = (2*cos(Pi/5))^2.
1, 1, 2, 5, 14, 42, 131, ... = A080937 and the convergent is 3.24697... = 2 + 2*cos(2*Pi/7) = (2*cos(Pi/7))^2.
1, 1, 2, 5, 14, 42, 132, 429, 1429, ... = A080938 and the convergent is 3.53208... = 2 + 2*cos(2*Pi/9) = (2*cos(Pi/9))^2. (End)
The characteristic polynomial for the N X N main submatrix M_N is Phi(N, x) = S(N, 2-x) - S(N-1, 2-x), with Chebyshev's S polynomial (see A049310) evaluated at 2-x. Proof by determinant expansion, to obtain the recurrence Phi(N, x) - (x-2)*Phi(N-1, x) - Phi(N-2, x), for N >= 2, and Phi(0, x) = 1 and Phi(1, x) = 1 - x, that is Phi(-1, x) = 1. The trace is tr(M_N) = 1 + 2^(N-1) = A000051(N-1), and Det(M_N) = 1. - Wolfdieter Lang, Mar 13 2020
The explicit form of the characteristic polynomial for the N X N main submatrix M_N is Phi(N, x) := Det(M_N - x*1_N) = Sum_{k=0..N} binomial(N+k, 2*k)*(-x)^k = Sum_{k=0..N} A085478(N, k)*(-x)^k, for N >= 0, with Phi(0, x) := 1. Proof from the recurrence given in the preceding comment. - Wolfdieter Lang, Mar 25 2020
For the proofs of the 2 X 2, 3 X 3 and 4 X 4 conjectures, see the comments in the respective A-numbers A001519, A080937 and A080938. - Wolfdieter Lang, Mar 30 2020
Replace the main diagonal 1,2,2,2,... of the matrix M with 1,0,0,0,...; 1,1,1,1,...; 1,3,3,3,...; 1,2,1,2,...; 1,2,3,4,...; 1,0,1,0...; and 1,1,0,0,1,1,0,0,.... Take powers of M and extract the upper left terms, resulting in respectively: A001405, A001006, A033321, A176677, A006789, A090344, and A007902. - Gary W. Adamson, Apr 12 2022
The statement that the upper left entry of M^n is a Catalan number is equivalent to Exercise 41 of R. Stanley, "Catalan Numbers." - Richard Stanley, Feb 28 2023
If the upper left 1 in matrix M is replaced with 3, taking powers of the resulting matrix and extracting the upper left terms apparently results in sequence A001700. - Gary W. Adamson, Apr 03 2023

Examples

			The matrix begins:
  1, 1, 0, 0, 0, ...
  1, 2, 1, 0, 0, ...
  0, 1, 2, 1, 0, ...
  0, 0, 1, 2, 1, ...
  0, 0, 0, 1, 2, ...
  ...
The first few antidiagonals are:
  1;
  1, 1;
  0, 2, 0;
  0, 1, 1, 0;
  0, 0, 2, 0, 0;
  0, 0, 1, 1, 0, 0;
  0, 0, 0, 2, 0, 0, 0;
  0, 0, 0, 1, 1, 0, 0, 0;
  0, 0, 0, 0, 2, 0, 0, 0, 0;
  0, 0, 0, 0, 1, 1, 0, 0, 0, 0;
  ...
Characteristic polynomial of the 3 X 3 matrix M_3: Phi(3, x) = 1 - 6*x + 5*x^2 - x^3, from {A085478(3, k)}_{k=0..3} = {1, 6, 5, 1}. - _Wolfdieter Lang_, Mar 25 2020
		

References

  • Richard P. Stanley, "Catalan Numbers", Cambridge University Press, 2015.

Crossrefs

Cf. A001333 (permanent of the matrix M).
Cf. A054142, A053123, A011973 (characteristic polynomials of submatrices of M).
Cf. A001700.

A330141 Triangle read by rows: binomial transform of triangle A330140.

Original entry on oeis.org

1, 1, 1, 1, 2, 2, 1, 3, 7, 4, 1, 4, 17, 20, 9, 1, 5, 36, 65, 60, 21, 1, 6, 72, 176, 249, 176, 51, 1, 7, 141, 434, 833, 889, 518, 127, 1, 8, 275, 1016, 2470, 3528, 3094, 1520, 323, 1, 9, 538, 2307, 6795, 12150, 14232, 10494, 4461, 835
Offset: 0

Author

Gary W. Adamson, Dec 02 2019

Keywords

Comments

Row sums = the binomial transform of the Catalan numbers, A007317.
Right border = the Motzkin numbers, A001006.

Examples

			First few rows of the triangle are:
  1;
  1,   1;
  1,   2,   2;
  1,   3,   7,   4;
  1,   4,  17,  20,   9;
  1,   5,  36,  65,  60,  21;
  1,   6,  72, 176, 249, 176,  51;
  1,   7, 141, 434, 833, 889, 518, 127;
  ...
		

Crossrefs