A mathematical model of tuberculosis (TB) transmission with children and adults groups: A fractional model

Fatmawati; Muhammad Altaf Khan; Ebenezer Bonyah; Zakia Hammouch; Endrik Mifta Shaiful; Fatmawati; Muhammad Altaf Khan; Ebenezer Bonyah; Zakia Hammouch; Endrik Mifta Shaiful

doi:10.3934/math.2020181

AIMS Mathematics

2020, Volume 5, Issue 4: 2813-2842. doi: 10.3934/math.2020181

Previous Article Next Article

Research article Special Issues

A mathematical model of tuberculosis (TB) transmission with children and adults groups: A fractional model

1.
Department of Mathematics, Faculty of Science and Technology, Universitas Airlangga Surabaya 60115, Indonesia
2.
Faculty of Natural and Agricultural Sciences, University of the Free State, South Africa
3.
Department of Mathematics Education, University of Education Winneba, Kumasi Campus, Kumasi, Ghana
4.
Departement de Mathématiques, Faculté des Sciences et Techniques Errachidia, Université Moulay Ismail, Morocco

Received: 01 December 2019 Accepted: 19 February 2020 Published: 17 March 2020
MSC : 34A34

We consider a novel fractional model to investigates the (tuberculosis) TB model dynamics with two age groups of human, that is, the children and the adults. First, we formulate the model and present the basic results associated to the model. Then, using the fractional operators, Caputo and the Atangana-Baleanu and obtain a generalized model. Further, we give a novel numerical approach for the solution of the fractional model and obtain their approximate solution. We show graphical results with various values of the fractional order. A comparison of the two operators are shown graphically. The results obtained through Atangana-Baleanu operator is flexible than that of Caputo derivative. The infection in tuberculosis (TB) infected people decreases fast when decreasing the fractional order.

Keywords:

Citation: Fatmawati, Muhammad Altaf Khan, Ebenezer Bonyah, Zakia Hammouch, Endrik Mifta Shaiful. A mathematical model of tuberculosis (TB) transmission with children and adults groups: A fractional model[J]. AIMS Mathematics, 2020, 5(4): 2813-2842. doi: 10.3934/math.2020181

Related Papers:

[1]	Ridha Dida, Hamid Boulares, Bahaaeldin Abdalla, Manar A. Alqudah, Thabet Abdeljawad . On positive solutions of fractional pantograph equations within function-dependent kernel Caputo derivatives. AIMS Mathematics, 2023, 8(10): 23032-23045. doi: 10.3934/math.20231172
[2]	Saeed M. Ali, Mohammed S. Abdo, Bhausaheb Sontakke, Kamal Shah, Thabet Abdeljawad . New results on a coupled system for second-order pantograph equations with $\mathcal{ABC}$ fractional derivatives. AIMS Mathematics, 2022, 7(10): 19520-19538. doi: 10.3934/math.20221071
[3]	Abdelkader Moumen, Hamid Boulares, Tariq Alraqad, Hicham Saber, Ekram E. Ali . Newly existence of solutions for pantograph a semipositone in $\Psi$ -Caputo sense. AIMS Mathematics, 2023, 8(6): 12830-12840. doi: 10.3934/math.2023646
[4]	Ayub Samadi, Chaiyod Kamthorncharoen, Sotiris K. Ntouyas, Jessada Tariboon . Mixed Erdélyi-Kober and Caputo fractional differential equations with nonlocal non-separated boundary conditions. AIMS Mathematics, 2024, 9(11): 32904-32920. doi: 10.3934/math.20241574
[5]	Hui Huang, Kaihong Zhao, Xiuduo Liu . On solvability of BVP for a coupled Hadamard fractional systems involving fractional derivative impulses. AIMS Mathematics, 2022, 7(10): 19221-19236. doi: 10.3934/math.20221055
[6]	Cuiying Li, Rui Wu, Ranzhuo Ma . Existence of solutions for Caputo fractional iterative equations under several boundary value conditions. AIMS Mathematics, 2023, 8(1): 317-339. doi: 10.3934/math.2023015
[7]	Mohamed Houas, Kirti Kaushik, Anoop Kumar, Aziz Khan, Thabet Abdeljawad . Existence and stability results of pantograph equation with three sequential fractional derivatives. AIMS Mathematics, 2023, 8(3): 5216-5232. doi: 10.3934/math.2023262
[8]	Ahmed M. A. El-Sayed, Wagdy G. El-Sayed, Kheria M. O. Msaik, Hanaa R. Ebead . Riemann-Liouville fractional-order pantograph differential equation constrained by nonlocal and weighted pantograph integral equations. AIMS Mathematics, 2025, 10(3): 4970-4991. doi: 10.3934/math.2025228
[9]	Isra Al-Shbeil, Abdelkader Benali, Houari Bouzid, Najla Aloraini . Existence of solutions for multi-point nonlinear differential system equations of fractional orders with integral boundary conditions. AIMS Mathematics, 2022, 7(10): 18142-18157. doi: 10.3934/math.2022998
[10]	Yujun Cui, Chunyu Liang, Yumei Zou . Existence and uniqueness of solutions for a class of fractional differential equation with lower-order derivative dependence. AIMS Mathematics, 2025, 10(2): 3797-3818. doi: 10.3934/math.2025176

Abstract

1. Introduction

The key to solving the general quadratic congruence equation is to solve the equation of the form $x^2\equiv a\mod p$ , where $a$ and $p$ are integers, $p > 0$ and $p$ is not divisible by $a$ . For relatively large $p$ , it is impractical to use the Euler criterion to distinguish whether the integer $a$ with $(a, p) = 1$ is quadratic residue of modulo $p$ . In order to study this issue, Legendre has proposed a new tool-Legendre's symbol.

Let $p$ be an odd prime, the quadratic character modulo $p$ is called the Legendre's symbol, which is defined as follows:

$\left(\frac{a}{p}\right)= \begin{cases}1, & \text { if } a \text { is a quadratic residue modulo } p; \\ -1, & \text { if } a \text { is a quadratic non-residue modulo } p ;\\ 0, & \text { if } p \mid a.\end{cases}$

The Legendre's symbol makes it easy for us to calculate the level of quadratic residues. The basic properties of Legendre's symbol can be found in any book on elementary number theory, such as ^[1,2,3].

The properties of Legendre's symbol and quadratic residues play an important role in number theory. Many scholars have studied them and achieved some important results. For examples, see the ^{[4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21]}.

One of the most representative properties of the Legendre's symbol is the quadratic reciprocal law:

Let $p$ and $q$ be two distinct odd primes. Then, (see Theorem 9.8 in ^[1] or Theorems 4–6 in ^[3])

$\left(\frac{p}{q}\right)\cdot\left(\frac{q}{p}\right) = (-1)^{\frac{(p-1)(q-1)}{4}}.$

For any odd prime $p$ with $p\equiv1\; \rm mod\; 4$ there exist two non-zero integers $\alpha\left(p\right)$ and $\beta\left(p\right)$ such that

$\begin{align*} p = \alpha^2\left(p\right)+\beta^2\left(p\right). \end{align*}$

(1)

In fact, the integers $\alpha\left(p\right)$ and $\beta\left(p\right)$ in the (1) can be expressed in terms of Legendre's symbol modulo $p$ (see Theorems 4–11 in ^[3])

$\alpha\left(p\right) = \frac{1}{2}\sum\limits_{a = 1}^{p-1}\left(\frac{a^3+a}{p}\right)\; \; \; \rm and\; \; \; \it \beta\left(\it p\right) = \frac{\rm1}{\rm 2}\sum\limits_{a = \rm1}^{p-\rm1}\left(\frac {\it a^{\rm3}+\it ra}{\it p}\right),$

where $r$ is any integer, and $\left(r, p\right) = 1$ , $\left(\frac{r}{p}\right) = -1$ , $\left(\frac{\ast}{p}\right) = \chi_2$ denote the Legendre's symbol modulo $p$ .

Noting that Legendre's symbol is a special kind of character. For research on character, Han ^[7] studied the sum of a special character $\chi\left(ma+\bar{a}\right)$ , for any integer $m$ with $\left(m, p\right) = 1$ , then

$\left|\sum\limits_{a = 1}^{p-1}\chi\left(ma+\bar{a}\right)\right|^2 = 2p+\left(\frac{m}{p}\right)\sum\limits_{a = 1}^{p-1}\chi(a)\sum\limits_{b = 1}^{p-1}\left(\frac{b(b-1)(a^2b-1)}{p}\right),$

which is a special case of a general polynomial character sums $\sum_{a = N+1}^{N+M}\chi\left(f\left(a\right)\right)$ , where $M$ and $N$ are any positive integers, and $f\left(x\right)$ is a polynomial.

In ^[8], Du and Li introduced a special character sums $C\left(\chi, m, n, c;p\right)$ in the following form:

$C\left(\chi, m, n, c;p\right) = \sum\limits_{a = 0}^{p-1}\sum\limits_{b = 0}^{p-1}\chi\left(a^2+na-b^2-nb+c\right)\cdot e\left(\frac{mb^2-ma^2}{p}\right),$

and studied the asymptotic properties of it. They obtained

$\begin{equation*} \label{test1} \begin{split} \sum\limits_{c = 1}^{p-1}\left|C\left(\chi, m, n, c;p\right)\right|^{2k} = \left\{ \begin{array}{l} &p^{2k+1}+\frac{k^2-3k-2}{2}\cdot p^{2k}+O\left(p^{2k-1}\right), \; \; \; \; \rm if\; \chi\ \rm is\ \rm the\ \rm Legendre\ \rm symbol\ \rm {modulo} \ \it p \rm; \\ &p^{2k+1}+\frac{k^2-3k-2}{2}\cdot p^{2k}+O\left(p^{2k-1/2}\right), \; \; \rm if\; \chi\ \rm is\ \rm a\ \rm complex\ \rm character\ \rm {modulo} \ \it p.\\ \end{array} \right. \end{split} \end{equation*}$

Recently, Yuan and Zhang ^[12] researched the question about the estimation of the mean value of high-powers for a special character sum modulo a prime, let $p$ be an odd prime with $p\equiv1\; \rm mod\; 6$ , then for any integer $k\geq0$ , they have the identity

$S_k\left(p\right) = \frac{1}{3}\cdot\left[d^k+\left(\frac{-d+9b}{2}\right)^k+\left(\frac{-d-9b}{2}\right)^k\right],$

where

$S_k\left(p\right) = \frac{1}{p-1}\sum\limits_{r = 1}^{p-1}A^k\left(r\right),$

$A\left(r\right) = 1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+r\bar{a}}{p}\right),$

and for any integer $r$ with $(r, p) = 1$ .

More relevant research on special character sums will not be repeated. Inspired by these papers, we have the question: If we replace the special character sums with Legendre's symbol, can we get good results on $p\equiv1\; \rm mod\; 4$ ?

We will convert $\beta\left(p\right)$ to another form based on the properties of complete residues

$\beta\left(p\right) = \frac{1}{2}\sum\limits_{a = 1}^{p-1}\left(\frac{a+n\bar{a}}{p}\right),$

where $\bar{a}$ is the inverse of $a$ modulo $p$ . That is, $\bar{a}$ satisfy the equation $x\cdot a\; \equiv1\mod p$ for any integer $a$ with $(a, p) = 1$ .

For any integer $k\geq0$ , $G\left(n\right)$ and $K_k(p)$ are defined as follows:

$G\left(n\right) = 1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right)\; \; \; \rm {and} \; \; \; \it K_{\it k}\left(\it p\right) = \frac{\rm 1}{\it p-\rm 1}\sum\limits_{\it n\rm = 1}^{\it p-\rm 1}\it G^{\it k}\left(\it n\right).$

In this paper, we will use the analytic methods and properties of the classical Gauss sums and Dirichlet character sums to study the computational problem of $K_k\left(p\right)$ for any positive integer $k$ , and give a linear recurrence formulas for $K_k\left(p\right)$ . That is, we will prove the following result.

Theorem 1. Let $p$ be an odd prime with $p\equiv1\; \rm mod\; 4$ , then we have

$K_k\left(p\right) = \left(4p+2\right)\cdot K_{k-2}\left(p\right)-8\left(2\alpha^2-p\right)\cdot K_{k-3}\left(p\right)+\left(16\alpha^4-16p\alpha^2+4p-1\right)\cdot K_{k-4}\left(p\right),$

for all integer $k\geq4$ with

$K_0\left(p\right) = 1, \; K_1\left(p\right) = 0, \; K_2\left(p\right) = 2p+1, \; K_3\left(p\right) = -3(4\alpha^2-2p),$

where

$\alpha = \alpha\left(p\right) = \sum\limits_{a = 1}^{\frac{p-1}{2}}\left(\frac{a+\bar{a}}{p}\right).$

Applying the properties of the linear recurrence sequence, we may immediately deduce the following corollaries.

Corollary 1. Let $p$ be an odd prime with $p\equiv1\; \rm mod\; 4$ . Then we have

$\frac{1}{p-1}\sum\limits_{n = 1}^{p-1}\frac{1}{1+\sum _{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right)} = \frac{16\alpha^2p-28\alpha^2-8p^2+14p} {16\alpha^4-16\alpha^2p+4p-1}.$

Corollary 2. Let $p$ be an odd prime with $p\equiv1\; \rm mod\; 4$ . Then we have

$\begin{align*} \frac{1}{p-1}\sum\limits_{n = 1}^{p-1}\sum\limits_{m = 0}^{p-1}\left(1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right)\right)\cdot e\left(\frac{nm^2}{p}\right) = -\sqrt{p}. \end{align*}$

Corollary 3. Let $p$ be an odd prime with $p\equiv1\; \rm mod\; 4$ . Then we have

$\frac{1}{p-1}\sum\limits_{n = 1}^{p-1}\sum\limits_{m = 0}^{p-1}\left[1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right)\right]^2\cdot e\left(\frac{nm^2}{p}\right) = \left(4\alpha^2-2p\right)\cdot\sqrt{p}.$

Corollary 4. Let $p$ be an odd prime with $p\equiv1\; \rm mod\; 8$ . Then we have

$\begin{align*} \sum\limits_{n = 1}^{p-1}\left(1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right)\right)\cdot \sum\limits_{m = 0}^{p-1}e\left(\frac{nm^4}{p}\right) = \sqrt{p}\left(-1+B(1)\right)-p, \end{align*}$

where

$B(1) = \sum\limits_{m = 0}^{p-1}e\left(\frac{m^4}{p}\right).$

If we consider such a sequence $F_k\left(p\right)$ as follows: Let $p$ be a prime with $p\equiv1\; \rm mod\; 8$ , $\chi_4$ be any fourth-order character modulo $p$ . For any integer $k\geq0$ , we define the $F_k\left(p\right)$ as

$F_k\left(p\right) = \sum\limits_{n = 1}^{p-1}\frac{1}{G^{k}(n)},$

we have

$\begin{align*} F_{k}\left(p\right) = &\frac{1}{16\alpha^4-16\alpha^2p+4p-1}F_{k-4}\left(p\right)-\frac{(4p+2)}{16\alpha^4-16\alpha^2p+4p-1}F_{k-2}\left(p\right)\\ &+\frac{4(4\alpha^2-2p)}{16\alpha^4-16\alpha^2p+4p-1}F_{k-1}\left(p\right). \end{align*}$

2. Some lemmas

Lemma 1. Let $p$ be an odd prime with $p\equiv1\; \rm mod\; 4$ . Then for any fourth-order character $\chi_4 \; \rm{mod} \; \it p$ , we have the identity

$\tau^2(\chi_4)+\tau^2(\bar{\chi_4}) = 2\sqrt{p}\cdot\alpha,$

where

$\tau(\chi_4) = \sum\limits_{a = 1}^{p-1}\chi_4(a)e\left(\frac{a}{p}\right)$

denotes the classical Gauss sums, $e(y) = e^{2\pi iy}, i^2 = -1$ , and $\alpha$ is the same as in the Theorem 1.

Proof. See Lemma 2.2 in ^[9].

Lemma 2. Let $p$ be an odd prime. Then for any non-principal character $\psi$ modulo $p$ , we have the identity

$\tau(\psi^2) = \frac{\psi^2(2)}{\tau(\chi_2)}\cdot\tau(\psi)\cdot\tau(\psi\chi_2),$

where $\chi_2 = \left(\frac{\ast}{p}\right)$ denotes the Legendre's symbol modulo $p$ .

Proof. See Lemma 2 in ^[12].

Lemma 3. Let $p$ be a prime with $p\equiv 1\; \rm mod\; 4$ , then for any integer $n$ with $\left(n, p\right) = 1$ and fourth-order character $\chi_{4}\; \rm mod \; \it p$ , we have the identity

$\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right) = -1-\chi_2(n)+\frac{1}{\sqrt{p}}\cdot(\chi_{4}(n)\cdot\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\cdot\tau^2(\chi_{4})).$

Proof. For any integer $a$ with $(a, p) = 1$ , we have the identity

$1+\chi_{4}(a)+\chi_2(a)+\bar{\chi_{4}}(a) = 4,$

if $a$ satisfies $a\equiv b^4\; \rm mod \; \it p$ for some integer $b$ with $(b, p) = 1$ and

$1+\chi_{4}(a)+\chi_2(a)+\bar{\chi_{4}}(a) = 0,$

otherwise. So from these and the properties of Gauss sums we have

$\begin{align*} \mathop\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right) = &\sum\limits_{a = 1}^{p-1}\left(\frac{a^2}{p}\right)\left(\frac{a^4+n}{p}\right)\nonumber\\ = &\sum\limits_{a = 1}^{p-1}\chi_2\left(a^4\right)\chi_2\left(a^4+n\right)\nonumber\\ = &\sum\limits_{a = 1}^{p-1}(1+\chi_{4}(a)+\chi_2(a)+\bar{\chi_{4}}(a))\cdot\chi_2(a)\cdot\chi_2(a+n)\nonumber\\ = &\sum\limits_{a = 1}^{p-1}(1+\chi_{4}(na)+\chi_2(na)+\bar{\chi_{4}}(na))\cdot\chi_2(na)\cdot\chi_2(na+n)\nonumber\\ = &\sum\limits_{a = 1}^{p-1}\chi_2(a)\chi_2(a+1)+\sum\limits_{a = 1}^{p-1}\chi_{4}(na)\chi_2(a)\chi_2(a+1) \end{align*}$

(2)

$\begin{align*} &+\sum\limits_{a = 1}^{p-1}\chi_2(na)\chi_2(a)\chi_2(a+1)+\sum\limits_{a = 1}^{p-1}\bar{\chi_{4}}(na)\chi_2(a)\chi_2(a+1)\nonumber\\ = &\sum\limits_{a = 1}^{p-1}\chi_2(1+\bar{a})+\sum\limits_{a = 1}^{p-1}\chi_{4}(na)\chi_2(a)\chi_2(a+1)\nonumber\\ &+\sum\limits_{a = 1}^{p-1}\chi_2(n)\chi_2(a+1)+\sum\limits_{a = 1}^{p-1}\bar{\chi_{4}}(na)\chi_2(a)\chi_2(a+1)\nonumber. \end{align*}$

Noting that for any non-principal character $\chi$ ,

$\sum\limits_{a = 1}^{p-1}\chi(a) = 0$

and

$\sum\limits_{a = 1}^{p-1}\chi(a)\chi(a+1) = \frac{1}{\tau(\bar{\chi})}\sum\limits_{b = 1}^{p-1}\sum\limits_{a = 1}^{p-1}\bar{\chi}(b)\chi(a)e\left(\frac{b(a+1)}{p}\right).$

Then we have

$\sum\limits_{a = 1}^{p-1}\chi_2(1+\bar{a}) = -1, \; \; \; \sum\limits_{a = 1}^{p-1}\chi_2(a+1) = -1,$

$\begin{align*} \sum\limits_{a = 1}^{p-1}\chi_{4}(a)\chi_2(a)\chi_2(a+1)& = \frac{1}{\tau(\chi_2)}\sum\limits_{b = 1}^{p-1}\sum\limits_{a = 1}^{p-1}\chi_2(b)\chi_{4}(a)\chi_2(a)e\left(\frac{b(a+1)}{p}\right)\nonumber\\ & = \frac{1}{\tau(\chi_2)}\sum\limits_{b = 1}^{p-1}\bar{\chi_{4}}(b)e\left(\frac{b}{p}\right)\sum\limits_{a = 1}^{p-1}\chi_{4}(ab)\chi_2(ab)e\left(\frac{ab}{p}\right) \end{align*}$

(3)

$\begin{align*} & = \frac{1}{\tau(\chi_2)}\cdot\tau(\bar{\chi_{4}})\cdot\tau(\chi_{4}\chi_2)\nonumber. \end{align*}$

For any non-principal character $\psi$ , from Lemma 2 we have

$\begin{align*} \tau(\psi^2) = \frac{\psi^2(2)}{\tau(\chi_2)}\cdot\tau(\psi)\cdot\tau(\psi\chi_2). \end{align*}$

(4)

Taking $\psi = \chi_{4}$ , note that

$\tau(\chi_2) = \sqrt{p}, \ \ \tau(\chi_{4})\cdot\tau(\bar{\chi_{4}}) = \chi_{4}(-1)\cdot p,$

from (3) and (4), we have

$\begin{align*} \sum\limits_{a = 1}^{p-1}\chi_{4}(a)\chi_2(a)\chi_2(a+1)& = \frac{\bar{\chi_{4}}^2(2)\cdot\tau(\chi_{4}^2)\cdot\tau(\chi_2)\cdot\tau(\bar{\chi_{4}})}{\tau(\chi_2)\cdot\tau(\chi_{4})}\nonumber\\ & = \frac{\chi_2(2)\cdot\tau(\chi_2)\cdot\tau^2(\bar{\chi_{4}})}{\tau(\chi_{4})\cdot\tau(\bar{\chi_{4}})}\nonumber\\ & = \frac{\chi_2(2)\cdot\sqrt{p}\cdot\tau^2(\bar{\chi_{4}})}{\chi_4(-1)\cdot p} \end{align*}$

(5)

$\begin{align*} & = \frac{\chi_2(2)\cdot\tau^2(\bar{\chi_{4}})}{\chi_4(-1)\cdot\sqrt{p}}.\nonumber \end{align*}$

Similarly, we also have

$\begin{align*} \sum\limits_{a = 1}^{p-1}\bar{\chi_{4}}(a)\chi_2(a)\chi_2(a+1) = \frac{\chi_2(2)\cdot\tau^2(\chi_{4})}{\chi_4(-1)\cdot\sqrt{p}}. \end{align*}$

(6)

Consider the quadratic character modulo $p$ , we have

$\begin{align*} \left(\frac{2}{p}\right) = \chi_2(2) = \begin{cases} \; \; 1, &\rm{ if\; \ p\equiv\pm1\; mod\; 8};\\ -1, &\rm{ if\; \ p\equiv\pm3\; mod\; 8}.\\ \end{cases} \end{align*}$

(7)

And when $p\equiv1\; \rm mod\; 8$ , we have $\chi_{4}(-1) = 1$ ; when $p\equiv5\; \rm mod\; 8$ , we have $\chi_{4}(-1) = -1.$ Combining (2) and (5)–(7) we can deduce that

$\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right) = -1-\chi_2(n)+\frac{1}{\sqrt{ p}}\cdot\left(\chi_{4}(n)\cdot\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\cdot\tau^2(\chi_{4})\right).$

This prove Lemma 3.

Lemma 4. Let $p$ be an odd prime with $p\equiv1\; \rm mod\; 4$ . Then for any integer $k\geq4$ and $n$ with $(n, p) = 1$ , we have the fourth-order linear recurrence formula

$G^k(n) = (4p+2)\cdot G^{k-2}(n)+8(p-2\alpha^2)\cdot G^{k-3}(n)+[(4\alpha^2-2p)^2-(2p-1)^2]\cdot G^{k-4}(n),$

where

$\alpha = \alpha(p) = \frac{1}{2}\sum\limits_{a = 1}^{p-1}\left(\frac{a^3+a}{p}\right) = \sum\limits_{a = 1}^{\frac{p-1}{2}}\left(\frac{a+\bar{a}}{p}\right),$

$\left(\frac{\ast} {p}\right) = \chi_2$ denotes the Legendre's symbol.

Proof. For $p\equiv1\; \rm mod\; 4$ , any integer $n$ with $\left(n, p\right) = 1$ , and fourth-order character $\chi_{4}$ modulo $p$ , we have the identity

$\chi_{4}^4(n) = \bar{\chi_{4}}^4(n) = \chi_0(n), \ \ \chi_{4}^2(n) = \chi_2(n),$

where $\chi_0$ denotes the principal character modulo $p$ .

According to Lemma 3,

$\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right) = -1-\chi_2(n)+\frac{1}{\sqrt{ p}}\cdot\left(\chi_{4}(n)\cdot\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\cdot\tau^2(\chi_{4})\right),$

$G(n) = 1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right).\quad\quad\quad\quad\quad\quad\quad\quad\ \$

We have

$\begin{align*} G(n) = -\chi_2(n)+\frac{1}{\sqrt{p}}\cdot\left(\chi_{4}(n)\cdot\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\cdot\tau^2(\chi_{4})\right), \end{align*}$

(8)

$\begin{align*} G^2(n) = &[-\chi_2(n)+\frac{1}{\sqrt{p}}\cdot(\chi_{4}(n)\cdot\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\cdot\tau^2(\chi_{4}))]^2\\ = &1-2\chi_2(n)\cdot\frac{1}{\sqrt{p}}\cdot\left(\chi_{4}(n)\cdot\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\cdot\tau^2(\chi_{4})\right)\\ &+\frac{1}{ p}\cdot\left(\chi_2(n)\cdot\tau^4(\bar{\chi_{4}})+\chi_2(n)\cdot\tau^4(\chi_{4})+2p^2\right)\\ = &1-2\chi_2(n)\cdot\frac{1}{\sqrt{p}}\cdot\left(\chi_{4}(n)\cdot\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\cdot\tau^2(\chi_{4})\right)\\ &+\frac{1}{p}\cdot\left(\chi_2(n)\cdot(\tau^4(\bar{\chi_{4}})+\tau^4(\chi_{4}))+2p^2\right). \end{align*}$

According to Lemma 1, we have

$\left(\tau^2(\chi_{4})+\tau^2(\bar{\chi_{4}})\right)^2 = \tau^4(\bar{\chi_{4}})+\tau^4(\chi_{4})+2p^2 = 4p\alpha^2.$

Therefore, we may immediately deduce

$\begin{align*} G^2(n) = &1-2(\chi_2(n)\cdot\left( G(n)+\chi_2(n)\right)\nonumber\\ &+\frac{1}{p}\left(\chi_2(n)\cdot(\tau^4(\bar{\chi_{4}})+\tau^4(\chi_{4}))+2p^2\right)\nonumber\\ = &1-2\chi_2(n)\cdot\left( G(n)+\chi_2(n)\right) \end{align*}$

(9)

$\begin{align*} &+\frac{1}{p}\cdot\left[\chi_2(n)((\tau^2(\bar{\chi_{4}})+\tau^2(\chi_{4}))^2-2p^2)+2p^2\right]\nonumber\\ = &2p-1-2\chi_2(n)\cdot G(n)+\left(4\alpha^2-2p\right)\cdot\chi_2(n)\nonumber, \end{align*}$

$\begin{align*} \; \; \; \; \; \; \; \; \; \; \; \; \; \; G^3(n) = &[-\chi_2(n)+\frac{1}{\sqrt{p}}\cdot(\chi_{4}(n)\cdot\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\cdot\tau^2(\chi_{4}))]^3\nonumber\\ = &\left(2p-1-2\chi_2(n)\cdot G(n)+\left(4\alpha^2-2p\right)\cdot\chi_2(n)\right)\cdot G(n) \end{align*}$

(10)

$\begin{align*} = &(4\alpha^2-2p)\chi_2(n)\cdot G(n)+(2p+3)G(n)-(4p-2)\chi_2(n) -2(4\alpha^2-2p)\nonumber \end{align*}$

and

$\begin{align*} \left[G^2(n)-(2p-1)\right]^2 = \left[\chi_2(n)\cdot(4\alpha^2-2p)-2\chi_2(n)\cdot G(n)\right]^2, \end{align*}$

which implies that

$\begin{align*} G^4(n) = (4p+2)\cdot G^{2}(n)+8(p-2\alpha^2)\cdot G(n)+[(4\alpha^2-2p)^2-(2p-1)^2]. \end{align*}$

(11)

So for any integer $k\geq4$ , from (8)–(11), we have the fourth-order linear recurrence formula

$\begin{align*} G^k(n)& = G^{k-4}(n)\cdot G^4(n)\\ & = (4p+2)\cdot G^{k-2}(n)+8(p-2\alpha^2)\cdot G^{k-3}(n)+[(4\alpha^2-2p)^2-(2p-1)^2]\cdot G^{k-4}(n). \end{align*}$

This proves Lemma 4.

3. Proof of the theorem

In this section, we will complete the proof of our theorem.

Let $p$ be any prime with $p\equiv 1\; \rm mod\; 4$ , then we have

$\begin{align*} K_0(p) = \frac{1}{p-1}\sum\limits_{n = 1}^{p-1}G^0(n) = \frac{p-1}{p-1} = 1.\; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \; \end{align*}$

(12)

$\begin{align*} K_1(p)& = \frac{1}{p-1}\sum\limits_{n = 1}^{p-1}G^1(n)\nonumber\\ & = \frac{1}{p-1}\sum\limits_{n = 1}^{p-1}\left(-\chi_2(n)+\frac{1}{\sqrt{p}}\cdot(\chi_{4}(n)\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\tau^2(\chi_{4}))\right)\\& = 0, \end{align*}$

(13)

$\begin{align*} \; \; \; \; \; K_2(p)& = \frac{1}{p-1}\sum\limits_{n = 1}^{p-1}G^2(n)\nonumber\\ & = \frac{1}{p-1}\sum\limits_{n = 1}^{p-1}\left(-\chi_2(n)+\frac{1}{\sqrt{p}}\cdot(\chi_{4}(n)\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\tau^2(\chi_{4}))\right)^2 \\ & = 2p+1,\end{align*}$

(14)

$\begin{align*} \; \; \; \; \; K_3(p)& = \frac{1}{p-1}\sum\limits_{n = 1}^{p-1}G^3(n)\nonumber\\ & = \frac{1}{p-1}\sum\limits_{n = 1}^{p-1}\left(-\chi_2(n)+\frac{1}{\sqrt{p}}\cdot(\chi_{4}(n)\tau^2(\bar{\chi_{4}})+\bar{\chi_{4}}(n)\tau^2(\chi_{4}))\right)^3 \\ & = -3(4\alpha^2-2p).\end{align*}$

(15)

It is clear that from Lemma 4, if $k\geq4$ , we have

$\begin{align*} \; \; \; \; \; K_k(p) = &\frac{1}{p-1}\sum\limits_{n = 1}^{p-1}G^k(n)\nonumber\\ = &(4p+2)\cdot K_{k-2}(p)-8\left(2\alpha^2-p\right)\cdot K_{k-3}(p) +(16\alpha^4-16p\alpha^2+4p-1)\cdot K_{k-4}(p). \end{align*}$

(16)

Now Theorem 1 follows (12)–(16). Obviously, using Theorem 1 to all negative integers, and that lead to Corollary 1.

This completes the proofs of our all results.

Some notes:

Note 1: In our theorem, know $n$ is an integer, and $(n, p) = 1$ . According to the properties of quadratic residual, $\chi_2(n) = \pm1$ , $\chi_4(n) = \pm1$ .

Note 2: In our theorem, we only discussed the case $p\equiv1\; \rm mod\; 8$ . If $p\equiv3\; \rm mod\; 4$ , then the result is trivial. In fact, in this case, for any integer $n$ with $(n, p) = 1$ , we have the identity

$\begin{align*} G(n)& = 1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+n\bar{a}^2}{p}\right) = 1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^4}{p}\right)\cdot\left(\frac{a^4+n}{p}\right)\\ & = 1+\sum\limits_{a = 1}^{p-1}\left(\frac{a}{p}\right)\cdot\left(\frac{a+n}{p}\right) = 1+\sum\limits_{a = 1}^{p-1}\left(\frac{a^2+na}{p}\right)\\ & = 1+\sum\limits_{a = 1}^{p-1}\left(\frac{1+n\bar{a}}{p}\right) = \sum\limits_{a = 0}^{p-1}\left(\frac{1+na}{p}\right) = 0. \end{align*}$

Thus, for all prime $p$ with $p\equiv3\; \rm mod\; 4$ and $k\geq1$ , we have $K_k(p) = 0$ .

4. Conclusions

The main result of this paper is Theorem 1. It gives an interesting computational formula for $K_k(p)$ with $p\equiv1\; \rm mod\; 4$ . That is, for any integer $k$ , we have the identity

$K_k(p) = (4p+2)\cdot K_{k-2}(p)-8\left(2\alpha^2-p\right)\cdot K_{k-3}(p)+(16\alpha^4-16p\alpha^2+4p-1)\cdot K_{k-4}(p).$

Thus, the problems of calculating a linear recurrence formula of one kind special character sums modulo a prime are given.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

The authors are grateful to the anonymous referee for very helpful and detailed comments.

This work is supported by the N.S.F. (11971381, 12371007) of China and Shaanxi Fundamental Science Research Project for Mathematics and Physics (22JSY007).

Conflict of interest

The authors declare no conflicts of interest.

References

[1]	World Health Organization, Anti-tuberculosis Drug Resistance in the World: Third Global Report, No. WHO / HTM / TB / 2004. 343, World Health Organization, 2004.
[2]	J. Cohen, Extensively drug-resistant TB gets foothold in South Africa, Science 313 (2006), 1554.
[3]	L. B. Reichman, J. H. Tanne, Timebomb: The Global Epidemic of Multi-Drug Resistant Tuberculosis, 2002.
[4]	Y. Zhou, K. Khan, Z. Feng, et al. Projection of tuberculosis incidence with increasing immigration trends, J. Theor. Biol., 254 (2008), 215–228.
[5]	P. Rodrigues, M. G. M. Gomes, C. Rebelo, Drug resistance in tuberculosis: a reinfection model, Theor. Popul. Biol., 71 (2007), 196–212.
[6]	N. Blaser, C. Zahnd, S. Hermans, et al. Tuberculosis in Cape Town: an age-structured transmission model, Epidemics, 14 (2016), 54–61.
[7]	C. P. Bhunu, W. Garira, Z. Mukandavire, et al. Tuberculosis transmission model with chemoprophylaxis and treatment, B. Math. Biol., 70 (2008), 1163–1191.
[8]	Centers for Disease Control and Prevention, TB in Children in the United States, CDC, 2014. Available from: https://www.cdc.gov/tb/topic/populations/tbinchildren/default.htm.
[9]	S. M. Blower, P. M. Small, P. C. Hopewell, Control strategies for tuberculosis epidemics: new models for old problems, Science, 273 (1996), 497–500.
[10]	Fatmawati, H. Tasman, An optimal treatment control of TB-HIV coinfection, International Journal of Mathematics and Mathematical Sciences, 2016 (2016).
[11]	R. I. Hickson, G. N. Mercer, K. M. Lokuge, A metapopulation model of tuberculosis transmission with a case study from high to low burden areas, PLoS One, 7 (2012).
[12]	R. M. G. J. Houben, T. Sumner, A. D. Grant, et al. Ability of preventive therapy to cure latent Mycobacterium tuberculosis infection in HIV-infected individuals in high-burden settings, P. Natl. A. Sci., 111 (2014), 5325–5330.
[13]	R. Kaplan, J. Caldwell, K. Middelkoop, et al. Impact of ART on TB case fatality stratified by CD4 count for HIV-positive TB patients in Cape Town, South Africa (2009-2011), J. Acq. Imm. Def., 66 (2014), 487–494.
[14]	Fatmawati, U. D. Purwati, F. Riyudha, et al. Optimal control of a discrete age-structured model for tuberculosis transmission, Heliyon, 6 (2020).
[15]	S. G, Samko, A. A, Kilbas, O. I. Marichev, Fractional Integrals and Derivatives: Theory and Applications, Gordon and Breach Science Publishers, 1993.
[16]	I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, California, USA, 1999.
[17]	T. Sardar, S. Rana, J. Chattopadhyay, A mathematical model of dengue transmission with memory, Commun. Nonlinear Sci., 22 (2015), 511–525.
[18]	J. Huo, H. Zhao, L. Zhu, The effect on backward bifurcation in a fractional orde HIV model, Nonlinear Analysis: Real World Applications, 26 (2015), 289–305.
[19]	M. Saeedian, M. Khalighi, N. Azimi-Tafreshi, et al. Memory effects on epidemic evolution: the susceptible-infected-recovered epidemic model, Phys. Rev. E., 95 (2017).
[20]	C. M. A. Pinto, A. R. M. Carvalho, The HIV/TB coinfection severity in the presence of TB multidrug resistant strains, Ecol. Complex., 32 (2017), 1–20.
[21]	Fatmawati, E. M. Shaiful, M. I. Utoyo, A fractional order model for HIV dynamics in a two-sex population, International Journal of Mathematics and Mathematical Sciences, 2018 (2018).
[22]	G. C. Wu, Z. G. Deng, D. Baleanu, et al. New variable-order fractional chaotic systems for fast image encryption, Chaos: An Interdisciplinary Journal of Nonlinear Science, 29 (2019).
[23]	G. C. Wu, T. Abdeljawad, J. Liu, et al. Mittag-Leffler stability analysis of fractional discrete-time neural networks via fixed point technique, Nonlinear Analysis: Modelling and Control, 24 (2019), 919–936.
[24]	G. C. Wu, D. Q. Zeng, D. Baleanu, Fractional impulisve differential equations: Exact solutions, integral equations and short memory case, Frac. Calc. Appl. Anal., 22 (2019), 180–192.
[25]	M. Itik, S. P. Banks, Chaos in a three-dimensional cancer model, Int. J. Bifurcat. Chaos, 20 (2010), 71–79.
[26]	M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Progr. Fract. Differ. Appl., 1 (2015), 1–13.
[27]	M. Caputo, M. Fabrizio, On the notion of fractional derivative and applications to the hysteresis phenomena, Meccanica, 52 (2017), 3043–3052.
[28]	T. Zhang, L. Xiong, Periodic motion for impulsive fractional functional differential equations with piecewise Caputo derivative, Appl. Math. Lett., 101 (2020), 106072.
[29]	A. Atangana, D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: theory and application to heat transfer model, Therm. Sci., 20 (2016).
[30]	A. Atangana, I. Koca, Chaos in a simple nonlinear system with Atangana-Baleanu derivatives with fractional order, Chaos, Solitons & Fractals, 89 (2016), 447–454.
[31]	A. Atangana, K. M. Owolabi, New numerical approach for fractional differential equations, Math. Model. Nat. Pheno., 13 (2018).
[32]	A. Atangana, J. F. Gomez-Aguilar, Decolonisation of fractional calculus rules: Breaking commutativity and associativity to capture more natural phenomena, The European Physical Journal Plus, 133 (2018).
[33]	M. A. Khan, S. Ullah, M. Farooq, A new fractional model for tuberculosis with relapse via Atangana-Baleanu derivative, Chaos, Solitons & Fractals, 116 (2018), 227–238.
[34]	E. Bonyah, Chaos in a 5-D hyperchaotic system with four wings in the light of non-local and non-singular fractional derivatives, Chaos, Solitons & Fractals, 116 (2018), 316–331.
[35]	K. Muhammad Altaf, A. Atangana, Dynamics of Ebola disease in the framework of different fractional derivatives, Entropy, 21 (2019).
[36]	R. Jan, M. A. Khan, P. Kumam, et al, Modeling the transmission of dengue infection through fractional derivatives, Chaos, Solitons & Fractals, 127 (2019), 189–216.
[37]	W. Wang, M. A. Khan, P. Kumam, et al. A comparison study of bank data in fractional calculus, Chaos, Solitons & Fractals, 126 (2019), 369–384.
[38]	M. A. Khan, The dynamics of a new chaotic system through the Caputo-Fabrizio and AtanaganBaleanu fractional operators, Adv. Mech. Eng., 11 (2019).
[39]	Fatmawati, M. A. Khan, M. Azizah, et al. A fractional model for the dynamics of competition between commercial and rural banks in Indonesia, Chaos, Solitons & Fractals, 122 (2019), 32–46.
[40]	S. Ullah, M. A. Khan, M. Farooq, et al. A fractional model for the dynamics of tuberculosis (TB) using Atangana-Baleanu derivative, Discrete Cont. Dyn. S, 13 (2019).
[41]	A. A. Velayati, Tuberculosis in children, International Journal of Mycobacteriology, 5 (2016).
[42]	C. Castillo-Chaves, B. Song, Dynamic models of tuberculosis and their applications, Math. Biosci. Eng., 1 (2004), 361–404.
[43]	Z. M. Odibat, N. T. Shawagfeh, Generalized Taylors formula, Appl. Math. Comput., 186 (2007), 286–293.
[44]	W. Lin, Global existence theory and chaos control of fractional differential equations, J. Math. Anal. Appl., 332 (2007), 709–726.
[45]	P. van den Driessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29–48.
[46]	K. Diethelm, N. J. Ford, A. D. Freed, A predictor-corrector approach for the numerical solution of fractional differential equations, Nonlinear Dynamics, 29 (2002), 3–22.
[47]	K. Diethelm, N. J. Ford, A. D. Freed, Detailed error analysis for a fractional Adams method, Numerical Algorithms, 36 (2004), 31–52.
[48]	C. P. Bhunu, Mathematical analysis of a three-strain tuberculosis transmission model, Appl. Math. Model., 35 (2011), 4647–4660.
[49]	S. Athithan, M. Ghosh, Optimal control of tuberculosis with case detection and treatment, World Journal of Modelling and Simulation, 11 (2015), 111–122.
[50]	Word Health Organization, Factsheet on the World Tuberculosis Report 2017, WHO, 2017. Available from: https://www.who.int/en/news-room/fact-sheets/detail/tuberculosis.
[51]	J. J. Tewa, S. Bowong, B. Mewoli, Mathematical analysis of two-patch model for the dynamical transmission of tuberculosis, Appl. Math. Model., 36 (2012), 2466–2485.
[52]	M. Toufik, A. Atangana, New numerical approximation of fractional derivative with non-local and non-singular kernel: application to chaotic models, The European Physical Journal Plus, 132 (2017), 444.
[53]	L. Xu, H. Hub, F. Qinc, Ultimate boundedness of impulsive fractional differential equations, Appl. Math. Lett., 62 (2016), 110–117.
[54]	L. Xu, J. Li, S. S. Ge, Impuls ivestabilization of fractional differential systems, ISA T., 70 (2017), 125–131.
[55]	L. Xu, X. Chu, H. Hu, Exponential ultimate boundedness of non-autonomous fractional differential systems with time delay and impulses, Appl. Math. Lett., 99 (2020), 106000.

This article has been cited by:

1.	Mounia Mouy, Hamid Boulares, Saleh Alshammari, Mohammad Alshammari, Yamina Laskri, Wael W. Mohammed, On Averaging Principle for Caputo–Hadamard Fractional Stochastic Differential Pantograph Equation, 2022, 7, 2504-3110, 31, 10.3390/fractalfract7010031
2.	Sabbavarapu Nageswara Rao, Manoj Singh, Ahmed Hussein Msmali, Abdullah Ali H. Ahmadini, Existence of Positive Solutions for a Coupled System of p-Laplacian Semipositone Hadmard Fractional BVP, 2023, 7, 2504-3110, 499, 10.3390/fractalfract7070499

Reader Comments

Your name:*

Email:*
© 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)