\documentclass[12pt,reqno]{article}

\usepackage[usenames]{color}
\usepackage{amssymb}
\usepackage{amsmath}
\usepackage{amsthm}
\usepackage{amsfonts}
\usepackage{amscd}
\usepackage{graphicx}
\usepackage{mathtools}

\usepackage[colorlinks=true,
linkcolor=webgreen,
filecolor=webbrown,
citecolor=webgreen]{hyperref}

\definecolor{webgreen}{rgb}{0,.5,0}
\definecolor{webbrown}{rgb}{.6,0,0}

\usepackage{color}
\usepackage{fullpage}
\usepackage{float}

\usepackage{psfig}
\usepackage{graphics}
\usepackage{latexsym}
\usepackage{epsf}
\usepackage{breakurl}

\setlength{\textwidth}{6.5in}
\setlength{\oddsidemargin}{.1in}
\setlength{\evensidemargin}{.1in}
\setlength{\topmargin}{-.1in}
\setlength{\textheight}{8.4in}

\newcommand{\seqnum}[1]{\href{https://oeis.org/#1}{\rm \underline{#1}}}
\newcommand\coolover[2]{\mathrlap{\smash{\overbrace{\phantom{%
    \begin{matrix} #2 \end{matrix}}}^{\mbox{$#1$}}}}#2}

\begin{document}

\begin{center}
\epsfxsize=4in
\leavevmode\epsffile{logo129.eps}
\end{center}

\theoremstyle{plain}
\newtheorem{theorem}{Theorem}
\newtheorem{corollary}[theorem]{Corollary}
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}

\theoremstyle{definition}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{example}[theorem]{Example}
\newtheorem{conjecture}[theorem]{Conjecture}

\theoremstyle{remark}
\newtheorem{remark}[theorem]{Remark}


\begin{center}
\vskip 1cm{\Large\bf 
Determinant of Three-Layer Toeplitz Matrices 
} 
\vskip 1cm
\large
Dmitry Efimov\\ 
Institute of Physics and Mathematics\\
Komi Science Centre UrD RAS\\
Kommunisticheskaya Street 24\\
167000 Syktyvkar\\
Russia\\
\href{mailto:dmefim@mail.ru}{\tt dmefim@mail.ru} \\
\end{center}

\vskip .2 in

\begin{abstract}
We obtain an explicit and efficient formula for the determinant of a three-layer Toeplitz matrix. 
We show that many well-known sequences, such as Jacobsthal numbers, generalized Fibonacci numbers, and $k$-Fibonacci numbers,
can be represented as sequences of determinants of three-layer Toeplitz matrices.
Further, we evaluate the spectrum for one of these matrices using the obtained formulae
and, as a consequence, discover some interesting factorizations of certain integer sequences
in terms of products of complex numbers, the imaginary parts of which are expressed using the tangent function.
\end{abstract}

\section{Introduction}
Owing to the Gaussian elimination method, no computational difficulty is encountered when calculating the determinant of an arbitrary matrix.
However, the problem of deriving new formulae for the determinants of various classes of matrices is yet to be solved;
such formulae for determinants can increase the speed of calculations, enable us to study the properties of the corresponding matrices more deeply, 
and expand the application potential.
Herein, we obtain explicit and efficient formulae for the determinants 
of {\itshape  three-layer Toeplitz matrices} $T=(t_{ij})$. Such matrices are given by
$$
t_{ij}=\left\{
 \begin{array}{ll}
   a,&\textrm{if } j-i\geq k;\\
	 b,&\textrm{if } -l<j-i<k;\\
   c,&\textrm{if } j-i\leq -l,
 \end{array}
\right.
$$
where $k,l\in \mathbb{Z}^+$ and $a$, $b$, and $c$ are elements of some field.
As the name suggests, one can clearly distinguish three diagonal layers in such matrices; each of these layers consists of the same elements. 
The middle layer contains the main diagonal.
If the elements of the side layers are the same, then we have a special case considered earlier in \cite{Efimov}. 
Herein, we use fairly elementary methods, namely Laplace expansion and combinations of various recurrence relations.

This paper is organized as follows. In the next section, 
we obtain explicit and efficient formulae for the determinants of three-layer Toeplitz matrices.
In the third section, we consider in detail the special case of such matrices
wherein the middle layer coincides with the main diagonal.
We show that many of the sequences in \cite{Sloane}
can be viewed as sequences of determinants of such matrices.
Several such sequences are combined in the table presented at the end of the paper.
As an application of the obtained formulae,
we calculate the spectrum for the three-layer Toeplitz matrices of a special type.
We also obtain interesting ratios that express elements of some integer sequences
in terms of products of complex numbers, the imaginary parts of which are expressed in terms of the tangent function.

\section{General formulae}
Let $n,k\in \mathbb{Z}^+$.
Consider a Toeplitz matrix $A=(a_{ij})$ of order $n$ with entries
$$
a_{ij}=\left\{
 \begin{array}{ll}
   a,&\textrm{if } j-i\geq k;\\
	 b,&\textrm{if } -1<j-i<k;\\
   c,&\textrm{if } j-i\leq -1,
 \end{array}
\right.
$$
where $a$, $b$, and $c$ are elements of some field.
\vspace{2mm}
\begin{equation}\label{14.04.21_1}
A=\left(
\begin{array}{cccccc}
\mathrlap{\smash{\overbrace{\phantom{\begin{matrix} b&\dots&b\ \end{matrix}}}^k}}b&\dots&b&a&\dots&a\\
c&\ddots&&\ddots&\ddots&\vdots\\
\vdots&\ddots&\ddots&&\ddots&a\\
\vdots&&\ddots&\ddots&&b\\
\vdots&&&\ddots&\ddots&\vdots\\
c&\dots&\dots&\dots&c&b
\end{array}
\right).
\end{equation}

\begin{theorem}\label{15.04.21_4}
Let $s$ be the quotient and $p$ the remainder when $n$ is divided by $k$. Then,
\begin{equation}\label{22.04.21_1}
 \det A=
 \begin{cases}
	   (b-c)^{n-s-\mathrm{sgn}(p)}\cfrac{a(b-c)^{s+\mathrm{sgn}(p)}-c(b-a)^{s+\mathrm{sgn}(p)}}{a-c},& a\not=c;\\
		 (b-c)^{n-1}[b+(s+\mathrm{sgn}(p)-1)c],& a=c,\phantom{\cfrac{u}{v}}
 \end{cases}
\end{equation}
where $\mathrm{sgn}(p)$ is the sign function (i.e., $\mathrm{sgn}(p)=0$ if $p=0$ and $\mathrm{sgn}(p)=1$ if $p>0$).
\end{theorem}

\begin{proof}
Let $d_n$ denote the determinant of a matrix $A$ of order $n$.
If $k\geq n$, then $A$ has the following form:
\begin{equation*}\label{14.04.21_5}
\left(
\begin{array}{ccccc}
b&\dots&\dots&\dots&b\\
c&\ddots&&&\vdots\\
\vdots&\ddots&\ddots&&\vdots\\
\vdots&&\ddots&\ddots&\vdots\\
c&\dots&\dots&c&b
\end{array}
\right).
\end{equation*}
Subtracting the penultimate row from the last one and expanding the determinant
along the last row, we have $d_n=(b-c)d_{n-1}$.
Applying this formula recursively, we obtain 
\begin{equation}\label{26.04.21_1}
 d_n=(b-c)^{n-1}b
\end{equation}
in this case. If $k=n$, then $s=1$ and $p=0$. If $k>n$, then $s=0$ and $p=n$. 
In both the cases, (\ref{26.04.21_1}) is consistent with the first and second lines of (\ref{22.04.21_1}).

Consider a complementary matrix of order $n$ of the following form:
$$
\left(
\begin{array}{ccccc|c}
b&&&&&a\\
c&\ddots&&T&&\vdots\\
\vdots&\ddots&\ddots&&&\vdots\\
\vdots&&\ddots&\ddots&&\vdots\\
c&\dots&\dots&c&b&a\\
\hline
c&\dots&\dots&\dots&c&a
\end{array}
\right).
$$
Here, the triangular area below the main diagonal is completely filled with elements $c$, 
the last column entirely consists of $a$,
and triangular area $T$ consists of arbitrary elements.
Let $f_n$ denote the determinant of such a matrix.
For $n=1$, we assume that the matrix consists of one element $a$ and that $f_1=a$.
If $n\geq 2$, then subtracting the penultimate row from the last one and expanding the determinant
along the last row, we obtain $f_n=(b-c)f_{n-1}$.
Applying this formula recursively, we obtain 
\begin{equation}\label{14.04.21_4}
 f_n=(b-c)^{n-1}a.
\end{equation}

Let $k<n$. 
Then, subtracting the penultimate row of $A$ from the last one, expanding the determinant
along the last row, and repeating this procedure $k$ times, we obtain
\begin{equation}\label{14.04.21_8}
 d_n=(b-c)^k f_{n-k}+(b-c)^{k-1}(b-a)d_{n-k}.
\end{equation}
By assumption, $n=ks+p$, $0\leq p< k$.
Applying (\ref{14.04.21_8}) recursively $s+\mathrm{sgn}(p)-1$ times, we obtain
$$
  d_n=(b-c)^{(s+\mathrm{sgn}(p)-1)(k-1)}(b-a)^{s+\mathrm{sgn}(p)-1} d_{k(1-\mathrm{sgn}(p))+p}+a\sum\limits_{i=1}^{s+\mathrm{sgn}(p)-1}(b-c)^{n-i}(b-a)^{i-1}.
$$
Because $k\geq k(1-\mathrm{sgn}(p))+p$, applying (\ref{26.04.21_1}) to $d_{k(1-\mathrm{sgn}(p))+p}$ yields
\begin{equation}\label{26.04.21_2}
  d_n=\sum\limits_{i=1}^{s+\mathrm{sgn}(p)}\left(\left\lfloor\cfrac{i}{s+\mathrm{sgn}(p)} \right\rfloor (b-a)+a\right)(b-c)^{n-i}(b-a)^{i-1}.
\end{equation}
Using the formula for the sum of a geometric progression, we find that (\ref{26.04.21_2}) is equivalent to (\ref{22.04.21_1}).
\end{proof}

Let $n,k,l\in \mathbb{Z}^+$, $1<l\leq k$.
Consider a Toeplitz matrix $B=(b_{ij})$ of order $n$ with entries
$$
b_{ij}=\left\{
 \begin{array}{ll}
   a,&\textrm{if } j-i\geq k;\\
	 b,&\textrm{if } -l<j-i<k;\\
   c,&\textrm{if } j-i\leq -l,
 \end{array}
\right.
$$
where $a, b, c$ are elements of an arbitrary field:
\vspace{5mm}
\begin{equation}\label{29.06_1}
B=
\hspace{-2mm} 
\begin{matrix}
 l\left\{\vphantom{\begin{matrix} b\\ \vdots\\ b \end{matrix}}\right. \\
 \vphantom{\ddots} \\ 
 \vphantom{\ddots}\\
 \vphantom{c}
\end{matrix}
\hspace{-1mm} 
\left(
\begin{array}{cccccc}
\mathrlap{\smash{\overbrace{\phantom{\begin{matrix} b&\dots&b\ \end{matrix}}}^k}}b&\dots&b&a&\dots&a\\
\vdots&\ddots&&\ddots&\ddots&\vdots\\
b&&\ddots&&\ddots&a\\
c&\ddots&&\ddots&&b\\
\vdots&\ddots&\ddots&&\ddots&\vdots\\
c&\dots&c&b&\dots&b
\end{array}
\right).
\end{equation}

\begin{theorem}\label{29.04.21_1}
Let $s$ be the quotient and $p$ the remainder when $n$ is divided by $k+l-1$. Then,
\begin{equation}\label{26.04.21_3}
 \det B=
 \begin{cases}
	   (b-a)^{(l-1)s}(b-c)^{(k-1)s}\cfrac{a(b-c)^{s+p}-c(b-a)^{s+p}}{(-1)^{(k-1)(l-1)s}(a-c)},& p\leq 1,\ a\not=c;\\
		 (-1)^{(k-1)(l-1)s}(b-c)^{n-1}[b+(s+p-1)c],& p\leq 1,\ a=c;\phantom{\cfrac{u}{v}}\\
     0,& p>1.\phantom{\cfrac{u}{v}}
 \end{cases}
\end{equation}
\end{theorem}

\begin{proof}
It is clear that if $k\geq n$, then matrix $B$ has at least two rows consisting entirely of elements $b$, and its determinant equals $0$.
In this case, $s=0$ and $p=n$, and the result is consistent with (\ref{26.04.21_3}).

Let $k<n$. 
Subtracting the first row from the second one and expanding the determinant
along the second row, we obtain
$$
\det B = (-1)^{k-1}(b-a)\det B',
$$ 
where $B'$ is the matrix in which the first and second rows only differ in the $(k+1)$-th column.
Repeating this procedure recursively $(l-1)$ times, we obtain
$$
\det B = (-1)^{(k-1)(l-1)}(b-a)^{l-1}\det B'',
$$
where
$$
B''=
\hspace{-2mm} 
\begin{matrix}
 k\left\{\vphantom{\begin{matrix} b\\ \ddots\\ b \end{matrix}}\right. \\
 \vphantom{b} \\ 
 \vphantom{\vdots}\\
 \vphantom{\ddots}\\
 \vphantom{\ddots}\\
 \vphantom{\ddots}\\
 \vphantom{c}
\end{matrix}
\hspace{-1mm} 
\left(
\begin{array}{ccc|cccccc}
b&\dots&b&&&&&&\\
&\ddots&\vdots&\ddots&&&a&&\\
c&&b&\dots&b&&&&\\
\hline
&&&b&\dots&b&&&a\\
&&&\vdots&\ddots&&\ddots&&\\
&c&&b&&\ddots&&\ddots&\\
&&&&\ddots&&\ddots&&b\\
&&&&&\ddots&&\ddots&\vdots\\
&&&c&&&b&\dots&b
\end{array}
\right).
$$
A part of matrix $B''$ located above the horizontal line contains $k$ rows,
and the first row contains exactly $k$ $b$'s.
The lower right part of the matrix is formed by the $n-k-l+1$ last rows and columns
of the initial matrix and clearly has the same form as that of the initial matrix. 

One can implement the algorithm recursively $s$ times with a shift of $k$ rows each time; 
i.e., at the next step, the $(k+1)$-th row of matrix $B''$ is subtracted from the $(k+2)$-th row, 
and the determinant is expanded along the $(k+2)$-th row; 
this procedure is continued recursively $(l-1)$ times.
As a result, we have
\begin{equation}\label{27.04.21_1}
  \det B = (-1)^{(k-1)(l-1)s}(b-a)^{(l-1)s}\det B''',
\end{equation}
where
$$
B'''=\left(
\begin{array}{ccccc|c}
b&\dots&b&&a&\\
&\ddots&&\ddots&&a\\
&&\ddots&&b&\\
&c&&\ddots&\vdots&\ddots\\
&&&&b&\dots\\
\hline
&&c&&&M
\end{array}
\right).
$$

The part of matrix $B'''$ located above the horizontal axis consists of $ks$ rows.
The number of elements $b$ in the first row is $k$.
Submatrix $M$ is formed by the intersection of the last $p=n-(k+l-1)s$ rows and columns of matrix $B$.
It is clear that if $p>1$, then matrix $M$ contains 
at least two rows consisting entirely of elements $b$.
Therefore, matrix $B'''$ also has two identical rows, and its determinant equals $0$.
Accordingly, the determinant of matrix $B$ in this case also equals $0$. 

If $p=0$, then $B'''$ is a matrix of the form (\ref{14.04.21_1}) and order $ks$.
By Theorem \ref{15.04.21_4}, we have
\begin{equation*}
 \det B'''=
 \begin{cases}
	   (b-c)^{(k-1)s}\cfrac{a(b-c)^s-c(b-a)^s}{a-c},& a\not=c;\\
		 (b-c)^{ks-1}(b+(s-1)c),& a=c,\phantom{\cfrac{u}{v}}
 \end{cases}
\end{equation*}
which in combination with (\ref{27.04.21_1}) agrees with (\ref{26.04.21_3}).

If $p=1$, then $B'''$ is a matrix of the form (\ref{14.04.21_1}) and order $ks+1$; therefore,
by Theorem \ref{15.04.21_4}, we have
\begin{equation*}
 \det B'''=
 \begin{cases}
	   (b-c)^{(k-1)s}\cfrac{a(b-c)^{s+1}-c(b-a)^{s+1}}{a-c},& a\not=c;\\
		 (b-c)^{ks}(b+sc),& a=c,\phantom{\cfrac{u}{v}}
 \end{cases}
\end{equation*}
which in combination with (\ref{27.04.21_1}) also agrees with (\ref{26.04.21_3}).
\end{proof}

\begin{remark}
 Theorems \ref{15.04.21_4} and \ref{29.04.21_1} allow us to compute 
 determinants of matrices (\ref{14.04.21_1}) and (\ref{29.06_1}) of order $n$ in $O(n)$ time.
\end{remark}

\section{Special case}
Consider in more detail  a three-layer Toeplitz matrix of order $n$ of the form
\begin{equation}\label{30.04.21_1}
C=\left(
\begin{array}{ccccc}
b&a&\dots&\dots&a\\
c&\ddots&\ddots&&\vdots\\
\vdots&\ddots&\ddots&\ddots&\vdots\\
\vdots&&\ddots&\ddots&a\\
c&\dots&\dots&c&b
\end{array}
\right).
\end{equation}
Here $k=l=1$. By Theorem \ref{15.04.21_4}, we have
$$
 \det C=
 \begin{cases}
	   \cfrac{a(b-c)^n-c(b-a)^n}{a-c},& a\not=c;\\
		 (b-c)^{n-1}(b+(n-1)c),& a=c.\phantom{\cfrac{u}{v}}
 \end{cases}
$$
Let $g_n$ denote the determinant of matrix $C$ of order $n$.
It is clear that the ordinary generating function of $(g_n)$ has the following form:
$$
   G(x)=\cfrac{(a+c-b)x+1}{(b-a)(b-c)x^2+(a+c-2b)x+1}.
$$
Accordingly, the following recurrence relation holds for sequence $(g_n)$: 
\begin{equation}\label{29.04.21_4}
   g_n=(2b-a-c)g_{n-1}-(b-a)(b-c)g_{n-2};\ g_0=1,\ g_1=b.
\end{equation}

Many sequences presented in \cite{Sloane} correspond to recurrence relation (\ref{29.04.21_4}) 
up to a shift by one or two elements.
For example, consider the so-called {\itshape metallic mean} \cite{Spinadel} of the following form:
$$
\sigma_k=\frac{k+\sqrt{k^2+4}}{2},\ \ \ k\in\mathbb{Z}^+.
$$
For $k=1$, we obtain the well-known {\itshape golden ratio}; for $k=2$, we obtain the {\itshape silver ratio}, and so on.
Now, we consider the following values as parameters $a$, $b$, and $c$:
$$
a=\sigma_k-k,\ b=0,\ c=-\sigma_k.
$$ 
Substituting these values into (\ref{29.04.21_4}), we obtain the following recurrence relation:
\begin{equation*}
   g_n=kg_{n-1}+g_{n-2};\ g_0=1,\ g_1=0,
\end{equation*}
which defines a sequence of the so-called {\itshape $k$-Fibonacci numbers} \cite{Falcon} shifted by one element to the right. 
For $k=1$, we obtain the Fibonacci sequence; for $k=2$, we obtain the sequence of  Pell numbers, and so on. 
Let $F_{k,n}$ denote the $n$-th member of the $k$-Fibonacci sequence. 
Thus the following proposition holds.
\begin{proposition}
\begin{equation*}
F_{k,n}=
\left|
\begin{array}{ccccc}
0&\sigma_k-k&\dots&\dots&\sigma_k-k\\
-\sigma_k&\ddots&\ddots&&\vdots\\
\vdots&\ddots&\ddots&\ddots&\vdots\\
\vdots&&\ddots&\ddots&\sigma_k-k\\
\mathrlap{\smash{\underbrace{\phantom{\begin{matrix} -\sigma_k&\dots&\dots&-\sigma_k&0&\sigma_k-k\ \end{matrix}}}_{n+1}}}-\sigma_k&\dots&\dots&-\sigma_k&0\\
\end{array}
\right|.
\end{equation*}
\end{proposition}
\vspace{5mm}

Some examples of other sequences from the 
{\it On-Line Encyclopedia of Integer Sequences}
(OEIS), which can also be considered
sequences of determinants of matrices of the form (\ref{30.04.21_1}), are presented in Table \ref{tab_20.03_1}.

Now, assume that $c=-a$; i.e., consider the matrices of the form
\begin{equation*}
D=\left(
\begin{array}{ccccc}
b&a&\dots&\dots&a\\
-a&\ddots&\ddots&&\vdots\\
\vdots&\ddots&\ddots&\ddots&\vdots\\
\vdots&&\ddots&\ddots&a\\
-a&\dots&\dots&-a&b
\end{array}
\right).
\end{equation*}
By Theorem \ref{15.04.21_4}, the determinant of matrix $D$ of order $n$ is given by
\begin{equation}\label{12.05.21_1}
 \det D=\cfrac{(b+a)^n+(b-a)^n}{2}.
\end{equation}

\begin{theorem}\label{12.05.21_3}
If $a,b\in {\mathbb R}$, then the spectrum of matrix $D$ of order $n$ has the following form:
\begin{equation}\label{13.05.21_1}
\sigma(D)=
\begin{cases}
 \left\{b+ia\tan{\cfrac{\pi k}{2m+1}},\ -m\leq k\leq m\right\},&n=2m+1;\\
 \left\{b+ia\tan{\cfrac{\pi (2k+1)}{4m}},\ -m\leq k\leq m-1\right\},&n=2m.
\end{cases}
\end{equation}
\end{theorem}

\begin{proof}
Without loss of generality, we may assume that $a>0$.
First, we find the spectrum of matrix $D$ for $b=0$.
Considering (\ref{12.05.21_1}), the characteristic equation $\det(D-\lambda I)=0$ in this case takes the following form:
$$
(a+\lambda)^n=(-1)^{n+1}(a-\lambda)^n.
$$
The numbers $a+\lambda$ and $a-\lambda$ must have the same modulus.
This is only possible when $\lambda$ is a purely imaginary number.
Let $\lambda=p i$, $p\in{\mathbb R}$; i.e., consider the equation
\begin{equation}\label{11.05.21_2}
  (a+pi)^n=(-1)^{n+1}(a-pi)^n.
\end{equation}
We denote the argument of $a+pi$ by $\phi$.
Because $a>0$, $\phi\in(-\pi/2, \pi/2)$.
Let $n=2m+1$. 
Equating arguments of the numbers on the right and left sides of equation (\ref{11.05.21_2}), we obtain $\phi=\pi k/n$, 
where $k$ takes integer values in the range $[-m,m]$.
Because $\tan{\phi}=p/a$, eigenvalues of $D$ in this case have the following form:
$$
 ia\tan{\cfrac{\pi k}{2m+1}},\ -m\leq k\leq m.
$$
Similarly, if $n=2m$, then equating the arguments of the numbers on the right and left sides of equation (\ref{11.05.21_2}),
we obtain $\phi=\pi (2k+1)/2n$, $k\in [-m,m-1]$, and accordingly, the eigenvalues of matrix $D$ have the following form:
$$
 ia\tan{\cfrac{\pi (2k+1)}{4m}},\ -m\leq k\leq m-1.
$$
Clearly, for an arbitrary $b$, the spectrum of matrix $D$ has the form (\ref{13.05.21_1}).
\end{proof}

As a consequence of (\ref{12.05.21_1}) and Theorem \ref{12.05.21_3}, we obtain the following interesting relation.
\begin{corollary}\label{12.05.21_2}
Let $a,b\in {\mathbb R}$. Then,
$$
\cfrac{(b+a)^n+(b-a)^n}{2}=
\begin{cases}
 \prod\limits_{k=-m}^m\left(b+ia\tan{\cfrac{\pi k}{2m+1}}\right),&n=2m+1;\\
 \prod\limits_{k=-m}^{m-1}\left(b+ia\tan{\cfrac{\pi (2k+1)}{4m}}\right),&n=2m.
\end{cases}
$$
\end{corollary}

\begin{example}[The sequence \seqnum{A007051}]
Let $a=1$, $b=2$. Then, by Corollary \ref{12.05.21_2}, we have
$$
\cfrac{3^n+1}{2}=
\begin{cases}
 \prod\limits_{k=-m}^m\left(2+i\tan{\cfrac{\pi k}{2m+1}}\right),&n=2m+1;\\
 \prod\limits_{k=-m}^{m-1}\left(2+i\tan{\cfrac{\pi (2k+1)}{4m}}\right),&n=2m.
\end{cases}
$$
For example,
\begin{eqnarray*}
 5 &=& (2-i\tan{\frac{\pi}{4}})\cdot (2+i\tan{\frac{\pi}{4}}),\\
 14 &=& 2\cdot(2-i\tan{\frac{\pi}{3}})\cdot (2+i\tan{\frac{\pi}{3}}),\\
 41 &=& (2-i\tan{\frac{3\pi}{8}})\cdot (2-i\tan{\frac{\pi}{8}})\cdot (2+i\tan{\frac{\pi}{8}})\cdot (2+i\tan{\frac{3\pi}{8}}),\\
 122 &=& 2\cdot (2-i\tan{\frac{2\pi}{5}})\cdot (2-i\tan{\frac{\pi}{5}})\cdot (2+i\tan{\frac{\pi}{5}})\cdot (2+i\tan{\frac{2\pi}{5}}),
\end{eqnarray*}
and so on.
\end{example}

\begin{table}[!ht]
\centering
{\renewcommand{\arraystretch}{0}%
\begin{tabular}{|c|c|c|c|c||l|}
 \hline\strut
 \small $a$&$b$&$c$&$k$&$l$&\hspace{20mm} $A$-number in the OEIS\\
 \hline
 \rule{0pt}{1.6pt}&&&&&\\
 \hline\strut
  $1$&\small $0$&$-2$&$1$&$1$&\seqnum{A078008} (Expansion of $(1-x)/((1+x)(1-2x))$)\\
 \hline\strut
  $1$&\small $0$&$-3$&$1$&$1$&\seqnum{A054878} (Number of closed walks of length $n$ ...)\\
 \hline\strut
  $1$& $2$& $-1$& $1$& $1$&\seqnum{A007051} ($a_n=(3^n+1)/2$)\\
 \hline\strut
  $2$&\small $0$&$-3$&$1$&$1$&\seqnum{A102901} ($a_n = a_{n-1} + 6a_{n-2}$, $a_0=1$, $a_1=0$)\\
 \hline\strut
  $2$&\small $1$&$-1$&$1$&$1$&\seqnum{A001045} (Jacobsthal numbers)\\
 \hline\strut
  $2$& $1$& $-2$& $1$& $1$&\seqnum{A046717} ($a_n=2a_{n-1}+3a_{n-2}$, $a_0=a_1=1$)\\
 \hline\strut
  $3$&\small $1$&$-2$&$1$&$1$&\seqnum{A015441} (Generalized Fibonacci numbers)\\
 \hline\strut
  $3$& $1$& $-3$& $1$& $1$&\seqnum{A003665} ($a_n=2^{n-1}(2^n + (-1)^n)$)\\
 \hline\strut
  $3$& $1$& $-4$& $1$& $1$&\seqnum{A320469} ($a_n = 3a_{n-1} + 10a_{n-2}$, $a_0=1$, $a_1=1$)\\
 \hline
\end{tabular}
}
\caption{Correspondence of sequences of determinants of  three-layer Toeplitz matrices to sequences from the OEIS}\label{tab_20.03_1}
\end{table}

\section{Acknowledgments}
The author would like to thank the anonymous referee for careful reading of the paper and helpful comments.

\begin{thebibliography}{9}

\bibitem{Efimov} D. B. Efimov, 
Determinants of generalized binary band matrices,
{\it Proceedings Komi Science Centre UrD RAS} \textbf{4} (28) (2016), 5--10.

\bibitem{Falcon} S. Falc\'on and \'A. Plaza,
On the Fibonacci $k$-numbers,
{\it Chaos Solitons Fractals} \textbf{32} (2007), 1615--1624. 

\bibitem{Spinadel} V. W. de Spinadel,
The metallic means family and renormalization group techniques,
{\it Proc. Steklov Inst. Math.} \textbf{6} (2000), 173--189. 

\bibitem{Sloane} N. J. A. Sloane et al., 
{\it The On-Line Encyclopedia of Integer Sequences}, 2021.  Available
at \url{https://oeis.org}.

\end{thebibliography}


\bigskip
\hrule
\bigskip

\noindent 2010 {\it Mathematics Subject Classification}:
Primary 15A15, Secondary 15B05, 15A18.

\noindent \emph{Keywords}: determinant, Toeplitz matrix, spectrum.
 
\bigskip
\hrule
\bigskip

\noindent(Concerned with sequences 
\seqnum{A000045},
\seqnum{A000129},
\seqnum{A001045},
\seqnum{A003665}, 
\seqnum{A007051},
\seqnum{A015441},
\seqnum{A046717},
\seqnum{A054878},
\seqnum{A078008},
\seqnum{A102901}, and
\seqnum{A320469}.)

\bigskip
\hrule
\bigskip

\vspace*{+.1in}
\noindent
Received June 10 2021;
revised version received  August 20 2021.
Published in {\it Journal of Integer Sequences}, 
October 25 2021.

\bigskip
\hrule
\bigskip

\noindent
Return to \htmladdnormallink{Journal of Integer Sequences home page}{http://www.cs.uwaterloo.ca/journals/JIS/}.
\vskip .1in

\end{document}

