The serendipity family of finite elements

Abstract: We give a new, simple, dimension-independent definition of the serendipity finite element family. The shape functions are the span of all monomials which are linear in at least s-r of the variables where s is the degree of the monomial or, equivalently, whose superlinear degree (total degree with respect to variables entering at least quadratically) is at most r. The degrees of freedom are given by moments of degree at most r-2d on each face of dimension d. We establish unisolvence and a geometric decomposition of the space.

\section{Introduction}
The serendipity family of finite element spaces are among the most popular finite element spaces
for parallelepiped meshes in two, and, to a lesser extent,
three dimensions.  For each such mesh and
each degree $r\ge 1$ they provide a finite element subspace with $C^0$ continuity which has
significantly smaller dimension than the more obvious alternative, the tensor product
Lagrange element family.
However, the serendipity elements
are rarely studied systematically, particularly in 3-D.  Usually only
the lowest degree examples are discussed, with the pattern for higher degrees not evident.
In this paper, we give a simple, but apparently new, definition of the serendipity
elements, by specifying in a dimension-independent fashion the space of shape functions
and a unisolvent set of degrees of freedom.

The serendipity finite element space $\S_r$ may be viewed as a reduction of
the space $\Q_r$, the
tensor product Lagrange finite element space of degree $r\ge1$.
The $\Q_r$ elements are certainly the simplest, and for many purposes the best,
$C^0$ finite elements on parallelepipeds.
They may be defined by specifying a space of polynomial shape functions,
$\Q_r(I^n)$, and a unisolvent set of degrees of freedom, for the unit cube $I^n$.
Here $I=[-1,1]$ is the
unit interval, $n\ge1$ is the space dimension, and
for $f\subset \mathbb{R}^n$ and $r\ge0$ the space $\Q_r(f)$ is defined as the restriction to $f$ of
the functions on $\mathbb{R}^n$ which are polynomial of degree at most $r$ in
each of the $n$ variables separately, so $\dim \Q_r(I^n)=(r+1)^n$.
In addition to the usual evaluation degrees of freedom associated  to each vertex, to each
face $f$ of $I^n$ of some dimension $d\ge 1$ are associated the degrees of freedom
\begin{equation}\label{dofq}
u \mapsto \int_f uq,\quad q\in\Q_{r-2}(f),
\end{equation}
(properly speaking, $q$ ranges over a basis of $\Q_{r-2}(f)$, but we
shall henceforth dispense with this distinction).
Since the number of faces of dimension
$d$ of $I^n$ is $2^{n-d}\binom{n}{d}$, and, by the binomial expansion,
$$
\sum_{d=0}^n 2^{n-d}\binom{n}{d}(r-1)^d = (r+1)^n,
$$
we see that the total number of degrees of freedom coincides with the dimension of $\Q_r(I^n)$.
The proof of unisolvence is straightforward, by using the degrees of freedom to show, in turn,
that $u$ vanishes on the faces of dimension $0,1,\ldots$.

The serendipity family in two dimensions is discussed in most finite element textbooks, especially
the lowest order cases, namely
the 4-node, 8-node, and 12-node rectangular elements.  The idea is to maintain the
same degrees of freedom as for $\Q_r$ on the boundary of the element, but to remove
interior degrees of freedom, and specify a correspondingly smaller shape function space
for which the smaller set of degrees of freedom is unisolvent.  In this way,
one obtains a space of lower dimension without sacrificing $C^0$ continuity and,
hopefully, without much loss of accuracy.
For $r=1$,
there are no interior degrees of freedom for $\Q_r$, so the space $\S_1$ is identical
to $\Q_1$. But the serendipity space $\S_2$ and $\S_3$ have only $8$ and $12$
degrees of freedom, respectively, compared to $9$ and $16$, respectively, for
$\Q_2$ and $\Q_3$.  A possible generalization to higher degree is to keep
only the $4r$ boundary degrees of freedom of $\Q_r$ and to seek $\S_r$ as a
subspace of $\Q_r$ of dimension $4r$ for which these degrees of freedom are
unisolvent.  This is easily accomplished by taking $\S_r$ to be the span
of the monomials $x^i$, $x^iy$, $y^i$, and $xy^i$, for $0\le i\le r+1$, ($4r$
monomials altogether, after accounting for duplicates).  The resulting finite element
space is referred to as the serendipity space in some of the literature,
e.g., \cite{MR1008473}.  However,
for $r>3$ this space does not contain the complete polynomial space $\P_r$, and
so does not achieve the same degree of approximation as $\Q_r$.  Therefore
the shape functions for the serendipity space in two dimensions is usually taken
to be the span of $\P_r(I^n)$ together with the above monomials, or equivalently,
$$
\S_r(I^n) = \P_r(I^n)+\spn[x^ry,xy^r].
$$
The degrees of freedom associated to the vertices and other faces of positive codimension
are taken to
be the same as for $\Q_r$, and the degrees of freedom in the interior of the element
can be taken as the moments $u\mapsto\int_{I^2} uq$, $q\in\P_{r-4}(I^2)$, resulting
in a unisolvent set.

This definition does not generalize in an obvious fashion to three (or more)
dimensions, but many texts discuss the
lowest order cases of serendipity elements in three dimensions: the 20-node brick, and possibly
the 32-node brick \cite{MR1008473, Kaliakin, MR0443377,MR1897985},
which have the same degrees of freedom as $\Q_r(I^3)$ on the boundary, $r=2,3$, but none in
the interior.  It is often remarked that the choice of shape function space is not
obvious, thus motivating the name ``serendipity.''
The pattern to extend these low degree cases to higher degree brick elements is not
evident and usually not discussed.  A notable exception is the text of
Szab\'o and Babu\v{s}ka \cite{MR1164869}, which defines the space of
serendipity polynomials on the three-dimensional
cube for all polynomial degrees, although without using the term serendipity and by an
approach quite different from that given here.  High-degree serendipity elements
on bricks have been used in the $p$-version of the finite element method \cite{MR1067945}.

In this paper we give a simple self-contained dimension-independent
definition of the serendipity family.  For general $n,r\ge 1$,
we define the polynomial space $\S_r(I^n)$ and the degrees of freedom associated to each
face of the $n$-cube and prove unisolvence.

\section{Shape functions and degrees of freedom}
\subsection{Shape functions}
We now give, for general dimension $n\ge1$ and general degree $r\ge 1$, a concise definition
of the space $\S_r(I^n)$ of shape functions for the serendipity finite element.
(As usual,
a monomial is said to be linear in some variable $x_i$ if it is divisible by $x_i$
but not $x_i^2$, and it is said to be superlinear if it is divisible by $x_i^2$.)
\begin{defn}
The serendipity space $\S_r(I^n)$ is the
span of all monomials in $n$ variables which are linear in at least $s-r$ of the
variables where $s$ is the degree of the monomial.
\end{defn}
\noindent
We may express this definition in an alternative form using the notion of the
\emph{superlinear degree}.
Define the superlinear degree of a monomial $p$, denoted $\deg_2 p$, to
be the total degree of $p$ with respect to variables which enter it superlinearly
(so, for example, $\deg_2 x^2 y z^3 = 5$) and define the superlinear
degree of a general polynomial
as the maximum of the superlinear degree of its monomials.
\begin{defnx}
The serendipity space $\S_r(I^n)$ is the
space of all polynomials in $n$ variables with superlinear degree at most $r$.
\end{defnx}
\noindent
It is easy to see that Definitions~2.1 and 2.1$'$ are equivalent, since
for a monomial $p$ which is linear in $l$ variables, $\deg p = \deg_2 p + l$.
Surprisingly, neither form of the definition seems to appear in the literature.

Since any monomial is, trivially, linear in at least $0$ variables, and no monomial in
$n$ variables is linear in more than $n$ variables, we have immediately from Definition~2.1
that $\P_r(I^n)\subset\S_r(I^n)\subset\P_{r+n}(I^n)$, where $\P_r(I^n)$ is the
space of polynomial functions of degree at most $r$ on $I^n$.  In fact, $\S_r(I^n)\subset\P_{r+n-1}(I^n)$, since the only monomial which is linear in all $n$ variables
is $x_1x_2\cdots x_n$, which is of degree $n$, not degree $r+n$.
In particular, the one-dimensional case is trivial: $\S_r(I)=\P_r(I)=\Q_r(I)$, for all $r\ge 1$.
The two-dimensional case is simple as well:
$\S_r(I^2)$ is spanned by $\P_r(I^2)$ and the two monomials $x_1x_2^r$ and $x_1^rx_2$ (these
two coincide when $r=1$).  Thus we recover the usual serendipity shape functions
in two dimensions.
In three dimensions, $\S_r(I^3)$ is obtained by adding to $\P_r(I^3)$ the span
of certain monomials of degrees $r+1$ and $r+2$, namely those of degree $r+1$ which are
linear in at least one of the $n$ variables, and those of degree $r+2$ which are linear
in at least two of them (there are three of these---$x_1x_2x_3^{r}$, $x_1x_2^{r}x_3$,
$x_1^{r}x_2x_3$---except for $r=1$, when all three coincide).

To calculate $\dim\S_r(I^n)$ we count the monomials in $n$ variables
with superlinear degree
at most $r$.  For any monomial $p$ in $x_1,\ldots,x_n$,
let $J\subset\range:=\{1,\ldots,n\}$ be the set of indices for which $x_i$ enters $p$ superlinearly
and let $d\ge0$ be the
cardinality of $J$.  Then
$$
p= \bigl(\prod_{j\in J} x_j^2\bigr)\x q \x \bigl(\prod_{i\in J^c} x_i^{a_i}\bigr),
$$
where $q$ is a monomial in the $d$ variables indexed by $J$ and each
$a_i$, $i\in J^c$ (the complement of $J$), equals either $0$ or $1$.  Note that $\deg_2p = 2d+\deg q$.
Thus we may uniquely specify a monomial in $n$ variables with superlinear degree
at most $r$
by choosing $d\ge 0$, choosing a set $J$ consisting of $d$ of the $n$ variables
(for which there are $\binom{n}{d}$ possibilities), choosing a monomial
of degree at most $r-2d$ in the $d$ variables ($\binom{r-d}{d}$
possibilities), and choosing the exponent to be
either $0$ or $1$ for the $n-d$ remaining indices
($2^{n-d}$ possibilities).  Thus
\begin{equation}\label{dimen}
\dim\S_r(I^n) = \sum_{d=0}^{\min(n,\lfloor r/2\rfloor)} 2^{n-d}\binom{n}{d}
\binom{r-d}{d}.
\end{equation}
Table~1 shows the dimension for small values of $n$ and $r$.
\begin{table}[t]
 \begin{tabular}{c|cccccccc}
 \multicolumn{4}{c}{} & \multicolumn{2}{c}{$r$}\\
 $n$ & 1 & 2 & 3 & 4 & 5 & 6 & 7 & 8\\
\hline
1 & 2 & 3& 4 & 5 & 6 & 7 & 8 & 9\\
2 & 4 & 8 & 12 & 17 & 23 & 30 & 38 & 47\\
3 & 8 & 20 & 32 & 50 & 74 & 105 & 144 & 192\\
4 & 16 & 48 & 80 & 136 & 216 & 328 & 480 & 681\\
5 & 32 & 112 & 192 & 352 & 592 & 952 & 1472 & 2202
 \end{tabular}
\bigskip
\caption{$\dim S_r(I^n)$}
\end{table}

\subsection{Degrees of freedom}
We complete the definition of the serendipity finite elements by specifying
a set of degrees of freedom and proving that they are unisolvent.  Let $f$ be
a face of $I^n$ of dimension $d\ge 0$.  Then the degrees of freedom associated
to $f$ are given by
\begin{equation}\label{serdof}
u \mapsto \int_f uq, \quad q\in \P_{r-2d}(f).
\end{equation}
Note that $\P_s(f)$ is defined to be the space of restrictions
to $f$ of $\P_s(I^n)$, so if $f$ is a vertex, then $\P_s(f)=\mathbb{R}$ for all $s\ge 0$.
In this case the integral is with respect to the counting measure,
so each vertex is assigned the evaluation degree of freedom.

In contrast to the degrees of freedom \eqref{dofq} of the tensor product
Lagrange family, for the serendipity family, the degrees of freedom on faces are
given in terms of the $\P_s(f)$ rather than $\Q_s(f)$.  For example, for $\S_6(I^2)$, there
are $6=\dim\P_{6-2\times 2}(I^2)$ degrees of freedom internal to
the square.  For $\S_6(I^3)$,
there are $6$ degrees of freedom on each face of the cube and one degree of freedom
internal to the cube.
Figure~\ref{fg:seren} shows degree of freedom diagrams for the spaces
$\S_r(I^n)$ for $r\le4$ and $n\le3$.
\begin{figure}[tb]
 \includegraphics[width=6in]{serendipity-fig.pdf}
\caption[]{Degrees of freedom for $\S_r(I^n)$ for $r=1,2,3,4$, $n=1,2,3$.
Dots indicate the number of degrees of freedom associated to each
face.}\label{fg:seren}
\end{figure}

\begin{thm}[Unisolvence] The degrees of freedom specified in \eqref{serdof} are unisolvent
 for $\S_r(I^n)$.
\end{thm}
\begin{proof}
 First we note that number of degrees of freedom equals the dimension
of the space.
Indeed, there are $2^{n-d}\binom{n}{d}$ faces of the cube of dimension $d$ (obtained
by fixing $n-d$ of the variables to $\pm 1$) and, for $f$ of dimension
$d$, $\dim\P_{r-2d}(f)=\binom{r-d}{d}$, so the total number of degrees of
freedom proposed is precisely $\dim\S_r(I^n)$ given by \eqref{dimen}.
It remains to show that if $u\in\S_r(I^n)$ and all the
quantities in \eqref{serdof} vanish, then $u$ vanishes.  We do this by induction
on $n$, the case $n=1$ being trivial.  Let $F$ be a face of $I^n$ of dimension $n-1$
and let $U=u|_F$.  The $U$ is a polynomial in $n-1$ variables with superlinear degree
at most $r$, i.e., $U\in\S_r(F)$.  Moreover, if $f$ is any face of $F$ of dimension
$d$, and $q\in\P_{r-2d}(f)$, then $\int_f Uq = \int_f uq=0$.  Therefore, by the
inductive hypothesis, $U$ vanishes.  In this way, we see that $u$ vanishes on all its
faces, i.e., whenever we fix some $x_i$ to $\pm 1$.  Therefore
$$
u=(1-x_1^2)\cdots(1-x_n^2)\, p
$$
for some polynomial $p$.  Note that $\deg_2 u = 2n + \deg p$.  Therefore $\deg p \le r-2n$,
and we make take $f=I^n$ and $q=p$ in \eqref{serdof} to find
$$
\int_{I^n}(1-x_1^2)\cdots(1-x_n^2) \, p^2 = 0.
$$
It follows that $p=0$.
\end{proof}

In the proof of unisolvence we established that the degrees of freedom associated to
a face and its subfaces determine the restriction of $u$ to the face.  This is important,
since it implies that
an assembled serendipity finite element function is continuous.

\section{Geometric decomposition}
In this section we give a geometric decomposition of $\S_r(I^n)$,
by which we mean a direct sum decomposition into subspaces associated to
the faces.  Such a decomposition can be used to derive explicit
local bases which are useful for the efficient implementation
of the elements, and also for insight.

First we introduce some notation.
Let $\Delta_d(I^n)$ denote the set of faces of the $n$-cube of dimension $d$
and $\Delta(I^n)=\bigcup_{d=0}^n\Delta_d(I^n)$ the set of all faces of all dimensions.
A face $f\in\Delta_d(I^n)$ is determined by the equations $x_j=c_j$ for $j\in J$
where $J \subset \range$ is a set of cardinality $n-d$ and each $c_j\in\{-1,1\}$.
We define the bubble function for $f$ as
$$
b_f = \prod_{j \in J^c} (1-x_j^2)\prod_{j \in J} (1+c_jx_j),
$$
the unique (up to constant multiple) nontrivial
polynomial of lowest degree vanishing on the $(n-1)$-dimensional faces of $I^n$ which
do not contain $f$.  Note that $b_f$ is strictly positive on
the relative interior of $f$ and vanishes on all faces of $I^n$ which do not contain $f$.
For $s\ge0$ we denote by $\P_s^f(I^n)$ the space of
polynomials of degree at most $s$ in the $d$ variables $x_j$, $j\in J^c$.  If $d=0$, then
$\P_s^f(I^n)$ is understood to be $\mathbb{R}$.

To a face $f$ of dimension $d$ we associate the space of polynomials
$V_f := \P_{r-2d}^f(I^n)b_f$. By definition,
any element $q\in V_f$ has the form $q=p\prod_{j \in J} (x_j+c_j)$ where $p\in\P_r^f(I^n)$,
and so $\deg_2 q \le r$.  Thus $V_f\subset\S_r(I^n)$.
The following theorem states that they do indeed form a geometric decomposition.

\begin{thm}[Geometric decomposition of the serendipity space]
Let $\S_r(I^n)$ denote the serendipity space of degree $r$, and for each $f\in\Delta_d(I^n)$,
let $V_f= \P_{r-2d}^f(I^n)b_f$.  Then
\begin{equation} \label{geo1}
\S_r(I^n) = \sum_{f\in\Delta(I^n)}V_f.
\end{equation}
Moreover, the sum is direct.
\end{thm}
\begin{proof}
Clearly $\dim V_f=\dim P_{r-2d}^f(I^n)=\binom{r-d}{d}$, and so \eqref{dimen} implies that
$$
\dim\S_r(I^n) = \sum_{f\in\Delta(I^n)}\dim V_f.
$$
Hence, it is sufficient to prove that if $p$ is a monomial with superlinear degree $\le r$,
then $p\in\sum_f V_f$.
Write $p=x_1^{\alpha_1}\cdots x_n^{\alpha_n}$ and let
\begin{equation*}
 J_0=\{\,i\in\range\,|\, \alpha_i=0\,\},\quad
 J_1=\{\,i\in\range\,|\, \alpha_i=1\,\},\quad
 J_2 = (J_0\cup J_1)^c,
\end{equation*}
denote the sets indexing the variables in which $p$ is constant, linear, and superlinear, respectively.
By assumption, $\sum_{j\in J_2}\alpha_j \le r$.
Now we expand
\begin{equation*}
x_j^{\alpha_j} =
 \begin{cases}
 \frac12(1+x_j)+\frac12(1-x_j), & j\in J_0, \\
\frac12(1+x_j)-\frac12(1-x_j), & j\in J_1, \\
\frac12(1+x_j) + \frac{(-1)^{\alpha_j}}2(1-x_j) + (1-x_j^2)r_j(x_j), & j\in J_2,
 \end{cases}
\end{equation*}
where $\deg r_j = \alpha_j-2$, and insert these into $p=x_1^{\alpha_1}\cdots x_n^{\alpha_n}$.
We find that $p$
is a linear combination of terms of the form
\begin{equation}\label{defq}
 q = \prod_{j\in J'} r_j(x_j)(1-x_j^2) \prod_{j\in J''\cup J_0 \cup J_1}(1\pm x_j),
\end{equation}
where $J'$ and $J''$ partition $J_2$.  Let $f$ be the face given by $1\mp x_j=0$, $j\in J''\cup J_0\cup J_1$,
where the signs are chosen opposite to those on \eqref{defq}.  Then $\dim f= \#J'$ and
$q = rb_f$ with $r\in \P_s^f(I^n)$, $s=\sum_{j\in J'}(\alpha_j-2)\le r-2\dim f$.  This shows
that $q\in V_f$ and $p \in \sum_f V_f$, as desired.
\end{proof}

\bibliographystyle{amsplain}


Envelope of intermediate lines of a plane curve

Abstract: For a pair of points in a smooth closed convex planar curve $\gamma$, its mid-line is the line containing its mid-point and the intersection point of the corresponding pair of tangent lines. It is well known that the envelope of the mid-lines ($EML$) is formed by the union of three affine invariants sets: Affine Envelope Symmetry Sets ($AESS$); Mid-Parallel Tangent Locus ($MPTL$) and Affine Evolute of $\gamma$. In this paper, we generalized these concepts by considering the envelope of the intermediate lines. For a pair of points of $\gamma$, its intermediate line is the line containing an intermediate point and the intersection point of the corresponding pair of tangent lines. Here, we present the envelope of intermediate lines ($EIL$) of the curve $\gamma$ and prove that this set is formed by three disconnected sets when the intermediate point is different from the mid-point: Affine Envelope of Intermediate Lines ($AEIL$); the curve $\gamma$ itself and the Intermediate-Parallel Tangent Locus ($IPTL$). When the intermediate point coincides with the mid-point, the $EIL$ coincides with the $EML$, and thus these sets are connected. Moreover, we introduce some standard techniques of singularity theory and use them to explain the local behavior of this set.

\section{Introduction}\label{intro}
Consider two different points $p_1$, $p_2$ of a smooth planar curve $\gamma$ which has non-parallel tangent lines at $p_1$ and $p_2$. The mid-line is the line that passes through the mid-point $M$ of $p_1$ and $p_2$ and the intersection of the tangent lines at $p_1$ and $p_2$. If these tangent lines are parallel, the mid-line is the line which passes through the point $M$ and is parallel to the both tangent lines. When $p_1=p_2$, the mid-line is just the affine normal at the point. The envelope of these mid-lines is an important affine invariant set associated with the curve, which has been studied by many authors  (see for instance \cite{giblin,Holtom,Warder,Bruce}). In \cite{ady1,ady2} the authors presented a generalization of these concepts for hypersurfaces in $\mathbb{R}^n$.

The {\it Envelope of Mid-Lines}, (denoted by $EML$), of the curve $\gamma$ is divided into three parts: The {\it Affine Envelope Symmetry Set} (denoted by $AESS$), the {\it Mid-Parallel Tangent Locus} (denoted by $MPTL$) and the {\it Affine Evolute}, which corresponded to the pairs of points $(p_1,p_2)$ with non-parallel, non-coincident parallel and coincident tangent lines, respectively.
The geometry of these curves has also been extensively studied by many mathematicians (see for example \cite{giblin},\cite{Holtom}).

Unsurprisingly, the envelope of the tangent lines to $\gamma$ contains -at least- $\gamma$ itself. It is natural to ask what lies between the envelope of tangents and the envelope of mid-lines. Consider two disjoint points $p_1$, $p_2$ of the curve $\gamma$ with non-parallel tangent lines. An intermediate point of the points $p_1$ and $p_2$ is the point $M_\alpha=(1-\alpha)p_1+\alpha p_2$ with $\alpha \in [0,1]$. The intermediate line is the line which passes through the intersection point of the tangent lines  at $p_1$ and $p_2$ and the point $M_{\alpha}$ (Figure \ref{fig:Inter} illustrates an intermediate line for the points $p_1$ and $p_2$).
\begin{figure}[tp]
	\includegraphics[scale=0.4]{reta_inter.png}
	\caption{Intermediate line of the points $p_1$ and $p_2$.}\label{fig:Inter}
\end{figure}
When $\alpha=0$ or $\alpha=1$ the intermediate line corresponds to the tangent lines in points $p_1$ and $p_2$ respectively. When $\alpha=1/2,$ the intermediate line corresponds to the mid-line.
Naturally, the {\it Envelope of Intermediate Lines} (denoted by $EIL$) is also divided into three subsets:  The \textit{Affine Envelope of Intermediate Lines} (denoted by $AEIL$), which corresponds to the pairs of points $(p_1,p_2)$ with non-parallel tangent lines; the \textit{Intermediate-Parallel Tangent Locus} (denoted by $IPTL$) is consonant with the pairs of disjoint points $(p_1,p_2)$ with parallel tangent lines; and another curve called \textit{Coincident Tangent Lines}, (denoted by $CTL$), which corresponds to the limit case of coincident points.
This article aims to explain how some basic techniques of Affine Geometry and Singularity Theory enable us to say exactly what happens with the envelope of intermediate lines when $\alpha$ varies in $[0,1]\setminus \{1/2\}$. Precisely speaking, we will obtain the accurate conditions for which the cusps appear and disappear for a fixed $\alpha$ and, crucially, explaining the contribution of inflexions.

We demonstrate that, unlike the case of $\alpha=1/2$ in which $AESS$ and $MPTL$ meet at their ordinary cusp singularities, when $\alpha\neq 1/2,$ $AEIL$ and $IPTL$ are disjointed.

The Affine Evolute of a planar curve is the limit of the $EML$, when $p_1=p_2$. In this paper we study the corresponding limit set $CTL$ and prove that this set is the curve $\gamma$ itself, when $\alpha\neq 1/2$.

The paper is organized as follows: In section \ref{sec:preliminaries}, we review some basic concepts of the Affine Differential Geometry of curves in $\mathbb{R}^2$. In section \ref{sec:AEIL}, we study the $AEIL$, giving necessary and sufficient conditions for a pair $(p_1,p_2)$ to contribute to this set.
In section \ref{sec:local}, we investigate the local behavior of the curves $AEIL$ and $IPTL$ and prove that the set $CTL$ is the proper curve $\gamma$. In this section, we also study non-oval smooth curves.

\section{Preliminaries}\label{sec:preliminaries}

In this section, we present the basic concepts of Affine Differential Geometry of planar smooth curves. We refer to \cite{Nomizu,Buchin} for more details on this subject.

Let $\gamma: [0,1] \longrightarrow \mathbb{R}^2$ be a $\mathcal{C}^{\infty}$ plane curve, without self-intersections or cusps, parametrized by $t$. The planar affine differential geometry is mainly focused on defining a new parametrization, $s$, which is an affine-invariant. The simplest affine-invariant parametrization $s$ for the curve $\gamma$ satisfies the following relation,
\begin{equation}\label{pcaa}
[ \gamma_s, \gamma_{ss}]=1,
\end{equation}

\noindent where $[\,,\,]$ is the notation for determinants.
When a curve satisfies equation \eqref{pcaa}, we say it this
parameterized by affine arc length parametrization $s$.

The vectors $\gamma_s$ and $\gamma_{ss}$ are the affine tangent and the affine normal of the curve $\gamma$, respectively.

The parameters $s$ and $t$ are related by
$$ [\gamma_t,\gamma_{tt}]=\left[\gamma_s\,s_t\,,\,
\gamma_{ss}\,s_t^2\,+\,\gamma_s\,s_{tt}\right]=s_t^3[\gamma_s,\gamma_{ss}]=s_t^3,$$
therefore,
$$ \dfrac{ds}{dt}= [\gamma_t,\gamma_{tt}]^{\frac{1}{3}}.$$
By differentiating the equation \eqref{pcaa}, we obtain
$$ [\gamma_s, \gamma_{sss}]=0 \Rightarrow
\gamma_{sss}+ \mu(s)\gamma_s=0,$$ for some $\mu(s) \in \mathbb{R}.$
The function $\mu(s)$ is called the affine curvature and is
the simplest non-trivial affine differential invariant, defining
up to an affine transformation. Notice that
$$ [\gamma_{s}, \gamma_{sss}]=0  \Rightarrow \gamma_{sss}=-\mu(s)\gamma_s,$$ therefore, we conclude that
$$\mu(s)=[\gamma_{ss}, \gamma_{sss}].$$

\begin{theorem}[\cite{Nomizu}] \label{thm1}
	Plane curves have constant affine curvature if and only if they are conic sections.
\end{theorem}

\begin{proposition}\label{Prop:Affine Normal}
	 Let $\gamma: \mathbb{R} \longrightarrow \mathbb{R}$ be a regular plane curve parametrized by an arbitrary parameter $t$. The affine normal $\xi(t)$ at $t$ is given by:
	\begin{eqnarray}\label{eq0}
	\xi(t)=\kappa^{-\dfrac{2}{3}}\gamma_{tt}-\dfrac{1}{3}\kappa_t\kappa^{-\dfrac{5}{3}}\gamma_t
	\end{eqnarray}
\end{proposition}
The affine curvature of a planar curve $\gamma$ parametrized by an arbitrary parameter is given in the next proposition.

\begin{proposition}\label{Prop:Affine Curvature}
	Let $\gamma$ be a smooth plane curve without inflection points pa\-ra\-me\-tri\-zed by an arbitrary parameter $t$. The affine curvature of $\gamma$ is given by
	\begin{eqnarray} \label{eq00}
	\mu = \dfrac{1}{9}\left(3\kappa\kappa_{tt}-5\kappa_t^2+9\kappa[\gamma_{tt},\gamma_{ttt}]\right)\kappa^{-\frac{8}{3}},
	\end{eqnarray}
	where $\kappa=[\gamma_t,\gamma_{tt}],$ denotes the Euclidean curvature of $\gamma$ and the indexes denote the derivative in respect to the given parameter.
\end{proposition}
\begin{proof}
	Note that $s_t=\kappa^{\frac{1}{3}}$ and also
	$\gamma_s=\gamma_t\kappa^{-\frac{1}{3}}$. Now, by calculating
	$\gamma_{ss}$, $\gamma_{sss}$ and using the fact that
	$\kappa_t=[\gamma_t,\gamma_{ttt}]$, one can conclude that
	$\mu=[\gamma_{ss},\gamma_{sss}].$
\end{proof}
Consider the plane curve $\gamma$ in the Monge form parametrization without Euclidean inflection at the origin, that is,

$$\gamma(t)=\left(t, \displaystyle\frac{1}{2}a_2t^2+ \cdots +\dfrac{1}{k!}a_kt^k+g(t)t^{k+1}
\right),$$ with $a_i \in \mathbb{R}$, $a_2\neq0$ and $g$ a smooth function. Using the Proposition \ref{Prop:Affine Curvature}, the affine curvature of $\gamma$ at the origin is
$$\mu(0)=\dfrac{3a_2a_4-5a_3^2}{9a_2^{\frac{8}{3}}}.$$

This means that the affine curvature depends on $4$-jet of the curve $\gamma$ (see \cite{giblin4} for further information on jet space).
Throughout this paper we deal with the envelope of a family of lines defined on the curve $\gamma$. To be more precise, we have the following definition.

\begin{definition} \label{env}
	The envelope, or discriminant, of a family $$F: \mathbb{R}^n \times \mathbb{R}^r \rightarrow \mathbb{R} $$ of $n$ parameters, is the set: $$\mathcal{D}_F = \left\{x \in \mathbb{R}^r  \, | \, \exists \, u \in \mathbb{R}^n, F(u,x)=\frac{\partial F}{\partial u_i}(u, x)=0, i=1, \dots, n  \right\}.$$
\end{definition}
Consider a smooth plane curve $\gamma$. Take the points $p_1, p_2 \in \gamma$ and let $\gamma_1(t)$ and $\gamma_2(s)$ be the local parameterizations around $p_1$ and $p_2$ with $l_1$ and $l_2$ as their tangent lines of $\gamma$ at $p_1$ and $p_2$, respectively. If $l_1$ and $l_2$ are concurrent then consider $R$ as their intersection point. For $0<\alpha<1$, let $M_\alpha$ be the intermediate point of the segment joining $p_1$ and $p_2$, i.e., $M_\alpha=(1-\alpha)\gamma_1(t) +\alpha\gamma_2(s)$.
\begin{definition}
	For a given oval smooth closed plane curve $\gamma$, the \textbf{intermediate line} of $\gamma$ is:
	\begin{enumerate}
		\item[$i.$] the line that passes through $R$ and $M_{\alpha}$ when $l_1$ and $l_2$ are concurrent with $t \neq s$;
		\item[$ii.$] the only line that passes through the point $M_{\alpha}$ and is parallel to the parallel tangent lines $l_1$ and $l_2$ with $t \neq s$;
		\item[$iii.$] the tangent line of $\gamma$ at $p_1$, when $t=s$. In the case of $\alpha=1/2$ (i.e. the mid-point), the intermediate line is the affine normal of $\gamma$ at $p_1$.
	\end{enumerate}
\end{definition}
We will use the co-normal maps to define the equation of the intermediate line .
\begin{definition}
	Let $\gamma$ be a smooth curve. For $p\in\gamma$, let $\nu=\nu(p)$ the linear functional, such that $$\nu(p)(\xi)=1 \ \ \textrm{and} \ \ \nu(p)(X)=0,$$ where $\xi$ is the affine normal of $\gamma$ and $X$ is a tangent vector of $\gamma.$ The dif\-fe\-ren\-ti\-a\-ble map $\nu$ is called the {\it co-normal map.}
\end{definition}
For further details about the properties of co-normal map, see for instance \cite{Nomizu}.

Let $\nu_i,$ with  $i=1, 2,$ be the co-normal maps at $p_i,$ i.e., $\nu_i=\nu(p_i)$ with $$\nu_i(\xi_i)=1 \mbox{ \space  and \space  } \nu_i(\gamma_i')=0,$$
where $\xi_i$ is the affine normal at the point $p_i$.
\begin{lemma}\label{eri}
	The equation of the intermediate line is given by
	$$(1-\alpha)\nu_2(C)\nu_1(X-M_\alpha)+\alpha\nu_1(C)\nu_2(X-M_\alpha)=0,$$ where $\alpha \in (0,1)$ and $C=p_2-p_1$ is the chord joining the points $p_1$ and $p_2$.
\end{lemma}
\begin{proof}
	Let $\nu$ be the functional that annul the director vector of the intermediate line. Then:  $$\nu=a\,\nu_1 + b\,\nu_2.$$
	Notice that:
	\begin{eqnarray*}
	\nu(R-M_\alpha)=0 & \Longleftrightarrow & (a\,\nu_1 + b\,\nu_2)((1-\alpha)C + R - p_2)=0  \\
	&  & a\nu_1((1-\alpha)C)+a\nu_1(R-p_2)+b\nu_2((1-\alpha)C)=0.
	\end{eqnarray*}

	In fact, $\nu_2(R-p_2)=0$, since $R-p_2$ is tangent to $\gamma_2$. As $R-p_2 = R-p_1 -C$, the last equation is equivalent to:

	\begin{align*}
	 a\nu_1((1-\alpha)C) +a\nu_1(R-p_1) -a\nu_1(C)+b\nu_2((1-\alpha)C) &=0\\
	 -a\nu_1(\alpha C) + b\nu_2((1-\alpha)C)&=0,
	\end{align*}
	because $R-p_1$ is tangent to $\gamma_1$.

	Therefore, taking $a=(1-\alpha)\nu_2(C)$ and $b=\alpha\nu_1(C)$, we obtain the result.
	\end{proof}

Consider the family $F:\gamma_1 \times \gamma_2 \times \mathbb{R}^2  \longrightarrow \mathbb{R}$ given by:
\begin{equation}\label{intermediate}
F(s,t,X)=(1-\alpha)\nu_2(C)\nu_1(X-M_\alpha)+\alpha\nu_1(C)\nu_2(X-M_\alpha),
\end{equation}
where $\alpha \in (0,1).$
For fixed points $p_1$ and $p_2$, and using the Lemma \ref{eri}, one can conclude that the equation of the intermediate line is $F(s,t,X)=0.$

According to the Definition \ref{env}, we intend to find the discriminant set $\mathcal{D}_F$ of the family $F$. More precisely,
\begin{align}\label{discriminant}
\mathcal{D}_F=\left\{X\in\mathbb{R}^2 \, | \, \exists \, (s,t); F=F_s=F_t=0 \, \textrm{at} \, (X,s,t)\right\}.
\end{align}
\section{Intermediate lines with non-parallel tangents}\label{sec:AEIL}

We begin with the following simple lemma:

\begin{lemma}\label{deco}
	The derivatives of co-normal maps $\nu_1$ and $\nu_2$ are given by:
	$$\nu_1^{\prime}=a\nu_1 + b\nu_2 \ \ \textrm{and} \ \ \nu_2^{\prime}=\bar{a}\nu_1 + \bar{b}\nu_2, $$
	where $a=\dfrac{\nu_1^{\prime}(\gamma_2^{\prime})}{\nu_1(\gamma_2^{\prime})}, b=\dfrac{\nu_1^{\prime}(\gamma_1^{\prime})}{\nu_2(\gamma_1^{\prime})}, \bar{a}=\dfrac{\nu_2^{\prime}(\gamma_2^{\prime})}{\nu_1(\gamma_2^{\prime})}$ and $\bar{b}=\dfrac{\nu_2^{\prime}(\gamma_1^{\prime})}{\nu_2(\gamma_1^{\prime})}.$
\end{lemma}

\begin{proof} Take the basis $\left\{\nu_1,\nu_2\right\}$ of the dual plane $\mathbb{R}_2.$ Thus, we can write the linear functional $\nu_1^{\prime}$ as a linear combination of the basis vectors, i.e., $\nu'_1=a\nu_1+b\nu_2.$ By applying $\nu'_1$ to $\gamma'_1$ we obtain $\nu'_1(\gamma'_1)= b\nu_2(\gamma'_1),$ which leads us to achieve
\[b=\dfrac{\nu_1^{\prime}(\gamma_1^{\prime})}{\nu_2(\gamma_1^{\prime})}.\]
The other coefficients are obtained analogously.
\end{proof}

The next result allows us to locate pairs of points of $\gamma$ that contribute to the envelope of intermediate lines, through which we can relate the parameters $s$ and $t$.

\begin{theorem}\label{pre}
	The necessary and sufficient condition for which the pair $(s, t)$ determines one point of the envelope of intermediate lines is

	\begin{eqnarray}\label{eq1}
	\nu_1(C)=-\lambda\nu_2(C),
	\end{eqnarray}
	 where
	$$\lambda=\left(\dfrac{1-\alpha}{\alpha}\right)^{\frac{1}{3}}.$$
\end{theorem}

\begin{proof} Since $\nu_1(\gamma'_1(t))=\nu_2(\gamma'_2(s))=0$, so the derivative of the Equation \eqref{intermediate} with respect to $t$ is
\begin{align*}
	F_t= & (1-\alpha)\nu_2(-\gamma'_1)\nu_1(X-M_\alpha)+(1-\alpha)\nu_2(C)(\nu'_1(X-M_\alpha))+ \\  & +\alpha\nu'_1(C)\nu_2(X-M_\alpha)-\alpha \nu_1(C)\nu_2((1-\alpha)\gamma'_1).
\end{align*}
Using Lemma \ref{deco}, we obtain:
	\begin{align}\label{eq2}
F_t = & -(1-\alpha)\nu_2(\gamma'_1)\nu_1(X-M_\alpha)+b\nu_2(C)\nu_2(X-M_\alpha) \nonumber \\  &-\alpha(1-\alpha)\nu_1(C)\nu_2(\gamma'_1)+aF.
	\end{align}
Similarly,
	\begin{align}\label{eq3}
	F_s = & \ \overline{a}\nu_1(C)\nu_1(X-M_\alpha)+\alpha\nu_1(\gamma'_2)\nu_2(X-M_\alpha) \nonumber \\ &- \alpha(1-\alpha)\nu_2(C)\nu_1(\gamma'_2)+\overline{b}F,
	\end{align}

By isolating $\nu_2(X-M_\alpha)$ from $F=0$ and replacing it in the equations $F_t=0$ and $F_s=0$, we obtain
\begin{eqnarray}\label{eq333}
	\nu_1(X-M_\alpha)=\dfrac{-\alpha^2\nu^2_1(C)\nu_2(\gamma'_1)     }{\alpha\nu_1(C)\nu_2(\gamma'_1) +b\nu^2_2(C)},
\end{eqnarray}
\begin{eqnarray}\label{eq444}
	\nu_1(X-M_\alpha)=\dfrac{\alpha(1-\alpha)\nu_1(C)\nu_2(C)\nu_1(\gamma'_2) }{\overline{a}\nu^2_1(C)-(1-\alpha)\nu_2(C)\nu_1(\gamma'_2)}.
\end{eqnarray}
	%

For the rest of the proof, we
 equate the equations \eqref{eq333} and \eqref{eq444}, and use the formulae of $b$ and $\bar{a}$ given in the Lemma \eqref{deco}.
	%

\end{proof}
\begin{remark}
	If $\gamma_1$ and $\gamma_2$  are parametrized by the affine arc length parametrization, then $\lambda=\beta.$ Thus, the condition given in the Theorem \ref{pre} is equivalent to
$$[\gamma_1^{\prime}+\lambda\gamma_2^{\prime}, \gamma_2-\gamma_1]=0,$$
where $\lambda=\left(\dfrac{1-\alpha}{\alpha}\right)^\frac{1}{3}.$
\end{remark}
\begin{theorem}\label{Theo:Solution of Discr}
Let $F$ be the family of the intermediate lines given in \eqref{intermediate}. The discriminant of the family $F$ admits a solution if and only the equation \eqref{pre} holds. Moreover, the solution of the system is given by
	\begin{equation}
	\displaystyle X-M_\alpha= \frac{\alpha\nu_1(C)}{\alpha\nu_2(\gamma_1')\nu_1(C)+ b\nu_2^2(C)  }\left(  (1-\alpha)\nu_2(C)\gamma_1'- \frac{\alpha\nu_1(C)\nu_2(\gamma_1')}{\nu_1(\gamma_2')}\gamma_2'       \right),

	\end{equation}
	where $\alpha \in (0,1).$
\end{theorem}
\begin{proof}
	Since  $\gamma'_1$ and $\gamma'_2$ are not parallels, we can write $X-M_\alpha= \overline{A}\gamma'_1 + \overline{B}\gamma'_2.$ Thus, by applying  $\nu_2$ the latter equation and isolating $\nu_2(X-M_\alpha)$ in the equation of intermediate line, we have:

	$$\overline{A}=\dfrac{\nu_2(X-M_\alpha)}{\nu_2(\gamma'_1)}=\dfrac{-(1-\alpha)\nu_2(C)\nu_1(X-M_\alpha)}{\alpha\nu_2(\gamma'_1)\nu_1(C)} $$

	Using the expression  \eqref{eq333} :
	\begin{align*}
	\overline{A} = & \dfrac{(1-\alpha)\nu_2(C)}{\alpha\nu_2(\gamma'_1)\nu_1(C)}\left( \dfrac{\alpha^2\nu^2_1(C)\nu_2(\gamma'_1)}{\alpha\nu_1(C)\nu_2(\gamma_1')+b\nu_2^2(C)}\right)\\
	= & \dfrac{\alpha(1-\alpha)\nu_1(C)\nu_2(C)}{\alpha\nu_1(C)\nu_2(\gamma'_1) + b\nu_2^2(C)}.
	\end{align*}

	Similarly,
	$$\overline{B}=\dfrac{-\alpha^2\nu_1^2(C)\nu_2(\gamma'_1)}{\nu_1(\gamma'_2)(\alpha\nu_2(\gamma'_1)\nu_1(C)+b\nu^2_2(C))}.$$

\end{proof}
\begin{remark}
	\
	\begin{enumerate}
		\item  If the curve $\gamma$ is parameterized by the affine arc length, the envelope of the intermediate lines is given by
	    \begin{align*}
		 X - M_\alpha & = \frac{ \alpha\lambda[\gamma_2^{\prime}, C ] [\gamma_2^{\prime}, \gamma_1^{\prime}]      }{     \alpha\lambda[\gamma_2^{\prime}, \gamma_1^{\prime}]^2 + [\gamma_2^{\prime}, C]}\left( (1-\alpha)\gamma_1^{\prime} -\alpha\lambda\gamma_2^{\prime}\right).
		\end{align*}
		\item Notice that when $\alpha=1/2$, we have the familiar set Affine Envelope Symmetry Set. More precisely,

		$$ X - M = \frac{1}{2}\left( \frac{ [\gamma_2^\prime,\gamma_1-\gamma_2] [\gamma_1^\prime,\gamma_2^{\prime}] }{     -2[\gamma_2^\prime,\gamma_1-\gamma_2] + [\gamma_1^\prime,\gamma_2^{\prime}]^2}\right) (\gamma_1^{\prime} -\gamma_2^{\prime}),$$
		for further details about this set, see \cite{giblin}.
	\end{enumerate}
\end{remark}
\section{Local structure of the envelope of Intermediate lines}\label{sec:local}
Let $\gamma$ be a $C^{\infty}$ closed plane curve. Also let $p_1=\gamma(t)$ and $p_2=\gamma(s)$ be two different points of $\gamma$ with $\kappa(p_1)\neq0$ and $\kappa^\prime(p_1)\neq0$, where $\kappa$ is the curvature at the point $p_1$. This means that $p_1$ is neither an inflection point nor a vertex point. Without loss of generality, we can assume that
\begin{equation}\label{P1andP2}
\begin{array}{ccl}
p_1 & = &(t,f(t))=\left(t,\frac{1}{2}t^2+a_3t^3+a_4 t^4+a_5 t^5 + O_6(t)\right),\\
p_2 & =& (s,g(s))=\left(s, b_0+b_1 s+b_2 s^2 +b_3s^3+b_4s^4+b_5 s^5 +O_6(s)\right),
\end{array}
\end{equation}
where $b_0>0$, $a_3>0$, $b_1^2-b_0>0$ and $O_6(t)$ (resp. $O_6(s)$) represents the terms of order grater than or equal to 6 with respect to the parameter $t$ (resp. $s$).
The tangent lines at the points $p_1$ and $p_2$ can be parallel or concurrent.

Let $F(x,y,t,s)=l_1(t,s) x + l_2(t,s) y+l_3(t,s)$ be the equation of the intermediate line passing through the point $M_{\alpha}$ given in (\ref{intermediate}). The condition of existence of solution for the envelope of the family of such lines is that the determinant of the following matrix be zero,
\begin{equation}\label{Matrix}
M=\left[
\begin{array}{ccc}
F_x(0,0,t,s) & F_y(0,0,t,s) & F(0,0,t,s) \\
F_{sx}(0,0,t,s) & F_{sy}(0,0,t,s) & F_s(0,0,t,s) \\
F_{tx}(0,0,t,s) & F_{ty}(0,0,t,s) & F_t(0,0,t,s)
\end{array}\right].
\end{equation}
The determinant of the matrix (\ref{Matrix}) is
\begin{align}\label{detM} \det(M)= \frac{1}{8} \alpha (\alpha -1) A(t,s) B(t,s),\end{align}
where $A(0,0)=0$ and $B(0,0)=8(2\alpha b_2-\alpha+1)b_0^3$ and we have
\begin{align}\label{A}
	A(t,s)= t-2b_2 s + 3 a_3 t^2-3b_3 s^2 + O_3(t,s),
\end{align}
where $O_3(t,s)$ denotes terms with respect to the parameters $t$ and $s$ with the order greater than or equal to three.
\subsection{Concurrent tangent lines}
If the tangent lines at the points $p_1$ and $p_2$ are not parallel, then the determinant of the matrix (\ref{Matrix}) given in (\ref{detM}) is zero, if  $B(t,s)=0$ (note that in this case $A(t,s)\neq 0$). Since $B(0,0)=8(2\alpha b_2-\alpha+1)b_0^3$ and $b_0\neq 0$. the expression $B(t,s)$, locally, has zero around the origin, if $b_2=\frac{\alpha-1}{2\alpha}$. Using this assumption, one can write the parameter $t$ as a function of $s$, due to the implicit function theorem. To be more precise, we have
\begin{align}\label{t}
	t = \frac{2\alpha b_0 b_3+b_1(\alpha-1)}{(\alpha-1)(2a_3b_0+b_1)} s +O_2(s).
\end{align}
\begin{remark}
	Note that as the point $p_1$ is not a vertex point, then, without loss of generality, we can assume that $b_0, b_1\,\,\text{and}\,\,a_3>0$. Therefore, $2a_3b_0+b_1 \neq 0$.
\end{remark}
\begin{theorem}\label{theo:non_Paralel}
	For a plane curve $\gamma$, suppose that the points $p_1$ and $p_2$ given in (\ref{P1andP2}) have non-parallel tangents. We have:
	\begin{itemize}
		\item[i)] If neither the point $p_1$ nor $p_2$ are inflection points, the condition of regularity of the envelope of the family of intermediate lines (i.e. the curve $AEIL$) is
		\begin{align*}
			\alpha & \neq \frac{b_0}{b_1^2},1,\\
			b_3 & \neq \frac{(\alpha-1)(-6\alpha a_3b_0b_1^2+4\alpha a_3b_0^2-3\alpha b_1^3-2a_3b_0^2)}{2\alpha b_0(6\alpha a_3b_0b_1+3\alpha b_1^2+2\alpha b_0-b_0)},
		\end{align*}
		\item[ii)] If the point $p_1$ is an inflection point, then the curve $AEIL$ is regular at $s=0$ if $b_3\neq 0$, equivalently, if $p_2$ is an inflection point at the origin.
 	\end{itemize}
\end{theorem}
\begin{proof}
	i) Let the points $p_1$ and $p_2$ given in \eqref{P1andP2} not be inflection points. After substituting the parameter $t$, given in \eqref{t}, in the envelope of intermediate lines \eqref{intermediate}, we reach the following parametrization
		\begin{align}\label{F non-parallel}
			F(t,\alpha) = \left(\tilde{a}_0+\tilde{a}_1 s+O(s^2), \tilde{b}_0+\tilde{b}_1 s + O(s^2)\right),
		\end{align}
		where
		\begin{align*}
			\tilde{a}_0 & = -\frac{\alpha b_1 b_0}{\alpha b_1^2-b_0},\\
			\tilde{a}_1 & = -\dfrac{b_0 \delta}{(b_1^2\alpha^2-(b_1^2+b_0)\alpha+b_0)(2a_3b_0+b_1)(b_1^2\alpha-b_0)},\\
			\tilde{b}_0 &= -\frac{\alpha  b_0^2}{\alpha b_1^2-b_0},\\
			\tilde{b}_1 & = -\dfrac{b_1 \delta}{(b_1^2\alpha^2-(b_1^2+b_0)\alpha+b_0)(2a_3b_0+b_1)(b_1^2\alpha-b_0)},
		\end{align*}
		and $\delta$ is a constant term. Note that, as we assume that $b_1^2-b_0>0$, the denominator of $\tilde{a}_1$ and $\tilde{b}_1$ are not nulls for any real value of $\alpha\neq b_0/b_1^2$ and $\alpha\neq1$. The envelope \eqref{F non-parallel} is singular at $s=0$ if $\delta=0$. This is observed if:
		\[
		b_3 = \frac{\alpha-1}{2\alpha b_0}\frac{(-6\alpha a_3b_0b_1^2+4\alpha a_3b_0^2-3\alpha b_1^3-2a_3b_0^2)}{(6\alpha a_3b_0b_1+3\alpha b_1^2+2\alpha b_0-b_0)}.
		\]
		This completes the proof of item (i).

		ii) Now suppose that the curve $\gamma$ is non-oval and has inflection point at the point $p_1$. Therefore, we can assume that the points $p_1$ and $p_2$ have the following local parametrizations: $p_1=(t,a_3 t^3+O(t^4))$ and $p_2=(s,b_0+b_1 s+b_2 s^2+O(s^3))$, with $b_0b_1\neq 0$. In this case, the component that annuls in determinant of the matrix $M$ given in \eqref{Matrix} becomes
		$$B(t,s)=16 \alpha b_0^3 b_2+O(t,s).$$
		Hence, $B(0,0)=0$ if and only if $b_2=0$. So, by using the Implicit Function Theorem, we obtain $t$ as function of $s$. To be more precise, we have
		\[
		t = \dfrac{\alpha b_3}{(\alpha-1)a_3}s+O(s^2).
		\]
		Therefore the parametrization of the envelope \eqref{intermediate} is as follow
		\[
		F(s,\alpha) = \left(-\frac{b_0}{b_1}+\hat{a} s+O(s^2) , \hat{b} s + O(s^2) \right),
		\]
		where
		\[
		\hat{a} \hat{b} = \dfrac{36\alpha b_3^2b_0^4}{b_1^5(\alpha-1)^2}.
		\]
		This means that $F(s,\alpha)$ is singular at the origin, if $b_3=0$. Equivalently, $F(s,\alpha)$ is regular at the origin, if $p_2$ is an inflection point (a schematic local behavior of $AEIL$ is given in Figure \ref{fig:InfNP}).
\end{proof}
\begin{figure}[h]
	\includegraphics[scale=0.4]{InfNP.pdf}
	\caption{The $AEIL$ (red curve) of a non-oval curve $\gamma$ (black curve), when $p_1$ and $p_2$ are inflection points.}\label{fig:InfNP}
\end{figure}
\subsection{Non-coincident parallels}
Suppose that the points $p_1$ and $p_2$, given in \eqref{P1andP2}, have paralell tangent lines. Hence, we have $b_1=0$. Also, the equation $A(t,s)=0$ given in \eqref{A}, gives the relation between the parameters $t$ and $s$ in this case.

Using the Implicit Function Theorem for the expression (\ref{A}), one can obtain $t$ as a function of $s$, as follows
\[ t = 2b_2 s + (-12a_3b_2^2+3b_3) s^2+ O_3(s).
\]

\begin{proposition}\label{prop:parallel}
	Let $\gamma$ be a smooth plane curve and $p_1$ and $p_2$, two points of $\gamma$, as given in (\ref{P1andP2}). If the tangent lines of $\gamma$ at the points $p_1$ and $p_2$ are parallel, then the $IPTL$ of the curve $\gamma$

	\begin{enumerate}
		\item passes through the point $M_\alpha(0,0)=(0,b_0\alpha)$.
		\item  is regular at $s=0$, if and only if $b_2\neq \dfrac{\alpha}{2(\alpha-1)}$.
		\item has an ordinary cusp singularity at $s=0$, if and only if
		$$	b_2 = \frac{\alpha}{2(\alpha-1)} \,\,\text{and}\,\, b_3\neq \left(\dfrac{\alpha}{\alpha-1}\right)^2 a_3.$$
		\item has a $(3,4)$-cusp singularity at $s=0$ if and only if
		$$	b_2 = \frac{\alpha}{2(\alpha-1)} \,\,,\,\, b_3 = \left(\dfrac{\alpha}{\alpha-1}\right)^2 a_3\,\,\,\text{and}\,\,\,
		b_4 \neq \left(\dfrac{\alpha}{\alpha-1}\right)^3 a_4.
		$$
	\end{enumerate}
\end{proposition}
\begin{proof}
	The proof comes from a straightforward calculation since a plane curve $\theta$ is regular at $t=0$ if $\theta'(0)\neq 0$ and has an ordinary cusp (resp. a $(3,4)$-cusp), singularity at the origin if $\theta'(0)=0$ and $[\theta^{(2)}(0),\theta^{(3)}(0)]\neq 0$ (resp. if $\theta'(0)=[\theta^{(2)}(0),\theta^{(3)}(0)]=0$ and $[\theta^{(3)}(0),\theta^{(4)}(0)]\neq 0$), where $[\cdot,\cdot]$ stands for the determinant of the $2\times 2$ matrix defined by $\mathbb{R}^2$ vectors.
\end{proof}

In the next proposition, we investigate the case when there is an inflection point. Let $p_1=(t,a_3 t^3+O(t^4))$ and $p_2$ is the same as in \eqref{P1andP2}.

\begin{proposition}\label{prop:parallel_p1_Inflection}
	By the above assumptions,
	if one of the points $p_1$ and $p_2$ or both of them have inflection at the origin, then the curve $IPTL$ passes through the point $M_{\alpha}$, and is regular at the origin, when $\alpha\in (0,1)$.

 Moreover, for $b_2\neq 0$, the curve $\gamma$ has an inflection of order $k$ at $p_1$, if and only if $IPTL$ has an inflection of order $k$ at the origin, for $k=1,2$.
\end{proposition}
\begin{proof}
	 Firstly, suppose that $p_2$ is not an inflection at $s=0$ i.e. $b_2\neq 0$.  Using the expression \eqref{A}, and by a straightforward calculation, one can write $s$ as a function of $t$. Therefore, we have the parametrization of the curve $IPTL$, as follows
	\[
	\left((1-\alpha) t + O(t^2), \alpha b_0 +a_3(1-\alpha) t^3+\frac{9\alpha a_3^2+4a_4 b_2(1-\alpha)}{b_2} t^4+O(t^5)\right).
	\]
	Now, suppose that the curve $\gamma$ has an inflection at the point $p_2$ as well. Therefore, using the relation $\nu_1(C)=-\beta\nu_2(C)$ where $\beta=\left((1-\alpha)/\alpha\right)^\frac{1}{3}$ (see Theorem \ref{pre}), we obtain:
	\[
	s= - \frac{a_3}{b_3}\left(\frac{\alpha-1}{\alpha}\right) t+O(t^2).
	\]
By replacing this in \eqref{intermediate}, we obtain the parametrization of the envelope of the intermediate lines
\[
F(t,\alpha)= \left(\frac{h_1(\alpha)}{2\alpha b_3} t+O(t^2) , \alpha b_0 + \frac{h_2(\alpha)}{2\alpha^2 b_3^2} t^3 + O(t^4) \right),
\]
where $h_1$ and $h_2$ are smooth functions with respect to $\alpha$ and their $1$-jets are:
\begin{align*}
	j^1 h_1(\alpha) = a_3-2(2a_3+b_3)\alpha\quad,\quad	j^1 h_2(\alpha) = a_3^3-3a_3^2(2a_3-b_3)\alpha.
\end{align*}
Figure \ref{fig:InfP} illustrates the $IPTL$ of the curve $\gamma$, when $p_1$ and $p_2$ are inflections with parallel tangents.
\end{proof}
\begin{figure}[tp]
	\includegraphics[scale=0.4]{InfP.pdf}
	\caption{The $IPTL$ (red curve) of the curve $\gamma$ (black curve) when $p_1$ and $p_2$ are inflections with parallel tangent lines.}\label{fig:InfP}
\end{figure}

\begin{remark}

\
\noindent {\bf (1)} Notice that when $\alpha=1/2$, the sets $AEIL=AESS$ and $IPTL=MPTL$ meet at their singularities (see Proposition 2.4.9 of \cite{Holtom}). For $\alpha\neq1/2$, these sets are disjoint.
\medskip

\noindent {\bf (2)} It is possible to determine the first $\alpha=\alpha_0$ for which the cusps appear ($\alpha$-born) in $IPTL$ or $AEIL$. Indeed, in the case of the $IPTL$ (resp. $AEIL$), one can find the relation between the parameters, as mentioned in this section. Then, according to the regularity conditions given in Proposition \ref{prop:parallel} (resp. Theorem \ref{theo:non_Paralel}), the $\alpha$-born can be found effortlessly.
\end{remark}

\subsection{Coincident tangent lines }
Here, we study the limit case of the envelope of intermediate lines for the smooth plane curve $\gamma$. This means that we desire to investigate the behavior of the intermediate lines when one point tends to the other point. For such, consider the local parametrization of the curve $\gamma$ at each point, i.e. let
$$\gamma_1(t)=(t, f(t)) \text{  \space and \space }  \gamma_2(s)=(s, f(s)). $$

Let  $C=(t-s, f(t)-f(s)), M_\alpha = \left((1-\alpha)t+\alpha s, (1-\alpha)f(t) + \alpha f(s)\right), N_1(t)=(-f^{\prime}(t), 1)$  and $N_2(s)=(-f^{\prime}(s),1),$ be the chord, the intermediate point and the Euclidean normal vectors at $\gamma_1$ and $\gamma_2$, respectively.

Note that  $\langle N_i, U \rangle=[\gamma_i^{\prime}, U]$, for $U \in \mathbb{R}^2$ and $i=1,2$, where $\langle\cdot,\cdot\rangle$ denotes the usual inner product of $\mathbb{R}^2$. Therefore we have
\begin{align}\label{con2}
\nu_1(U)=\frac{[\gamma_1^{\prime}, U]}{[\gamma_1^{\prime}, \gamma_1^{\prime\prime}]^{\frac{1}{3}}}=\frac{\langle N_1, U \rangle}{[\gamma_1^{\prime}, \gamma_1^{\prime\prime}]^{\frac{1}{3}}}\text{\space and \space} \nu_2(U)=\frac{[\gamma_2^{\prime}, U]}{[\gamma_2^{\prime}, \gamma_2^{\prime\prime}]^{\frac{1}{3}}}=\frac{\langle N_2, U \rangle}{[\gamma_2^{\prime}, \gamma_2^{\prime\prime}]^{\frac{1}{3}}}.
\end{align}

Thus, if we denote $\langle N_i, U \rangle = N_i(U)$, then the equation of the intermediate line becomes
$$(1-\alpha)N_2(C)N_1(X-M_\alpha)+\alpha N_1(C)N_2(X-M_\alpha)=0.  $$

From Theorem \ref{pre}, we obtain

$$\frac{N_1(C)  }{(f^{\prime\prime}(t_1))^{\frac{1}{3 }}}=-\lambda\frac{N_2(C)  }{(f^{\prime\prime}(t_2))^{\frac{1}{3 }}} \Rightarrow \frac{N_1(C)}{N_2(C)}=-\left(\frac{(1-\alpha)f^{\prime\prime}(t_1)}{\alpha f^{\prime\prime}(t_2)} \right)^{\frac{1}{3}}.   $$

Thus, the angular coefficient $A=A(s,t)$ of the intermediate line is given by

\begin{align*}
	  A(s,t)= \frac{(1-\alpha)f^{\prime}(t)(\alpha f^{\prime\prime}(s))^{\frac{1}{3}} - \alpha f^{\prime}(s)((1-\alpha)f^{\prime\prime}(t))^{\frac{1}{3}}                    }{(1-\alpha)(\alpha f^{\prime\prime}(s))^{\frac{1}{3}}-   \alpha((1-\alpha)f^{\prime\prime}(t))^{\frac{1}{3}}                        }.
\end{align*}
Therefore, for $\alpha=1/2$ we have
\begin{align*}
\lim_{s\longrightarrow t}A(s,t)= \frac{    \frac{1}{3}f'(t)(f''(t))^{-\frac{2}{3}}f'''(t) - f''(t)^\frac{4}{3} }{\frac{1}{3}(f''(t))^{-\frac{2}{3}}f'''(t)},
\end{align*}
which is the affine normal of $\gamma$ at the point $p_1$.

For $\alpha\neq 1/2$, we have:
\begin{align*}
\lim_{s\longrightarrow t}A(s,t)= \frac{f'(t)(f''(t))^\frac{1}{3}((1-\alpha)\alpha^\frac{1}{3} -\alpha(1-\alpha)^\frac{1}{3})    }{(f''(t))^\frac{1}{3}((1-\alpha)\alpha^\frac{1}{3} -\alpha(1-\alpha)^\frac{1}{3})}=f'(t).
\end{align*}
Therefore, we have the following result:
\begin{theorem}
	Let $\gamma$ be a smooth plane curve, $p_1=\gamma_1(t)=(t,f(t))$ and $p_2=\gamma_2(s)=(s,f(s))$ two points of $\gamma$.
	When the parameter $s$ tends to $t$, then the intermediate line ($\alpha\neq 1/2$) tends to the tangent line of $\gamma$, at the point $p_1$.
\end{theorem}

\begin{figure}[htp]
\includegraphics[scale=0.4]{example.pdf}
\caption{A schematic example of Envelope of Intermediate Lines of the bean curve $\gamma(t)=\left(0.1\cos(2\pi t)+\cos(\pi t),0.1\sin(2\pi t+1)+\sin(\pi t)\right)$. The left figure illustrates the envelope when $\alpha=0.5$ where $MPTL$ (red), $AESS$ (blue) and Affine Evolute (green) appear. The right figure illustrates the envelope for $\alpha=0.6$, where $IPTL$ (red) and $AEIL$ (blue) appear.}
\end{figure}

\subsection{Singularity theory approach}
 Our object in this section is not to study a single value of $\alpha$, but to study what happens to the envelope, as  $\alpha$ varies in $(0,1)\setminus\{1/2\}$.

Here, we shall state the results from the singularity theory, which allow us to make precise statements about the way in which the envelope of the family of intermediate lines evolves as $\alpha$ changes.
\begin{definition}[\cite{Giblin3}]
	For $X=(x,y,\alpha)$ in the discriminant surface \eqref{discriminant}, the function $f(t) = F(X, t)$ has
	\begin{enumerate}
\item type  $A_2$ at $t = t_0$ if $f^\prime(t_0) = f^{\prime\prime}(t_0) = 0, f^{\prime\prime\prime}(t_0) \neq 0$, \\
\item type $A_3$ at $t = t_0$ if $f^\prime(t_0) =  f^{\prime\prime}(t_0) =f^{\prime\prime\prime}(t_0)=0, f^{(4)}(t_0) \neq 0$.
	\end{enumerate}
\end{definition}
\begin{remark}
	The conditions given in Propositions \ref{prop:parallel},\ref{prop:parallel_p1_Inflection} and Theorem \ref{theo:non_Paralel} are the precise conditions for which the envelope has $A_2$ or $A_3$ singularity at the origin.
\end{remark}
The additional required condition, besides $A_2$ or $A_3$ singularity, is presented bellow.

\begin{definition}[\cite{giblin4}]
	Suppose that $t=t(s)$ and $F=F_s=0$ at $(X_0, s_0)=(x_0,y_0,\alpha_0,s_0)$. Also suppose that $f(s) = F(X_0, s)$ has an $A_r$-singularity at $s_0$. Consider the partial derivatives $F_{x}$, $F_{y}$ and $F_{\alpha}$ with respect to the parameters $x_i$, evaluated at $X_0$,
	their Taylor polynomials $\mathcal{T}_i$ up to degree $r-1$, expanded about $s_0$ (so these have $r$ terms). The family $F(X, s)$ is called a versal unfolding of $f$ at $s_0$ if the $\mathcal{T}_i$ span a vector space of dimension $r$. Thus, if the coefficients in the $\mathcal{T}_i$ are placed as the columns of an $(r-1) \times 3$ matrix, the rank
	is $r-1$. Obviously, this is possible only for $r \leq 3$.
\end{definition}
In the next result, we determine the versality condition of the family of the Affine Envelope of Intermediate Lines (i.e. $AEIL$) for an $A_2$-singularity at the origin. Analogously, one can find the versality conditions for the other cases.
\begin{theorem} Let $\gamma$ be a closed plane curve, $p_1=\gamma(t)$ and $p_2=\gamma(s)$ two distintc points of $\gamma$ with concurrent tangent lines which have local parametrizations given in \eqref{P1andP2}. Consider $F$ as the family of intermediate lines of $\gamma$, given in \eqref{intermediate} . Suppose that $f(s)=F(x_0,y_0,\alpha_0,s)$ has an $A_2$-singularity at the origin. Then, $F$ is versal, if and only if
		\[
		a_3=-\frac{1}{6}\frac{5\alpha-1}{\alpha^2b_1}.
		\]

\end{theorem}
\begin{proof} Using the relation between the parameters given in \eqref{t} and the conditions of $A_2$ singularity, as mentioned in Theorem \ref{theo:non_Paralel}, we have
$$F(x,y,\alpha,s)=\alpha^2b_1^2(-\alpha^2b_1^2-b_1x+y)+(1-\alpha)\alpha b_1^2(-\alpha^2b_1^2+y) + O(s).$$

Notice that

\begin{eqnarray*}
j^1\left(\dfrac{\partial F}{\partial x}(X_0,s)\right)(0) & = & m_{11}s  \\
j^1\left(\dfrac{\partial F}{\partial y}(X_0,s)\right)(0) & = & m_{21}s \\
j^1\left(\dfrac{\partial F}{\partial \alpha}(X_0,s)\right)(0) & = & m_{31}s,  \\
\end{eqnarray*}
where $X_0=(x_0,y_0,\alpha_0)=(-\alpha b_1,0,\alpha)$ and

\begin{eqnarray*}
m_{11} & = & -\dfrac{\alpha b_1^2(6\alpha^2a_3b_1+5\alpha-1)}{2(3\alpha a_3b_1+\alpha+1)} \\
m_{12} & = & \dfrac{b_1(6\alpha^2a_3b_1+5\alpha-1)}{2(3\alpha a_3b_1+\alpha+1)} \\
m_{13} & = & \dfrac{\alpha b_1^3(36\alpha^3a_3^2b_1^2(\alpha-1)+12\alpha^3a_3b_1+9\alpha^2a_3b_1-24\alpha a_3b_1+\alpha^2+2\alpha-5)}{2(3\alpha a_3b_1+\alpha+1)}.
\end{eqnarray*}

The matrix $M=(m_{11},m_{12}, m_{13})$ has rank 1, if and only if $a_3=-\frac{1}{6}\frac{5\alpha-1}{\alpha^2b_1}.$
\end{proof}

Now let $X$ denote the parameters $(x,y,\alpha)$ and $X_0$, a fixed value of $X$. Suppose that the family of intermediate lines, $F$, be a versal family of $f(t)=F(X_0,t)$ where $f(t)$ has type $A_k$, $k=2$ or $3$. It is true that the discriminant of $F$, i.e. $\mathcal{D}_F$ given in (\ref{discriminant}), in a small neighborhood  of $X_0$, is locally diffeomorphic to the standard discriminant of a family of an $A_k$-singularity (details about these concepts are found in various places such as \cite{giblin4}). Figure \ref{fig:versal} illustrates the standard transitions of co-dimension one, which are called lips, beaks and swallowtail transitions.

\begin{figure}[htp]
	\includegraphics[scale=0.5]{versal.pdf}
	\caption{The standard transitions of cusps, beaks, lips (up left to right, respectively) and swallowtail (down) of a versal family $G=G(x_1,x_2,x_3,t)$ of a function of type $A_k$-singularity, $k=2$ or $3$. This Figure is given in \cite{Giblin3} page 7.}\label{fig:versal}
\end{figure}
\medskip
\newpage
\bibliographystyle{amsplain}


On proportionally modular numerical semigroups that are generated by arithmetic progressions

Abstract: A numerical semigroup is a submonoid of ${\mathbb Z}_{\ge 0}$ whose complement in ${\mathbb Z}_{\ge 0}$ is finite. For any set of positive integers $a,b,c$, the numerical semigroup $S(a,b,c)$ formed by the set of solutions of the inequality $ax \bmod{b} \le cx$ is said to be proportionally modular. For any interval $[\alpha,\beta]$, $S\big([\alpha,\beta]\big)$ is the submonoid of ${\mathbb Z}_{\ge 0}$ obtained by intersecting the submonoid of ${\mathbb Q}_{\ge 0}$ generated by $[\alpha,\beta]$ with ${\mathbb Z}_{\ge 0}$.For the numerical semigroup $S$ generated by a given arithmetic progression, we characterize $a,b,c$ and $\alpha,\beta$ such that both $S(a,b,c)$ and $S\big([\alpha,\beta]\big)$ equal $S$.

\section{Introduction} \label{intro}
\vskip 10pt

A numerical semigroup $S$ is a submonoid of ${\mathbb Z}_{\ge 0}$ whose complement ${\mathbb Z}_{\ge 0} \setminus S$ is finite. For the complement to be finite, it is necessary and sufficient that $\gcd(S)=1$. For a given subset $A$ of positive integers, we write
\[ \langle A \rangle = \big\{ a_1x_1+\cdots+a_kx_k: a_i \in A, x_i \in {\mathbb Z}_{\ge 0}, k \in \mathbb N \big\} \:\bigcup\: \{0\}. \]
Note that $\langle A \rangle$ is a submonoid of ${\mathbb Z}_{\ge 0}$, and that $S=\langle A \rangle$ is a numerical semigroup if and only if $\gcd(A)=1$.

For a numerical semigroup $S$, the complement ${\mathbb Z}_{\ge 0} \setminus S$ is denoted by $\text{G}(S)$, and called the gap set of $S$. The largest element in $\text{G}(S)$ is called the Frobenius number of $S$, and denoted by $\text{F}(S)$. The cardinality of $\text{G}(S)$ is called the genus of $S$, and denoted by $g(S)$.

We say that $A$ is a set of generators of the semigroup $S$, or that the semigroup $S$ is generated by the set $A$, when $S=\langle A \rangle$. Further, $A$ is a minimal set of generators for $S$ if $A$ is a set of generators of $S$ and no proper subset of $A$ generates $S$. Every semigroup has a unique minimal set of generators. The embedding dimension $e(S)$ of $S$ is the size of the minimal set of generators.

Given a subset $A$ of ${\mathbb Q}_{\ge 0}$, we denote by $\langle A \rangle$ the submonoid of ${\mathbb Q}_{\ge 0}$ generated by $A$:
\[ \langle A \rangle = \big\{ a_1x_1+\cdots+a_kx_k: a_i \in A, x_i \in {\mathbb Z}_{\ge 0}, k \in \mathbb N \big\} \:\bigcup\: \{0\}. \]
Then $S(A)=\langle A \rangle \cap {\mathbb Z}_{\ge 0}$ is a submonoid of ${\mathbb Z}_{\ge 0}$. We say that $S(A)$ is the numerical semigroup associated to $A$. When $A=I$ is an interval, there is a simple way to describe $\langle A \rangle$.
\vskip 5pt

\begin{prop} {\bf (\cite[Lemma 1, pp. 282]{RGGU03})} \label{basic1}
Let $x \in {\mathbb Q}^+$ and let $I$ be an interval. Then $x \in \langle I \rangle$ if and only if there exists a positive integer $n$ such that $\frac{x}{n} \in I$.
\end{prop}
\vskip 5pt

\noindent Toms \cite{Tom03} in an attempt to solve problems in classification theory in $C^{\star}$-algebras defined Toms decomposable numerical semigroups. The concept of proportionally modular numerical semigroups was introduced by Rosales {\it et al\/} in \cite{RGU05}. Rosales \& Garcia-Sanchez \cite{RG-S08} showed that a numerical semigroup is Toms decomposable if and only if it is the intersection of finitely many proportionally modular numerical semigroups.
\vskip 5pt

\noindent A proportionally modular Diophantine inequality is an expression of the form $ax \bmod{b} \le cx$, where $a,b,c$ are positive integers. The set of integer solutions of a proportionally modular Diophantine inequality form a numerical semigroup. For $a,b,c \in \mathbb N$, define
\[ S(a,b,c) = \{ x \in {\mathbb Z}_{\ge 0}: ax \bmod{b} \le cx \}. \]
A numerical semigroup of this form is called proportionally modular. Since the inequality $ax \bmod{b} \le cx$ has the same set of integer solutions as $(a \bmod{b})x \bmod{b} \le cx$, we may assume that $a \in \{1,\ldots,b\}$. We note that $S(b,b,c)={\mathbb Z}_{\ge 0}$.
\vskip 5pt

\noindent Given a proportionally modular numerical semigroup, there exists a dual proportionally modular numerical semigroup.

\begin{prop} {\bf (\cite[Proposition 1, pp. 416]{DR06})} \label{basic2}
Let $a,b,c \in \mathbb N$ with $c<a<b$. Then
\[ S(a,b,c) = S(b+c-a,b,c). \]
\end{prop}
\vskip 5pt

We note that $S\big([\alpha,\beta]\big)={\mathbb Z}_{\ge 0}$ if and only if $\alpha \le \frac{1}{n} \le \beta$ for some $n \in \mathbb N$, and that $S(a,b,c)={\mathbb Z}_{\ge 0}$ if and only if $a \bmod{b} \le c$. The following two Propositions give a connection between $S(a,b,c)$ and $S\big([\alpha,\beta]\big)$.
\vskip 5pt

\begin{prop} {\bf (\cite[Lemma 1, pp. 454]{RGU08})} \label{basic3}
Let $a,b,c \in \mathbb N$ with $c<a$. Then
\[ S(a,b,c) = S\left(\left[\textstyle\frac{b}{a},\textstyle\frac{b}{a-c}\right]\right). \]
\end{prop}

\begin{prop} {\bf (\cite[Lemma 1, pp. 454]{RGU08})} \label{basic4}
Let $\frac{p}{q},\frac{r}{s} \in \mathbb Q$ with $\frac{p}{q}<\frac{r}{s}$. Then
\[ S\left(\left[\textstyle\frac{p}{q},\textstyle\frac{r}{s}\right]\right) = S(qr,pr,qr-ps). \]
\end{prop}
\vskip 5pt

Let $a,b,c \in \mathbb N$ with $c<a<b$. Using Propositions \ref{basic2} and \ref{basic3}, there exist intervals $[\alpha,\beta]$ and $[{\alpha}^{\prime},{\beta}^{\prime}]$, with rational endpoints, and with $\alpha>1$ and ${\alpha}^{\prime}>1$, such that
\[ S\left([\alpha,\beta]\right) = S\left([{\alpha}^{\prime},{\beta}^{\prime}]\right). \]
Moreover, it easily follows from Propositions \ref{basic2} and \ref{basic3} that
\begin{equation} \label{basic5}
{\alpha}^{\prime} = \frac{\beta}{\beta-1}, \quad {\beta}^{\prime} = \frac{\alpha}{\alpha-1}.
\end{equation}
\vskip 5pt

\noindent B\'{e}zout sequences are closely connected to the study of proportionately modular numerical semigroups; see \cite{RGU08}. A B\'{e}zout sequence is an increasing finite sequence of rational numbers $\left\{\frac{p_1}{q_1},\ldots,\frac{p_n}{q_n}\right\}$ with $\gcd(p_k,q_k)=1$ satisfying $p_{k+1} q_k-p_k q_{k+1}=1$ for $k \in \{1,\ldots,n-1\}$.  A B\'{e}zout sequence is said to be {\em proper} if $p_{k+\ell} q_k-p_{k+\ell} q_k>1$ for each $k$ and each $\ell>1$. As a consequence, we have
\begin{equation} \label{convex}
p_1 > \cdots > p_r < p_{r+1} < \cdots < p_n
\end{equation}
for some $r \in \{1,\ldots,n\}$; see \cite[Corollary 18, pp. 459]{RGU08}.

From the proof of Lemma 22 and Theorem 23 in \cite{RGU08}, given an interval $I=[\alpha,\beta]$, $\alpha,\beta \in \mathbb Q$ and $S=S(I)$, there exists a permutation of the generators $a_1,\ldots,a_e$ of $S$ and positive integers $b_1,\ldots,b_e$ such that
\begin{equation} \label{Bez1}
\left\{\frac{a_n}{b_n}\right\}_{n \ge 1} \:\text{forms a proper B\'{e}zout sequence, and}\:\alpha \le \frac{a_1}{b_1} < \frac{a_e}{b_e} \le \beta.
\end{equation}
\vskip 5pt

\noindent Given rational numbers $p_1/q_1$ and $p_n/q_n$ in reduced form, there exists a unique proper B\'{e}zout sequence with first term $p_1/q_1$ and last term $p_n/q_n$. This unique B\'{e}zout sequence has been determined algorithmically; see \cite[Algorithm 3.5]{BR09}.
\vskip 5pt

The following Proposition gives a characterization of proportionally modular numerical semigroups in terms of its minimal set of generators.
\vskip 5pt

\begin{prop} {\bf (\cite[Theorem 31, pp. 463]{RGU08})} \label{Bez2}
A numerical semigroup $S$ is proportionally modular if and only if there exists a permutation $a_1,\ldots,a_e$ of its minimal set of generators such that the following two conditions hold:
\begin{itemize}
\item[{\rm (i)}]
$\gcd(a_k,a_{k+1})=1$ for $1 \le k \le e-1$, and \\[-16pt]
\item[{\rm (ii)}]
$a_k \mid (a_{k-1}+a_{k+1})$ for $2 \le k \le e-1$.
\end{itemize}
\end{prop}
\vskip 5pt

\noindent Certain aspects of the numerical semigroup generated by an arithmetic progression $a,a+d,a+2d,\ldots,a+(k-1)d$ has been studied by many authors; see \cite{Bat58, Bra42, Gra73, NW72, Rob56, Tri94}. In particular, the gap set $\text{G}(S)$ and the Frobenius number $\text{F}(S)$ are determined by the following result in \cite{Tri94}. It can be shown that the embedding dimension $e(S)=\min\{a,k\}$.

\begin{prop} {\bf (\cite[Theorem 1, pp. 780]{Tri94})} \label{T}
Let $a,d,k$ be positive integers, with $\gcd(a,d)=1$ and $k \ge 2$. Let $S=\langle a,a+d,a+2d,\ldots,a+(k-1)d \rangle$. Then
\begin{itemize}
\item[{\rm (i)}]
$\text{G}(S) = \left\{ ax+dy: 0 \le x \le \left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor, 1 \le y \le a-1 \right\}$. \\[-12pt]
\item[{\rm (ii)}]
$\text{F}(S) = \max \text{G}(S) = a\left\lfloor \textstyle\frac{a-2}{k-1} \right\rfloor+d(a-1)$.
\end{itemize}
\end{prop}
\vskip 5pt

\noindent Given a numerical semigroup $S$, it is natural to ask for a characterization of positive integers $a,b,c$ for which $S(a,b,c)=S$, and also for a characterization of intervals $I=[\alpha,\beta]$ for which $S\big([\alpha,\beta]\big)=S$. We note that Propositions \ref{basic2} and \ref{basic3} connect the numerical semigroups $S(a,b,c)$ and $S\big([\alpha,\beta]\big)$, so characterizing one characterizes the other. Let $S=S\big(AP(a,d;k)\big)=\langle a,a+d,a+2d,\ldots,a+(k-1)d \rangle$, $k \ge 2$. Since $S\big(AP(a,d;k)\big)=S\big(AP(a,d;a)\big)$ whenever $k>a$, we may assume, without loss of generality, that $k \le a$.

In this paper, we characterize $a,b,c$ and $\alpha,\beta$ such that both $S(a,b,c)$ and $S\big([\alpha,\beta]\big)$ equal the numerical semigroup $S=S\big(AP(a,d;k)\big)$ generated by any arithmetic progression. The two characterizations appear as Theorems \ref{I} and \ref{II}, and are based on Propositions \ref{max}, \ref{min}, and \ref{link}. Delgado \& Rosales \cite{DR06} provide an algorithmic formula to determine $\text{F}\big(S(a,b,c)\big)$ that does not lead to an explicit formula. The determination of an explicit formula for $\text{F}\big(S(a,b,c)\big)$ therefore remains an open problem. Our characterization of $a,b,c$ for which $S=S\big(AP(a,d;k)\big)$ yields an explicit formula in these special cases since the problem of determining the Frobenius number for semigroups generated by arithmetic progressions is well known, as mentioned above. We note in passing that the algorithm in \cite[Algorithm 3.5]{BR09} to compute the unique proper B\'{e}zout sequence connecting two given reduced rational numbers $p_1/q_1$ and $p_n/q_n$ relies on knowing these rational numbers. Since there is no obvious relationship between $p_1/q_1$ and $p_n/q_n$, and the numerical semigroup $S$, these endpoints of the proper B\'{e}zout sequence has to be determined first. This means that the algorithm in \cite[Algorithm 3.5]{BR09} is not directly useful in resolving the problem.

\section{Main Results} \label{main}
\vskip 10pt

\noindent In this section, we consider an inverse problem related to two types of numerical semigroups that were introduced in Section \ref{intro} -- numerical semigroups generated by intervals $[\alpha,\beta]$ and proportionally modular numerical semigroups $S(a,b,c)$. Propositions \ref{basic2} and \ref{basic3} in Section \ref{intro} connect these two types of numerical semigroups.

In order that the semigroup generated by an arithmetic progression $a,a+d,a+2d,\ldots,a+(k-1)d$ be a numerical semigroup, it is necessary and sufficient that $\gcd(a,d)=1$. Henceforth we assume $\gcd(a,d)=1$, and write $S\big(AP(a,d;k)\big)=\langle a,a+d,a+2d,\ldots,a+(k-1)d \rangle$. The set $\{a,a+d,a+2d,\ldots,a+(k-1)d\}$ is the minimal set of generators for $S\big(AP(a,d;k)\big)$ when $2 \le k \le a$, whereas the set $\{a,a+d,a+2d,\ldots,a+(a-1)d\}$ is the minimal set of generators for $S\big(AP(a,d;k)\big)$ for $k>a$. Throughout this section, we may therefore assume without loss of generality that $2 \le k \le a$.
\vskip 5pt

\noindent We illustrate our main result, Theorem \ref{I}, with the following two examples. In Example \ref{SG1} we determine intervals $I$ for which $S(I)$ is the ordinary numerical semigroup $\{0,5,\rightarrow\}$, and in Example \ref{SG2} intervals $I$ for which $S(I)$ is the numerical semigroup $\langle 5,7 \rangle$ with embedding dimention two. These are both special cases of our result in Theorem \ref{I}.
\vskip 5pt

\begin{eg} \label{SG1}
We find $I=[\alpha,\beta]$, $\alpha>1$, for which $S(I)=\{0,5,\rightarrow\}$. From the argument leading up to eqn.~\eqref{Bez1}, we have $\{a_1,\ldots,a_5\}=\{5,\ldots,9\}$, giving rise to the following four B\'{e}zout sequences:
\begin{itemize}
\item $\frac{5}{1}$, $\frac{6}{1}$, $\frac{7}{1}$, $\frac{8}{1}$, $\frac{9}{1}$.
\item $\frac{9}{2}$, $\frac{5}{1}$, $\frac{6}{1}$, $\frac{7}{1}$, $\frac{8}{1}$.
\item $\frac{9}{8}$, $\frac{8}{7}$, $\frac{7}{6}$, $\frac{6}{5}$, $\frac{5}{4}$.
\item $\frac{8}{7}$, $\frac{7}{6}$, $\frac{6}{5}$, $\frac{5}{4}$, $\frac{9}{7}$.
\end{itemize}
This leads to
\[ \alpha \le \tfrac{5}{1} < \tfrac{9}{1} \le \beta \;\;\text{or}\;\; \alpha \le \tfrac{9}{8} < \tfrac{5}{4} \le \beta \;\;\text{or}\;\; \alpha \le \tfrac{9}{2} < \tfrac{8}{1} \le \beta \;\;\text{or}\;\; \alpha \le \tfrac{8}{7} < \tfrac{9}{7} \le \beta. \]
Note that the pairs of intervals $[5,9]$, $[\tfrac{9}{8},\tfrac{5}{4}]$ and $[\tfrac{9}{2},8]$, $[\tfrac{8}{7},\tfrac{9}{7}]$ are connected by eqn.~\eqref{basic5}.
\vskip 5pt

\noindent From Proposition \ref{basic1}, $\frac{g}{n} \notin I$ whenever $g \in \text{G}(S)=\{1,2,3,4\}$ and $n \in \mathbb N$. This results in
\[ \alpha \in (4,5], \beta \in [9,\infty) \;\;\text{or}\;\; \alpha \in \left(1,\tfrac{9}{8}\right], \beta \in \left[\tfrac{5}{4},\tfrac{4}{3}\right) \;\;\text{or}\;\;
    \alpha \in \left(4,\tfrac{9}{2}\right], \beta \in [8,\infty) \;\;\text{or}\;\; \alpha \in \left(1,\tfrac{8}{7}\right], \beta \in \left[\tfrac{9}{7},\tfrac{4}{3}\right).
\]
Note that the pairs of intervals $(4,5]$, $\left[\tfrac{5}{4},\tfrac{4}{3}\right)$, $[9,\infty)$, $\left(1,\tfrac{9}{8}\right]$,  $\left(4,\tfrac{9}{2}\right]$, $\left[\tfrac{9}{7},\tfrac{4}{3}\right)$, and $[8,\infty)$, $\left(1,\tfrac{8}{7}\right]$ are connected by eqn.~\eqref{basic5}.

\begin{figure}[h!]
\incfig{intex1}
\caption{For $I$ such that $S(I)=\{0,5,\rightarrow\}$, $[5,9] \subset I \subset (4,\infty)$ or $[\frac{9}{2},8] \subset I \subset (4,\infty)$ or $[\frac{9}{7},\frac{5}{4}] \subset I \subset (1,\frac{4}{3})$ or $[\frac{8}{7},\frac{9}{7}] \subset I \subset (1,\frac{4}{3})$.}
\end{figure}
\end{eg}
\vskip 5pt

\begin{eg} \label{SG2}
From Proposition \ref{Bez2}, it follows that every numerical semigroup with embedding dimension $2$ is proportionally modular. We find $I=[\alpha,\beta]$, $\alpha>1$, for which $S(I)=\langle 5,7 \rangle$. From the argument leading up to eqn.~\eqref{Bez1}, we have $\{a_1,a_2\}=\{5,7\}$, so that $1<\frac{5}{b_1}<\frac{7}{b_2}$ and $7b_1-5b_2=1$ or $1<\frac{7}{b_1}<\frac{5}{b_2}$ and $5b_1-7b_2=1$. This leads to
\[ \alpha \le \tfrac{5}{3} < \tfrac{7}{4} \le \beta \;\;\text{or}\;\; \alpha \le \tfrac{7}{3} < \tfrac{5}{2} \le \beta. \]
Note that the intervals $[\tfrac{5}{3},\frac{7}{4}]$ and $[\tfrac{7}{3},\tfrac{5}{2}]$ are connected by eqn.~\eqref{basic5}.
\vskip 5pt

\noindent From Proposition \ref{basic1}, $\frac{g}{n} \notin I$ whenever $g \in \text{G}(S)=\{1,2,3,4,6,8,9,11,13,16,18,23\}$ and $n \in \mathbb N$. This results in
\[ \alpha \in \left(\tfrac{23}{14},\tfrac{5}{3}\right], \beta \in \left[\tfrac{7}{4},\tfrac{23}{13}\right) \quad \text{or} \quad
    \alpha \in \left(\tfrac{23}{10},\tfrac{7}{3}\right], \beta \in \left[\tfrac{5}{2},\tfrac{23}{9}\right).
\]
Note that the pairs of intervals $\left(\tfrac{23}{14},\tfrac{5}{3}\right]$, $\left[\tfrac{5}{2},\tfrac{23}{9}\right)$ and $\left[\tfrac{7}{4},\tfrac{23}{13}\right)$, $\left(\tfrac{23}{10},\tfrac{7}{3}\right]$ are connected by eqn.~\eqref{basic5}.

\begin{figure}[h!]
\incfig{intex2}
\caption{For $I$ such that $S(I)=\langle 5,7 \rangle$, $[\frac{5}{3},\frac{7}{4}] \subset I \subset (\frac{23}{14},\frac{23}{13})$ or $[\frac{7}{3},\frac{5}{2}] \subset I \subset (\frac{23}{10},\frac{23}{9})$.}
\end{figure}
\end{eg}
\vskip 5pt

The main result in this section is the characterization of $\alpha,\beta$ such that $S\big([\alpha,\beta]\big)=S\big(AP(a,d;k)\big)$ in Theorem \ref{I}. Essential to our proof of Theorem \ref{I} are three results listed as Propositions \ref{max}, \ref{min}, and \ref{link}. We also characterize $a,b,c$ such that $S(a,b,c)=S\big(AP(a,d;k)\big)$, via Proposition \ref{basic2} and Theorem \ref{I}. We also state, without proof, the two special cases $S\big(AP(a,1;k)\big)$ and $S\big(AP(a,1;a)\big)$.
\vskip 5pt

\begin{prop} \label{max}
Let $a,d,k$ be positive integers, with $\gcd(a,d)=1$ and $2 \le k \le a$. Let $S\big(AP(a,d;k)\big)=\langle a,a+d,a+2d,\ldots,a+(k-1)d \rangle$. If $dd_a=1+q_a a$, with $1 \le d_a<a$ and $0 \le q_a<d$, and $\lambda=\left\lfloor\frac{a-2}{k-1}\right\rfloor$, then
\[ \max \left\{ \frac{g}{\left\lceil \frac{d_a g}{a} \right\rceil}: g \in \text{G}(S) \right\} =  \frac{\lambda a+d(a-1)}{\lambda d_a+q_a(a-1)+1} =
     \frac{\text{F}(S)}{\left\lceil \frac{d_a\text{F}(S)}{a} \right\rceil},
\]
with the maximum achieved at $g=\lambda a+d(a-1)$.
\end{prop}

\begin{Pf}
If $g \in \text{G}(S)$, then $g=ax+dy$, with $0 \le x \le \left\lfloor \frac{y-1}{k-1} \right\rfloor$ and $1 \le y \le a-1$. Let $d_a$ be such that $dd_a \equiv 1\pmod{a}$, $d_a>0$, and write $dd_a=q_a a+1$, $q_a \in {\mathbb Z}_{\ge 0}$. Then
\[ \frac{d_a g}{a} = \frac{d_a(ax+dy)}{a} = d_a x + q_a y + \frac{y}{a}. \]
Since $1 \le y<a$,
\[ \frac{g}{\left\lceil \textstyle\frac{d_a g}{a} \right\rceil} = \frac{ax+dy}{d_a x+q_a y+1} \stackrel{\mbox{\tiny{def}}}{=} f_{\ell}(x,y). \]
We show that
\begin{equation} \label{lbasic1}
\max \left\{ f_{\ell}(x,y): 0 \le x \le \left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor, 1 \le y \le a-1 \right\} = f_{\ell}\left( \left\lfloor \textstyle\frac{a-2}{k-1} \right\rfloor, a-1\right) = \frac{\text{F}(S)}{\left\lceil \frac{d_a\text{F}(S)}{a} \right\rceil}.
\end{equation}
\vskip 5pt

\noindent Fix $y \in \{1,\ldots,a-1\}$. For $x \in \left\{0,1,2,\ldots,\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor -1\right\}$, the numerator of
\[ f_{\ell}(x+1,y) - f_{\ell}(x,y) = \frac{a(x+1)+dy}{d_a(x+1)+q_a y+1} - \frac{ax+dy}{d_a x+q_a y+1} \]
is
\[ a(d_a x+q_a y+1)-d_a(ax+dy) = a(q_a y+1)-y(q_a a+1) =  a-y > 0. \]
Therefore,
\begin{equation} \label{lbasic2}
\max \left\{ f_{\ell}(x,y): 0 \le x \le \left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor, 1 \le y \le a-1 \right\} = \max \left\{ f_{\ell}\left(\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor,y\right): 1 \le y \le a-1 \right\}.
\end{equation}
\vskip 5pt

\noindent For $y \in \{1,\ldots,a-2\}$, the numerator of
\[ f_{\ell}\left(\left\lfloor \textstyle\frac{y}{k-1} \right\rfloor,y+1\right) - f_{\ell}\left(\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor,y\right) =
    \frac{a\left\lfloor \textstyle\frac{y}{k-1} \right\rfloor+d(y+1)}{d_a\left\lfloor \textstyle\frac{y}{k-1} \right\rfloor+q_a(y+1)+1} -
    \frac{a\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor+dy}{d_a\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor+q_a y+1}
\]
is
\[ d\left(d_a\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor+q_a y+1\right) - q_a\left(a\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor+dy\right) =
    \left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor + d > 0 \:\: \text{if}\:\: (k-1) \nmid y,
\]
and
\[ (a+d)\left(d_a\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor+q_a y+1\right) - (d_a+q_a)\left(a\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor+dy\right) =
    \left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor - y + a + d > 0 \:\: \text{if}\:\: (k-1) \mid y.
\]
Therefore,
\begin{equation} \label{lbasic3}
\max \left\{ f_{\ell}\left(\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor,y\right): 1 \le y \le a-1 \right\} = f_{\ell}\left(\left\lfloor \textstyle\frac{a-2}{k-1} \right\rfloor,a-1\right).
\end{equation}
The claim in eqn.~\eqref{lbasic1} follows from eqn.~\eqref{lbasic2} and eqn.~\eqref{lbasic3}.
\end{Pf}

\begin{prop} \label{min}
Let $a,d,k$ be positive integers, with $\gcd(a,d)=1$ and $2 \le k \le a$. Let $S\big(AP(a,d;k)\big)=\langle a,a+d,a+2d,\ldots,a+(k-1)d \rangle$, and set $b=a+(k-1)d$. If $dd_b=1+q_b b$, with $1 \le d_b<b$ and $0 \le q_b<d$, and $a-2=\lambda(k-1)+\mu$, $0 \le \mu < k-1$, then
\[ \min \left\{ \frac{g}{\left\lfloor \frac{d_b g}{b}\right\rfloor}: g \in \text{G}(S), d_b g \ge b \right\} = \frac{\lambda b+d}{\lambda d_b+q_b} =
    \frac{\text{F}(S)-\mu d}{\left\lfloor \frac{d_b\big(\text{F}(S)-\mu d\big)}{b} \right\rfloor}, \]
with the minimum achieved at $g=\lambda a+d\big(\lambda(k-1)+1\big)=\lambda b+d$.
\end{prop}

\begin{Pf}
If $g \in \text{G}(S)$, then $g=ax+dy$, with $0 \le x \le \left\lfloor \frac{y-1}{k-1} \right\rfloor$ and $1 \le y \le a-1$. Let $d_b$ be such that $dd_b \equiv 1\pmod{b}$, $d_b>0$, and write $dd_b=q_b b+1$, $q_b \in {\mathbb Z}_{\ge 0}$. Then
\begin{eqnarray*}
\frac{d_b g}{b} & = & \frac{d_b x\big(b-(k-1)d\big)+d_b dy}{b} \\
& = & d_b x + \frac{d_b d}{b}\big(y-(k-1)x\big) \\
& = & d_b x + q_b\big(y-(k-1)x\big) + \frac{y-(k-1)x}{b}.
\end{eqnarray*}
Since $0<y-(k-1)x<a<b$, we have
\[ \frac{g}{\left\lfloor \frac{d_b g}{b}\right\rfloor} = \frac{ax+dy}{d_b x+q_b\big(y-(k-1)x\big)} \stackrel{\mbox{\tiny{def}}}{=} f_u(x,y). \]
We show that
\begin{equation} \label{ubasic1}
\min \left\{ f_u(x,y): 0 \le x \le \left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor, 1 \le y \le a-1 \right\} = f_u\left( \lambda, \lambda(k-1)+1 \right)
= \frac{\text{F}(S)-\mu d}{\left\lfloor \textstyle\frac{d^{-1}\big(\text{F}(S)-\mu d\big)}{b} \right\rfloor}.
\end{equation}
\vskip 5pt

\noindent Fix $y \in \{1,\ldots,a-1\}$. For $x \in \left\{0,1,2,\ldots,\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor -1\right\}$, the numerator of
\[f_u(x+1,y) - f_u(x,y) = \frac{a(x+1)+dy}{d_b(x+1)+q_b\big(y-(k-1)(x+1)\big)} - \frac{ax+dy}{d_b x+q_b\big(y-(k-1)x\big)} \]
is
\begin{eqnarray*}
a\left(d_b x+q_b\big(y-(k-1)x\big)\right)-\big(d_b-q_b(k-1)\big)(ax+dy) = y\left( aq_b-d(d_b-q_b(k-1)\big) \right) \\
= y \left( q_b\big(a+(k-1)d\big) - (q_b b+1) \right) = -y < 0.
\end{eqnarray*}
Therefore,
\begin{equation} \label{ubasic2}
\min \left\{ f_u(x,y): 0 \le x \le \left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor, 1 \le y \le a-1 \right\} = \min \left\{ f_u\left(\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor,y\right): 1 \le y \le a-1 \right\}.
\end{equation}
\vskip 5pt

\noindent Since $y-(k-1)\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor=1+\big(y-1 \!\!\pmod{k-1}\big)$, we may write
\begin{equation} \label{ubasic3}
f_u\left(\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor,y\right) = \frac{a\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor+dy}{d_b\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor + q_b\left(1+\big(y-1 \!\!\!\!\pmod{k-1}\big)\right)}.
\end{equation}
\vskip 5pt

\noindent Write $y-1=q(k-1)+r$, $0 \le r<k-1$. Then the minimum on the right-side of eqn.~\eqref{ubasic2} can be written as
\begin{equation} \label{ubasic4}
\min \left\{ \frac{qa+dy}{qd_b+q_b(r+1)}: 0 \le q \le \lambda -1, 0 \le r < k-1 \;\text{or}\; q=\lambda, 0 \le r \le \mu \right\}
\end{equation}
in view of eqn.~\eqref{ubasic3}.
\vskip 5pt

\noindent Fix $q$. For $0 \le q<\lambda$, let $0 \le r < k-1$, and for $q=\lambda$, let $0 \le r \le \mu$. The numerator of
\begin{eqnarray*}
f_u\left(\left\lfloor \textstyle\frac{y}{k-1} \right\rfloor,y+1\right) - f_u\left(\left\lfloor \textstyle\frac{y-1}{k-1} \right\rfloor,y\right) \\
= \frac{qa+d(y+1)}{qd_b+q_b(r+2)} - \frac{qa+dy}{qd_b+q_b(r+1)} = \frac{d\big(qd_b+q_b(r+1)\big) - q_b(qa+dy)}{\big(qd_b+q_b(r+1)\big) \big(qd_b+q_b(r+2) \big) }
\end{eqnarray*}
is
\[ q(q_b b-q_b a+1)-q_b d\big(y-(r+1)\big) = q + q\big(q_b(k-1)d\big) - q_b dq(k-1) = q > 0. \]
\vskip 5pt

\noindent So for fixed $q$, the minimum is achieved when $r=0$, and we now have
\begin{eqnarray} \label{ubasic5}
\min \left\{ \frac{qa+dy}{qd_b+q_b(r+1)}: 0 \le q \le \lambda -1, 0 \le r < k-1 \;\text{or}\; q=\lambda, 0 \le r \le \mu \right\} \nonumber \\
= \min \left\{ \frac{qb+d}{qd_b+q_b}: 0 \le q \le \lambda \right\}.
\end{eqnarray}
\vskip 5pt

\noindent Finally for $q \in \{0,1,2,\ldots,\lambda -1 \}$, the numerator of
\[ \frac{(q+1)b+d}{(q+1)d_b+q_b} - \frac{qb+d}{qd_b+q_b} = \frac{b(qd_b+q_b) - d_b(qb+d)}{(qd_b+q_b) \big((q+1)d_b+q_b\big)} \]
is
\[ q_b b - (q_b b+1) < 0. \]
\vskip 5pt

\noindent Therefore, the minimum is achieved at $q=\lambda$, and we have
\begin{equation} \label{ubasic6}
\min \left\{ \frac{qb+d}{qd_b+q_b}: 0 \le q \le \lambda \right\} = \frac{\lambda b+d}{\lambda d_b+q_b} = f_u\left( \lambda, \lambda(k-1)+1 \right).
\end{equation}
The claim in eqn.~\eqref{ubasic1} follows from eqns.~\eqref{ubasic2} - ~\eqref{ubasic6}.
\end{Pf}

\begin{prop} \label{link}
Let $a,d$ be positive integers, with $\gcd(a,d)=1$ and let $b=a+(a-1)d$. Let $S\big(AP(a,d;a)\big)=\langle a,a+d,a+2d,\ldots,a+(a-1)d\rangle$, and let
\[ dd_a = 1+q_a a, \quad dd_b = 1+q_b b, \quad aa_b \equiv 1 \!\!\!\!\pmod{b}, \quad dd^{\prime} \equiv 1 \!\!\!\!\pmod{b-d}, \]
where
\[ 1 \le d_a < a, \quad  0 \le q_a < d, \quad 1 \le d_b < b, \quad  0 \le q_b < d, \quad 1 \le a_b < b, \quad 1 \le d^{\prime} < b-d. \]
\begin{itemize}
\item[{\rm (i)}]
In terms of $d_a$ and $q_a$, we have
\[ q_b = q_a, \quad d_b = d_a+(a-1)q_a, \quad d^{\prime} = d_a+(a-2)q_a, \quad a_b = (a-1)q_a+d_a+1. \]
\item[{\rm (ii)}]
\[ \left\lceil\frac{d_a g}{a}\right\rceil = \left\lceil\frac{a_b g}{b}\right\rceil \;\;\text{and}\;\; \left\lfloor\frac{d_b g}{b}\right\rfloor = \left\lfloor \frac{d^{\prime} g}{b-d}
    \right\rfloor.
\]
for each $g \in \text{G}(S)$.
\end{itemize}
\end{prop}

\begin{Pf}
If $g \in \text{G}(S)$, then $g=dy$, with $1 \le y \le a-1$. From Propositions \ref{max} and \ref{min},
\begin{equation} \label{linkeq}
\left\lceil\frac{d_a g}{a}\right\rceil = q_a y + 1 \:\:\text{and}\:\: \left\lfloor\frac{d_b g}{b}\right\rfloor = q_b y.
\end{equation}

Multiplying both sides of $dd_a-q_a a=1$ by $a-1$, then adding and subtracting $ad_a$, and rearranging, we get
\[ a\big( (a-1)q_a + d_a + 1 \big) - bd_a = 1. \]
Thus, $a\big((a-1)q_a+d_a+1\big) \equiv 1\pmod{b}$ and $1 \le (a-1)q_a+d_a+1<b$, so that
\begin{equation}\label{a_b}
a_b = (a-1)q_a + d_a  + 1.
\end{equation}
Therefore
\[ \frac{a_b g}{b} = \frac{\big((a-1)q_a+d_a+1\big)dy}{b} = \frac{(b-a)q_a y+\big(1+q_a a\big)y+dy}{b} = q_a y + \frac{(d+1)y}{b}. \]
Since $0 \le (d+1)y \le (d+1)(a-1) < b$, we have $\left\lceil\frac{a_b g}{b}\right\rceil=q_a y+1=\left\lceil\frac{d_a g}{a}\right\rceil$ by eqn.~\eqref{linkeq}. This proves the first part of the Proposition.
\vskip 5pt

For $j \in \{0,\ldots,a-1\}$, we have $d\big(d_a+jq_a\big)-(a+jd)q_a=dd_a-q_a a=1$. Thus, $d\big(d_a+jq_a\big) \equiv 1\pmod{a+jd}$, and since $1 \le d_a+jq_a<a+jd$, we have $d_a+jq_a$ equals $d^{-1}\bmod{a+jd}$ for each $j \in \{0,\ldots,a-1\}$. In particular,
\begin{equation} \label{d_b}
d_b = d_a+(a-1)q_a.
\end{equation}
Multiplying both sides of eqn.~\eqref{d_b} by $d$ we get
\[ 1+q_b b = dd_b = d\big(d_a+(a-1)q_a\big) = (1+q_a a) + (a-1)dq_a = 1+q_a b. \]
Since $d^{\prime}$ is $d^{-1}\bmod{a+jd}$ for $j=a-2$, we have
\begin{equation}\label{q_b}
d^{\prime} = d_a + (a-2)q_a \:\:\text{and}\:\: q_a = q_b.
\end{equation}
Eqns.~\eqref{a_b}, \eqref{d_b}, \eqref{q_b} prove part (i).
\vskip 5pt

We now have
\[ \frac{d^{\prime} g}{b-d} = \frac{\big((a-2)q_a+d_a\big)dy}{b-d} = \frac{\big((b-d)-a\big)q_a y+\big(1+q_a a\big)y}{b-d} = q_b y + \frac{y}{b-d}. \]
Since $0 \le y < b-d$, we have $\left\lfloor\frac{d^{\prime} g}{b-d}\right\rfloor=q_b y=\left\lfloor\frac{d_b g}{b}\right\rfloor$ by eqn.~\eqref{linkeq}. This proves part (ii).
\end{Pf}
\vskip 5pt

\noindent Let $S=S\big(AP(a,d;k)\big)$ denote the numerical semigroup generated by $\text{AP}(a,d;k)$. Let $a,b,c$ be any positive integers. Since $S(a,b,c)=S(a\bmod{b},b,c)$ and $S(a,b,c)={\mathbb Z}_{\ge 0}$ for $c \ge a$, we may assume without loss of generality that $c<a<b$. Propositions \ref{basic3} and \ref{basic4} provide a one-one correspondence between $S(a,b,c)$ and $S([\alpha,\beta])$. Via this correspondence, each interval $[\alpha,\beta]$, with $\alpha \ge 0$, corresponds to a unique triple $(a,b,c)$, with $a \in \mathbb N$. The assumption $c<a<b$ on the positive integers $a,b,c$ gives rise to intervals $[\alpha,\beta]$ with $\alpha>1$. In Theorem \ref{I}, we determine all intervals $[\alpha,\beta]$ with $\alpha>1$ such that $S([\alpha,\beta])=S$. We then use Proposition \ref{basic3} to determine all triples $(a,b,c)$ with $c<a<b$ such that $S(a,b,c)=S$. To extend the result of Theorem \ref{I} to all intervals $[\alpha,\beta]$ with $\alpha \ge 0$, one may combine Proposition \ref{basic3} and the fact that $S(a,b,c)=S(a^{\prime},b,c)$ for $a \equiv a^{\prime}\pmod{b}$.

We characterize intervals $[\alpha,\beta]$, $\alpha>1$ in Theorem \ref{I} by calculating certain B\'{e}zout sequences related to $\text{AP}(a,d;k)$. Bullejos \& Rosales \cite{BR09} provide an algorithm to determine the B\'{e}zout sequence connecting any two rational numbers $p_1/q_1$ and $p_2/q_2$. However, it remains to determine $p_1/q_1$ and $p_2/q_2$ in order to characterize $[\alpha,\beta]$. We give a direct method to explicitly determine the desired B\'{e}zout sequences.
\vskip 5pt

\begin{thm} \label{I}
Let $a,d,k$ be positive integers, with $\gcd(a,d)=1$, $2 \le k \le a$, $b=a+(k-1)d$, and let $\lambda=\left\lfloor\frac{a-2}{k-1}\right\rfloor$. Define $dd_a=1+q_a a$,
where $1 \le d_a < a$ and $0 \le q_a < d$.
Let $S\big(AP(a,d;k)\big)=\langle a,a+d,a+2d,\ldots,a+(k-1)d \rangle$. Then for $\alpha,\beta \in [1,\infty) \cap \mathbb Q$,
\[ S\big([\alpha,\beta]\big) = S\big(AP(a,d;k)\big) \]
if and only if
\[ \frac{\lambda a+d(a-1)}{\lambda d_a+q_a(a-1)+1} < \alpha \le \frac{a}{d_a}, \quad \frac{b}{d_a+(a-1)q_a} \le \beta < \frac{\lambda b+d}{\lambda\big(d_a+(a-1)q_a\big)+q_a}, \]
or
\begin{eqnarray*}
\frac{\lambda b+d}{\lambda\big(b-d_a-(a-1)q_a\big)+(d-q_a)} < \alpha \le \frac{b}{b-d_a-(a-1)q_a},  \\
\frac{a}{a-d_a} \le \beta < \frac{\lambda a+d(a-1)}{\lambda(a-d_a)+(d-q_a)(a-1)-1}.
\end{eqnarray*}
\noindent Moreover, for $k=a$, we also have
\[ \frac{d(a-1)}{q_a(a-1)+1} < \alpha \le \frac{b}{(a-1)q_a+d_a+1}, \quad \frac{b-d}{(a-2)q_a+d_a} \le \beta < \frac{d}{q_a}, \]
or
\[ \frac{d}{d-q_a} < \alpha \le \frac{b-d}{b-d-(a-2)q_a-d_a}, \quad \frac{b}{b-(a-1)q_a-d_a-1} \le \beta < \frac{d(a-1)}{(d-q_a)(a-1)-1}, \]
\noindent We adopt the convention $\frac{1}{0}=\infty$.
\end{thm}

\begin{Pf}
Throughout this proof, we define
\[ dd_b = 1+q_b b, \quad aa_b \equiv 1 \!\!\!\!\pmod{b}, \quad dd^{\prime} \equiv 1 \!\!\!\!\pmod{b-d}, \]
with $1 \le d_b < b$, $0 \le q_b < d$, $1 \le a_b<b$, and $1 \le d^{\prime}<b-d$.
\vskip 5pt

Let $\alpha,\beta \in [1,\infty) \cap \mathbb Q$ be such that $S(I)=S(\text{AP}(a,d;k))=S$, with $I=[\alpha,\beta]$. By Proposition \ref{basic1}, it is enough to show that $\frac{i}{n} \notin I$ when $i \in \text{G}(S)$ and $n \in \mathbb N$, and $\frac{j}{n} \in I$ when $j \in \text{AP}(a,d;k)$ for some $n=n(j) \in \mathbb N$ that depends on $j$. By eqn.~\eqref{Bez1}, there exists a proper B\'{e}zout sequence $\left\{\frac{a_1}{b_1},\ldots,\frac{a_k}{b_k}\right\}$ such that $\alpha \le \frac{a_1}{b_1}<\frac{a_k}{b_k} \le \beta$ and $a_1,\ldots,a_k$ is a permutation of $\text{AP}(a,d;k)$. Then there exists $r \in \{1,\ldots,k\}$ such that
\[ a_1 > \cdots > a_r < a_{r+1} < \cdots < a_k \]
by eqn.~\eqref{convex}. Therefore, $a_r=a$ and $a+d=\min\{a_{r-1},a_{r+1}\}$, unless $r=1$ or $k$.

If $r=1$, then $a_1,\ldots,a_k$ is the increasing sequence $a,a+d,\ldots,a+(k-1)d$. Thus, $a_i=a+(i-1)d$, $1 \le i \le k$. The sequence of $b_i$'s satisfy $a_{i+1}b_i-a_ib_{i+1}=1$, so that $b_i \equiv a_{i+1}^{-1} \pmod{a_i} \equiv d^{-1} \pmod{a_i}$ for $1 \le i \le k-1$. From $a_kb_{k-1}-a_{k-1}b_k=1$ we have $b_k \equiv -(a_{k-1})^{-1} \pmod{a_k} \equiv d^{-1} \pmod{a_k}$. Thus, $b_i \equiv d^{-1} \pmod{a_i}$ for each $i \in \{1,\ldots,k\}$. Since $b_i,d^{-1}\pmod{a_i}$ belong to $\{1,\ldots,a_i-1\}$, we have $b_i=d^{-1}\pmod{a_i}$, for each $i \in \{1,\ldots,k\}$. So $\alpha,\beta$ must satisfy
\begin{equation} \label{r=1}
\alpha \le \frac{a_1}{b_1} = \frac{a}{d_a} < \frac{a_k}{b_k} = \frac{b}{d_b} \le \beta.
\end{equation}
In order that $\frac{i}{n} \notin I$ whenever $i \in \text{G}(S)$ and $n \in \mathbb N$, we must have either $\frac{i}{n}<\alpha$ or $\frac{i}{n}>\beta$. From eqn.~\eqref{r=1}, this amounts to $\frac{i}{n}<\frac{a}{d_a}$ or $\frac{i}{n}>\frac{b}{d_b}$, or to $n>\frac{id_a}{a}$ or $n<\frac{id_b}{b}$. Therefore, $\alpha$ is greater than the maximum of $i/\left\lceil\frac{id_a}{a}\right\rceil$ as $i$ runs through values in $\text{G}(S)$ and $\beta$ is less than the minimum of $i/\left\lfloor\frac{id_b}{b}\right\rfloor$ as $i$ runs through values in $\text{G}(S)$. Therefore, by Propositions \ref{max} and \ref{min}
\begin{equation} \label{r=1bounds}
\alpha > \frac{\lambda a+d(a-1)}{\lambda d_a+q_a(a-1)+1} \;\;\text{and}\;\; \beta < \frac{\lambda b+d}{\lambda d_b+q_b}.
\end{equation}
This gives the first of the two ranges for $\alpha$ and $\beta$ in the theorem using Proposition \ref{link}, part (i).
\vskip 5pt

If $r=k$, then $a_1,\ldots,a_k$ is the decreasing sequence $a+(k-1)d,a+(k-2)d,\ldots,a$.  Thus, $a_i=a+(k-i)d$, $1 \le i \le k$. Arguing as above, we see that $b_i=a_i-\left(d^{-1} \pmod{a_i}\right)$ for each $i \in \{1,\ldots,k\}$. So $\alpha,\beta$ must satisfy
\begin{equation} \label{r=k}
\alpha \le \frac{a_1}{b_1} = \frac{b}{b-d_b} < \frac{a_k}{b_k} = \frac{a}{a-d_a} \le \beta.
\end{equation}
In order that $\frac{i}{n} \notin I$ whenever $i \in \text{G}(S)$ and $n \in \mathbb N$, we must have either $\frac{i}{n}<\alpha$ or $\frac{i}{n}>\beta$. From eqn.~\eqref{r=k}, this amounts to $\frac{i}{n}<\frac{b}{b-d_b}$ or $\frac{i}{n}>\frac{a}{a-d_a}$, or to $n>\frac{i(b-d_b)}{b}$ or $n<\frac{i(a-d_a)}{a}$. Therefore, $\alpha$ is greater than the maximum of $i/\left\lceil\frac{i(b-d_b)}{b}\right\rceil$ as $i$ runs through values in $\text{G}(S)$ and $\beta$ is less than the minimum of $i/\left\lfloor\frac{i(a-d_a)}{a}\right\rfloor$ as $i$ runs through values in $\text{G}(S)$. Since
\[ \frac{m}{\left\lceil\frac{(b-d_b)m}{b}\right\rceil} < \frac{n}{\left\lceil\frac{(b-d_b)n}{b}\right\rceil}  \Longleftrightarrow
    \frac{n}{\left\lfloor\frac{d_b n}{b}\right\rfloor} < \frac{m}{\left\lfloor\frac{d_b m}{b}\right\rfloor}
\]
and
\[ \frac{m}{\left\lfloor\frac{(a-d_a)m}{a}\right\rfloor} < \frac{n}{\left\lfloor\frac{(a-d_a)n}{a}\right\rfloor}  \Longleftrightarrow
    \frac{n}{\left\lceil\frac{d_a n}{a}\right\rceil} < \frac{m}{\left\lceil\frac{d_a m}{a}\right\rceil}
\]
for every $m,n \in \mathbb N$, we have
\begin{equation} \label{r=kbounds}
\alpha > \frac{\lambda b+d}{\lambda(b-d_b)+(d-q_b)} \;\;\text{and}\;\; \beta < \frac{\lambda a+d(a-1)}{\lambda(a-d_a)+(d-q_a)(a-1)-1},
\end{equation}
by Propositions \ref{max} and \ref{min}. This gives the second of the two ranges for $\alpha$ and $\beta$ in the theorem using Proposition \ref{link}, part (i).
\vskip 5pt

Henceforth, assume $r \ne 1,k$. Two cases arise: (i) $a_{r+1}=a+d$, and (ii) $a_{r-1}=a+d$. These cases arise only if $k=a$, as we show below.
\vskip 5pt

\noindent {\sc Case} (i) From $a_r \mid (a_{r-1}+a_{r+1})$ we have $a_{r-1} \equiv -a_{r+1} \equiv -d\pmod{a}$. If $a_{r-1}=a+id$, with $0 \le i \le k-1 \le a-1$, then $i \equiv -1\pmod{a}$. This is possible only if $k=a$, and then $a_{r-1}=a+(a-1)d=b$. From eqn.~\eqref{convex}, the sequence $a_1,\ldots,a_k$ must be $b=a+(a-1)d,a,a+d,\ldots,a+(a-2)d$. The corresponding sequence $b_2,\ldots,b_k$ is the same as in the case $r=1$ above. Moreover, $ab_1-\big(a+(a-1)d\big)b_2=1$, so that $b_1 \equiv a_b \pmod{b}$. Since $b_1,a_b$ belong to $\{1,\ldots,b-1\}$, we have $b_1=a_b$.
So $\alpha,\beta$ must satisfy
\begin{equation} \label{r=1,k}
\alpha \le \frac{a_1}{b_1} = \frac{b}{a_b} < \frac{a_k}{b_k} = \frac{b-d}{d^{\prime}} \le \beta.
\end{equation}
In order that $\frac{i}{n} \notin I$ whenever $i \in \text{G}(S)$ and $n \in \mathbb N$, we must have either $\frac{i}{n}<\alpha$ or $\frac{i}{n}>\beta$. From eqn.~\eqref{r=1,k}, this amounts to $\frac{i}{n}<\frac{b}{a_b}$ or $\frac{i}{n}>\frac{b-d}{d^{\prime}}$, or to $n>\frac{i a_b}{b}$ or $n<\frac{i d^{\prime}}{b-d}$. Therefore, $\alpha$ is greater than the maximum of $i/\left\lceil\frac{i a_b}{b}\right\rceil$ as $i$ runs through values in $\text{G}(S)$ and $\beta$ is less than the minimum of $i/\left\lfloor\frac{i d^{\prime}}{b-d}\right\rfloor$ as $i$ runs through values in $\text{G}(S)$. By Proposition \ref{link}, part (ii), the bounds for $\alpha$ and $\beta$ are as given by eqn.~\eqref{r=1bounds} with $\lambda=0$. This gives the first of the two ranges for $\alpha$ and $\beta$ in the additional case $k=a$ in the theorem using Proposition \ref{link}, part (i).
\vskip 5pt

\noindent {\sc Case} (ii) From $a_r \mid (a_{r-1}+a_{r+1})$ we have $a_{r+1} \equiv -a_{r-1} \equiv -d\pmod{a}$. Arguing as in Case (i), the sequence $a_1,\ldots,a_k$ must be $a+(a-2)d,a+(a-3)d,\ldots,a,a+(a-1)d=b$, and this is possible only if $k=a$. The corresponding sequence $b_1,\ldots,b_{k-1}$ is the same as in the case $r=k$ above, and $b_k=b-\left(a_b \pmod{b}\right)$. So $\alpha,\beta$ must satisfy
\begin{equation} \label{r=k}
\alpha \le \frac{a_1}{b_1} = \frac{b-d}{b-d-d^{\prime}} < \frac{a_k}{b_k} = \frac{b}{b-a_b} \le \beta.
\end{equation}
In order that $\frac{i}{n} \notin I$ whenever $i \in \text{G}(S)$ and $n \in \mathbb N$, we must have either $\frac{i}{n}<\alpha$ or $\frac{i}{n}>\beta$. From eqn.~\eqref{r=k}, this amounts to $\frac{i}{n}<\frac{b-d}{b-d-d^{\prime}}$ or $\frac{i}{n}>\frac{b}{b-a_b}$, or to $n>\frac{i(b-d-d^{\prime})}{b-d}$ or $n<\frac{i(b-a_b)}{b}$. Therefore, $\alpha$ is greater than the maximum of $i/\left\lceil\frac{i(b-d-d^{\prime})}{b-d}\right\rceil$ as $i$ runs through values in $\text{G}(S)$ and $\beta$ is less than the minimum of $i/\left\lfloor\frac{i(b-a_b)}{b}\right\rfloor$ as $i$ runs through values in $\text{G}(S)$. Since
\[ \frac{m}{\left\lceil\frac{(b-d-d^{\prime})m}{b-d}\right\rceil} < \frac{n}{\left\lceil\frac{(b-d-d^{\prime})n}{b-d}\right\rceil}  \Longleftrightarrow
     \frac{n}{\left\lfloor\frac{d^{\prime}n}{b-d}\right\rfloor} < \frac{m}{\left\lfloor\frac{d^{\prime}m}{b-d}\right\rfloor}
\]
and
\[ \frac{m}{\left\lfloor\frac{(b-a_b)m}{b}\right\rfloor} < \frac{n}{\left\lfloor\frac{(b-a_b)n}{b}\right\rfloor}  \Longleftrightarrow
     \frac{n}{\left\lceil\frac{a_b n}{b}\right\rceil} < \frac{m}{\left\lceil\frac{a_b m}{b}\right\rceil}
\]
for every $m,n \in \mathbb N$, we have the bounds for $\alpha$ and $\beta$ as given by eqn.~\eqref{r=kbounds} with $\lambda=0$ by Proposition \ref{link}, part (ii). This gives the second of the two ranges for $\alpha$ and $\beta$ in the additional case $k=a$ in the theorem using Proposition \ref{link}, part (i).
\end{Pf}
\vskip 5pt

\begin{cor} \label{spcase1}
Let $a,k$ be positive integers, with $2 \le k \le a$, $b=a+k-1$, and $\lambda=\left\lfloor\frac{a-2}{k-1}\right\rfloor$. Let $S\big(AP(a,1;k)\big) = \langle a,a+1,a+2,\ldots,a+k-1 \rangle$. Then for $\alpha,\beta \in [1,\infty) \cap \mathbb Q$,
\[ S\big([\alpha,\beta]\big) = S\big(AP(a,1;k)\big) \]
if and only if
\[ a-\frac{1}{\lambda+1} < \alpha \le a, \quad b \le \beta < b+\frac{1}{\lambda} \]
or
\[ 1+\frac{\lambda}{\lambda(b-1)+1} < \alpha \le 1+\frac{1}{b-1}, \quad 1+\frac{1}{a-1} \le \beta < 1+\frac{\lambda+1}{(\lambda+1)(a-1)-1}. \]
\noindent Moreover, for $k=a$, we also have
\[ a-1 < \alpha \le \frac{b}{2}, \quad b-1 \le \beta, \]
or
\[ 1 < \alpha \le 1+\frac{1}{b-2}, \quad 1+\frac{2}{b-2} \le \beta < 1+\frac{1}{a-2}. \]
\noindent We adopt the convention $\frac{1}{0}=\infty$.
\end{cor}
\vskip 5pt

\begin{rem} \label{spcase1_k=a}
For any positive integer $a>1$ and for $\alpha,\beta \in [1,\infty) \cap \mathbb Q$,
\[ S\big([\alpha,\beta]\big) = S\big(AP(a,1;a)\big) = \big\{0,a,\to\big\} \]
is the special case $\lambda=0$ in Corollary \ref{spcase1}. We adopt the convention $\frac{1}{0}=\infty$.
\end{rem}
\vskip 5pt

\begin{thm} \label{II}
Let $a,d,k$ be positive integers, with $\gcd(a,d)=1$, $2 \le k \le a$, $b=a+(k-1)d$, and let $\lambda=\left\lfloor\frac{a-2}{k-1}\right\rfloor$. Define $dd_a=1+q_a a$,
where $1 \le d_a < a$ and $0 \le q_a < d$. Let $S\big(AP(a,d;k)\big)=\langle a,a+d,a+2d,\ldots,a+(k-1)d \rangle$. Then for positive integers $p,q,m$ with
$q<p<m$,
\[ S\big(p,m,q\big) = S\big(AP(a,d;k)\big) \]
if and only if
\begin{eqnarray*}
\frac{md_a}{a} \le p < \frac{m\big(\lambda d_a+q_a(a-1)+1\big)}{\lambda a+d(a-1)}, \\
p + \frac{m\big((a-1)q_a+d_a\big)}{b} \le q < p - \frac{m\left(\lambda\big((a-1)q_a+d_a\big)+q_a\right)}{\lambda b+d}
\end{eqnarray*}
or
\begin{eqnarray*}
m - \frac{m\big((a-1)q_a+d_a\big)}{b} \le p < m - \frac{m\left(\lambda\big((a-1)q_a+d_a\big)+q_a\right)}{\lambda b+d}, \\
p - m + \frac{md_a}{a} \le q < p - m - \frac{m\big((\lambda d_a+q_a)(a-1)-1\big)}{\lambda a+d(a-1)}.
\end{eqnarray*}
\noindent Moreover, for $k=a$, we also have
\[ \frac{m\big((a-1)q_a+d_a+1\big)}{b} \le p < \frac{m\big(q_a(a-1)+1\big)}{d(a-1)}, \quad p - \frac{m\big((a-2)q_a+d_a\big)}{b-d} \le q < p - \frac{mq_a}{d}, \]
or
\begin{eqnarray*}
m - \frac{m\big((a-2)q_a+d_a\big)}{b-d} \le p < m - \frac{mq_a}{d}, \\
p - m + \frac{m\big((a-1)q_a+d_a+1\big)}{b} \le q < p - m + \frac{m\big(q_a(a-1)+1\big)}{d(a-1)}.
\end{eqnarray*}
\end{thm}

\begin{Pf}
This follows directly from Proposition \ref{basic2} and Theorem \ref{I}. If $S\big(p,m,q\big)=S\big([\alpha,\beta]\big)$, then $\frac{m}{p}=\alpha$ and                $\frac{m}{p-q}=\beta$, so that $p=\frac{m}{\alpha}$ and $p-q=\frac{m}{\beta}$. The latter implies $q=p-\frac{m}{\beta}$. The bounds on $\alpha,\beta$ in Theorem \ref{I} translate to bounds on $p,q$ in terms of $m$ in Theorem \ref{II}.
\end{Pf}
\vskip 5pt

\begin{cor} \label{spcase2}
Let $a,k$ be positive integers, with $2 \le k \le a$, $b=a+k-1$, and $\lambda=\left\lfloor\frac{a-2}{k-1}\right\rfloor$. Let $S\big(AP(a,1;k)\big) = \langle a,a+1,a+2,\ldots,a+k-1 \rangle$.  Then for positive integers $p,q,m$ with $q<p<m$,
\[ S\big(p,m,q\big) = S\big(AP(a,1;k)\big) \]
if and only if
\[ \frac{m}{a} \le p < \frac{m\big(\lambda+1\big)}{(\lambda+1)a-1}, \quad p + \frac{m}{b} \le q < p - \frac{m\lambda}{\lambda b+1} \]
or
\[ m - \frac{m}{b} \le p < m - \frac{m\lambda}{\lambda b+1}, \quad p - m + \frac{m}{a} \le q < p - m - \frac{m\big(\lambda(a-1)-1\big)}{(\lambda+1)a-1)}. \]
\noindent Moreover, for $k=a$, we also have
\[ \frac{2m}{b} \le p < \frac{m}{a-1}, \quad p - \frac{m}{b-1} \le q < p, \]
or
\[ m - \frac{m}{b-1} \le p < m, \quad p - m + \frac{2m}{b} \le q < p - m + \frac{m}{a-1}. \]
\end{cor}
\vskip 5pt

\begin{rem} \label{spcase2_k=a}
For any positive integer $a>1$ and positive integers $p,q,m$ with $q<p<m$,
\[ S\big(p,m,q\big) = S\big(AP(a,1;a)\big) = \big\{0,a,\to\big\} \]
is the special case $\lambda=0$ in Corollary \ref{spcase2}.
\end{rem}
\vskip 20pt