Identities on q-Harmonic Numbers

Mengxiao Zhou, Haitao Jin^*, Huanhuan Zheng

School of Science, Tianjin University of Technology and Education, Tianjin, China.

DOI: 10.4236/jamp.2024.125111                

Abstract

With the help of the classical Abel’s lemma on summation by parts and algorithm of q-hypergeometric summations, we deal with the summation, which can be written as multiplication of a q-hypergeometric term and q-harmonic numbers. This enables us to construct and prove identities on q-harmonic numbers. Several examples are also given.

Keywords

Harmonic Numbers, q-Zeilberger Algorithm, Abel’s Lemma

1. Introduction

Harmonic numbers are a class of famous sequences in combinatorics, number theory and computer science. Many properties could be found in the literature. Recall that for positive integers m and n, the classical generalized harmonic numbers are given by

$H_n^{\left(m\right)}={\displaystyle \underset{k=1}{\overset{n}{\sum }}\frac{1}{k^m}}.$ (1.1)

For convenience, we have $H_0^{\left(m\right)}=0$ . As usually, $H_n^{\left(1\right)}=H_n$ .

In recent years, many identities involving the generalized harmonic numbers have been established [1] - [8] . For example, in [2] , the author reconsidered the following type of sums

$S\left(m\right)={\displaystyle \underset{k=1}{\overset{n}{\sum }}\left(\begin{array}{c}n\\ k\end{array}\right){\left(-1\right)}^k{k}^m{H}_k}.$ (1.2)

With the help of binomial transformations and difference operators, he reproved some well-known identities such as

${\displaystyle \underset{k=1}{\overset{n}{\sum }}\left(\begin{array}{c}n\\ k\end{array}\right){\left(-1\right)}^k{H}_k=-\frac{1}{n}},n>0,$ (1.3)

${\displaystyle \underset{k=1}{\overset{n}{\sum }}\left(\begin{array}{l}n\\ k\end{array}\right){\left(-1\right)}^{k}k{H}_K=\frac{1}{n-1},n>1},$ (1.4)

${\displaystyle \underset{k=1}{\overset{n}{\sum }}\left(\begin{array}{l}n\\ k\end{array}\right){\left(-1\right)}^k{k}^2{H}_k=\frac{-n}{\left(n-1\right)\left(n-2\right)}},n>2.$ (1.5)

q-analogues of combinatorial identities always appear in other branches such as physics [9] . In this paper, we focus on the q-analogues of identities involving harmonic numbers. The well-known two kinds of q-harmonic numbers are defined as follows.

Definition 1. For positive integers n, two kinds of the q-harmonic numbers are defined as

$H_q\left(n\right)={\displaystyle \underset{k=1}{\overset{n}{\sum }}\frac{1}{\left[k\right]}},{\tilde{H}}_q\left(n\right)={\displaystyle \underset{k=1}{\overset{n}{\sum }}\frac{q^k}{\left[k\right]}},$ (1.6)

where $\left[k\right]=\frac{1-q^k}{1-q}$ is the q-integer. We also set $H_q\left(0\right)={\tilde{H}}_q\left(0\right)=0$ .

Note that when $q\to 1$ , $\left[k\right]\to k$ , thus they are q-analogues of classical harmonic numbers.

Definition 2. For non-negative integers $n,k$ and a complex number q with $k\le n$ , $\left|q\right|<1$ , the q-binomial coefficients are defined as

$\left[\begin{array}{c}n\\ k\end{array}\right]=\frac{{\left(q;q\right)}_n}{{\left(q;q\right)}_k{\left(q;q\right)}_{n-k}},$ (1.7)

where ${\left(a;q\right)}_n=\left(1-a\right)\left(1-aq\right)\cdots \left(1-a{q}^{n-1}\right)$ is the well-known q-factorial.

Note that we have ${\left(a;q\right)}_0=1$ and $\left[\begin{array}{c}n\\ k\end{array}\right]=0$ when $k>n$ or $k<0$ .

Mansour et al. [10] discovered the following elegant identity by using partial fraction decomposition

${\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^{n-k}{q}^{\left(\begin{array}{c}k-n\\ 2\end{array}\right)}\left[\begin{array}{c}n+k\\ k\end{array}\right]\left[\begin{array}{c}n\\ k\end{array}\right]H_q\left(k\right)={\displaystyle \underset{k=1}{\overset{n}{\sum }}{q}^{k^2}\frac{\left[2k\right]}{{\left[k\right]}^2}}},n\ge 1,$ (1.8)

which could be known as a q-analogue of the following identity

${\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^{n-k}\left(\begin{array}{c}n+k\\ k\end{array}\right)\left(\begin{array}{c}n\\ k\end{array}\right)H_k=2H_n},n\ge 1.$ (1.9)

We should point out that there exist some typos in this result. The correct version is given in Section 3.

Wei and Gu [11] considered the q-analogues of the following type of sums

$W\left(m\right)={\displaystyle \underset{\text{k}=0}{\overset{n}{\sum }}{\left(\begin{array}{c}n\\ k\end{array}\right)}^m\left\{1+m\left(n-2k\right)H_k\right\}}.$ (1.10)

Note that when $-2\le m\le 6$ and $k\ne 0$ , the “closed form” of $W\left(m\right)$ have been known. By using Watson’s q-Whipple transformation, they discovered the q-analogues of these identities. For example,

$\begin{array}{l}{\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left[\begin{array}{c}n\\ k\end{array}\right]}^{-1}{q}^{-k}\left\{\left[2k-n\right]\left(2k\left(1-q\right)+{\tilde{H}}_q\left(k\right)\right)+\left(1+2k\right)q^{2k-n}-2k\right\}}\\ =\frac{\left[1+n\right]}{q^{n+1}}{\tilde{H}}_q\left(n+1\right),\end{array}$ (1.11)

${\displaystyle \underset{k=0}{\overset{n}{\sum }}\left[\begin{array}{c}n\\ k\end{array}\right]q^{k\left(k-n-1\right)}\left\{-\left[2k-n\right]{\tilde{H}}_q\left(k\right)+\left(1+k\right)q^{2k-n}-k\right\}=1},$ (1.12)

which reduce to

$W\left(-1\right)={\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(\begin{array}{c}n\\ k\end{array}\right)}^{-1}\left\{1-\left(n-2k\right)H_k\right\}}=\left(1+n\right)H_{n+1},$ (1.13)

$W\left(1\right)={\displaystyle \underset{k=0}{\overset{n}{\sum }}\left(\begin{array}{c}n\\ k\end{array}\right)}\left\{1+\left(n-2k\right)H_k\right\}=1.$ (1.14)

Chen, Hou and Jin [12] proposed the Abel-Zeilberger algorithms to prove combinatorial identities on non-hypergeometric terms. In this paper, we observe that summations on q-harmonic numbers fall in the scope of a simplified q-version of this method. This enables us to prove and establish many identities on q-harmonic numbers. For more detail of the q-analogues identities, see [13] .

2. Our Method

The classical Abel’s lemma on summation by parts is as follows.

Lemma 1. For any two sequences $\left\{a_n\right\}$ , $\left\{b_n\right\}$ , we have

${\displaystyle \underset{k=m}{\overset{n-1}{\sum }}\left(a_{k+1}-a_k\right)b_k}={\displaystyle \underset{k=m}{\overset{n-1}{\sum }}{a}_{k+1}\left(b_k-b_{k+1}\right)+a_n{b}_n-a_m{b}_m}.$ (2.1)

Using the difference operator, Abel’s lemma can be rewritten as:

${\displaystyle \underset{k=m}{\overset{n-1}{\sum }}{b}_k\Delta a_k}=-{\displaystyle \underset{k=m}{\overset{n-1}{\sum }}{a}_{k+1}\Delta b_k+a_n{b}_n-a_m{b}_m}.$ (2.2)

Given a q-hypergeometric term $F\left(n,k\right)$ , namely, $\frac{F\left(n+1,k\right)}{F\left(n,k\right)},\frac{F\left(n,k+1\right)}{F\left(n,k\right)}$

are both rational functions in $q^n$ and $q^k$ , the q-Zeilberger algorithm tries to find polynomials $a_0\left(q^n\right),\cdots ,a_d\left(q^n\right)$ and rational functions $R\left(q^n,q^k\right)$ such that

$\begin{array}{l}{a}_0\left(q^n\right)F\left(n,k\right)+a_1\left(q^n\right)F\left(n+1,k\right)+\cdots +a_d\left(q^n\right)F\left(n+d,k\right)\\ =\Delta R\left(q^n,q^k\right)F\left(n,k\right).\end{array}$ (2.3)

Then summing this skew-recurrence relation over k, one can find a recurrence relation of the sum

$S\left(n\right)={\displaystyle \underset{k}{\sum }F}\left(n,k\right).$ (2.4)

For more detail of hypergeometric algorithms, see [14] .

Let us consider the sum $S\left(n\right)={\displaystyle \underset{k}{\sum }F\left(n,k\right)b_k}$ , where $F\left(n,k\right)$ is q-hypergeometric term, b_k is a sequence satisfying $\Delta b_k$ is q-hypergeometric. The method can be described as follows.

Step 1. Applying the q-Zeilberger algorithm to $F\left(n,k\right)$ , we find the skew-recurrence relation

$a_0\left(q^n\right)F\left(n,k\right)+a_1\left(q^n\right)F\left(n+1,k\right)+\cdots +a_d\left(q^n\right)F\left(n+d,k\right)=\Delta G\left(n,k\right),$ (2.5)

where $G\left(n,k\right)=R\left(q^n,q^k\right)F\left(n,k\right)$ .

Step 2. Multiplying both sides of this relation by b_k and summing over k, we have

$a_0\left(q^n\right)S\left(n\right)+a_1\left(q^n\right)S\left(n+1\right)+\cdots +a_d\left(q^n\right)S\left(n+d\right)={\displaystyle \underset{k}{\sum }\left(\Delta G\left(n,k\right)\right)b_k}.$ (2.6)

Step 3. Applying Abel’s lemma to the right-hand side, we will transform it to a hypergeometric sum since $\Delta b_k$ is q-hypergeometric.

$\begin{array}{l}{a}_0\left(q^n\right)S\left(n\right)+a_1\left(q^n\right)S\left(n+1\right)+\cdots +a_d\left(q^n\right)S\left(n+d\right)\\ =-{\displaystyle \underset{k}{\sum }G\left(n,k+1\right)\Delta }{b}_k+W\left(n\right).\end{array}$ (2.7)

Note that $W\left(n\right)=0$ in many cases.

Step 4. Denote by $T\left(n\right)={\displaystyle \underset{k}{\sum }G\left(n,k+1\right)\Delta b_k}$ . If a closed form for $T\left(n\right)$ could be found by using the q-Zeilberger algorithm, we thus obtain a recurrence relation for the original sum $S\left(n\right)$ .

Remark. The Maple package for the q-Zeilberger algorithms can also be found in [14] . Here we use the package APCI, which is written by Hou Q. H. [15] . After loading the package, we input the command qZeil (F, n, k, q, “cert”) and then get the skew-recurrence relation (2.5).

Let us give an example to illustrate this method.

Example 1. For positive integer n, we have

${\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^k{q}^{\left(\begin{array}{c}k+1\\ 2\end{array}\right)-n\left(k+1\right)}\left[\begin{array}{c}n\\ k\end{array}\right]{\tilde{H}}_q\left(k\right)}=-\frac{1}{\left[n\right]},$ (2.8)

${\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^k{q}^{\left(\begin{array}{c}k\\ 2\end{array}\right)}\left[\begin{array}{c}n\\ k\end{array}\right]H_q\left(k\right)}=-\frac{1}{\left[n\right]}.$ (2.9)

The first result appeared in [16] and the second appeared in [17] . Clearly, both of them could be seen as q-analogue of

${\displaystyle \underset{k=1}{\overset{n}{\sum }}\left(\begin{array}{c}n\\ k\end{array}\right){\left(-1\right)}^k{H}_k}=-\frac{1}{n},n>0.$ (2.10)

Proof. For brevity, we only prove the first identity. Denote the left-hand side sum by $S\left(n\right)$ and let

$F\left(n,k\right)={\left(-1\right)}^k{q}^{\left(\begin{array}{c}k+1\\ 2\end{array}\right)-n\left(k+1\right)}\left[\begin{array}{c}n\\ k\end{array}\right].$ (2.11)

By the q-Zeilberger algorithm, we find that

$F\left(n,k\right)=G\left(n,k+1\right)-G\left(n,k\right),$ (2.12)

where $G\left(n,k\right)=-\frac{q^n\left(1-q^k\right)}{\left(1-q^n\right)q^k}F\left(n,k\right)$ . Then with the help of Abel’s lemma, we

obtain the following formula by multiplying both sides of the above equation by ${\tilde{H}}_q\left(k\right)$ and summing over k from 0 to +∞

$S\left(n\right)={\displaystyle \underset{k=0}{\overset{n}{\sum }}\Delta G\left(n,k\right){\tilde{H}}_q\left(k\right)}=-{\displaystyle \underset{k=0}{\overset{+\infty }{\sum }}\frac{G\left(n,k+1\right)}{1-q^{k+1}}\left(1-q\right)q^{k+1}}.$ (2.13)

Denote the new sum by $T\left(n\right)$ and $F_1\left(n,k\right)=\frac{G\left(n,k+1\right)}{1-q^{k+1}}\left(1-q\right)q^{k+1}$ , by the q-Zeilberger algorithm, we further find $F_1\left(n,k\right)=G_1\left(n,k+1\right)-G_1\left(n,k\right)$ , where $G_1\left(n,k\right)=-\frac{q^n\left(1-q^{k+1}\right)}{\left(1-q^n\right)q^{k+1}}{F}_1\left(n,k\right)$ .

By summing over k from 0 to ∞, we obtain

$T\left(n\right)=-G_1\left(n,0\right)=\frac{1-q}{1-q^n}=\frac{1}{\left[n\right]}.$ (2.14)

Thus, we finally have $S\left(n\right)=-T\left(n\right)=-\frac{1}{\left[n\right]}$ .

3. Applications

By using the above method, we can construct or prove many identities on q-harmonic numbers. Here, we only give two examples. The first is the correct version of identity

${\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^{n-k}{q}^{\left(\begin{array}{c}k-n\\ 2\end{array}\right)}\left[\begin{array}{c}n+k\\ k\end{array}\right]\left[\begin{array}{c}n\\ k\end{array}\right]H_q\left(k\right)}={\displaystyle \underset{k=1}{\overset{n}{\sum }}{q}^{k^2}\frac{\left[2k\right]}{{\left[k\right]}^2}},n\ge 1.$ (3.1)

Theorem 1. For positive integer n, we have

${\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^{n-k}{q}^{\left(\begin{array}{c}k-n\\ 2\end{array}\right)}\left[\begin{array}{c}n+k\\ k\end{array}\right]\left[\begin{array}{c}n\\ k\end{array}\right]H_q\left(k\right)}={\displaystyle \underset{k=1}{\overset{n}{\sum }}\frac{\left[2k\right]}{{\left[k\right]}^2}}.$ (3.2)

Proof. Denote the left-hand side sum by $S\left(n\right)$ and let $F\left(n,k\right)={\left(-1\right)}^{n-k}{q}^{\left(\begin{array}{c}k-n\\ 2\end{array}\right)}\left[\begin{array}{c}n+k\\ k\end{array}\right]\left[\begin{array}{c}n\\ k\end{array}\right]$ . By the q-Zeilberger algorithm, we find that

$F\left(n,k\right)-F\left(n+1,k\right)=G\left(n,k+1\right)-G\left(n,k\right),$ (3.3)

where $G\left(n,k\right)=\frac{q^{n+1}{\left(1-q^k\right)}^2\left(1+q^{n+1}\right)}{\left(1-q^{n+1}\right)\left(q^{\text{n+}1}-q^k\right)}F\left(n,k\right)$ . Then with the help of Abel’s

lemma, we obtain the following formula by multiplying both sides of Equation (3.3) by $H_q\left(k\right)$ and summing over k from 0 to +∞,

$S\left(n\right)-S\left(n+1\right)={\displaystyle \underset{k=0}{\overset{n}{\sum }}\Delta G\left(n,k\right)H_q\left(k\right)}=-{\displaystyle \underset{k=0}{\overset{+\infty }{\sum }}\frac{G\left(n,k+1\right)}{1-q^{k+1}}}\left(1-q\right).$ (3.4)

Denote the new sum by $T\left(n\right)$ and $F_1\left(n,k\right)=\frac{G\left(n,k+1\right)}{1-q^{k+1}}\left(1-q\right)$ , by the q-Zeilberger algorithm, we further find

$\begin{array}{l}\left(1+q^{n+2}\right)\left(q^{n+1}-1\right)F_1\left(n,k\right)-\left(1+q^{n+1}\right)\left(q^{n+2}-1\right)F_1\left(n+1,k\right)\\ =G_1\left(n,k+1\right)-G_1\left(n,k\right),\end{array}$ (3.5)

where $G_1\left(n,k\right)=-\frac{\left(1+q^{n+2}\right)\left(q^k-1\right)\left(q^{k+1}-1\right)q^{n+1}\left(q^{2n+3}-1\right)}{\left(q^{n+1}-q^k\right)\left(q^{n+2}-1\right)}{F}_1\left(n,k\right)$ . By summing over k from 0 to ∞, we obtain

$\left(1+q^{n+2}\right)\left(q^{n+1}-1\right)T\left(n\right)-\left(1+q^{n+1}\right)\left(q^{n+2}-1\right)T\left(n+1\right)=0,$ (3.6)

Taking the initial value $T\left(1\right)=G\left(1,0\right)+G\left(1,1\right)=\frac{1+q^2}{1+q}$ into account, we obtain

$S\left(n\right)-S\left(n+1\right)=-T\left(n\right)=-\frac{1+q^{n+1}}{1-q^{n+1}}\left(1-q\right).$ (3.7)

Furthermore, $S\left(1\right)=1+q$ , we thus have

$S\left(n\right)=S\left(1\right)+{\displaystyle \underset{k=2}{\overset{n}{\sum }}\frac{1-q^{2k}}{{\left(1-q^k\right)}^2}\left(1-q\right)}={\displaystyle \underset{k=1}{\overset{n}{\sum }}\frac{\left[2k\right]}{{\left[k\right]}^2}}.$ (3.8)

The second example is the following q-analogues of (1.4) and (1.5).

Theorem 2. For positive integer n, we have

${\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^k{q}^{\left(\begin{array}{c}k+1\\ 2\end{array}\right)-n\left(k+1\right)}\left[k\right]\left[\begin{array}{c}n\\ k\end{array}\right]{\tilde{H}}_q\left(k\right)}=\frac{1}{q\left[n-1\right]},\left(n>1\right),$ (3.9)

${\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^k{q}^{\left(\begin{array}{c}k+1\\ 2\end{array}\right)-n\left(k+1\right)}{\left[k\right]}^2\left[\begin{array}{c}n\\ k\end{array}\right]{\tilde{H}}_q\left(k\right)}=-\frac{\left[n\right]}{q^2\left[n-1\right]\left[n-2\right]},\left(n>2\right).$ (3.10)

Proof. We only prove the first identity. Denote the summation by $S\left(n\right)$ and

$F\left(n,k\right)={\left(-1\right)}^k{q}^{\left(\begin{array}{c}k+1\\ 2\end{array}\right)-n\left(k+1\right)}\left[k\right]\left[\begin{array}{c}n\\ k\end{array}\right].$ (3.11)

Then, by the q-Zeilberger algorithm, we know

$F\left(n,k\right)=G\left(n,k+1\right)-G\left(n,k\right),$ (3.12)

where

$G\left(n,k\right)=-\frac{q^n\left(q-q^k\right)}{\left(q-q^n\right)q^k}F\left(n,k\right).$ (3.13)

Noting that the summation range of the series is equivalent to summing from 0 to ∞, we multiply both sides of Equation (3.13) by ${\tilde{H}}_q\left(k\right)$ and sum over k from 0 to ∞. Using Abel’s lemma, we get (with the reminder term being 0).

$\begin{array}{c}S\left(n\right)={\displaystyle \underset{k=0}{\overset{n}{\sum }}{\left(-1\right)}^k{q}^{\left(\begin{array}{c}k+1\\ 2\end{array}\right)-n\left(k+1\right)}\left[k\right]\left[\begin{array}{c}n\\ k\end{array}\right]{\tilde{H}}_q\left(k\right)}\\ \text{=}{\displaystyle \underset{k=0}{\overset{n}{\sum }}\Delta G\left(n,k\right){\tilde{H}}_q\left(k\right)}\\ \text{=}-{\displaystyle \underset{k=0}{\overset{+\infty }{\sum }}\frac{G\left(n,k+1\right)}{1-q^{k+1}}\left(1-q\right)q^{k+1}}\\ \text{=}-T\left(n\right).\end{array}$ (3.14)

Let

$F_1\left(n,k\right)=\frac{G\left(n,k+1\right)}{1-q^{k+1}}\left(1-q\right)q^{k+1}.$ (3.15)

Furthermore, by the q-Zeilberger algorithm, we know

$F_1\left(n,k\right)=G_1\left(n,k+1\right)-G_1\left(n,k\right),$ (3.16)

where

$G_1\left(n,k\right)=\frac{\left(q^{k+n+1}+q^2-q^{k+1}-2q^{n+1}+q^n\right)q^n\left(q^{k+1}-1\right)}{\left(-1+q^k\right)\left(q-q^n\right)q^k\left(q^n-1\right)q^2}{F}_1\left(n,k\right).$ (3.17)

Taking the sum over k from 0 to ∞ for Equation (3.17), we obtain:

$T\left(n\right)=-G_1\left(n,0\right)=-\frac{1}{q\left[n-1\right]}.$ (3.18)

Thus, we have:

$S\left(n\right)=-T\left(n\right)=\frac{1}{q\left[n-1\right]}\text{,}\left(n>1\right).$ (3.19)

4. Conclusion

Sometimes, finding q-analogues for given combinatorial identities may be challenging. In this paper, we used the ideal of the Abel-Zeilberger algorithm to construct and prove several q-analogues of identities involving harmonic numbers. In particular, we obtained the q-analogues identities for the classical kind of

${\displaystyle \underset{k=1}{\overset{n}{\sum }}\left(\begin{array}{c}n\\ k\end{array}\right){\left(-1\right)}^k{k}^m{H}_k}$ when $m=0,1,2$ . For $m\ge 3$ , we can also construct the corresponding identities. In fact, we have verified almost all identities on q-harmonic

numbers appeared in the references. We also point out that some examples of the method have been included in Zheng huanhuan’s master thesis.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Jin, H.T. and Du, D.K. (2015) Abel’s Lemma and Identities on Harmonic Numbers. Integers, 15, A22.
[2]	Boyadzhiev, K.N. (2014) Binomial Transform and the Backward Difference. Mathematics, 13, 43-63.
[3]	Chen, K.W. and Chen, Y.H. (2020) Infinite Series Containing Generalized Harmonic Functions. Notes Number Theory and Discrete Mathematics, 26, 85-104. https://doi.org/10.7546/nntdm.2020.26.2.85-104
[4]	Chu, W.C. and Donno, L.D. (2005) Hypergeometric Series and Harmonic Number Identities. Advances in Applied Mathematics, 34, 123-137. https://doi.org/10.1016/j.aam.2004.05.003
[5]	Frontczak, R. (2021) Binomial Sums with Skew-Harmonic Numbers. Palestine Journal of Mathematics, 10, 756-763.
[6]	Guo, D.W. (2022) Some Combinatorial Identities Concerning Harmonic Numbers and Binomial Coefficients. Discrete Mathematics Letters, 8, 41-48. https://doi.org/10.47443/dml.2021.0090
[7]	Kargın, L. and Can M. (2020) Harmonic Number Identities via Polynomials with r-Lah Coefficients. Comptes Rendus Mathématique, 358, 535-550. https://doi.org/10.5802/crmath.53
[8]	Sun, Z.W. (2023) Series with Summands Involving Higher Harmonic Numbers. arXiv:2210.07238v8.
[9]	Andrews, G.E. (1986) q-Series: Their Development and Application in Analysis, Number Theory, Combinatorics, Physics, and Computer Algebra. CBMS Regional Conference Series in Mathematics, Vol. 66. https://doi.org/10.1090/cbms/066
[10]	Mansour, T., Mansour, M. and Song, C.W. (2012) q-Analogs of Identities Involving Harmonic Numbers and Binomial Coefficient. Applications and Applied Mathematics, 7, 22-36.
[11]	Wei, C.A. and Gu, Q.P. (2010) q-Generalizations of a Family of Harmonic Number Identities. Advances in Applied Mathematics, 45, 24-27. https://doi.org/10.1016/j.aam.2009.11.007
[12]	Chen, W.Y.C., Hou, Q.H. and Jin, H.T. (2011) The Abel-Zeilberger Algorithm. The Electronic Journal of Combinatorics, 18, 17. https://doi.org/10.37236/2013
[13]	Zriaa, S. and Mouçouf, M. (2023) Algebraic Identities on q-Harmonic Numbers and q-Binomial Coefficient. arXiv.2301.13747v1.
[14]	Petkovsěk, M., Wilf, H.S. and Zeilberger, D. (1996) A = B. A.K. Peters Ltd., New York. https://doi.org/10.1201/9781439864500
[15]	Chen, W.Y.C., Hou, Q.H. and Mu, Y.P. (2012) The Extended Zeilberger Algorithm with Parameters. Journal of Symbolic Computation, 47, 643-654. https://doi.org/10.1016/j.jsc.2011.12.024
[16]	Xu, A.M. (2014) On a General q-Identity. The Electronic Journal of Combinatorics, 21, 2.28. https://doi.org/10.37236/3962
[17]	Chu, W.C. and Wang, X.Y. (2019) Harmonic Number Sums and q-Analogues. International Journal of Computer Mathematics: Computer Systems Theory, 4, 48-56. https://doi.org/10.1080/23799927.2019.1570974