Next Article in Journal
Application of a Multiplier Transformation to Libera Integral Operator Associated with Generalized Distribution
Previous Article in Journal
Direct Training via Backpropagation for Ultra-Low-Latency Spiking Neural Networks with Multi-Threshold
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Repdigits as Sums of Four Tribonacci Numbers

School of Science, University of Science and Technology Liaoning, Anshan 114051, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Symmetry 2022, 14(9), 1931; https://doi.org/10.3390/sym14091931
Submission received: 30 August 2022 / Revised: 11 September 2022 / Accepted: 13 September 2022 / Published: 16 September 2022
(This article belongs to the Topic Engineering Mathematics)

Abstract

:
In this paper, we show that 66666 is the largest repdigit expressible as the sum of four tribonacci numbers. We used Binet’s formula, Baker’s theory, and a reduction method during the proving procedure. We also used the periodic properties of tribonacci number modulo 9 to deal with three individual cases.

1. Introduction

A palindromic number is a number that has reflectional symmetry across a vertical axis. A repdigit is a palindromic number N that only has one distinct digit when it is represented in base 10. That is, N has the form
d 10 l 1 9 ,
for some positive integers d and l with 1 d 9 and l 1 .
The problem of finding all repdigits that are perfect powers was posed by Obláth [1] and settled in 1999 by Bugeaud and Mignotte [2]. In 2000, Luca [3] proved that the largest repdigits in Fibonacci and Lucas sequences were F 10 = 55 and L 5 = 11 . Afterward, finding repdigits in recurrence sequences has drawn much attention in the literature. For instance, Dıaz-Alvarado and Luca [4] proved that the largest Fibonacci number that can be represented as a sum of two repdigits is F 20 = 6765 = 6666 + 99 . Luca [5] answered the question of finding all repdigits as sums of three Fibonacci numbers. Luca, Normenyo, and Togbé [6] obtained all repdigits expressible as sums of three Lucas numbers. They also solved a similar problem involving three Pell numbers in [7]. Ddamulira [8] identified all repdigits as sums of three balancing numbers. In [9], the authors confirmed a conjecture by Luca [5]. More precisely, they determined all repdigits expressible as sums of four Fibonacci or Lucas numbers. Afterward, in [10], the authors obtained analogous results for Pell numbers. Keskin and Erduvan [11] tackled the same problem with four balancing numbers.
Recently, Ddamulira [12] found all repdigits as sums of three Padovan numbers. As far as we know, this is the only reference involving repdigits that are expressible as sums of more than two numbers in a high-order recurrence sequence. Compared to the second-order recurrence sequence, it is difficult to solve a similar question with a high-order recurrence sequence. Therefore, it is interesting to find all repdigits that are sums of four numbers in some other third-order recurrence sequence.
In this paper, we investigate the presence of repdigits as sums of four tribonacci numbers. More precisely, we prove the following theorem.
Theorem 1.
All non-negative integer solutions m 1 , m 2 , m 3 , m 4 , N , d , of the Diophantine equation
N = T m 1 + T m 2 + T m 3 + T m 4 = d 10 l 1 9 with d { 1 , , 9 }
have
N 4 , 5 , 6 , 7 , 8 , 9 , 11 , 22 , 33 , 44 , 55 , 66 , 77 , 88 , 99 , 111 , 333 , 555 , 666 , 999 , 2222 , 3333 , 66666 .
Our paper is organized as follows. In the following section, we recall some important results that are useful for the proof of our main theorem. We use them in Section 3 to prove Theorem 1. During the proof, first, we use Baker’s method several times to obtain a bound of m 1 , which is too large to conduct a brute force search. We then apply the reduction method of de Weger several times to find a very low bound for m 1 , which enables us to run a simple computer program in Mathematica to find the small solutions. It is worth mentioning that the reduction method is invalid in three cases during the computations. For these cases, we use periodic properties of { T n } ( mod 9 ) to reach the contradictions.

2. Auxiliary Results

We use a definition of tribonacci numbers (see [13]) with little difference from the common one that starts from 0 , 0 , 1 . Let T n n 0 be the tribonacci sequence satisfying the recurrence relation T n + 3 = T n + 2 + T n + 1 + T n with initial conditions T 0 = 0 and T 1 = T 2 = 1 . The first few terms of this sequence are
0 , 1 , 1 , 2 , 4 , 7 , 13 , 24 , 44 , 81 , 149 , 274 , 504 , 927 , 1705 , 3136 , 10609 , 19513 .
Its characteristic equation, z 3 z 2 z 1 = 0 , has one real root α and two complex roots β and γ = β ¯ . In 1982, Spickerman [14] found the following “Binet-like” formula:
T n = a α n + b β n + c γ n , for all n 0 ,
where
a = 1 α 2 + 4 α 1 , b = 1 β 2 + 4 β 1 , c = 1 γ 2 + 4 γ 1 = b ¯ .
Numerically, we have 1.83 α 1.84 , 0.73 | β | = | γ | = α 1 2 0.74 , 0.33 a 0.34 and 0.25 | b | = | c | 0.26 . It follows that the complex conjugate roots β and γ have little influence on the right side of Equation (2), setting
e n : = T n a α n = b β n + c γ n then | e ( n ) | 1 α n 2
holds for all n 1 . In addition, it is known that:
α n 2 T n α n 1 holds for all n 1 .
Let η ( 1 ) = η be an algebraic number of degree d with a minimal primitive polynomial
f ( X ) = a 0 i = 1 d ( X η ( i ) ) Z [ X ] ,
where the positive integer a 0 is the leading coefficient and η ( i ) ( i = 2 , , d ) are the conjugates of η . The logarithmic height of η is given by
h η = 1 d log | a 0 | + i = 1 d log max η i , 1 .
In particular, if η = p / q is a rational number with gcd ( p , q ) = 1 and q > 0 , then h ( η ) = log max { | p | , q } . The following properties of h ( η ) will be used in the next section.
Lemma 1
([15]). Let a and b be algebraic numbers. Then,
h a ± b h a + h b + log 2 , h a b ± h a + h b , h a r = | r | h a ( r Z ) .
We recall a variation of a result of Matveev [16] due to Bugeaud, Mignotte, and Siksek [17], which helps us to give an upper bound of m 1 .
Lemma 2
([17], Theorem 9.4). Let η 1 , , η t be positive real algebraic numbers in a real algebraic number field K R of degree D K , b 1 , , b t be nonzero integers and assume that
Λ : = η 1 b 1 η t b t 1
is nonzero. Then
log | Λ | > 1.4 × 30 t + 3 × t 4.5 × D K 2 1 + log D K 1 + log B A 1 A 2 A t ,
where
B max | b 1 | , , | b t | ,
and
A i max D K h η i , | log η i | , 0.16 , for all i = 1 , , t .
To reduce the bound we obtain from Lemma 2, we introduce a modified version of the Baker and Davenport reduction method that appears in [18]. Let v 1 , v 2 , β R , and let x 1 , x 2 Z be unknown. Let
Λ = β + x 1 v 1 + x 2 v 2 .
Let c , δ be positive constants. Set X = max { | x 1 | , | x 2 | } . Let X 0 be a (large) positive constant. Assume that
| Λ | < c · exp ( δ · Y ) ,
X X 0 .
When β = 0 in (6), we have
Λ = x 1 v 1 + x 2 v 2 .
Put v = v 1 / v 2 . Let the continued fraction expansion of v be given by
[ a 0 , a 1 , a 2 , ] ,
and let the k-th convergent of v be p k / q k for k = 0 , 1 , 2 , . We may assume without loss of generality that | v 1 | < | v 2 | and that x 1 > 0 . We have the following result.
Lemma 3
([18], Lemma 3.1). (1)   If (7) and (8) hold for x 1 , x 2 with X X * , then ( x 2 , x 1 ) = ( p k , q k ) for an index k that satisfies
k 1 + log ( 1 + X 0 5 ) log ( 1 + 5 2 ) : = Y 0 .
Moreover, the partial quotient a k + 1 satisfies
a k + 1 > 2 + | v 2 | exp ( δ q k ) c q k .
(2)   If for some k with q k X * , we have
a k + 1 > | v 2 | exp ( δ q k ) c q k ,
then (7) holds for ( x 2 , x 1 ) = ( p k , q k ) .
Lemma 4
([18], Lemma 3.2). Let
A = max 0 k Y 0 a k + 1 .
If (7) and (8) hold for x 1 , x 2 and β = 0 , then
Y < 1 δ log c ( A + 2 ) X 0 | v 2 | .
When β v 1 v 2 0 in (6), put v = v 1 / v 2 and ψ = β / v 2 . Then we have
Λ v 2 = ψ x 1 v + x 2 .
Let p / q be a convergent of v with q > X 0 . For a real number x, we define
x = min { | x n | , n Z }
be the distance from x to the nearest integer. We have the following Davenport lemma.
Lemma 5
([18], Lemma 3.3). Suppose that
q ψ > 2 X 0 q .
Then, the solution of (7) and (8) satisfy
Y < 1 δ log q 2 c | v 2 | X 0 .
Finally, the following result of Le [19] will help us to deal with some inequalities involving logarithms.
Lemma 6
([19]). Let f ( y ) R [ x ] be a polynomial with degree n and f ( m ) ( y ) be its m-th derivative. If there is a real number x 0 satisfying
x 0 > max 0 , f log x 0 , f ( 1 ) log x 0 , , f ( n ) log x 0 ,
then x f ( log x ) > 0 when x x 0 .

3. Proof of Theorem 1

3.1. Bounding the Variables

We assume that m 1 m 2 m 3 m 4 0 . From (1) and (4), we have
α m 1 2 T m 1 T m 1 + T m 2 + T m 3 + T m 4 = d 10 1 9 10
and
10 1 d 10 1 9 T m 1 + T m 2 + T m 3 + T m 4 4 T m 1 < α m 1 + 2 ,
where we use 4 < α 3 . Thus,
( m 1 2 ) log α log 10
and
1 ( m 1 + 2 ) log α log 10 .
Since
1 6 < log α log 10 = 0.264649 < 1 3 ,
we have
m 1 2 6 < < m 1 + 5 3
by (15) and (16). Therefore, if m 1 330 , then 1 111 . Running a Mathematica program in the range 0 m 4 m 3 m 2 m 1 330 , 1 d 9 , and 1 111 , we find no other solutions except those listed in Theorem 1.
From now on, we assume that m 1 > 330 . By using (3), Equation (1) can be written as
a α m 1 + e ( m 1 ) + a α m 2 + e ( m 2 ) + a α m 3 + e ( m 3 ) = d 10 1 9 .
We then consider (17) in four different cases as follows.

3.2. Case 1

Identity (17) is equivalent to
a α m 1 d · 10 9 = d 9 a ( α m 2 + α m 3 + α m 4 ) e ( m 1 ) e ( m 2 ) e ( m 3 ) e ( m 4 ) .
Taking the absolute value on both sides of (18), estimate the size of the right-hand side of the equality. We have
a α m 1 d · 10 9 d 9 + a ( α m 2 + α m 3 + α m 4 ) + | e ( m 1 ) | + | e ( m 2 ) | + | e ( m 3 ) | + | e ( m 4 ) | < 1 + 3 a α m 2 + 4 α m 4 / 2 < 3 a + 3 a α m 2 + 12 a α ( 2 m 2 m 4 ) / 2 < 18 a α m 2 .
To apply Lemma 2. We divide (through (19)) by a α m 1 to obtain
10 · α m 1 d 9 a 1 < 18 α m 1 m 2 .
Then, let
Λ 1 : = 10 · α m 1 d 9 a 1 .
It is sufficient to check that Λ 1 0 . Suppose that Λ 1 = 0 , then we have
10 · d 9 = a α m 1 .
Now, we consider the Q -automorphism σ of the Galois extension Q α , β over Q given by σ α : = β and σ β : = α . Since Λ 1 = 0 , we have σ Λ 1 = 0 . Thus, conjugating the relation (22) under σ , and taking absolute values on both sides, we have
10 · d 9 = | σ a α m 1 | = b | β | m 1 < | b | < 1 2 ,
which is false for 1 and d 1 . Therefore, Λ 1 0 .
Then, we apply Lemma 2 with the data
η 1 : = 10 , η 2 : = α , η 3 : = d 9 a , b 1 : = , b 2 : = m 1 , b 3 : = 1 , t : = 3 .
The minimal polynomial of a is 44 x 3 2 x 1 and has roots a , b , c . Since | b | = | c | < | a | < 1 , we have
h a = 1 3 log 44 .
Let K : = Q ( α ) . Then D K = 3 because η 1 , η 2 , η 3 Q α . Since max , m 1 , 1 m 1 , we take B : = m 1 . Further, the minimal polynomial of α over Z is x 3 x 2 x 1 has roots α , β , γ with 1.83 α 1.84 and | β | = | γ | < 1 . Thus, we have h α = 1 3 log α . Since h 10 = log 10 and
h η 3 h d + h 9 + h a 4 log 3 + 1 3 log 44 < 6 log 3 ,
we take A 1 : = 3 log 10 , A 2 : = log α and A 3 : = 18 log 3 . Then, from Lemma 2, (21) is bounded below by
log | Λ 1 | > 1.4 × 30 6 × 3 4.5 × 3 2 1 + log 3 1 + log m 1 3 log 10 log α 18 log 3 > 7.39 × 10 14 log m 1 log α .
Combining the above inequality with (20), we have
m 1 m 2 7.41 × 10 14 log m 1 .

3.3. Case 2

Identity (17) is also equivalent to
a ( α m 1 + α m 2 ) d · 10 9 = d 9 a ( α m 3 + α m 4 ) e ( m 1 ) e ( m 2 ) e ( m 3 ) e ( m 4 ) .
Thus, it follows that
a ( α m 1 + α m 2 ) d · 10 9 d 9 + a ( α m 3 + α m 4 ) + | e ( m 1 ) | + | e ( m 2 ) | + | e ( m 3 ) | + | e ( m 4 ) | < 1 + 2 a α m 3 + 4 α m 4 / 2 < 3 a + 2 a α m 3 + 12 a α ( 2 m 3 m 4 ) / 2 < 17 a α m 3 .
As before, we divide both sides of (23) by a α m 1 + α m 2 to have
10 α m 2 d 9 a ( 1 + α m 1 m 2 ) 1 < 17 α m 3 m 1 1 + α m 2 m 1 ,
which yields
10 α m 2 d 9 a 1 + α m 1 m 2 1 < 17 α m 1 m 3 .
Let
Λ 2 : = 10 α m 2 d 9 a 1 + α m 1 m 2 1 .
Again, we need to check that Λ 2 0 . Suppose that Λ 2 = 0 , then
a α m 1 + α m 2 = 10 · d 9 .
We consider the Q -automorphism σ of the Galois extension Q α , β over Q given by σ α : = β and σ β : = α . We have σ Λ 2 = 0 because Λ 2 = 0 . Thus, conjugating the relation (26) under σ , and taking the absolute values on both sides, we have
10 · d 9 = | σ a α m 1 + α m 2 | = | b | | β | m 1 + | β | m 2 < 2 | b | < 1 ,
which is false for 1 and d 1 . Therefore, Λ 2 0 . Hence, we apply Lemma 2 with the data
η 1 : = 10 , η 2 : = α , η 3 : = d 9 a 1 + α m 1 m 2 , b 1 : = , b 2 : = m 2 , b 3 : = 1 , t : = 3 .
Let K : = Q ( α ) . Then D K = 3 because η 1 , η 2 , η 3 Q α . Since max , m 1 , 1 m 1 , we take B : = m 1 . Since
h η 3 h d + h 9 + h a + h 1 + α m 1 m 2 4 log 3 + 1 3 log 44 + m 1 m 2 log α + log 2 < 4.52 × 10 14 log m 1 ,
we can take A 1 : = 3 log 10 , A 2 : = log α and A 3 : = 1.36 × 10 15 log m 1 . So, Lemma 2 reveals that (25) is bounded below by
log | Λ 2 | > 1.4 × 30 6 × 3 4.5 × 3 2 1 + log 3 1 + log m 1 3 log 10 log α × 1.36 × 10 15 log m 1 > 5.09 × 10 28 log m 1 2 log α .
Combining the above inequality with (24), we have
m 1 m 3 5.11 × 10 28 log m 1 2 .

3.4. Case 3

Rewriting (17) as below:
a ( α m 1 + a α m 2 + α m 3 ) d · 10 9 = d 9 a α m 4 e ( m 1 ) e ( m 2 ) e ( m 3 ) e ( m 4 ) .
Thus, it follows that
a ( α m 1 + α m 2 + a α m 3 ) d · 10 9 d 9 + a α m 4 + | e ( m 1 ) | + | e ( m 2 ) | + | e ( m 3 ) | + | e ( m 4 ) | < 1 + a α m 4 + 4 α m 4 2 < 3 a + a α m 4 + 12 a α m 4 2 < 16 a α m 4 ,
Dividing through (27) by a α m 1 + α m 2 + α m 3 , we have
10 α m 3 d 9 a ( 1 + α m 1 m 3 + α m 2 m 3 ) 1 < 16 α m 4 m 1 1 + α m 2 m 1 + α m 3 m 1 .
This means that
10 α m 3 d 9 a 1 + α m 1 m 3 + α m 2 m 3 1 < 16 α m 1 m 4 .
Thus, we put
Λ 3 : = 10 α m 3 d 9 a 1 + α m 1 m 3 + α m 2 m 3 1 .
Then, we need to verify that Λ 3 0 . Suppose that Λ 3 = 0 , then we have
a α m 1 + α m 2 + α m 3 = 10 · d 9 .
To see that this is not true, we again consider the Q -automorphism σ of the Galois extension Q α , β over Q given by σ α : = β and σ β : = α . Since Λ 3 = 0 , we have σ Λ 3 : = 0 . Thus, conjugating the relation (30) under σ , and taking absolute values on both sides, we have
10 · d 9 = | σ a α m 1 + α m 2 α m 3 | = | b | | β | m 1 + | β | m 2 + | β | m 3 < 3 | b | < 3 2 ,
which is false for 1 and d 1 . Therefore, Λ 3 0 . So, we apply Lemma 2 with the data
η 1 : = 10 , η 2 : = α , η 3 : = d 9 a 1 + α m 1 m 3 + α m 2 m 3 , b 1 : = , b 2 : = m 3 , b 3 : = 1 , t : = 3 .
Since η 1 , η 2 , η 3 Q ( α ) , we take the field K : = Q α with degree D K : = 3 . Since max , m 1 , 1 m 1 , we take B : = m 1 . Further,
h η 3 h d + h 9 + h a + h 1 + α m 1 m 3 + α m 2 m 3 4 log 3 + 1 3 log 44 + m 1 m 3 log α + m 2 m 3 log α + 2 log 2 4 log 3 + 1 3 log 44 + m 1 m 2 + 2 m 2 m 3 log α + 2 log 2 < 9.36 × 10 28 log m 1 2 .
Thus, we can take A 1 : = 3 log 10 , A 2 : = log α , and A 3 : = 2.81 × 10 29 log m 1 2 . According to Lemma 2, (29) is bounded below by
log | Λ 2 | > 1.4 × 30 6 × 3 4.5 × 3 2 1 + log 3 1 + log m 1 3 log 10 log α × 2.81 × 10 29 log m 1 2 > 1.05 × 10 43 ( log m 1 ) 3 log α .
By comparing the above inequality with the right-hand side of (28), we have
m 1 m 4 1.07 × 10 43 ( log m 1 ) 3 .

3.5. Case 4

Equation (17) is equivalent to
a ( α 1 m + α 2 m + α 3 m + α 4 m ) d · 10 9 = 1 9 e ( m 1 ) e ( m 2 ) e ( m 3 ) e ( m 4 ) .
Taking the absolute value on both sides, (31) shows that
a ( α 1 m + α 2 m + α 3 m + α 4 m ) d · 10 9 1 9 + | e ( m 1 ) | + | e ( m 2 ) | + | e ( m 3 ) | + | e ( m 4 ) | < 1 + 4 α m 4 2 < 3 a + 12 a α m 4 2 < 15 a .
Dividing both sides by a ( α m 1 + α m 2 + α m 3 + α m 4 ) , we have
10 α m 4 d 9 a ( α m 1 m 4 + α m 2 m 4 + α m 3 m 4 ) 1 < 15 a α m 1 a ( α m 2 m 1 + α m 3 m 1 α m 4 m 1 + 1 ) ,
which implies
10 α m 4 d 9 a ( α m 1 m 4 + α m 2 m 4 + α m 3 m 4 ) 1 < 15 α m 1 .
Put
Λ 4 : = 10 α m 4 d 9 a ( α m 1 m 4 + α m 2 m 4 + α m 3 m 4 ) 1 .
As before, one can justify that Λ 4 0 . Suppose that Λ 4 = 0 , then we have
a ( α m 1 + α m 2 + α m 3 + α m 4 ) = 10 · d 9 .
To see that this is not true, we consider the Q -automorphism σ of the Galois extension Q ( α , β ) over Q given by σ ( α ) : = β and σ ( β ) : = α . Now, since Λ 4 = 0 , we have σ ( Λ 4 ) : = 0 . Thus, conjugating the relation (34) under σ , and taking absolute values on both sides, we have
10 · d 9 = | σ a α m 1 + α m 2 α m 3 + α m 4 | = | b | | β | m 1 + | β | m 2 + | β | m 3 + | β | m 4 < 4 | b | < 2 ,
which is false for 2 and d 1 . Therefore, Λ 4 0 .
Thus, we apply Lemma 2 with the data
η 1 : = 10 , η 2 : = α , η 3 : = d 9 a ( 1 + α m 1 m 4 + α m 2 m 4 + α m 3 m 4 ) , b 1 : = , b 2 : = m 4 , b 3 : = 1 , t : = 3 .
where η 1 , η 2 , η 3 Q ( α ) . We can take the field K : = Q ( α ) with degree D K : = 3 . As max , m 1 , 1 m 1 , we can also take B : = m 1 . Further,
h ( η 3 ) h ( d ) + h ( 9 ) + h ( a ) + h ( 1 + α m 1 m 4 + α m 2 m 4 + α m 3 m 4 ) 4 log 3 + 1 3 log 17 + ( m 1 m 4 ) log α + ( m 2 m 4 ) log α + ( m 3 m 4 ) log α + 3 log 2 4 log 3 + 1 3 log 17 + ( ( m 1 m 2 ) + 2 ( m 2 m 3 ) log α + 3 ( m 3 m 4 ) ) log α + 3 log 2 < 3.92 × 10 43 ( log m 1 ) 3 .
Thus, we choose A 1 : = 3 log 10 , A 2 : = log α and A 3 : = 1.18 × 10 44 ( log m 1 ) 3 . From Lemma 2, we obtain
log | Λ 4 | > 1.4 × 30 6 × 3 4.5 × 3 2 ( 1 + log 3 ) ( 1 + log m 1 ) ( 3 log 10 ) ( log α ) × 1.18 × 10 44 ( log m 1 ) 3 > 4.41 × 10 57 ( log m 1 ) 4 log α ,
which, combined with (33), gives us
m 1 4.43 × 10 57 ( log m 1 ) 4 .
From Lemma 6, we can choose f ( x ) : = 4.43 × 10 57 x 4 . Then, we obtain m 1 2.42 × 10 66 . We record the above as the following Lemma.
Lemma 7.
Let ( m 1 , m 2 , m 3 , m 4 , N , d , ) be the nonnegative integer solutions to the Diophantine Equation (1) with m 1 m 2 m 3 m 4 0 , 1 d 9 , and 2 . Then we have
< m 1 2.42 × 10 66 .

4. Reducing The Bounds

The bounds obtained in Lemma 7 are too large to carry out meaningful computations on the computer. Thus, we need to reduce these bounds. To do so, we return to (20), (24), (28) and (33) apply Lemma 5 via the following procedure.

4.1. Step 1

First, let
τ 1 : = log 10 m 1 log α + log d 9 a , 1 d 9 .
For technical reasons, we assume that m 1 m 2 20 for the moment and go to (19). We will obtain a bound of m 1 m 2 larger than 20. Thus, we can get rid of this condition in both cases. Note that e τ 1 1 = Λ 1 0 . Thus, τ 1 0 . If τ 1 > 0 , then
0 < τ 1 < e τ 1 1 = | Λ 1 | < 18 α m 1 m 2 .
If τ 1 < 0 , by (13), we have
e τ 1 1 = 9 a α m 1 d 10 1 < 9 · 18 a α m 2 d 10 α m 2 + 7 α m 1 2 1 α m 1 m 2 9 1 α 11 < 1 2 .
Hence, e τ 1 < 2 . Thus, we have
0 < | τ 1 | < e | τ 1 | 1 = e | τ 1 | Λ 1 < 36 α m 1 m 2 .
Therefore, in both cases, we have
0 < τ 1 = log 10 m 1 log α + log d 9 a < 36 α m 1 m 2 .
This means that
| τ 1 | < 36 α m 2 m 1 < α m 2 m 1 + 6 < α 6.1 exp ( 0.609 ( m 1 m 2 ) ) ,
with
X = max { m 1 , n } = m 1 2.42 × 10 66 .
Dividing through (35) by log 10 , we have
τ 1 log 10 = log ( d 9 a ) log 10 m 1 log α log 10 + n .
Then, we put
c 1 : = α 6.1 , δ : = 0.609 , X 0 : = 2.42 × 10 66 , ψ : = log ( d 9 a ) log 10 , Y : = m 1 m 2 , v : = log α log 10 , v 1 : = log α , v 2 : = log 10 , β : = log d 9 a .
We now apply Lemma 5 on (36). A quick computer search in Mathematica reveals that the convergent
p 130 q 130 = 1779234883646329125716285138173060634490565368549796513064476646863 6722987436594440887072037743863558091187499203687945561784124798326
of τ is such that q 130 > X 0 . Therefore, we find that q = q 134 satisfies the hypothesis of Lemma 5 for d = 1 , , 9 . Applying Lemma 5, we have m 1 m 2 273 .

4.2. Step 2

Next, we put
τ 2 : = log 10 m 2 log α + log d 9 a 1 + α m 1 m 2 , 1 d 9 .
For technical reasons, as before, we assume that m 1 m 3 20 for the moment and go to (23). We will obtain a bound of m 1 m 3 larger than 20. Thus, we can get rid of this condition in both cases. Note that e τ 2 1 = Λ 2 0 . Thus, τ 2 0 . If τ 2 > 0 , then
0 < τ 2 < e τ 2 1 = | Λ 2 | < 17 α m 1 m 3 .
If τ 2 < 0 , we have
e τ 2 1 = 9 a α m 1 + α m 2 d 10 1 < 9 · 17 a α m 3 d 10 α m 3 + 7 α m 1 2 1 α m 1 m 3 9 1 α 11 < 1 2 .
Hence, e τ 2 < 2 . Thus, we have
0 < | τ 2 | < e | τ 2 | 1 = e τ 2 Λ 2 < 34 α m 1 m 3 .
Therefore, in both cases, we have
0 < | τ 2 | = log 10 m 2 log α + log d 9 a 1 + α m 1 m 2 < 36 α m 1 m 3 .
This means that
τ 2 < α m 3 m 1 + 6 < α 6.1 exp ( 0.609 ( m 1 m 3 ) ) .
Dividing through (37) by log 10 , we have
τ 2 log 10 = 1 log 10 log d 9 a ( α m 1 m 2 + 1 ) m 2 log α log 10 + .
Thus, we put
c 1 : = α 6.1 , δ : = 0.609 , X 0 : = 2.42 × 10 66 , Y : = m 1 m 3 , v : = log α log 10 , v 1 : = log α , v 2 : = log 10 , ψ : = 1 log 10 log d 9 ( α m 1 m 2 + 1 ) , β : = log d 9 ( α m 1 m 2 + 1 ) .
We now apply Lemma 5 on (38). We found that q = q 138 satisfies the hypothesis of Lemma 5 for d = 1 , , 9 . Thus, we have m 1 m 3 298 . Hence, m 3 32 by the assumption m 1 330 .

4.3. Step 3

Now, we put
τ 3 : = log 10 m 3 log α + log d 9 a 1 + α m 1 m 3 + α m 2 m 3 , 1 d 9 .
For technical reasons, we assume that m 1 m 4 20 for the moment and go to (27). We will obtain a bound of m 1 m 4 larger than 20. Thus, we can remove this condition in both cases. Note that e τ 3 1 = Λ 3 0 . Thus, τ 3 0 . If τ 3 > 0 , then
0 < τ 3 < e τ 3 1 = | Λ 3 | < 16 α m 4 m 1 .
If τ 3 < 0 , we have
e τ 3 1 = 9 a α 1 m + α 2 m + α 3 m d 10 1 < 9 · 16 a α m 4 d 10 α m 4 + 7 α m 1 2 1 α m 1 m 4 9 1 α 11 < 1 2 .
Hence, e τ 3 < 2 . Thus, we have
0 < | τ 3 | < e | τ 3 | 1 = e τ 3 Λ 3 < 32 α m 1 m 4 .
Therefore, in both cases, we have
0 < | τ 3 | = log 10 m 3 log α + log d 9 a 1 + α m 1 m 3 + α m 2 m 3 < 36 α m 1 m 4 .
This means that
τ 3 < α m 4 m 1 + 6 < α 6.1 exp ( 0.609 ( m 1 m 4 ) ) .
Dividing through (39) by log 10 , we have
τ 3 log 10 = 1 log 10 log d 9 a ( α m 1 m 3 + α m 2 m 3 + 1 ) m 3 log α log 10 +
Thus, we can take
c 1 : = α 6.1 , δ : = 0.609 , X 0 : = 4.1 × 10 66 , Y : = m 1 m 4 , v : = log α log 10 , v 1 : = log α , v 2 : = log 10 , ψ : = 1 log 10 log d 9 a ( α m 1 m 3 + α m 2 m 3 + 1 ) , β : = log d 9 a ( α m 1 m 3 + α m 2 m 3 + 1 ) .
We find that q = q 143 satisfies the hypothesis of Lemma 5 for 1 d 9 , 0 m 2 m 3 m 1 m 3 298 . Applying Lemma 5, we have m 1 m 4 313 and, hence, m 4 17 .

4.4. Step 4

Lastly, we put
τ 4 : = log 10 m 4 log α + log d 9 a 1 + α m 1 m 4 + α m 2 m 4 + α m 3 m 4 1 d 9 .
We use the original assumption that m 1 > 330 and go to (32). Note that e τ 4 1 = Λ 4 0 . Thus, τ 4 0 . If τ 4 > 0 , then
0 < τ 4 < e τ 4 1 = | Λ 4 | < 15 α m 1 .
If τ 4 < 0 , we have
e τ 4 1 = 9 a α m 1 + α m 2 + α m 3 + α m 4 d 10 1 < 9 · 15 a d 10 α 7 α m 1 2 1 α m 1 9 1 α 11 < 1 2 .
Hence, e τ 4 < 2 . Thus, we have
0 < | τ 4 | < e | τ 4 | 1 = e τ 4 Λ 4 < 30 α m 1 .
Therefore, in both cases, we have
0 < | τ 4 | = log 10 m 4 log α + log d 9 a 1 + α m 1 m 4 + α m 2 m 4 + α m 3 m 4 < 36 α m 1 .
This means that
τ 4 < α m 1 + 6 < α 6.1 exp ( 0.609 m 1 ) .
Dividing through (40) by log 10 , we have
τ 4 log 10 = 1 log 10 log d 9 a ( α m 1 m 4 + α m 2 m 4 + α m 3 m 4 + 1 ) m 4 log α log 10 +
Thus, we can take
c 1 : = α 6.1 , δ : = 0.609 , X 0 : = 4.1 × 10 66 , Y : = m 1 , v : = log α log 10 , v 1 = log α , v 2 = log 10 , ψ : = 1 log 10 log d 9 a ( α m 1 m 4 + α m 2 m 4 + α m 3 m 4 + 1 ) , β : = log d 9 a ( α m 1 m 4 + α m 2 m 4 + α m 3 m 4 + 1 ) .
We find that q = q 147 satisfies the hypothesis of Lemma 5 for 1 d 9 , 0 m 3 m 4 m 2 m 4 m 1 m 4 313 except for three special cases ( d , m 1 m 4 , m 2 m 4 , m 3 m 4 ) = ( 9 , 3 , 0 , 1 ) , ( 9 , 6 , 5 , 1 ) , ( 9 , 6 , 4 , 4 ) . Applying Lemma 5, we have m 1 322 , which contradicts the assumption that m 1 > 330 . Now, we consider the three cases
( d , m 1 m 4 , m 2 m 4 , m 3 m 4 ) = ( 9 , 3 , 1 , 0 ) , ( 9 , 6 , 5 , 1 ) , ( 9 , 6 , 4 , 4 ) .
Obviously, d ( 10 m 1 ) / 9 0 ( mod 9 ) when d = 9 . It is easy to see that the period of tribonacci numbers modulo 9 is 39. Since
T n + T n + T n + 1 + T n + 3 ¬ 0 ( mod 9 ) , T n + T n + 1 + T n + 5 + T n + 6 ¬ 0 ( mod 9 ) , T n + T n + 4 + T n + 4 + T n + 6 ¬ 0 ( mod 9 )
for 1 n 39 , there is no solution to Equation (1) in the above cases.  □

5. Conclusions

In this paper, we completely solved the diophantine Equation (1). More precisely, we found that 66666 is the largest repdigit expressible as the sum of four tribonacci numbers. Our method is based on Baker’s method. We first gave a larger upper bound of m 1 . Then, the reduction method reduced such a bound to an applicable one. During the reduction procedure, the periodic properties of { T n } ( mod 9 ) were used to deal with three individual cases.
It is worth mentioning that our method could be applied in b-repdigits. For each b, we may give a large bound. It would be a challenge to find all solutions for every b. We may not find a unified bound of b. It likely has infinite solutions.

Author Contributions

Conceptualization, Y.Z. and P.Y.; validation, S.Z. and K.Z.; writing—original draft preparation, Y.Z. and S.Z.; writing—review and editing, P.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Foundation of Liaoning Educational Committee, project 2019LNJC08.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors express their sincere gratitude to the referees for their valuable comments and advice to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Obláth, R. Une proprété des puissances parfaites. Mathesis 1956, 65, 356–364. [Google Scholar]
  2. Bugeaud, Y.; Mignotte, M. On integers with identical digits. Mathematika 1999, 46, 411–417. [Google Scholar] [CrossRef]
  3. Luca, F. Fibonacci and Lucas Numbers with Only One Distinct Digit. Port. Math. 2000, 57, 242–254. [Google Scholar]
  4. Dıaz-Alvarado, S.; Luca, F. Fibonacci numbers which are sums of two repdigits. In Proceedings of the XIVth International Conference on Fibonacci numbers and their applications, Morelia, Mexico, 5 July 2011; pp. 5–9. [Google Scholar]
  5. Luca, F. Repdigits as sums of three Fibonacci numbers. Math. Commun. 2012, 17, 1–11. [Google Scholar]
  6. Luca, F.; Normenyo, B.V.; Togbé, A. Repdigits as sums of three Lucas numbers. Colloq. Math. 2019, 156, 255–265. [Google Scholar] [CrossRef]
  7. Normenyo, B.V.; Luca, F.; Togbé, A. Repdigits as sums of three Pell numbers. Period. Math. Hungar. 2018, 77, 318–328. [Google Scholar] [CrossRef]
  8. Ddamulira, M. Repdigits as sums of three balancing numbers. Math. Slovaca 2020, 70, 557–566. [Google Scholar] [CrossRef]
  9. Normenyo, B.V.; Luca, F.; Togbé, A. Repdigits as Sums of Four Fibonacci or Lucas Numbers. J. Integer Seq. 2018, 21, Article 18.7.7. [Google Scholar]
  10. Luca, F.; Normenyo, B.V.; Togbé, A. Repdigits as sums of four Pell numbers. Bol. Soc. Mat. Mex. 2019, 25, 249–266. [Google Scholar] [CrossRef]
  11. Keskin, R.; Erduvan, F. Repdigits in the base b as sums of four balancing numbers. Math. Bohem. 2021, 146, 55–68. [Google Scholar] [CrossRef]
  12. Ddamulira, M. Repdigits as sums of three Padovan numbers. Bol. Soc. Mat. Mex. 2020, 26, 247–261. [Google Scholar] [CrossRef] [PubMed]
  13. Tribonacci Number. Available online: https://mathworld.wolfram.com/TribonacciNumber.html (accessed on 10 September 2022).
  14. Spickerman, W.R.; Joyner, R.N. Binet’s formula for the recursive sequence of order k. Fibonacci Quart. 1984, 22, 327–331. [Google Scholar]
  15. Bombieri, E.; Gubler, W. Heights in Diophantine Geometry; Cambridge University Press: Cambridge, UK, 2006. [Google Scholar]
  16. Matveev, E.M. An explicit lower bound for ahomogeneous rational linear form in logarithms of algebraic numbers. II. Izv. Math. 2000, 64, 1217–1269. [Google Scholar] [CrossRef]
  17. Bugeaud, Y.; Mignotte, M.; Siksek, S. Classical and modular approaches to exponential Diophantine equations I. Fibonacci and Lucas perfect powers. Ann. Math. 2006, 163, 969–1018. [Google Scholar] [CrossRef] [Green Version]
  18. de Weger, B.M.M. Algorithms for Diophantine Equations; Centrum voor Wiskunde en Informatica: Amsterdam, The Netherlands, 1989; pp. 38–40. [Google Scholar]
  19. Le, M.H. An Inequality Concerning the Gel’fond-Baker Method. J. Guangdong Educ. Inst. 2008, 28, 11–12. [Google Scholar]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, Y.; Yang, P.; Zhang, S.; Zhang, K. Repdigits as Sums of Four Tribonacci Numbers. Symmetry 2022, 14, 1931. https://doi.org/10.3390/sym14091931

AMA Style

Zhou Y, Yang P, Zhang S, Zhang K. Repdigits as Sums of Four Tribonacci Numbers. Symmetry. 2022; 14(9):1931. https://doi.org/10.3390/sym14091931

Chicago/Turabian Style

Zhou, Yuetong, Peng Yang, Shaonan Zhang, and Kaiqiang Zhang. 2022. "Repdigits as Sums of Four Tribonacci Numbers" Symmetry 14, no. 9: 1931. https://doi.org/10.3390/sym14091931

APA Style

Zhou, Y., Yang, P., Zhang, S., & Zhang, K. (2022). Repdigits as Sums of Four Tribonacci Numbers. Symmetry, 14(9), 1931. https://doi.org/10.3390/sym14091931

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop