• 検索結果がありません。

[email protected] DepartmentofMathematicsKasetsartUniversityBangkok10900Thailand VichianLaohakosol [email protected] GraduateSchoolofScienceandTechnologyHirosakiUniversityHirosaki,036-8561Japan TakaoKomatsu s OntheSumofReciprocalsofNumbersSatisfy

N/A
N/A
Protected

Academic year: 2022

シェア "[email protected] DepartmentofMathematicsKasetsartUniversityBangkok10900Thailand VichianLaohakosol [email protected] GraduateSchoolofScienceandTechnologyHirosakiUniversityHirosaki,036-8561Japan TakaoKomatsu s OntheSumofReciprocalsofNumbersSatisfy"

Copied!
9
0
0

読み込み中.... (全文を見る)

全文

(1)

23 11

Article 10.5.8

Journal of Integer Sequences, Vol. 13 (2010),

2 3 6 1

47

On the Sum of Reciprocals of Numbers Satisfying a Recurrence Relation of Order s

Takao Komatsu

1

Graduate School of Science and Technology Hirosaki University

Hirosaki, 036-8561 Japan

[email protected]

Vichian Laohakosol

2

Department of Mathematics

Kasetsart University Bangkok 10900

Thailand [email protected]

Abstract We discuss the partial infinite sum P

k=nu−s

k for some positive integern, where uk satisfies a recurrence relation of orders,un=aun−1+un−2+· · ·+un−s(n≥s), with initial valuesu0≥0,uk∈N(0≤k≤s−1), where aand s(≥2) are positive integers.

Ifa= 1, s= 2, andu0 = 0, u1 = 1, then uk =Fk is the k-th Fibonacci number. Our results include some extensions of Ohtsuka and Nakamura. We also consider continued fraction expansions that include such infinite sums.

1The first author was supported in part by the Grant-in-Aid for Scientific Research (C) (No. 22540005), the Japan Society for the Promotion of Science.

2 The second author was supported in part by the Commission on Higher Education and the Thailand Research Fund RTA5180005.

(2)

1 Introduction

The so-calledFibonacci zeta function is defined by ζF(s) =

X

n=1

1 Fns,

where Fn is the n-th Fibonacci number satisfying the recurrence formula Fn =Fn−1+Fn−2 (n ≥2), F0 = 0, F1 = 1.

Ohtsuka and Nakamura [8] studied the partial infinite sums of reciprocal Fibonacci num- bers P

k≥n1/Fks. They gave an explicit formula for the integer part of (P

k≥nFk−1)1 and (P

k≥nFk2)−1. Holliday and Komatsu [2] generalized these results to the cases of Gn and Hn, satisfyingGn =aGn−1+Gn−2 (n≥2) with G0 = 0 and G1 = 1, andHn =Hn−1+Hn−2

(n ≥ 2) with H0 = c and H1 = 1, where a ≥ 1 and c ≥ 0 are integers. In this paper we shall not consider the integer part, but the nearest integer function of (P

k≥nuk1)−1, where {uk}k≥0 is a sequence of non-negative integers satisfying a linear recurrence formula of the type

un=aun−1+un−2+· · ·+un−s,

where a and s (≥ 2) are positive integers. Here, k · k denotes the nearest integer3, namely, kxk=⌊x+ 1/2⌋. Our main result is the following:

Theorem 1. Let {un}n≥0 be an integer sequence satisfying the recurrence formula

un =aun−1+un−2+· · ·+un−s (n≥s) (1) with initial conditions

u0 ≥0, uk ∈N (0≤k ≤s−1), (2)

where a and s(≥2) are positive integers. Then there is a positive integer n0 such that

X

k=n

1 uk

!−1

=un−un−1 (n≥n0). (3)

If a = 1, u0 = 0, u1 = u2 = 1, u3 = 2, · · ·, us−1 = 2s−3, then the un’s are generalized Fibonacci numbers (sometimes called “Fibonaccis-step numbers” [1]). Ifs = 2, thenuk =Fk

are Fibonacci numbers, while if s= 3, then uk =Tk are Tribonacci numbers.

We need the following two lemmas in order to prove this theorem.

Lemma 2. Let a, s∈N, s≥2 and let

f(x) = xs−axs−1−xs−2− · · · −x−1.

Then

3In other contexts, this notation is sometimes used for the distance from the nearest integer.

(3)

(a) f(x) has exactly one positive simple zero α∈R with a < α < a+ 1;

(b) the remaining s−1 zeros of f(x) lie within the unit circle in the complex plane.

Proof. The case where a= 1 can be found in [7], so we assume from now on that a≥2.

By Descarte’s rule of signs, we see thatf(x) has at most one positive real zero. Since f(a)<0 and f(a+ 1)>0, its unique positive real zero, sayα, satisfies (2≤) a < α < a+ 1.

Since multiple roots are counted separately by Descarte’s rule again, part (a) is proved.

Observe from part (a) that

for real x > α,we have f(x)>0, (4) while for real 0< x < α,we have f(x)<0. (5) Next, let

g(x) = (x−1)f(x) =xs+1−(a+ 1)xs+ (a−1)xs−1+ 1.

Observe further that

for real x > α,we have g(x)>0, (6) while for real 1< x < α,we have g(x)<0. (7) To prove part (b), we proceed by establishing several claims.

Claim 1. f(x) has no complex zero z1 with |z1|> α.

Proof of Claim 1. If 0 =f(z1) = z1s−azs−11 −z1s−2− · · · −z1−1, then

|z1|s≤a|z1|s−1+|z1|s−2+· · ·+|z1|+ 1, which implies that f(|z1|)≤0, contradicting (4).

Claim 2. f(x) has no complex zero z2 with 1<|z2|< α.

Proof of Claim 2. Iff(z2) = 0, then 0 =g(z2) = z2s+1−(a+ 1)z2s+ (a−1)z2s−1+ 1 and so (a+ 1)|z2|s ≤ |z2|s+1+ (a−1)|z2|s−1+ 1,

i.e., g(|z2|)≥0, contradicting (7).

Claim 3. f(x) has no complex zero z3 6=α, with either |z3|=α or|z3|= 1.

Proof of Claim 3. Iff(z3) = 0, then

0 =g(z3) = z3s+1−(a+ 1)z3s+ (a−1)z3s−1+ 1 (8) so that

(a+ 1)|z3|s ≤ |z3|s+1+ (a−1)|z3|s−1+ 1. (9)

(4)

If |z3| = α or |z3| = 1, then g(|z3|) = 0 and so (9) must be an equality. Then the two conditions z3s+1 and z3s−1 ∈ R, or z3s+1 = −(a−1)z3s−1 follow from two applications of the fact that

|z1+z2|=|z1|+|z2| ⇐⇒ z1

z2

∈R≥0 (forz2 6= 0) and from

z1+z2 ∈R∧ z1

z2

∈R

=⇒ (z1, z2 ∈R∨z1 =−z2).

• If z3s+1 and z3s−1 ∈ R, then (8) shows that z3s ∈ R, which in turn forces z3 ∈ R. Thus, z3 = ± α or z3 = ± 1. The possibility z3 = α is ruled out by the hypothesis, and the possibilityz3 = 1 is ruled out by (5). To rule out the remaining two possibilities of negative zeros, consider

g(−x) =

(−xs+1−(a+ 1)xs−(a−1)xs−1+ 1, if s is even;

xs+1+ (a+ 1)xs+ (a−1)xs−1+ 1, if s is odd. (10) By Descarte’s rule of signs applied to g(−x), we deduce that g(x) and so also f(x), has at most one real negative zero if s is even and has no real negative zeros if s is odd. When s is even, since f(0) = −1, f(−1) = a > 0, should f(x) have a real negative zero, such zero must lie in the interval (−1,0) and so can neither be −α nor −1.

• Ifz3s+1 =−(a−1)z3s−1, then (8) gives z3s = a+11 . Thus, either 2s≤as <|α|s=|z3|s = 1

a+ 1 ≤ 1

3 or 1 =|z3|s = 1

a+ 1 ≤ 1 3. Both possibilities are untenable and Claim 3 is proved.

Part (b) now follows from Claims 1–3.

We shall keep the notation of Lemma 2throughout the rest of the paper.

Lemma 3. Lets≥2and let{un}n≥0 be an integer sequence satisfying the recurrence formula (1) and the initial conditions (2). Then there are real numbers c >0, d >1, andα > a such that

un=cαn+O(d−n) (n → ∞). (11)

Proof. Let α1 = α, α2, . . . , αt with |αj| <1 (2 ≤ t ≤ s) be the distinct roots of f(x), and letrj for j = 2,3, . . . , t denote the multiplicity of the root αj. Then the expansion formula (11) follows from the shape of un, which is given by

un=cαn+

t

X

j=2

Pj(n)αnJ, where

Pj(x)∈R[x], degPj =rJ −1, 1 +r2+r3+· · ·+rt=s .

(5)

Proof of Theorem 1. Applying Lemma3 and the expansion formula 1

1±ǫ = 1∓ǫ+O(ǫ2) = 1 +O(ǫ) (ǫ→0), we have

1 uk

= 1

k+O(d−k) = 1

k(1 +O((αd)−k)))

= 1

k(1 +O((αd)−k))) = 1

k +O((α2d)−k), Since

X

k=n

1 uk

= 1 c

X

k=n

1 αk +O

X

k=n

2d)−k

!

= α

c(α−1)α−n+O((α2d)−n), we obtain

X

k=n

1 uk

!−1

= α−1

α cαn+O(d−n)

=un−un−1+O(d−n).

Theorem 1 follows by choosing n ≥n0 sufficient large so that the modulus of the last error term becomes less than 1/2.

2 Related results

The following results are similarly obtained. Here,n1, n2,n3,n4 andn5 are positive integers depending only on a.

(6)

Theorem 4.

X

k=n

1 u2k

!−1

=u2n−u2n−2 (n≥n1). (12)

X

k=n

1 u2k−1

!−1

=u2n−1−u2n−3 (n ≥n2). (13)

X

k=n

(−1)k uk

!1

= (−1)n(un+un−1) (n ≥n3). (14)

X

k=n

(−1)k u2k

!−1

= (−1)n(u2n+u2n−2) (n≥n4). (15)

X

k=n

(−1)k u2k−1

!−1

= (−1)n(u2n−1+u2n−3) (n ≥n5). (16)

Proof. We shall prove only (14). The other identities are proved similarly. By (11) we get

X

k=n

(−1)k uk

=

X

k=n

(−1)kk+O(d−k)

=

X

k=n

(−1)k

k 1 +O((αd)−k)

= α

c(−α)n(α+ 1) +O (−α2d)−n .

By taking its reciprocal, we have

X

k=n

(−1)k uk

!1

= c(−α)n(α+ 1)

α 1 +O((αd)−n)

= (−1)n(cαn+cαn−1) +O((−d)−n)

= (−1)n(un+un−1) +O(d−n).

The identity (14) follows by choosing n ≥ n3 sufficiently large so that the modulus of the last error term becomes less than 1/2.

3 The sum of reciprocal Tribonacci numbers

The so-called Tribonacci numbers Tn ([6, Ch. 46], [9, sequence A000073], [3]) are defined by Tn=Tn−1+Tn−2+Tn−3 (n≥3), T0 = 0, T1 =T2 = 1.

(7)

By setting a = 1, s = 3 and uk = Tk (k ≥ 0) in Theorem 1 and Theorem 4, we get some identities about the partial Tribonacci zeta functions. Numerical evidences imply that identities hold for smaller positive integers n, as indicated in the identities. The detailed explanations for small n can be seen in [5].

Corollary 5.

X

k=n

1 Tk

!−1

=Tn−Tn−1 (n≥1). (17)

X

k=n

1 T2k

!−1

=T2n−T2n−2 (n≥1). (18)

X

k=n

1 T2k−1

!−1

=T2n−1−T2n−3 (n ≥2). (19)

X

k=n

(−1)k Tk

!−1

= (−1)n(Tn+Tn−1) (n≥2). (20)

X

k=n

(−1)k T2k

!−1

= (−1)n(T2n+T2n−2) (n≥1). (21)

X

k=n

(−1)k T2k−1

!−1

= (−1)n(T2n−1+T2n−3) (n≥2). (22)

4 Continued fraction expansion of generalized m-step zeta functions

The first author [4] studied several continued fraction expansions of some types of Fibonacci zeta functionsζF(s) :=P

n=1Fn−s and Lucas zeta functions inζL(s) := P

n=1L−sn , whereLn

is the n-th Lucas number defined by

Ln=Ln−1+Ln−2 (n ≥2) L0 = 2, L1 = 1.

A continued fraction expansion of the generalizedm-step zeta functions defined byζu(m)(s) :=

P

n=1u−sn , where

un =aun−1+un−2 +· · ·+un−m (n≥m)

(8)

with initial positive integral values uk (0≤k≤m−1), is given by

ζu(m)(s) = 1

us1− u21s

us1+us2 − u2s2

us2+us3− u23s

us3+us4 − ...− u2n−1s usn−1+usn− · · ·

.

DefineAn(respectively, Bn) as the numerator (respectively, denominator) of thenth conver- gent of the continued fraction expansion given for ζu(m)(s):

An

Bn

= 1

us1− u2s1

us1+us2− u22s

us2 +us3− u23s

us3+us4− ...− u2n−1s usn−1 +usn

.

Hence {Aν}ν≥0 and {Bν}ν≥0 satisfy the following recurrence formulas.

Aν = (usν−1+usν)Aν−1−u2ν−1s Aν−2 (ν ≥2), A0 = 0, A1 = 1;

Bν = (usν−1+usν)Bν−1−u2ν−1s Bν−2 (ν ≥2), B0 = 1, B1 =us1

In fact,Aν and Bν can be expressed explicitly as follows.

Lemma 6. For n = 1,2, . . .

An = (u1u2· · ·un)s

n

X

ν=1

1

usν, Bn= (u1u2· · ·un)s.

Proof. By induction we have Bn= (u1u2· · ·un)s. Thus, An=Bn

n

X

ν=1

1

usν = (u1u2· · ·un)s

n

X

ν=1

1 usν .

Theorem 1 provides us with interesting information about the nearest integer of the reciprocal ofζu(m)(s)−An/Bn.

(9)

References

[1] Eric Weisstein, World of Mathematics, published electronically at http://mathworld.wolfram.com/Fibonaccin-StepNumber.html.

[2] S. H. Holliday and T. Komatsu, On the sum of reciprocal generalized Fibonacci numbers, submitted.

[3] E. Kili¸c, Tribonacci sequences with certain indices and their sums, Ars Combin. 86 (2008), 13–22.

[4] T. Komatsu, On continued fraction expansions of Fibonacci and Lucas Dirichlet series, Fibonacci Quart. 46/47 (2008/2009), 268–278.

[5] T. Komatsu, On the sum of reciprocal Tribonacci numbers,Ars Combin., to appear.

[6] T. Koshy,Fibonacci and Lucas Numbers with Applications, John Wiley & Sons, 2001.

[7] M. D. Miller, On generalized Fibonacci numbers, Amer. Math. Monthly 78 (1971), 1108–1109.

[8] H. Ohtsuka and S. Nakamura, On the sum of reciprocal Fibonacci numbers,Fibonacci Quart. 46/47 (2008/2009), 153–159.

[9] N. J. A. Sloane, (2010), The On-Line Encyclopedia of Integer Sequences, published electronically at http://www.research.att.com/~njas/sequences/.

2010 Mathematics Subject Classification: Primary 11A55; Secondary 11B39.

Keywords: Fibonacci numbers, recurrence relations of s-th order, partial infinite sum.

(Concerned with sequenceA000073.)

Received January 20 2010; revised version received May 19 2010. Published in Journal of Integer Sequences, May 20 2010.

Return to Journal of Integer Sequences home page.

参照

関連したドキュメント

van Lint raised the problem whether the number 120 is the unique (positive) integer n for which the set { 1, 3, 8, 120 } constitutes a solution for Diophantus’ problem above; in

Zeleke, On proofs of certain combinatorial identities, to appear in Proceedings of the 11th International Conference on Fibonacci Numbers and Their Applica- tions.

Zeleke, On proofs of certain combinatorial identities, to appear in Proceedings of the 11th International Conference on Fibonacci Numbers and Their Applica- tions..

The method employed to prove indecomposability of the elements of the Martin boundary of the Young lattice can not be applied to Young-Fibonacci lattice, since the K 0 -functor ring

Leonhard Euler (1707−1783) named the equation after John Pell by mis- take, studied the infinite continued fractions and proved that a finally periodic continued fraction describes

Some algorithmic problems for holonomic distributions 16:10 – 17:00 Hikosaburo Komatsu (Univ. of Tokyo). History of mathematics of the world due

Since any density may be approximated in L 1 ( R ) by densities with finite total variation, our approach is no less general than the Fourier method. Section 3 contains some

Next, we present a slightly different continued fraction expansion based on the Euclidean algorithm to the nearest integer, namely the centred continued fraction ex- pansion. The