Discrete Hardy-type Inequalities
This paper studies the Hardy-type inequalities on the discrete intervals. The first result is the variational formulas of the optimal constants. Using these formulas, one may obtain an approximating procedure and the known basic estimates of the optimal constants. The second result, which is the main innovation of this paper, is about the factor of basic upper estimates. An improved factor is presented, which is smaller than the known one and is best possible. Some comparison results are included for comparing the optimal constants on different intervals.
\section{Introduction}

For given two constants $p$ and $q$ with $1< p\leqslant q< \infty$, two positive sequences $\mathbf{u}$ and $\mathbf{v}$ on a discrete interval $[1, N]:= \{1, 2, \dots, N \}$ with $N\leqslant +\infty$, here is the discrete Hardy-type inequality: 
\begin{equation}\label{Hardy}
\left[\sum_{n= 1}^N u_n \left(\sum_{i=1}^n x_i \right)^q \right]^{1/q} \leqslant A \left(\sum_{n=1}^N v_n x_n^p \right)^{1/p},
\end{equation}
where $\mathbf{x}$ is an arbitrary non-negative sequence on $[1, N]$. For saving notations, the constant $A$ is assumed to be optimal. 
 
The purpose of this paper is two-fold. First, we give some variational formulas of the optimal constants. The primary applications of the variational formulas are the approximating procedure and the basic estimates. It is necessary to review the advance of the basic estimates in recent research, cf. \cite{Bliss, Chen4, Chen3, Maz'ya, Opic}. In continuous case, the following result, due to B. Opic \rf{Opic}{Theorem 1.14} and V. G. Maz'ya \rf{Maz'ya}{Theorem 1, pp. 42-43}, is well known
\bg{equation}\label{basic}
B \leqslant A \leqslant \tilde{k}_{q, p} B, 
\end{equation}
where $B$ is a quantity described by $N$, $p$, $q$, $\mathbf{u}$ and $\mathbf{v}$, the factor $\tilde{k}_{q, p}$ is a constant defined by $p$ and $q$:
\bg{equation} \label{tilde k_qp}
\tilde{k}_{q,p} = \left(1+ \frac{q}{p^*} \right)^{1/q} \left(1+ \frac{p^*}{q} \right)^{1/p^*}, 
\end{equation}
here $p^*$ is the conjugate number of $p$, i.e. $1/p + 1/{p^*}= 1$. In particular, $\tilde{k}_{p,p}= p^{1/p}(p^*)^{1/{p^*}}$. Afterwards, Chen \rf{Chen3}{Theorem 2.1} gets the same conclusions through the variational formulas of the optimal constants. Furthermore, there is an approximating procedure \rf{Chen3}{Theorem 2.2} based on the variational formulas, which can improve the estimates of the optimal constants step by step. In discrete context, when $p=q$, Chen, Wang and Zhang \cite{Chen2} arrive the corresponding variational formulas and basic estimates which are similar to (\ref{basic}), of course, $B$ must be adjusted appropriately in discrete case. When $p\ne q$, Mao \rf{Mao}{Proposition A.1} gets the similar result, but the factor of basic upper estimates is $p^{1/q} (p^*)^{1/p^*}$, which is a little coarser than $\tilde{k}_{q,p}$. Our first destination is to show the corresponding variational formulas in discrete context with the condition of $p\neq q$. Later on, as applications of these formulas, we obtain the basic estimates and the approximating procedure. Overall, these results can be regarded as an extension of the studies in continuous context \cite{Chen3}. 

Second, we study the upper bounds of the basic estimates of the optimal constants in discrete case. Our result is the factor $\tilde{k}_{q, p}$ in (\ref{basic}) can be improved to $k_{q, p}$: 
\bg{equation} \label{k_qp}
k_{q,p} = \left(\frac{r}{B(\frac{1}{r}, \frac{q- 1}{r})} \right)^{1/p- 1/q}, 
\end{equation}
where $B(a, b)= \int_0^1 x^{a-1} (1- x)^{b- 1} \text{\rm d} x$ is the Beta function and $r= q/p -1$. Moreover, our result shows that the factor is best possible and is consistent with the result of continuous case. In continuous case, the improvement has been worked out, cf. \rf{Bennett3}{Theorem 8}, \rf{Manakov}{Theorem 2}, and \rf{Kufner}{pp. 45-47}. The key is the result of Bliss \cite{Bliss}, which gives an integral inequality that the optimal constants can be attained. However, the analogue of the conclusion in the discrete context is nontrivial, as mentioned in \rf{Bennett3}{page 170, two lines above (61)}, \textquotedblleft I have been unable to prove the discrete analogue of Theorem 8\textquotedblright (here the last result is the continuous case). We are lucky to be able to prove this conclusion which constitutes the second part of this paper. When $p=q$, it is well known that the factor $\tilde{k}_{q,p}$ is sharp (see for instance \rf{Hardy_book}{Theorem 326 and 327}). Note that if we allow $q\rightarrow p$, by the identity 
$$
\lim_{r\rightarrow 0^{+}} B^r \left(\frac{a}{r}, \frac{b}{r} \right) = \frac{a^a b^b}{(a+ b)^{a+ b}}, \qquad (a,b >0),
$$
we have
\bg{equation}\label{k_pp}
k_{p,p}= \tilde{k}_{p,p}= p^{1/p}(p^*)^{1/{p^*}}. 
\end{equation}
It means our improved factor is consistent with the original one when $p= q$. Thus, our main results are devoted to the case of $p< q$.

It is a long time for the research about Hardy-type inequalities, which are one of the major themes of harmonic analysis and represent useful tools e.g. in the theory and practice of differential equations, in the theory of approximation etc. From probabilistic consideration, these inequalities are important tools to study the convergence rate of the corresponding processes. These are the origin and motive of this study. 

This paper is organized as follows. The rest of this section, we give some notations and definitions, then illustrate the main results. In Section 2 and Section 3, we prove the conclusions on discrete half line (i.e. $N= \infty$). The case of finite interval (i.e. $N< \infty$) will be handled unitedly in the final section, which gives some comparison results for the optimal constants and their basic estimates on different intervals. 

For the simplicity of illustration, we need some notations. Let $\hat{v}_i = v_i^{1- p^*}, 1\leqslant i\leqslant N$.  For any sequence $\mathbf{x}$ on $[1, N]$, define an operator $H$: 
\bg{equation}\label{H}
H\mathbf{x}(n)= 
\begin{cases}
0, & n= 0, \\
\sum_{i=1}^n x_i, & n=1, 2, \cdots, N,
\end{cases}
\end{equation}
which means the partial summation of $\mathbf{x}$. The following notations are used frequently: 
$$
\alpha}  \def\bz{\beta \wedge \bz = \min \{ \alpha}  \def\bz{\beta, \bz \}, \qquad \alpha}  \def\bz{\beta \vee \bz = \max \{ \alpha}  \def\bz{\beta, \bz \}.
$$

Set
$$
{\scr A}}   \def\bq{{\scr B}[1, N]= \{\mathbf{x}: \text{$x_1>0$, and $x_i\geqslant 0$ for $i=2, \dots, N$} \}.
$$
In this article, we use the convention $1/0 = \infty$. Then the optimal constant $A$ can be denoted by the following variational formula: 
\begin{equation}\label{A}
A = \sup_{\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, N]} \frac{\left[\sum_{n= 1}^N u_n \left(\sum_{i=1}^n x_i \right)^q \right]^{1/q}}{\left(\sum_{n= 1}^N v_n x_n^p \right)^{1/p}} = \sup_{\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, N]} \frac{\| H\mathbf{x} \|_{l^q(u)}}{\| \mathbf{x} \|_{l^p(v)}}, 
\end{equation}
where $\| \mathbf{x} \|_{l^q(u)}= \left[ \sum_{n=1}^{N} u_n x_n^q \right]^{1/q}$ and similarly to $\| \mathbf{x} \|_{l^p(v)}$. 

For upper estimates, define the single summation operator $I^*$ and the double summation operator $I\!I^*$ as:
\begin{align}
I_n^*(\mathbf{x}) &= \frac{\hat{v}_n}{x_n} \left(\sum_{i=n}^N u_i (H\mathbf{x}(i))^{q/p^*} \right)^{p^*/q}, \label{I*} \\
I\!I_n^*(\mathbf{x}) &= \frac{1}{H\mathbf{x}(n)} \sum_{i=1}^n \hat{v}_i \left(\sum_{j=i}^N u_j (H\mathbf{x}(j))^{q/p^*} \right)^{p^*/q}, \label{II*}
\end{align}
with domain ${\scr A}}   \def\bq{{\scr B}[1, N]$. 

For lower estimates, there are some differences:
\begin{align}
I_n (\mathbf{x})&= \frac{\hat{v}_n}{x_n} \left(\sum_{i=n}^N u_i (H\mathbf{x}(i))^{q- 1} \right)^{p^*- 1}, \label{I} \\
I\!I_n (\mathbf{x})&= \frac{1}{H \mathbf{x} (n)} \sum_{i=1}^n \hat{v}_i \left(\sum_{j=i}^N u_j (H\mathbf{x}(j))^{q- 1} \right)^{p^*- 1}. \label{II} \\
\end{align}
It is easy to see that $I\!I^*= I\!I$ and $I^*= I$ when $p= q$. To avoid the non-summability problem, the domain of $I$ and $I\!I$ have to be modified to:
$$
{\scr A}}   \def\bq{{\scr B}_0 [1, N]= \left\{\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, N]: \sum_{i=1}^N v_i x_i^p < \infty \right\}.
$$

With these notations, we can give the main conclusions of the variational formulas of the optimal constants.

\nnd\begin{thm1}\label{Var}
{\cms
The optimal constant $A$ in the Hardy-type inequality {\rm (\ref{Hardy})} satisfies

(i) upper estimates:
\begin{equation}
A\leqslant \inf_{\mathbf{x}\in {\scr A}}   \def\bq{{\scr B}[1, N]} \left(\sup_{n\in [1, N]} I\!I_n^*(\mathbf{x}) \right)^{1/p^*} = \inf_{\mathbf{x}\in {\scr A}}   \def\bq{{\scr B}[1, N]} \left(\sup_{n\in [1, N]} I_n^*(\mathbf{x}) \right) ^{1/p^*}.
\end{equation}

(ii) lower estimates:
\begin{equation}
\aligned
A&\geqslant \sup_{\mathbf{x}\in {\scr A}}   \def\bq{{\scr B}_0 [1, N]} \|\mathbf{x}\|_{l^p(v)}^{p/q -1} \left(\inf_{n\in [1,N]} I\!I_n (\mathbf{x})\right)^{{(p-1)}/q} \\
&\geqslant \sup_{\mathbf{x}\in {\scr A}}   \def\bq{{\scr B}_0 [1, N]} \|\mathbf{x}\|_{l^p(v)}^{p/q -1} \left(\inf_{n\in [1,N]} I_n (\mathbf{x})\right)^{{(p-1)}/q}. \\
\endaligned
\end{equation}
}
\end{thm1}

Using Theorem \ref{Var} on a appropriate test function, we can obtain the basic estimates (\ref{basic}). To be specific, that is:

\nnd\begin{crl1}\label{basic estimates}
{\cms
The inequality {\rm (\ref{Hardy})} holds for every $\mathbf{x}\in {\scr A}}   \def\bq{{\scr B}[1, N]$ if and only if
$B< \infty$, where
\bg{equation}\label{B}
B= \sup_{n\in [1, N]} \left( \sum_{i=1}^n \hat{v}_i \right)^{1/p^*} \left(\sum_{j=n}^N u_j \right)^{1/q}.
\end{equation}
Moreover, we have
$$
B\leqslant A\leqslant \tilde{k}_{q,p}B,
$$
where $\tilde{k}_{q,p}$ is defined as (\ref{tilde k_qp}), which is independent of $\mathbf{u}$, $\mathbf{v}$ and $N$. 
}
\end{crl1}

Roughly speaking, the conclusion of Corollary \ref{basic estimates} is from the first iteration through an appropriate test function. Moreover, we can improve the estimates step by step through multiple iterations on this test function. The following corollary is based on this idea, which is of great significance to numerical computation.

\nnd\begin{crl1}\label{approximating}
{\cms
(i) Define
$$
\aligned
x^{(1)}_n &= (H \mathbf{\hat{v}} (n))^\alpha}  \def\bz{\beta - (H \mathbf{\hat{v}} (n- 1))^\alpha}  \def\bz{\beta, \\
x^{(m+ 1)}_n &= \hat{v}_n \left(\sum_{i= n}^N u_i (H \mathbf{x}^{(m)} (i))^{q/{p^*}} \right)^{{p^*}/q}, 
\endaligned
$$
where $\alpha}  \def\bz{\beta= q/(p^*+ q)$. If $\dz_1:=  \sup_{n\in [1, N]} \left( I\!I_n (\mathbf{x}^{(1)}) \right)^{1/p^*}< \infty$, define a sequence as: 
\bg{equation}
\dz_m= \sup_{n\in [1, N]} \left( I\!I^*_n(\mathbf{x}^{(m)})\right)^{1/p^*}, \quad m=2, 3,  \dots. 
\end{equation}
Otherwise, define $\dz_m \equiv \infty$. Then $\dz_m$ is a non-increasing sequence (denote $\dz_\infty$ be the limit of $\dz_m$) and we have
$$
A\leqslant \dz_\infty \leqslant \cdots \leqslant \dz_1 \leqslant \tilde{k}_{q,p} B. 
$$

(ii) Fix $k\in [1, N]$, define
$$
\aligned
y^{(k, 1)}_n &=
\begin{cases}
\hat{v}_n, & 1\leqslant n\leqslant k, \\
0, & n> k, 
\end{cases} \\
y^{(k, m+1)}_n &= \hat{v}_n \left(\sum_{i=n}^N \left(H \mathbf{y}^{(k, m)} (i) \right)^{q- 1} \right)^{p^*- 1},
\endaligned
$$
and define a sequence as
$$
\aligned
\widetilde{\dz}_m&= \sup_{k\in [1, N]} \big\| \mathbf{y}^{(k, m)} \big\|_{l^p (v)}^{p/q -1} \left(\inf_{n\in [1, N]} I\!I_n \left( \mathbf{y}^{(k, m)} \right) \right)^{{(p- 1)}/q}, \\
\overline{\dz}_m&= \sup_{k\in [1, N]} \frac{\left[\sum_{n= 1}^N u_n \left( H \mathbf{y}^{(k, m)} (n) \right)^q \right]^{1/q}}{\left[\sum_{n=1}^N v_n \left(y_n^{(k, m)} \right)^p \right]^{1/p}}. \\
\endaligned
$$
Then we have $A\geqslant \widetilde{\dz}_m \vee \overline{\dz}_m$ for all $m\geqslant 1$.
}
\end{crl1}

Another main result of this paper is about the factor in (\ref{basic}). Just like the continuous case, the factor of the basic upper estimates can be improved. Furthermore, we can prove the improved factor is best possible. This result is described in detail below.
\nnd\begin{thm1}\label{kqpB}
{\cms
The basic upper estimates can be improved to: 
\begin{equation}
A\leqslant k_{q, p}B, 
\end{equation}
where $k_{q,p}$ is defined as (\ref{k_qp}). In particular, when $N= \infty$ and $\sum_{i=1}^\infty \hat{v}_i = \infty$, the factor $k_{q,p}$ is sharp. 
}
\end{thm1}

\section{Proof of Theorem \ref{Var} }

In this section, we always assume $N= \infty$, and the case of the finite interval will be discussed in Section \ref{Interval}. 

The first proposition is about the property of the sequences which reach the equality case of (\ref{Hardy}). This result is useless in the proof of Theorem \ref{Var}, but is necessary to study the property of the optimal constants. 
 
\nnd\begin{prp1}\label{decreasing}
{\cms
If the optimal constant $A$, appearing in (\ref{Hardy}), is attained by some non-negative sequence $\mathbf{x}$. Define $w_n = \hat{v}_n^{-1} x_n$, then the sequence $\mathbf{w}$ is decreasing. 
}
\end{prp1}

\medskip \noindent {\bf Proof}. 
With the definition of $\mathbf{w}$, we can rewrite the Hardy- type inequalities (\ref{Hardy}) as:
\bg{equation}\label{Hardy*}
\left[ \sum_{n=1}^\infty u_n \left(\sum_{i=1}^n \hat{v}_i w_i \right)^q \right]^{1/q}\leqslant A \left(\sum_{n=1}^\infty \hat{v}_n w_n^p \right)^{1/p},
\end{equation}
and $A$ is attained at $\mathbf{w}$.

The idea in the remainder of this proof is from Bennett \rf{Bennett1}{Section 3}. Using reduction to absurdity, assume there exist integers $i$ and $j$ with $1\leqslant i< j< \infty$, but $w_i< w_j$. We can construct a new sequence $\mathbf{w'}$ from $\mathbf{w}$ as
\begin{equation}\label{decreasing1}
w'_i= w'_j= w_0,
\end{equation}
where $w_0$ satisfies 
\begin{equation}\label{decreasing2}
(\hat{v}_i+ \hat{v}_j) w_0^p= \hat{v}_i w_i^p+ \hat{v}_j w_j^p.
\end{equation}

On the one hand, from (\ref{decreasing2}), the right side of (\ref{Hardy*}) is unchanged when $\mathbf{w}$ is replaced by $\mathbf{w'}$. On the other hand, since $p> 1$ and (\ref{decreasing2}), we have
\bg{equation}\label{decreasing3}
w_i< w_0< w_j.
\end{equation}
\bg{equation}\label{decreasing4}
\hat{v}_i w_i + \hat{v}_j w_j < (\hat{v}_i + \hat{v}_j) w_0.
\end{equation}
Combine (\ref{decreasing3}) and (\ref{decreasing4}), we obtain 
$$
\sum_{k=1}^n \hat{v}_k w_k < \sum_{k=1}^n \hat{v}_k w'_k, \qquad \forall n \geqslant 1. 
$$
It means the left side of (\ref{Hardy*}), increases strictly when $\mathbf{w}$ is replaced by $\mathbf{w'}$. Hence $A$ is not attained at $\mathbf{w}$, which is a contradiction.
\quad $\square$ \medskip

\noindent {\bf Proof of Theorem \ref{Var}.} The briefing of the proof of Theorem \ref{Var} is given as follows.

(a) First, we need to verify the relation of the single summation operator $I^*$ and the double summation operator $I\!I^*$:
$$
\inf_{\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, \infty)} \left[\sup_{n\in [1, \infty)} I\!I_n^*(\mathbf{x}) \right]^{1/p^*} = \inf_{\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, \infty)} \left[\sup_{n\in [1, \infty)} I_n^*(\mathbf{x}) \right]^{1/p^*}.
$$

For any $\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, \infty)$, as an application of the proportional property, we get
$$
\aligned
\sup_{n\in [1, \infty)} I\!I_n^*(\mathbf{x}) &= \sup_{n\in [1, \infty)} \frac{1}{H \mathbf{x} (n)} \left[ \sum_{i= 1}^{n} \hat{v}_i \left(\sum_{j=i}^{\infty} u_j (H \mathbf{x} (j))^{q/{p^*}} \right)^{{p^*}/q}\right] \\
&\leqslant \sup_{n\in [1, \infty)} \frac{1}{x_n} \left[ \hat{v}_n \left(\sum_{i=n}^{\infty} u_i (H \mathbf{x} (i))^{q/{p^*}} \right)^{{p^*}/q}\right] \\
&= \sup_{n\in [1, \infty)} I_n^*(\mathbf{x}).
\endaligned
$$
Hence, we have
$$
\inf_{\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, \infty)} \left[\sup_{n\in [1, \infty)} I\!I_n^*(\mathbf{x}) \right]^{1/p^*} \leqslant\inf_{\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, \infty)} \left[\sup_{n\in [1, \infty)} I_n^*(\mathbf{x}) \right]^{1/p^*}.
$$

On the other hand, for any $\mathbf{x}\in {\scr A}}   \def\bq{{\scr B}[1, \infty)$, define
$$
y_n= \hat{v}_n \left(\sum_{i= n}^{\infty} u_i (H \mathbf{x} (i))^{q/{p^*}} \right)^{{p^*}/q}.
$$
Obviously, $y_n> 0$ on $[1, \infty)$, then $\mathbf{y} \in {\scr A}}   \def\bq{{\scr B}[1, \infty)$. Again, using the proportional property, we have
$$
\aligned
\sup_{n\in [1, \infty)} I_n^*(\mathbf{y}) &= \sup_{n\in [1, \infty)} \left[ \frac{\sum_{i=n}^\infty u_i (H \mathbf{y} (i))^{q/{p^*}}}{\sum_{i=n}^\infty u_i (H \mathbf{x} (i))^{q/{p^*}}} \right]^{{p^*}/q} \\
&\leqslant \sup_{n\in [1, \infty)} \frac{1}{H \mathbf{x} (n)} \sum_{i= 1}^n \hat{v}_i \left(\sum_{j= i}^\infty u_j (H \mathbf{x} (j))^{q/{p^*}} \right)^{{p^*}/q} \\
&= \sup_{n\in [1, \infty)} I\!I_n^*(\mathbf{x}). 
\endaligned
$$
Since $\mathbf{x}$ is arbitrary, we obtain the conclusion we need.

(b) The next step is to show the upper estimates of the optimal constants. Assume $A$ is attained at a non-negative sequence $\mathbf{a}$. For each positive sequence $\mathbf{h}$, as an application of the H\"older inequality and the H\"older-Minkowski inequality, we have
$$
\aligned
\sum_{n= 1}^{\infty} u_n (H \mathbf{a} (n))^q &= \sum_{n= 1}^{\infty} u_n \left(\sum_{i= 1}^n a_i v_i^{1/p} h_i^{-1} v_i^{-1/p} h_i \right)^q \\
&\leqslant \sum_{n= 1}^{\infty} u_n \left(\sum_{i= 1}^n a_i^p v_i h_i^{-p} \right)^{q/p} \left(\sum_{k= 1}^n v_k^{-{p^*}/{p}} h_k^{p^*} \right)^{q/{p^*}} \\
&\leqslant \left\{ \sum_{n= 1}^{\infty} a_n^p v_n h_n^{-p} \left[\sum_{i= n}^{\infty} u_i \left(\sum_{k= 1}^i \hat{v}_k h_k^{p^*} \right)^{q/{p^*}} \right]^{p/q} \right\}^{q/p}.
\endaligned
$$
At the last step, we use the H\"older-Minkowski inequality, which needs the condition $p< q$. In particular, when $p= q$, it is Fubini theorem. Now, making a power $1/q$, we get
\begin{align}\label{==1}
\left[\sum_{n= 1}^{\infty} u_n (H \mathbf{a} (n))^q \right]^{1/q} &\leqslant \left\{ \sum_{n= 1}^{\infty} a_n^p v_n h_n^{-p} \left[\sum_{i= n}^{\infty} u_i \left(\sum_{k= 1}^i \hat{v}_k h_k^{p^*} \right)^{q/{p^*}} \right]^{p/q} \right\}^{1/p} \notag \\
&\leqslant \sup_{n\in [1, \infty)} \left[ \frac{1}{h_n^q} \sum_{i=n}^{\infty} u_i \left(\sum_{k= 1}^i \hat{v}_k h_k^{p^*} \right)^{q/{p^*}} \right]^{1/q} \left(\sum_{j= 1}^{\infty} v_j a_j^p \right)^{1/p}.
\end{align}
For any $\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, \infty)$, let 
$$
h_n= \left( \sum_{i=n}^{\infty} u_i (H \mathbf{x} (i))^{q/{p^*}} \right)^{1/q}, 
$$
by the proportional property, we have
$$
\aligned
\sup_{n\in [1, \infty)} & \left[\frac{1}{h_n^q} \sum_{i=n}^{\infty} u_i \left(\sum_{j= 1}^i \hat{v}_j h_j^{p^*} \right)^{q/{p^*}} \right]^{1/q} \\
&\leqslant \sup_{n\in [1, \infty)} \left[\frac{1}{H \mathbf{x} (n)} \cdot \sum_{i= 1}^n \hat{v}_i \left(\sum_{j= i}^{\infty} u_j (H \mathbf{x} (j))^{q/p^*} \right)^{p^*/q} \right]^{1/p^*} \\
&= \sup_{n\in [1, \infty)} {I\!I_n^*(\mathbf{x})}^{1/p^*}.
\endaligned
$$
Inserting this formula into (\ref{==1}), we obtain
$$
A = \frac{\left( \sum_{n= 1}^\infty u_n (H \mathbf{a} (n))^q \right)^{1/q}}{\left( \sum_{n= 1}^\infty v_n a_n^p \right)^{1/p}}\leqslant \sup_{n\in [1, \infty)} {I\!I_n^*(\mathbf{x})}^{1/p^*}.
$$
Since $\mathbf{x}$ is arbitrary, it follows that
$$
A \leqslant \inf_{\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, \infty)} \sup_{n \in [1, \infty)} {I\!I_n^*(\mathbf{x})}^{1/p^*}. 
$$

(c) For the lower estimates, again, we consider the relation of $I$ and $I\!I$ first. For any $\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}_0 [1, \infty)$, we need to show:
$$
\inf_{n\in [1, \infty)} I_n (\mathbf{x}) \leqslant \inf_{n\in [1, \infty)} I\!I_n (\mathbf{x}). 
$$
In fact, with the help of proportional property, an argument similar to the one used in (a) can show this result. Next, we should show the variational formulas of $A$.

Since the summability of sequence $\mathbf{x}$, without loss of generality, we may assume $\sum_{i= 1}^\infty v_i x_i^p =1$. Hence our next step is to proof
$$
\sup_{\mathbf{x} \in \tilde{{\scr A}}   \def\bq{{\scr B}}[1, \infty) } \inf_{n\in [1, \infty)} I\!I_n(\mathbf{x})^{(p- 1)/q} \leqslant A,
$$
where $\tilde{{\scr A}}   \def\bq{{\scr B}}[1, \infty)= \{ \mathbf{x} \in {\scr A}}   \def\bq{{\scr B}_0 [1, \infty): \sum_{i=1}^\infty v_i x_i^p = 1 \}$.

We would begin with the classical variational formulas (\ref{A}) of the optimal constants. For any $\mathbf{x} \in \tilde{{\scr A}}   \def\bq{{\scr B}}[1, \infty)$, define
$$
y_n= \hat{v}_n \left(\sum_{i= n}^\infty u_i (H \mathbf{x} (i))^{q-1} \right)^{p^*- 1},
$$
then we have
\bg{equation}\label{==2}
A \geqslant \frac{\| H \mathbf{y} \|_{l^q (u)}}{\| \mathbf{y} \|_{l^p(v)}} = \frac{\left[\sum_{n= 1}^\infty u_n (H \mathbf{y} (n))^{q} \right]^{1/q}}{\left[\sum_{n= 1}^\infty \hat{v}_n \left(\sum_{i=n}^\infty u_i (H \mathbf{x} (i))^{q-1} \right)^{p^*} \right]^{1/p}}.
\end{equation}
Consider the denominator of (\ref{==2}), according to Fubini theorem and the definition of $\mathbf{y}$, we obtain
\begin{align}\label{==3}
\sum_{i=1}^\infty & y_i \left(\sum_{j= i}^\infty u_j (H \mathbf{x} (j))^{q-1} \right) = \sum_{j=1}^\infty u_j (H \mathbf{y} (j)) (H \mathbf{x} (j))^{q-1}\notag \\
&\leqslant \left[\sum_{j= 1}^\infty u_j (H \mathbf{y} (j))^{q/p} (H \mathbf{x} (j))^{q/{p^*}} \right]^{p/q} \left[\sum_{j= 1}^\infty u_j (H \mathbf{x} (j))^q \right]^{(q-p)/q}.
\end{align}
The last step is based on the H\"older inequality, which needs the condition $p< q$. Moreover, since $\sum_{i= 1}^\infty v_i x_i^p = 1$ we have
\bg{equation}\label{==4}
\left[\sum_{j= 1}^\infty u_j (H \mathbf{x} (j))^q \right]^{(q-p)/q} \leqslant A^{q- p}.
\end{equation}
Combining (\ref{==2}), (\ref{==3}), (\ref{==4}) and using the proportional property, we obtain
$$
A \geqslant \left[\frac{\sum_{i=1}^\infty u_i (H \mathbf{y} (i))^q}{\sum_{i=1}^\infty u_i (H \mathbf{y} (i))^{q/p} (H \mathbf{x} (i))^{q/{p^*}}} \right]^{p/{q^2}} \geqslant \inf_{n \in [1, \infty)} \left[ \frac{H \mathbf{y} (n)}{H \mathbf{x} (n)} \right]^{(p-1)/q}.
$$
By the definition of $\mathbf{y}$, we get
$$
\inf_{n\in [1, \infty)} I\!I_n(\mathbf{x})^{(p-1)/q} \leqslant A.
$$
Since $\mathbf{x}$ is arbitrary, we obtain the variational formulas of $A$. The proof is completed in the case $N= \infty$.
\quad $\square$ \medskip

\noindent {\bf Proof of Corollary \ref{basic estimates}.} With the help of Theorem \ref{Var}, we can obtain the basic estimates if we choose an appropriate test function.

We consider the upper estimates first. Before the proof, we need some preparations. Given an increasing positive sequence $\mathbf{\Phi}    \def\qqz{\Psi}$ on $[1, N]$, for any $n\in [1, N- 1]$ and $0< \alpha}  \def\bz{\beta< 1$, we assert
\begin{align}\label{==5}
\sum_{i=n+ 1}^{N} \left[\left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_n}\right)^\alpha}  \def\bz{\beta- \left(\frac{\Phi}    \def\qqz{\Psi_{i- 1}}{\Phi}    \def\qqz{\Psi_n}\right)^\alpha}  \def\bz{\beta \right] \left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_n}\right)^{-1} \leqslant \frac{\alpha}  \def\bz{\beta}{1- \alpha}  \def\bz{\beta}.
\end{align}
In fact, it can be proved by induction. Assume $n=N- 1$. Let $y= \Phi}    \def\qqz{\Psi_N/ \Phi}    \def\qqz{\Psi_{N- 1}$, then $y\geqslant 1$ (since $\mathbf{\Phi}    \def\qqz{\Psi}$ is increasing). Through simple calculations, we know the function
$$
f(x)= (x^\alpha}  \def\bz{\beta- 1) x^{-1}, \qquad x\geqslant 1,
$$
reaches the maximum when $x= \left(\frac{1}{1- \alpha}  \def\bz{\beta}\right)^{1/\alpha}  \def\bz{\beta}$. Hence
$$
\aligned
\left[\left(\frac{\Phi}    \def\qqz{\Psi_N}{\Phi}    \def\qqz{\Psi_{N- 1}}\right)^\alpha}  \def\bz{\beta- 1 \right] \left(\frac{\Phi}    \def\qqz{\Psi_N}{\Phi}    \def\qqz{\Psi_{N- 1}}\right)^{-1} &= (y^\alpha}  \def\bz{\beta- 1) y^{-1} \\
&\leqslant (1- \alpha}  \def\bz{\beta)^{1/ \alpha}  \def\bz{\beta}\left(\frac{\alpha}  \def\bz{\beta}{1- \alpha}  \def\bz{\beta}\right) \\
&\leqslant \frac{\alpha}  \def\bz{\beta}{1- \alpha}  \def\bz{\beta}.
\endaligned
$$
For any $1\leqslant m\leqslant N- 1$, assume the inequality (\ref{==5}) is true when $n= m$. Now consider $n=m- 1$. Let $y= \Phi}    \def\qqz{\Psi_{m}/ \Phi}    \def\qqz{\Psi_{m- 1}$, then $y\geqslant 1$. By the assumption, we have
$$
\aligned
\sum_{i= m}^N & \left[\left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_{m- 1}}\right)^\alpha}  \def\bz{\beta - \left(\frac{\Phi}    \def\qqz{\Psi_{i- 1}}{\Phi}    \def\qqz{\Psi_{m- 1}}\right)^\alpha}  \def\bz{\beta \right] \left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_{m- 1}}\right)^{-1} \\
&= \left(\frac{\Phi}    \def\qqz{\Psi_{m}}{\Phi}    \def\qqz{\Psi_{m- 1}}\right)^{\alpha}  \def\bz{\beta- 1} \sum_{i=m}^{N} \left[\left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_{m}}\right)^\alpha}  \def\bz{\beta - \left(\frac{\Phi}    \def\qqz{\Psi_{i- 1}}{\Phi}    \def\qqz{\Psi_{m}}\right)^\alpha}  \def\bz{\beta \right] \left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_{m}}\right)^{-1} \\
&\leqslant \left(\frac{\Phi}    \def\qqz{\Psi_{m}}{\Phi}    \def\qqz{\Psi_{m- 1}}\right)^{\alpha}  \def\bz{\beta- 1} \left[\frac{\alpha}  \def\bz{\beta}{1- \alpha}  \def\bz{\beta}+ 1- \left(\frac{\Phi}    \def\qqz{\Psi_{m- 1}}{\Phi}    \def\qqz{\Psi_{m}}\right)^\alpha}  \def\bz{\beta \right] \\
&= \frac{1}{1- \alpha}  \def\bz{\beta} y^{\alpha}  \def\bz{\beta- 1}- y^{-1}.
\endaligned
$$
Again, by simple calculations, we know the function 
$$
f(x)= \frac{1}{1- \alpha}  \def\bz{\beta} x^{\alpha}  \def\bz{\beta- 1}- x^{-1}, \qquad x\geqslant 1,
$$
reachs the maximum $\dfrac{\alpha}  \def\bz{\beta}{1- \alpha}  \def\bz{\beta}$ when $x= 1$. Hence
$$
\sum_{i= m}^N \left[\left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_{m- 1}}\right)^\alpha}  \def\bz{\beta - \left(\frac{\Phi}    \def\qqz{\Psi_{i- 1}}{\Phi}    \def\qqz{\Psi_{m- 1}}\right)^\alpha}  \def\bz{\beta \right] \left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_{m- 1}}\right)^{-1} \leqslant \frac{\alpha}  \def\bz{\beta}{1- \alpha}  \def\bz{\beta}.
$$
By induction, we prove this assertion. We should notice that the right side of (\ref{==5}) is independent of $N$, then (\ref{==5}) also be true when $N \rightarrow \infty$.

Let $\alpha}  \def\bz{\beta\in (0, 1)$ be an undetermined parameter. Having the inequality (\ref{==5}) in hand, we can prove:
\bg{equation}\label{==6}
\left(\sum_{i= n}^\infty u_i (H \mathbf{\hat{v}} (i))^{\alpha}  \def\bz{\beta q/p^*} \right)^{1/q}\leqslant B (H \mathbf{\hat{v}} (n))^{(\alpha}  \def\bz{\beta- 1)/p^*} \left(\frac{1}{1- \alpha}  \def\bz{\beta}\right)^{1/q}, \quad 1\leqslant n< \infty.
\end{equation}
With the definition of $B$, for any $n$, we have $(\sum_{i=n}^\infty u_i)^{1/q}\leqslant B (H \mathbf{\hat{v}} (n))^{-1/{p^*}}$. Let $\Phi}    \def\qqz{\Psi_n= (H \mathbf{\hat{v}} (n))^{q/{p^*}}$, summation by parts, we have
$$\aligned
\sum_{i= n}^\infty u_i \Phi}    \def\qqz{\Psi_i^\alpha}  \def\bz{\beta &= \Phi}    \def\qqz{\Psi_n^\alpha}  \def\bz{\beta \left(\sum_{i=n}^\infty u_i \right) + \sum_{i=n+ 1}^\infty \left(\Phi}    \def\qqz{\Psi_i^\alpha}  \def\bz{\beta- \Phi}    \def\qqz{\Psi_{i- 1}^\alpha}  \def\bz{\beta \right) \left(\sum_{j=i}^\infty u_j \right) \\
&\leqslant B^q \Phi}    \def\qqz{\Psi_n^{\alpha}  \def\bz{\beta- 1} + B^q \sum_{i=n+ 1}^\infty \left(\Phi}    \def\qqz{\Psi_i^\alpha}  \def\bz{\beta- \Phi}    \def\qqz{\Psi_{i- 1}^\alpha}  \def\bz{\beta \right) \Phi}    \def\qqz{\Psi_i^{-1} \\
&= B^q \Phi}    \def\qqz{\Psi_n^{\alpha}  \def\bz{\beta- 1} \left\{ 1+ \sum_{i=n+ 1}^\infty \left[\left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_n}\right)^\alpha}  \def\bz{\beta- \left(\frac{\Phi}    \def\qqz{\Psi_{i- 1}}{\Phi}    \def\qqz{\Psi_n}\right)^\alpha}  \def\bz{\beta \right] \left(\frac{\Phi}    \def\qqz{\Psi_i}{\Phi}    \def\qqz{\Psi_n}\right)^{-1} \right\} \\
&\leqslant B^q \Phi}    \def\qqz{\Psi_n^{\alpha}  \def\bz{\beta- 1} \left(\frac{1}{1- \alpha}  \def\bz{\beta} \right).
\endaligned$$
Making a power $1/q$, we obtain the inequality (\ref{==6}).

Now, we are ready to prove the upper bounds of the basic estimates. For any $n\in[1, \infty)$, Let 
$x_n = (H \mathbf{\hat{v}} (n))^{\alpha}  \def\bz{\beta} - (H \mathbf{\hat{v}} (n- 1))^{\alpha}  \def\bz{\beta}$, we have
$$
\aligned
{I_n^*(\mathbf{x})}^{1/p^*} &= \left[\frac{\hat{v}_n}{(H \mathbf{\hat{v}} (n))^{\alpha}  \def\bz{\beta} - (H \mathbf{\hat{v}} (n- 1))^{\alpha}  \def\bz{\beta}} \left(\sum_{i= n}^\infty u_i (H \mathbf{\hat{v}} (i))^{\alpha}  \def\bz{\beta q /p^*} \right)^{{p^*}/q} \right]^{1/{p^*}} \\
&\leqslant \left[\frac{1}{\alpha}  \def\bz{\beta (H \mathbf{\hat{v}} (n))^{\alpha}  \def\bz{\beta- 1} } \left(\sum_{i= n}^\infty u_i (H \mathbf{\hat{v}} (i))^{\alpha}  \def\bz{\beta q/ p^*} \right)^{{p^*}/q} \right]^{1/{p^*}} \\
&= \alpha}  \def\bz{\beta^{-1/{p^*}} \left(H \mathbf{\hat{v}} (n) \right)^{(1- \alpha}  \def\bz{\beta)/{p^*}} \left(\sum_{i= n}^\infty u_i (H \mathbf{\hat{v}} (i))^{\alpha}  \def\bz{\beta q/{p^*}} \right)^{1/q} \\
&\leqslant B \alpha}  \def\bz{\beta^{-1/{p^*}} (1- \alpha}  \def\bz{\beta)^{-1/q}.
\endaligned
$$
An easy calculation shows that the function
$$
f(x)= x^{-1/{p^*}} (1- x)^{-1/q}, \qquad 0< x<1, 
$$
reaches the maximum
$$
\tilde{k}_{q, p}= \left(1+ \frac{q}{p^*}\right)^{1/q} \left(1+ \frac{p^*}{q}\right)^{1/p^*}
$$
when $x= \dfrac{q}{p^*+ q}$. Hence we take $\alpha}  \def\bz{\beta= \dfrac{q}{p^*+ q}$, then we have
$$
\sup_{n\in [1, \infty)} \left( I_n^* (\mathbf{x}) \right)^{1/p^*}\leqslant \tilde{k}_{q, p} B.
$$
By Theorem \ref{Var}, we get the basic upper estimates.

The basic lower estimates are more straightforward. For any $n\in [1, \infty)$, we can choose a test sequence as 
$$
x_i^{(n)}= 
\begin{cases}
\hat{v}_i, & 1\leqslant i\leqslant n, \\
0, & n< i< \infty. 
\end{cases}
$$
It is obvious that $\mathbf{x}^{(n)} \in {\scr A}}   \def\bq{{\scr B}_0 [1, \infty)$, then by (\ref{A}) we have
$$
\aligned
A&\geqslant \sup_{n\in [1, \infty)} \frac{\left[\sum_{i= 1}^\infty u_i (H \mathbf{x}^{(n)} (i))^q \right]^{1/q}}{\left[ \sum_{i= 1}^\infty v_i \left(x_i^{(n)} \right)^p \right]^{1/p}} \\
&= \sup_{n\in [1, \infty)} \left(\sum_{i= 1}^n \hat{v}_i \right)^{1/p^*} \left[ \left(\sum_{i= 1}^n \hat{v}_i \right)^{-q} \left( \sum_{i= 1}^{n- 1} u_i \left(H \mathbf{\hat{v}} (i) \right)^q \right) + \sum_{i= n}^\infty u_i \right]^{1/q} \\
&\geqslant B.
\endaligned
$$
This completes the proof of Corollary \ref{basic estimates}.
\quad $\square$ \medskip

\medskip
\noindent {\bf Proof of Corollary \ref{approximating}.} By the proportional property, we can obtain the monotonicity of \{$\dz_n$\}. The approximating sequence $\{ \dz_n \}$ comes from the upper estimates of the variational formula, and $\{\widetilde{\dz}_n\}$ comes from the lower one. These results are the simple applications of Theorem \ref{Var}. The sequence $\{\overline{\dz}_n\}$ is the straightforward application of the classical variational formula (\ref{A}).
\quad $\square$ \medskip

\section{Proof of Theorem \ref{kqpB}}

Again, we assume $N= \infty$, and the case of finite interval will be discussed in Section \ref{Interval}. We begin with the following well-known lemma. 

\nnd\begin{lmm1}\label{compare}
{\cms
Let $\mathbf{a}$, $\mathbf{b}$ be sequences with non-negative entries. If
$$
\sum_{k=i}^\infty a_k \leqslant \sum_{k=i}^\infty b_k, \qquad (\forall ~ i=1,2, \cdots)
$$
then for any increasing non-negative sequence $\mathbf{c}$, we have
$$
\sum_{k=1}^\infty a_k c_k \leqslant \sum_{k=1}^\infty b_k c_k.
$$
}
\end{lmm1}

\medskip \noindent {\bf Proof}. 
Set $c_0= 0$. Summation by parts, we have
$$
\aligned
\sum_{k=1}^\infty a_k c_k &= \sum_{k=1}^\infty \left(\sum_{i=1}^k (c_i- c_{i-1}) \right) a_k = \sum_{i=1}^\infty \left(\sum_{k=i}^\infty a_k \right) (c_i- c_{i-1})  \\
&\leqslant \sum_{i=1}^\infty \left(\sum_{k=i}^\infty b_k \right) (c_i- c_{i-1})= \sum_{k=1}^\infty b_k c_k.
\endaligned
$$ 
This completes the proof of Lemma \ref{compare}. 
\quad $\square$ \medskip

The conclusion of Lemma \ref{compare} is about increasing sequences, analogously, there is a conclusion corresponding to the decreasing sequences, cf. \rf{Bennett2}{Lemma 1}. The following lemma is due to Bliss \cite{Bliss}, which gives a special Hardy-type inequality in the continuous case.
\nnd\begin{lmm1}\label{Bliss_lemma}
{\cms
For any non-negative real function $f(x)$, we have
\bg{equation}
\left(\int_0^\infty \frac{1}{x^{q- r}} \left(\int_0^x f(t) \text{\rm d} t \right)^q \text{\rm d} x \right)^{1/q}\leqslant k_{q,p} \left(\frac{p^*}{q} \right)^{1/q} \left(\int_0^\infty f^p(x) \text{\rm d} x \right)^{1/p}, 
\end{equation}
where $r= q/p -1$ and $k_{q,p}$ is the optimal constant, which is defined as (\ref{k_qp}). Moreover, the optimal constant is attained when 
$$
f(x)= \frac{c}{(d \cdot x^r + 1)^{(r+ 1)/r}}, 
$$
where $c$ and $d$ are non-negative constants. 
}
\end{lmm1}

\medskip
\noindent { {\bf Proof of Theorem \ref{kqpB}}.  By Corollary \ref{basic estimates}, it is obvious that $A= \infty$ if $B= \infty$. To avoid this trivial case, we assume $B< \infty$. 

(a) First we consider the case that $H \mathbf{\hat{v}} (\infty)= \lim_{n \rightarrow \infty} H \mathbf{\hat{v}} (n) =\infty$. Similar to Proposition \ref{decreasing}, we can rewrite the Hardy-type inequalities (\ref{Hardy}) as:
$$
\sum_{n=1}^\infty u_n \left(\sum_{i=1}^n \hat{v}_i x_i \right)^q \leqslant A^q \left(\sum_{n=1}^\infty \hat{v}_n x_n^p \right)^{q/p}.
$$

Define sequence $\mathbf{\tilde{u}}$
\bg{equation}\label{tilde_u}
\tilde{u}_n = B^q \left( (H \mathbf{\hat{v}} (n))^{-q/p^*}- (H \mathbf{\hat{v}} (n+1))^{-q/p^*} \right), \qquad n\geqslant 1. 
\end{equation}
By direct summation and $H \mathbf{\hat{v}} (\infty)= \infty$, we have 
$$
\sum_{i= n}^\infty \tilde{u}_i = B^q \left(H \mathbf{\hat{v}} (n)^{-q/p^*} - H \mathbf{\hat{v}} (\infty)^{-q/p^*} \right)= B^q H \mathbf{\hat{v}} (n)^{-q/p^*} \geqslant \sum_{i= n}^\infty u_i. 
$$
Applying Lemma \ref{compare}, for any non-negative sequence $\mathbf{x}$, we obtain
\bg{equation}\label{kqp1}
\sum_{n= 1}^\infty u_n \left(\sum_{i=1}^n \hat{v}_i x_i \right)^q \leqslant \sum_{n= 1}^\infty \tilde{u}_n \left(\sum_{i=1}^n \hat{v}_i x_i \right)^q. 
\end{equation}

The next is to show
$$
\sum_{n= 1}^\infty \tilde{u}_n \left(\sum_{i=1}^n \hat{v}_i x_i \right)^q \leqslant A^q \left(\sum_{n=1}^\infty  \hat{v}_n x_n^p \right)^{q/p}. 
$$
In order to use Lemma \ref{Bliss_lemma}, we should construct a function which connects summation with integration. Defined function $f: [0, \infty) \rightarrow [0, \infty)$
\bg{equation}\label{prp1}
f(x)= 
\begin{cases}
x_n, & H \mathbf{\hat{v}} (n-1) \leqslant x< H \mathbf{\hat{v}} (n), \\
0, & x\geqslant \sup_n H \mathbf{\hat{v}} (n).
\end{cases}
\end{equation}
It is clear that 
\bg{equation}\label{prp2}
\sum_{i=1}^\infty \hat{v}_i x_i^p = \int_{0}^\infty f^p(x) \text{\rm d} x, 
\end{equation}
and
\bg{equation}\label{prp3}
\sum_{i=1}^n \hat{v}_i x_i\leqslant \int_0^\alpha}  \def\bz{\beta f(t) \text{\rm d} t,
\end{equation}
where $H \mathbf{\hat{v}} (n) \leqslant \alpha}  \def\bz{\beta< H \mathbf{\hat{v}} (n+ 1)$. 

For convenience, write $\tilde{u}_0= 0, \hat{v}_0= 0$. Applying (\ref{prp3}), Lemma \ref{Bliss_lemma} and (\ref{prp2}), we see that
\begin{align}
\allowdisplaybreaks
\sum_{n=0}^\infty \tilde{u}_n \left(\sum_{k=0}^n \hat{v}_k x_k \right)^q &= \sum_{n=0}^\infty B^q \left( (H \mathbf{\hat{v}} (n))^{-q/p^*}- (H \mathbf{\hat{v}} (n+ 1))^{-q/p^*} \right) \left(\sum_{k= 0}^n \hat{v}_k x_k \right)^q \notag \\
&= \sum_{n=0}^\infty B^q \left(\frac{q}{p^*} \int_{H \mathbf{\hat{v}} (n)}^{H \mathbf{\hat{v}} (n+ 1)} x^{-q/p^* - 1} \text{\rm d} x\right) \left(\sum_{k= 0}^n \hat{v}_k x_k \right)^q \notag \\
&\leqslant \frac{q}{p^*} B^q \left(\sum_{n=0}^\infty \int_{H \mathbf{\hat{v}} (n)}^{H \mathbf{\hat{v}} (n+ 1)} x^{-q/p^* - 1} \left(\int_0^x f(t) \text{\rm d} t\right)^q \text{\rm d} x \right) \notag \\
&= \frac{q}{p^*} B^q \int_0^\infty x^{r- q} \left(\int_0^x f(t) \text{\rm d} t \right)^q \text{\rm d} x \quad (r=q/p - 1)\notag \\
&\leqslant B^q k_{q, p}^q \left(\int_0^\infty f^p(x) \text{\rm d} x \right)^{q/p} \notag \\
&= B^q k_{q, p}^q \left(\sum_{i=0}^\infty \hat{v}_i x_i^p \right)^{q/p}. \notag
\end{align}
By the definition of the optimal constant, we have $A\leqslant k_{q, p} B$.

(b) To show the factor of the basic upper estimate is best possible, we attempt to mimic the extremal function in Lemma \ref{Bliss_lemma}: 
$$
f(x)= \frac{c}{\left(dx^r+ 1\right)^{(r+1)/r}} \quad \text{and} \quad \int_0^x f(t) \text{\rm d} t = \frac{cx}{\left(dx^r +1 \right)^{1/r}},
$$
where $c$ and $d$ are arbitrary positive constants. We start with this form, set
$$
u_n= n^{-q/p^*}- (n+1)^{-q/p^*}, \qquad v_n \equiv 1,
$$
and 
$$
x_n = \frac{c n}{(n^r+ d)^{1/r}}- \frac{c (n-1)}{((n-1)^r+ d)^{1/r}}. 
$$
Obviously, the form of $\mathbf{x}$ comes from the difference of $\int_0^x f(t) \text{\rm d} t$. In this case, we have $B= 1$. Here we are free to choose $c$ and $d$, however, no matter what choices, there is some loss of precision between integrals and series. But by direct calculation, we find that this loss becomes negligible when $c/d \rightarrow 0$. Without loss of generality, we choose $c= 1$ and let $d$ be a positive and large enough real number. Next, we calculate the left and the right side of (\ref{Hardy}). 

The calculation of the right side of (\ref{Hardy}) is direct. By the definition of $\mathbf{x}$, we have
\begin{align}\label{right}
\sum_{n= 1}^\infty x_n^p &= \sum_{n=1}^\infty \left[\frac{n}{(n^r + d)^{1/r}} -\frac{n- 1}{((n- 1)^r + d)^{1/r}} \right]^p \notag \\
&= \sum_{n=1}^\infty \left[\int_{n- 1}^n \frac{d}{(x^r+ d)^{1/r +1}} \text{\rm d} x \right]^p \notag \\
&\leqslant \sum_{n= 1}^\infty \int_{n- 1}^n \left(\frac{d}{(x^r+ d)^{1/r +1}} \right)^p \text{\rm d} x \notag \\
&= r^{-1} d^{(1- p)/r} B\left(\frac{1}{r}, \frac{q- 1}{r}\right).
\end{align}

The left side of (\ref{Hardy}) is difficult. First, we assert there is a large enough integer $N$ such that 
\bg{equation}\label{left1}
\int_N^{\infty} x^{-q/{p^*}- 1} \left[ \frac{x^q}{(x^r+ d)^{q/r}} \right] \text{\rm d} x \leqslant \int_1^{\infty} (x+ 1)^{-q/{p^*}- 1} \left[ \frac{x^q}{(x^r+ d)^{q/r}} \right] \text{\rm d} x. 
\end{equation}
In fact, we have
$$
\int_N^{\infty}  x^{-q/{p^*}- 1} \left[ \frac{x^q}{(x^r+ d)^{q/r}} \right] \text{\rm d} x \leqslant \int_N^{\infty} x^{-q/{p^*}- 1} \text{\rm d} x = \frac{p^*}{q} N^{-q/{p^*}}. 
$$
The existence of $N$ is obvious since the left side of (\ref{left1}) decreases to $0$ as $N \uparrow \infty$. Fix this sufficiently large integer $N$, then the left side of (\ref{left1}) is calculable. Using the integral transform $s^{-1}= d^{-1} x^r +1$, we have 
\bg{equation}\label{left2}
\int_N^{\infty} x^{-q/{p^*}- 1} \left[ \frac{x^q}{(x^r+ d)^{q/r}} \right] \text{\rm d} x = r^{-1} d^{- \frac{q}{rp^*}} B\left(\frac{1+r}{r}, \frac{q- r- 1}{r}, \frac{d}{N^r+ d}\right),
\end{equation}
where $B(a, b, x)$ is the incomplete Beta function: 
$$
B(a, b, x)= \int_{0}^{x} s^{a- 1} (1- s)^{b- 1} \text{\rm d} s.
$$
Applying the mean value theorem, (\ref{left1}) and (\ref{left2}), we have
\begin{align}\label{left3}
\int_1^\infty & \left[x^{-q/{p^*}}- (x+ 1)^{-q/{p^*}} \right] \frac{x^q}{(x^r+ d)^{q/r}} \text{\rm d} x \notag \\
& \geqslant \int_1^\infty \frac{q}{p^*} (x+ 1)^{-q/{p^*}- 1} \left[\frac{x^q}{(x^r+ d)^{q/r}} \right] \text{\rm d} x \notag \\
& \geqslant d^{- \frac{q}{rp^*}} \left(\frac{q}{p^*} \right) r^{-1} B\left(\frac{1+r}{r}, \frac{q- r- 1}{r}, \frac{d}{N^r+ d}\right). 
\end{align}

Now, it's very easy to calculate the optimal constants. Using the relation
$$
B(a+ 1, b- 1)= \frac{a}{b- 1} B(a, b),
$$
it follows from (\ref{A}), (\ref{right}) and (\ref{left3}) that
$$
\aligned
A^q &\geqslant \left[\sum_{n=1}^{\infty} \left(n^{-q/p^*}- (n+ 1)^{-q/p^*} \right) (H \mathbf{x} (n))^q \right] \left( \sum_{n=1}^\infty x_n^p \right)^{- q/p} \\
&\geqslant \left[\int_{1}^{\infty} \left[x^{-q/p^*} - (x+ 1)^{-q/p^*} \right] \frac{x^q}{(x^r+ d)^{\frac{q}{r}}} \text{\rm d} x \right] \left( \sum_{n=1}^\infty x_n^p \right)^{- q/p} \\
&\geqslant \left(\frac{q}{p^*}\right) r^{q/p - 1} \cdot B\left(\frac{1+ r}{r}, \frac{q- 1- r}{r}, \frac{d}{N^r+ d} \right) \cdot B\left(\frac{1}{r}, \frac{q- 1}{r}\right)^{- q/p} \\
&\rightarrow k_{q,p}^q \qquad \text{(as $d \rightarrow \infty$)}. 
\endaligned
$$
Hence, the factor of basic upper estimate is best possible. 

(c) The final step is to remove condition $H \mathbf{\hat{v}} (\infty)= \infty$. We use Proposition \ref{C. estimates} of section \ref{Interval}. Fix $N_0 < \infty$. For given $\mathbf{u}$ and $\mathbf{v}$ on $[1, \infty)$, we define $\mathbf{u}^{N_0}$ and $\mathbf{v}^{N_0}$ to be the restriction of $\mathbf{u}$ and $\mathbf{v}$ on $[1, N_0]$. Then define 
$$
\overline{u}_n =
\begin{cases}
u_n^{N_0}, & 1\leqslant n \leqslant N_0, \\
0, & n> N_0, 
\end{cases}
$$
and 
$$
\overline{v}_n =
\begin{cases}
v_n^{N_0}, & 1\leqslant n\leqslant N_0, \\
1, & n> N_0.
\end{cases}
$$
Obviously, we have $H \mathbf{\overline{v}} (\infty) = \infty$. Applying the result of (a), we have 
$$
A(\mathbf{\overline{u}}, \mathbf{\overline{v}}) \leqslant k_{q, p} B(\mathbf{\overline{u}}, \mathbf{\overline{v}}).
$$
By Proposition \ref{C. estimates}, we get
$$
A(\mathbf{u}^{N_0}, \mathbf{v}^{N_0}) \leqslant k_{q, p} B(\mathbf{u}^{N_0}, \mathbf{v}^{N_0}).
$$
The assertion follows by letting $N_0\rightarrow \infty$. This completes the proof of Theorem \ref{kqpB} in the case $N= \infty$. 
\quad $\square$ \medskip

Review part (a) of the proof of Theorem (\ref{kqpB}), when 
\bg{equation}\label{v}
N= \infty, \qquad H \mathbf{\hat{v}} (\infty) =\infty,
\end{equation}
we give the method to construct $\mathbf{u}$ from $\mathbf{v}$ such that the Hardy-type inequalities (\ref{Hardy})  hold with these $\mathbf{u}$ and $\mathbf{v}$. The part (b) show the optimal constant reaches the upper bound of the basic estimate. It means that the basic upper estimate with the improved factor $k_{q,p}$ holds for a large class of $(\mathbf{u}, \mathbf{v})$. The original idea of this construction is from Chen \rf{Chen3}{Proposition 4.5}. To distinguish it from Theorem \ref{kqpB}, we give the following proposition. 

\nnd\begin{prp1}\label{Bliss_prp}
{\cms
For any positive sequences $\mathbf{v}$ and constant $0< C< \infty$, the discrete Hardy-type inequalities (\ref{Hardy}) hold on $[1, \infty)$ with 
\bg{equation}
\tilde{u}_n = C^q \left( (H \mathbf{\hat{v}} (n))^{-q/p^*}- (H \mathbf{\hat{v}} (n+1))^{-q/p^*} \right), \qquad n\geqslant 1, 
\end{equation}
and its optimal constant $A$ satisfies 
\bg{equation}
A \leqslant k_{q, p} C, 
\end{equation}
where $k_{q,p}$ is defined as (\ref{k_qp}). Moreover, when $N$ and $\mathbf{\hat{v}}$ satisfy (\ref{v}), the upper bound is sharp with $C= B$.
}
\end{prp1}

\section{Hardy-type Inequalities on Interval}\label{Interval}
In this section, we study the comparison results of the optimal constants and the basic estimates on different intervals. In continuous case, the corresponding comparison results have been done by Chen \rf{Chen3}{Appendix}. With these results, we can get the complete proofs of Theorem \ref{Var} and Theorem \ref{kqpB}. 

Before specific discussion, we need some notations. Fix two natural numbers $N$ and $N'$ with $N< N'$. Given two positive sequences $\mathbf{u}$ and $\mathbf{v}$ on $[1, N]$, we can extend them to $[1, N']$ as follows:
\bg{equation}\label{u'}
u'_i =
\begin{cases}
u_i, & 1\leqslant i\leqslant N, \\
0, & N< i\leqslant N';
\end{cases}
\end{equation}
\bg{equation}\label{v'}
v'_i =
\begin{cases}
v_i, & 1\leqslant i\leqslant N, \\
\#, & N< i\leqslant N',
\end{cases}
\end{equation}
where $\#$ means arbitrary positive numbers. Denote $A_N(\mathbf{u}, \mathbf{v})$ be the optimal constant of the Hardy-type inequalities (\ref{Hardy}) in the interval $[1, N]$ with sequences $\mathbf{u}$ and $\mathbf{v}$, and similar for $B_N(\mathbf{u}, \mathbf{v})$. 

The first result is a comparison for the optimal constants on different intervals. 

\nnd\begin{prp1}\label{C. constant}
{\cms
Given two positive sequences $\mathbf{u}'$ and $\mathbf{v}'$ on $[1, N']$. Use $\mathbf{u}$ and $\mathbf{v}$ to denote their restrictions to $[1, N]$. Then we have $A_N(\mathbf{u}, \mathbf{v}) \uparrow A_{N'}(\mathbf{u}', \mathbf{v}')$ as $N \uparrow {N'} \leqslant \infty$. 

In particular, if the inequality (\ref{Hardy}) holds on $[1, N']$, then it also holds with the same constant $A_{N'}(\mathbf{u}', \mathbf{v}')$ on $[1, N]$. 
}
\end{prp1}

\medskip \noindent {\bf Proof}. 
(a) Given a non-negative sequence $\mathbf{x}$ on $[1, N]$, we can extend to $[1, N']$ by setting 
\bg{equation}\label{extension of x}
x'_i =
\begin{cases}
x_i & 1\leqslant i\leqslant N, \\
0 & N< i\leqslant N'.
\end{cases} 
\end{equation}
Then we have 
$$
\aligned
\left[ \sum_{n=1}^N u_n \left(H \mathbf{x} (n) \right)^q \right]^{1/q} &= \left[ \sum_{n=1}^{N'} u'_n \left(H \mathbf{x'} (n) \right)^q \right]^{1/q} \\ 
&\leqslant A_{N'}(\mathbf{u}', \mathbf{v}') \left[ \sum_{n=1}^{N'} v'_n {x'}_n^p \right]^{1/p} \\
&= A_{N'}(\mathbf{u}', \mathbf{v}') \left[ \sum_{n=1}^{N} v_n {x}_n^p \right]^{1/p}.
\endaligned
$$
It means that $A_N(\mathbf{u}, \mathbf{v}) \leqslant A_{N'}(\mathbf{u}', \mathbf{v}')$. 

(b) Our next goal is to show the convergence. First we consider the case that $\sum_{n=1}^{N'} u'_n = \infty$. Clearly, in this case we have $N'= \infty$ and $A_{N'}(\mathbf{u}', \mathbf{v}')= \infty$. Besides, restricting to $[1, n]$ and choosing $\mathbf{x}= (1, 0, \dots, 0)$, we obtain 
$$
A_n(\mathbf{u}, \mathbf{v}) \geqslant \left(\sum_{i=1}^n u_i \right)^{1/q} v_1^{- 1/p} \rightarrow \infty, \qquad \text{as $n\rightarrow \infty$}.
$$
Hence the convergence holds in this case. 

(c) Let $\sum_{n=1}^{N'} u'_n < \infty$. For every non-negative sequence $\mathbf{x}$ on $[1, N']$ with $\sum_{n = 1}^{N'} v'_n x_n^p < \infty$, we get
$$
\frac{\left[\sum_{n= 1}^N u_n (H \mathbf{x} (n))^q \right]^{1/q}}{\left(\sum_{n= 1}^N v_n x_n^p \right)^{1/p}} \rightarrow \frac{\left[\sum_{n= 1}^{N'} u'_n (H \mathbf{x} (n))^q \right]^{1/q}}{ \left( \sum_{n= 1}^{N'} v'_n x_n^p \right)^{1/p}} \leqslant A_{N'}(\mathbf{u}', \mathbf{v}'),
$$
as $N \uparrow N'$. With (\ref{A}), for every $\varepsilon} \def\xz{\xi> 0$, we can choose a sequence $\mathbf{x}$ such that
$$
A_{N'}(\mathbf{u}', \mathbf{v}') \leqslant \frac{\left[\sum_{n= 1}^{N'} u'_n (H \mathbf{x} (n))^q \right]^{1/q}}{\left(\sum_{n= 1}^{N'} v'_n x_n^p \right)^{1/p}} + \varepsilon} \def\xz{\xi.
$$
Then we can choose $N$ closed to $N'$ such that
$$
\frac{\left[\sum_{n= 1}^{N'} u'_n (H \mathbf{x} (n))^q \right]^{1/q}}{\left(\sum_{n= 1}^{N'} v'_n x_n^p \right)^{1/p}} \leqslant \frac{\left[\sum_{n= 1}^N u_n (H \mathbf{x} (n))^q \right]^{1/q}}{\left(\sum_{n= 1}^N v_n x_n^p \right)^{1/p}} + \varepsilon} \def\xz{\xi.
$$
Hence, we get
$$
A_N(\mathbf{u}, \mathbf{v}) \leqslant A_{N'}(\mathbf{u}', \mathbf{v}') \leqslant \frac{\left[\sum_{n= 1}^N u_n (H \mathbf{x} (n))^q \right]^{1/q}}{\left(\sum_{n= 1}^N v_n x_n^p \right)^{1/p}} + 2\varepsilon} \def\xz{\xi \leqslant A_N(\mathbf{u}, \mathbf{v}) + 2\varepsilon} \def\xz{\xi.
$$
It means that the convergence holds.
\quad $\square$ \medskip

The following result is about the factor in the basic estimates. 

\nnd\begin{prp1}\label{C. estimates}
{\cms
Given two positive sequences $\mathbf{u}$ and $\mathbf{v}$ on $[1, N]$, $\mathbf{u}'$ and $\mathbf{v}'$, defined by (\ref{u'}) and (\ref{v'}), are the extensions on $[1, N']$. Suppose that $A_{N'}(\mathbf{u}', \mathbf{v}')\leqslant k B_{N'}(\mathbf{u}', \mathbf{v}')$ for a universal constant $k$, then we have $A_N(\mathbf{u}, \mathbf{v})\leqslant k B_N(\mathbf{u}, \mathbf{v})$. 
}
\end{prp1}

\medskip \noindent {\bf Proof}. 
Given a sequence $\mathbf{x}$ on $[1, N]$, we can extend it from $[1, N]$ to $[1, N']$ by (\ref{extension of x}). Then we have
$$
\aligned
\left( \sum_{n=1}^N u_n (H \mathbf{x} (n))^q \right)^{1/q} &= \left( \sum_{n=1}^{N'} u'_n (H \mathbf{x}' (n))^q \right)^{1/q} \\
&\leqslant A_{N'}(\mathbf{u}', \mathbf{v}') \left( \sum_{n= 1}^{N'} v'_n {x'}_n^p \right)^{1/p} \\
&\leqslant k B_{N'}(\mathbf{u}', \mathbf{v}') \left( \sum_{n= 1}^{N'} v'_n {x'}_n^p \right)^{1/p} \\
&= k B_{N'}(\mathbf{u}', \mathbf{v}') \left( \sum_{n= 1}^{N} v_n x_n^p \right)^{1/p}. 
\endaligned
$$
With the definition of the extensions (\ref{u'}) and (\ref{v'}), we can easily check that 
$$
B_{N'}(\mathbf{u}', \mathbf{v}')= B_{N}(\mathbf{u}, \mathbf{v}).
$$
It follows that 
$$
\left( \sum_{n=1}^N u_n (H \mathbf{x} (n))^q \right)^{1/q} \leqslant k B_{N}(\mathbf{u}, \mathbf{v}) \left( \sum_{n= 1}^{N} v_n x_n^p \right)^{1/p}. 
$$
Hence $A_N(\mathbf{u}, \mathbf{v})\leqslant k B_N(\mathbf{u}, \mathbf{v})$ as required. 
\quad $\square$ \medskip

With the help of Proposition \ref{C. constant} and Proposition \ref{C. estimates}, we know the variational formulas of the optimal constants, the basic estimates and the improved factor of the basic upper estimates are true when $N< \infty$. So far, we complete the proofs of our main results. 

The following result gives an opposite view of Proposition \ref{C. constant}: from some local sub-intervals to the whole interval. It gives us an approximating procedure for the unbounded interval. 

\nnd\begin{prp1}\label{C. interval}
{\cms
Given two positive sequences $\mathbf{u}$ and $\mathbf{v}$ on $[1, N]$, extend them to $[1, N']$ by (\ref{u'}) and (\ref{v'}). Then we have $A_N (\mathbf{u}, \mathbf{v}) = A_{N'} (\mathbf{u}', \mathbf{v}')$. 
}
\end{prp1}

\medskip \noindent {\bf Proof}. 
For any sequence $\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, N]$, let $\mathbf{x}'$ be the extension of $\mathbf{x}$ from $[1, N]$ to $[1, N']$ by (\ref{extension of x}). Obviously, we have $\mathbf{x}' \in {\scr A}}   \def\bq{{\scr B}[1, N']$. The inequalities in $[1, N']$ are
$$
\| H \mathbf{x}' \|_{l^q(u')} \leqslant A_{N'} (\mathbf{u}', \mathbf{v}') \| \mathbf{x}' \|_{l^p(v')}.
$$
With (\ref{u'}), (\ref{v'}) and (\ref{extension of x}), it follows that
$$
\| H \mathbf{x} \|_{l^q(u)} \leqslant A_{N'} (\mathbf{u}', \mathbf{v}') \| \mathbf{x} \|_{l^p(v)}.  
$$
Because $\mathbf{x}$ is arbitrary, it implies that $A_N (\mathbf{u}, \mathbf{v}) \leqslant A_{N'} (\mathbf{u}', \mathbf{v}')$. 

Conversely, for any $\mathbf{x} \in {\scr A}}   \def\bq{{\scr B}[1, N']$, we have
$$
\aligned
\left(\sum_{n=1}^{N'} {u'}_n (H \mathbf{x} (n))^q \right)^{1/q} &= \left(\sum_{n=1}^{N} u_n (H \mathbf{x} (n))^q \right)^{1/q} \\
&\leqslant A_N (\mathbf{u}, \mathbf{v}) \left( \sum_{n=1}^N v_n \mathbf{x}_n^p \right)^{1/p} \\
&\leqslant A_N (\mathbf{u}, \mathbf{v}) \left( \sum_{n=1}^{N'} v'_n \mathbf{x}_n^p \right)^{1/p}. 
\endaligned
$$
This implies that $A_{N'} (\mathbf{u}', \mathbf{v}') \leqslant A_N (\mathbf{u}, \mathbf{v})$ and then the equality holds.
\quad $\square$ \medskip

\section{Examples}\label{examples}

As mentioned in introduction, Hardy-type inequalities play important role in probability theory. The first example is from birth-death processes which is standard having constant birth and death rates, cf. \rf{Chen5}{Example 5.3}. We present this example to illustrate the power of out results. 
\nnd\begin{xmp1}\label{example1}
{\cms
Let $p=q=2$ and $N= \infty$. For $n \geqslant 1$, let $u_n = \gamma}  \def\kz{\kappa^n$, $v_n= b \gamma}  \def\kz{\kappa^n$, where $\gamma}  \def\kz{\kappa$ and $b$ are constants with $\gamma}  \def\kz{\kappa< 1$ and $b> 0$. Then
$$
B < \widetilde{\dz}_1= \overline{\dz}_1 < A = \dz_1 < 2B,
$$
where $B= \dfrac{1}{\sqrt{b} (1- \gamma}  \def\kz{\kappa)}$, $\widetilde{\dz}_1= \overline{\dz}_1= \dfrac{\sqrt{1+ \gamma}  \def\kz{\kappa}}{\sqrt{b} (1- \gamma}  \def\kz{\kappa)}$, $A = \dz_1= \dfrac{1}{\sqrt{b}(1- \sqrt{\gamma}  \def\kz{\kappa})}$. Moreover, the optimal constant is attained at sequence
$$
a_n= \gamma}  \def\kz{\kappa^{(-n +1)/2} \left[n - (n- 1) \gamma}  \def\kz{\kappa^{1/2} \right], \qquad n\geqslant 1.
$$
}
\end{xmp1}

\medskip \noindent {\bf Proof}.  (a) First, $B$ is easy to calculate. By the definition, we have
$$
B= \sup_{n\in [1, \infty)}  \left(\sum_{i=1}^n b^{-1} \gamma}  \def\kz{\kappa^{- i} \right)^{1/2} \left(\sum_{j= n}^{\infty} \gamma}  \def\kz{\kappa^{j}\right)^{1/2} = \frac{1}{\sqrt{b} (1- \gamma}  \def\kz{\kappa)}. 
$$
Next, by (\ref{k_pp}), we have $k_{2, 2}= 2$. By Corollary \ref{basic estimates}, we obtain the basic estimates of the optimal constants: 
\bg{equation}\label{EX13}
\frac{1}{\sqrt{b} (1- \gamma}  \def\kz{\kappa)} \leqslant A \leqslant \frac{2}{\sqrt{b} (1- \gamma}  \def\kz{\kappa)}. 
\end{equation}

(b) To compute $\dz_1$, we use Corollary \ref{approximating}. Let 
$$
x_n^{(1)}= (H \mathbf{\hat{v}} (n))^{1/2}- (H \mathbf{\hat{v}} (n- 1))^{1/2}, \qquad n\geqslant 1,
$$
and then
$$
H \mathbf{x}^{(1)} (n) = (H \mathbf{\hat{v}} (n))^{1/2} = \left[\frac{\gamma}  \def\kz{\kappa^{-n}- 1}{b (1- \gamma}  \def\kz{\kappa)}\right]^{1/2}. 
$$
For convenience, we use $\fz_n= \gamma}  \def\kz{\kappa^{-n} -1$ in the following. By direct computations, we have
\begin{align}\label{EX11}
I\!I_n^* \left(\mathbf{x}^{(1)}\right) &= \frac{1}{H \mathbf{x}^{(1)} (n)} \sum_{i=1}^n \hat{v}_i \left(\sum_{j=i}^\infty u_j \left(H \mathbf{x}^{(1)} (j) \right) \right) \notag \\
&= \frac{b^{-3/2}}{(1- \gamma}  \def\kz{\kappa)^{1/2}} \frac{1}{H \mathbf{x}^{(1)} (n)} \sum_{i=1}^n \gamma}  \def\kz{\kappa^{-i}  \left(\sum_{j=i}^\infty \gamma}  \def\kz{\kappa^j \fz_j^{1/2} \right) \notag \\
&= \frac{\fz_n^{-1/2}}{b (1- \gamma}  \def\kz{\kappa)} \left[\sum_{j=1}^n \gamma}  \def\kz{\kappa^{j} \fz_j^{3/2} + \fz_n \sum_{j=n+ 1}^\infty \gamma}  \def\kz{\kappa^{j} \fz_j^{1/2}\right]. 
\end{align}
At the last step, we exchange the order of summation. From (\ref{EX11}), it is easy to check that $I\!I_n^* \left(\mathbf{x}^{(1)}\right)$ reaches the maximum when $n \rightarrow \infty$. Hence, by L'Hospital's rule, we obtain
$$
\aligned
\dz_{1}^2 &= \sup_{n\in [1, \infty)} I\!I_n^* \left(\mathbf{x}^{(1)} \right) \\
&= \frac{1}{b (1- \gamma}  \def\kz{\kappa)} \left[ \lim_{n\rightarrow \infty} \fz_n^{-1/2} \sum_{j=1}^n \gamma}  \def\kz{\kappa^{j} \fz_j^{3/2} + \lim_{n\rightarrow \infty} \fz_n^{1/2} \sum_{j=n+ 1}^\infty \gamma}  \def\kz{\kappa^{j} \fz_j^{1/2} \right] \\
&= \frac{1}{b (1- \gamma}  \def\kz{\kappa)} \left[ \frac{1}{1- \sqrt{\gamma}  \def\kz{\kappa} } + \frac{\sqrt{\gamma}  \def\kz{\kappa}}{1- \sqrt{\gamma}  \def\kz{\kappa}} \right] \\
&= \frac{1}{b (1- \sqrt{\gamma}  \def\kz{\kappa})^2}. 
\endaligned
$$

(c) Similarly, we use Corollary \ref{approximating} to compute $\overline{\dz}_1$ and $\widetilde{\dz}_1$. Fix $k> 0$, let 
$$
y_n^{(k, 1)}= 
\begin{cases}
b^{-1} \gamma}  \def\kz{\kappa^{-n}, & n\leqslant k, \\
0, & n>k,
\end{cases}
$$
and then 
$$
H \mathbf{y}^{(k, 1)} (n) = \frac{\gamma}  \def\kz{\kappa^{-(n \wedge k)} - 1}{b (1- \gamma}  \def\kz{\kappa)}= \frac{\fz_{n \wedge k}}{b (1- \gamma}  \def\kz{\kappa)}.
$$
By lots of tedious calculations, we have 
$$
\aligned
I\!I_n \left( \mathbf{y}^{(k, 1)} \right) &= \frac{1}{H \mathbf{y}^{(k, 1)} (n)} \sum_{i=1}^n \hat{v}_i \left(\sum_{j= i}^\infty u_j \left(H \mathbf{y}^{(k, 1)} \right) \right) \\
&= \frac{1}{b \fz_{n \wedge k}} \sum_{i=1}^n \gamma}  \def\kz{\kappa^{-i} \left(\sum_{j= i}^\infty \gamma}  \def\kz{\kappa^{j} \fz_{j \wedge k} \right) \\
&= \frac{1}{b (1- \gamma}  \def\kz{\kappa)} \left[\frac{1+ \gamma}  \def\kz{\kappa}{1- \gamma}  \def\kz{\kappa} - \frac{2 (n \wedge k)}{\fz_{n \wedge k}} + (k- n)\vee 0 - \frac{\gamma}  \def\kz{\kappa^{k+ 1} \fz_{(n- k) \vee 0}}{1- \gamma}  \def\kz{\kappa}  \right].
\endaligned
$$
Next, note that $I\!I_n \left( \mathbf{y}^{(k, 1)} \right)$ reaches the minimum when $n= k$, and then
$$
\aligned
\widetilde{\dz}_1^2 &= \sup_{k\in [1, \infty)} \inf_{n\in [1, \infty)} I\!I_n \left(\mathbf{y}^{(k, 1)} \right) \\
&= \sup_{k\in [1, \infty)} \frac{1}{b (1- \gamma}  \def\kz{\kappa)} \left(\frac{1+ \gamma}  \def\kz{\kappa}{1- \gamma}  \def\kz{\kappa} - \frac{2 k}{\fz_{k}} \right) \\
&= \frac{1}{b (1- \gamma}  \def\kz{\kappa)} \lim_{k\rightarrow \infty} \left(\frac{1+ \gamma}  \def\kz{\kappa}{1- \gamma}  \def\kz{\kappa} - \frac{2 k}{\fz_{k}} \right) \\
&= \frac{1+ \gamma}  \def\kz{\kappa}{b (1- \gamma}  \def\kz{\kappa)^2}.
\endaligned
$$

Now, we consider $\overline{\dz}_1$. Since
$$
\sum_{n=1}^\infty v_n \left(y_n^{(k, 1)} \right)^2 = \frac{\fz_k}{b(1- \gamma}  \def\kz{\kappa)}, 
$$
and 
$$
\sum_{n=1}^\infty u_n \left( H \mathbf{y}^{(k, 1)} (n) \right)^2 = \frac{1}{b^2 (1- \gamma}  \def\kz{\kappa)^2} \left[\sum_{n= 1}^k \gamma}  \def\kz{\kappa^n \fz_n^2 + \frac{\gamma}  \def\kz{\kappa^{k+ 1} \fz_k^2}{1- \gamma}  \def\kz{\kappa} \right], 
$$
we have 
$$
\aligned
\overline{\dz}_1^2&= \sup_{k\in [1, \infty)} \frac{1}{b(1- \gamma}  \def\kz{\kappa)} \left[\fz_k^{-1} \sum_{n=1}^k \gamma}  \def\kz{\kappa^n \fz_n^2 + \frac{\gamma}  \def\kz{\kappa- \gamma}  \def\kz{\kappa^{k+ 1}}{1- \gamma}  \def\kz{\kappa}  \right] \\
&= \frac{1}{b(1- \gamma}  \def\kz{\kappa)} \left[ \lim_{k\rightarrow \infty} \fz_k^{-1} \sum_{n=1}^k \gamma}  \def\kz{\kappa^n \fz_n^2 + \frac{\gamma}  \def\kz{\kappa }{1- \gamma}  \def\kz{\kappa}  \right] = \frac{1+ \gamma}  \def\kz{\kappa}{b (1- \gamma}  \def\kz{\kappa)^2}.
\endaligned
$$
In the last step, the L'Hospital's rule is used to calculate the limitation of $k$. 

(d) So far, by Corollary \ref{approximating}, we obtain the estimates of the optimal constants,  which is more precise than the basic estimates (\ref{EX13}) 
\bg{equation}\label{EX12}
\frac{\sqrt{1+ \gamma}  \def\kz{\kappa}}{\sqrt{b} (1- \gamma}  \def\kz{\kappa)} \leqslant A \leqslant \frac{1}{\sqrt{b} (1- \sqrt{\gamma}  \def\kz{\kappa})}. 
\end{equation}

In fact, the optimal constant can be accurately calculated. Let $a_n= \gamma}  \def\kz{\kappa^{(-n +1)/2} \left[n - (n- 1) \gamma}  \def\kz{\kappa^{1/2} \right] (n\geqslant 1)$, then
$$
H \mathbf{a} (n)= n \gamma}  \def\kz{\kappa^{(-n+ 1)/2}. 
$$
Here we want to use $\mathbf{a}$ instead of $\mathbf{y}^{(k, 1)}$ to get the lower estimates. However, it is easy to check that $\mathbf{a}$ is non-summability. It means that Theorem \ref{Var} is invalid. By the classical variational formula (\ref{A}) and the L'Hospital's rule, we have
$$
\aligned
A^2 &\geqslant \frac{\sum_{n= 1}^\infty \gamma}  \def\kz{\kappa^n a_n^2 }{\sum_{n=1}^\infty b \gamma}  \def\kz{\kappa^n \left(a_n - a_{n-1} \right)^2} \\
&= b^{-1} \lim_{n\rightarrow \infty} \frac{n^2}{\left[n- (n- 1) \gamma}  \def\kz{\kappa^{1/2} \right]^2} \\
&= \frac{1}{b (1- \sqrt{\gamma}  \def\kz{\kappa})^2}.
\endaligned
$$
As a consequence, we obtain $A = \dfrac{1}{\sqrt{b} (1- \sqrt{\gamma}  \def\kz{\kappa})}$.
\quad $\square$ \medskip

To distinguish the first example, the second one is about the nonlinear situation $p \neq q$, which is from proof of Theorem \ref{kqpB}. The optimal constant is clear in this example.  
\nnd\begin{xmp1}\label{example2}
{\cms
Let $p \neq q$ and $N= \infty$. For $n\geqslant 1$, let $u_n= n^{-q/{p^*}}- (n+ 1)^{-q/{p^*}}$, $v_n \equiv 1$. Then

(1) The optimal constant is $A= k_{q, p}$, which is attained at sequence $\mathbf{x}$: 
$$
x_n= \frac{cn}{(n^r+ d)^{1/r}} - \frac{c(n- 1)}{((n-1)^r+ d)^{1/r}}, \qquad n \geqslant 1,
$$
where $r=q/p - 1$, $k_{q, p}$ is defined as (\ref{k_qp}), $c$ and $d$ are arbitrary positive constants. 

(2) The basic estimates and the approximating procedure are
$$
B \leqslant \overline{\dz}_1 \vee \widetilde{\dz}_1 \leqslant A = k_{q, p} B \leqslant \dz_1, 
$$ 
where $B=1$, $\overline{\dz}_1 \geqslant 1$, $\widetilde{\dz}_1 \geqslant 1$ and $\dz_1 \leqslant \left(1+ \dfrac{q}{p^*} \right)^{1/q+ 1/p^*}$.
}
\end{xmp1}

\medskip \noindent {\bf Proof}.  The first part has been done in Theorem \ref{kqpB}. The remainder of this proof is to compute $\dz_1$, $\overline{\dz}_1$ and $\widetilde{\dz}_1$. 

To compute $\dz_1$, let
$$
x_n^{(1)}= n^{q/(p^*+ q)} - (n- 1)^{q/(p^*+ q)}, n\geqslant 1, 
$$
then we have 
$$
\aligned
I\!I_n^* \left(\mathbf{x}^{(1)}\right) &= n^{- \frac{q}{p^*+q}} \sum_{i= 1}^n \left\{\sum_{j=i}^{\infty} \left[j^{-\frac{q}{p^*}}- (j+ 1)^{- \frac{q}{p^*}} \right] j^{\frac{q^2}{p^*(p^*+ q)}} \right\}^{p^*/q} \\
&\leqslant n^{- \frac{q}{p^*+q}} \sum_{i= 1}^n \left\{ \left(\frac{q}{p^*} \right) \sum_{j=i}^{\infty} \int_j^{j+1} x^{\frac{q^2}{p^*(p^*+ q)}- \frac{q}{p^*}- 1} \text{\rm d} x \right\}^{p^*/q} \\
&= n^{- \frac{q}{p^*+q}} \sum_{i= 1}^n \left[ \left(\frac{p^*+ q}{p^*} \right) i^{- \frac{q}{p^*+ q}} \right]^{p^*/q} \\
&\leqslant n^{- \frac{q}{p^*+q}} \left(\frac{p^*+ q}{p^*} \right)^{p^*/q} \left(1+ \int_1^n x^{- \frac{p^*}{p^*+ q}} \text{\rm d} x\right) \\
&= \left(1+ \frac{q}{p^*} \right)^{p^*/q + 1}. \\
\endaligned
$$
Therefore, we obtain
$$
\dz_1 = \sup_{n\in [1, \infty)} \left[ I\!I_n^* \left(\mathbf{x}^{(1)} \right)\right]^{1/p^*} \leqslant \left(1+ \frac{q}{p^*} \right)^{1/q + 1/p^*}.
$$

To compute $\overline{\dz}_1$ and $\widetilde{\dz}_1$, let
$$
y_n^{(k, 1)} =
\begin{cases}
1, & n\leqslant k, \\
0, & n> k. 
\end{cases}
$$
Obviously, we have $H \mathbf{y}^{(k, 1)} (n) = n \wedge k$, $\|\mathbf{y}^{(k, 1)} \|_{l^p(v)} = k^{1/p}$ and 
$$
\aligned
\|H \mathbf{y}^{(k, 1)} \|_{l^q(u)} &= \left[ \sum_{n=1}^\infty u_n (n \wedge k)^q \right]^{1/q} \\
&= \left[ \sum_{n=1}^{k- 1} \left(n^{-\frac{q}{p^*}}- (n+1)^{-\frac{q}{p^*}} \right) n^q + k^{q/p} \right]^{1/q}.
\endaligned
$$
Hence, we obtain
$$
\aligned
\overline{\dz}_1 &= \sup_{k\in [1, \infty)} \frac{\|H \mathbf{y}^{(k, 1)} \|_{l^q(u)}}{\|\mathbf{y}^{(k, 1)} \|_{l^p(v)}} \\
&= \sup_{k\in [1, \infty)} \left[k^{-q/p} \sum_{n=1}^{k- 1} \left(n^{-q/p^*}- (n+1)^{-q/p^*} \right) n^q + 1 \right]^{1/q} \\ 
&\geqslant 1.
\endaligned
$$
Now, we consider $\widetilde{\dz}_1$. With directly calculating, we have
$$
\aligned
I\!I_n & \left(\mathbf{y}^{(k, 1)}\right) = \inf_{n\in [1, \infty)} \frac{1}{n\wedge k} \sum_{i=1}^n \left[\sum_{j=i}^{\infty} u_j \left(j \wedge k\right)^{q- 1} \right]^{p^*- 1}\\
&= \frac{1}{n \wedge k} \sum_{i=1}^{k \wedge n} \left[k^{q/p - 1} + \sum_{j=i}^{k-1} u_j j^{q-1} \right]^{p^*- 1} + \mathbbm{1}_{\{n> k\}} k^{\frac{q- p}{p- 1}} \sum_{i=k +1}^n i^{-q/p}. 
\endaligned
$$
Obviously, $I\!I_n \left(\mathbf{y}^{(k, 1)}\right)$ is increasing when $n\geqslant k$. It means that $I\!I_n \left(\mathbf{y}^{(k, 1)}\right)$ reaches its minimum at $n\in [1, k]$. Thus, we obtain
$$
\aligned
\inf_{n\in [1, \infty)} I\!I_n \left(\mathbf{y}^{(k, 1)}\right) &= \inf_{n\leqslant k} \frac{1}{n} \sum_{i=1}^{n} \left[k^{q/p - 1} + \sum_{j=i}^{k-1} u_j j^{q-1} \right]^{p^*- 1} \\
&\geqslant \inf_{n\leqslant k} \left[k^{q/p - 1} + \sum_{j=n}^{k-1} u_j j^{q-1} \right]^{p^*- 1} \\
&= k^{(q/p - 1)(p^*- 1)}.
\endaligned
$$
Therefore, we obtain
$$
\aligned
\widetilde{\dz}_1 &= \sup_{k\in [1, \infty)} k^{1/q- 1/p} \left( \inf_{n\in [1, \infty)} I\!I_n \left(\mathbf{y}^{(k, 1)}\right) \right)^{(p- 1)/q} \\
&\geqslant \sup_{k\in [1, \infty)} k^{1/q- 1/p} \left[k^{(q/p - 1)(p^*- 1)} \right]^{(p- 1)/q}= 1. \qquad \square
\endaligned
$$

\noindent {\bf Acknowledgements} $\quad$ This paper is based on the series of studies of my supervisor Prof. M. F. Chen. Heartfelt thanks are given to my supervisor for his careful guidance and helpful suggestions. Thanks are also given to Prof. Y. H. Mao, Prof. F. Y. Wang and Prof. Y. H. Zhang for their comments and suggestions, which lead to lots of improvements of this paper. 

The research is supported by NSFC (Grant No. 11131003) and by the ``985'' project from the Ministry of Education in China.


2014-06-24T02:11:17.000
1406
1406.1984
en
https://arxiv.org/abs/1406.1984
This paper studies the Hardy-type inequalities on the discrete intervals. The first result is the variational formulas of the optimal constants. Using these formulas, one may obtain an approximating procedure and the known basic estimates of the optimal constants. The second result, which is the main innovation of this paper, is about the factor of basic upper estimates. An improved factor is presented, which is smaller than the known one and is best possible. Some comparison results are included for comparing the optimal constants on different intervals.
Functional Analysis (math.FA)
Discrete Hardy-type Inequalities
Qwen/Qwen-72B
1. YES
2. YES