16 6
16 6
16 6
VOLUME 16
laBr
NUMBER 6
CONTENTS
On the Density of the Image Sets of Certain
Arithmetic FunctionsIII. . . . .
.
.Rosalind
Guaraldo 481
Evaluation of Sums of Convolved Powers Using
Stirling and Eulerian Numbers.
H. W. Gould 488
b-kdlc Numbers in Pascalfs Triangle Modulo b. ......Heiko Harborth 497
Divisibility Properties of Polynomials in Pascal's
Triangle*..... V. E. Hoggatt3 Jr.3 and Marjorie Bieknell-Johnson
501
The Golden Section in the Earliest Notated
Western Music
.
Paul Larson 513
On Fibonacci Numbers Which Are Powers. .
Neville
Bobbins 515
Primes, Powers, and Partitions.
B. de La Rosa 518
On Odd Perfect Numbers.
G. L. Cohen 523
A Simple Continued Fraction Represents a Mediant
Nest of Intervals.
Irving Adler 527
Folded Sequences and Bode!s Problem
....17. E. Greig 530
Fibonacci Numbers in Coin Tossing
Sequences.
.Mark Finkelstein
and Robert Whitley 539
Strong Divisibility Sequences with Nonzero
Initial Term
Clark Kimberling
541
Minimum Periods Modulo n for Bernoulli Numbers.....
W. Herget 544
The Rank-Vector of a Partition. .
Hansraj Gupta 548
The Andrews Formula for Fibonacci Numbers.
.....Hansraj
Gupta 552
Operational Formulas for Unusual Fibonacci Series.....H. W. Gould 555
A Figurate Number Curiosity: Every Integer is a
Quadratic Function of a Figurate Number...
Harvey J. Hindin 561
Elementary Problems and Solutions.........Edited by A. P. Hillman 562
Advanced Problems and Solutions
Edited by Raymond E. Whitney 566
DECEMBER
1978
Volume Index
570
^ e Fibonacci Quarterly
THE OFFICIAL JOURNAL OF THE FIBONACCI ASSOCIATION
DEVOTED TO THE STUDY
OF INTEGERS WITH SPECIAL
PROPERTIES
EDITOR
Verner E. Eoggatt3
Jr.
CO-EDITOR
Gerald E. Bergum
EDITORIAL BOARD
H. L. Alder
Marjorie Bicknell-Johnson
Paul F. Byrd
L. Carlitz
H. W. Gould
A. P. Hillman
David A. Klarner
Leonard Klosinski
Donald E. Knuth
C. T. Long
M. N. S. Swamy
D. E. Thoro
Maxey Brooke
Bro. A. Brousseau
Calvin D. Crabill
T. A. Davis
A. F. Horadam
Dov Jarden
L. H. Lange
James Maxwell
Sister M. DeSales
McNabb
John Mitchem
D. W. Robinson
Lloyd Walker
Charles H. Wall
The California Mathematics Council
ROSALIND GUARALDO
College,
Brooklyn,
Francis
1.
NY 11201
INTRODUCTION
S(w) = ^2f(dd9
j),
j-i
^(n) = n +
Sin),
p)
^P = {a:|a: = Tin)
C= {x\x
for any
T{n)
n).
/ ( 0 , j ) = 0, j = 1 , 2 , . . .
fid,
j ) = O(Q\) u n i f o r m l y i n j , i . e . ,
sup {fid,
j), 0 < d < j } = o ( j ! )
t h e n t h e d e n s i t y of <R e x i s t s .
VKOOfa: We f i r s t show t h a t
(2.2)
Didkl,
dk\
+ P ) = ( * ) , 0 <_r k\
- 1.
x = dk\ + Y^ djjl
J
Clearly, T(x)
= T(y)
and
k-l
y = ^! + ]C ^P'!
=i
j= i
if and only if
481
If
482
= dk_
= = dk_
= dk_
[Dec.
= = dk_
= 0).
(d + l)kl
D(dkl9
- 1) = D(k\
- 1).
fe_+i
and ^ Z?(dfe! - 1 ) .
+ 1 ) ! - 1)
c? = 0
Lojnmci 2.4:
There exists an integer k0 such that for all k >_ k0 the sets
ft(0, k\ - 1) ,fi(/c!, 2k\ - 1) , . . . ,fi(fefe!, (k + 1) ! - 1) are pairwise disjoint,
except possibly for adjacent pairs.
Vsioofi: The maximum value in Q,(dkl , (d + l)kl
- 1) is at most (d + )k\ 1 + kMk, where Mk- max {f(d9 j) , 1 . J ,fc},and the minimum value in Q((d +
2)kl , (6? + 3)fc! - 1) is at least (d + 2)/c!. By assumption (b) , there exists
krQ such that f(d9 j) < j!/2, for all j >L ?c0f, and there exists kQ _> kr0 such
that f(d> j) < jl/2
- k0'Mfc(r , for all j >. k0 , where Mfer = max {f(d9
j) | 1 <_ J <.
/CQ}. Therefore, if k >_k09 we have
E^>
+ E
J
= k
+ 1
j=k
J !/2 <
'
^!-
+l
, (d + 2)fc! - l ) | , 0 <_
W + Dkl - 1) - e,
J=o
where ^ depends on the number of elements that the sets 7(0, Zc! - 1 ) , Q(kl ,
2fc! - 1 ) , ..., ft(/c&!, (A: + 1) ! - l) have in common, we obtain
(2.5)
Z?((fe + 1 ) ! - 1) = (k + DD(kl
fc-i
- 1) - ] T X d . j k .
i =o
fc-i
Let 4 k = (&! - l ) / f c ! and Ek = ] T A d > k /(fe + 1) ! , k >_kQ.
d=o
Therefore,
Then 2 . 5 becomes
1978]
k+l
Ak
~
A
k-1
483
~ek
~
=
~ek-l
- Av = - e k
k
so Ak + 1 - Ak^ = -J2eJ>
obtain
k-1
i
- e * ' ^k + i
= A
k0 ~ J2
j = k0
J'
R e p l a c i n g k + 1 by k, we
j=k0
^ =i, o - > , .
(2.6)
J =
' "
fc-i
Clearly, 1/&! 4fe 1 and ]T E^ = Ak^ - Ak <_ Ak <_ 1.
Thus, ]P e^ is a
L = AK - ^ e j .
Ak = L + ]T j.,fc^ k Q .
^=L+o(l).
Multiplying both sides of this equation by kl and using the definition of the
Aki we obtain
(2.9)
Using 2.3, 2.4, 2.9, and the definition of the A!s and the e's, we have
D dkl
d-l
- 1) = X ^ ^ ! , ^ ^
kl
d-l
" D " XA ^
(3=0
O= 0
d-l
= ^
(Lfc! + o{k\))
+ o((fc +
l)lek)
e=0
D(dk\ - 1) = dklL
+ o((fc + 1)!).
- 1) +
j-o
D(dkkl,
dkkl +dk_i(k
- 1) I + ) - , wherefiis the number of elements that
the sets Q(0, dkkl -1) and Q,(dkkl , dkkl + ^ ^ ( k - l ) ! + ) have in common.
Hence, if n is sufficiently large, then, by using 2.2, 2.10, and the definition of the A's, we obtain
D(n) = dkk\L
+ D{dk_1(k-l)l
+ o{(k + l)l)
o((k+l)\).
k8k
[Dec.
- 1)! +
o(kl)
- 1)!L + o{(k
+ 1)!).
o(k\)9
D(n)/n
2.12:
is not
" ^ e j = Ak ~ Z ec < f o r a 1 1 ^ ^ o w n e r e ^ o i s
Ja
fe
J = ^o
defined as in Lemma 2.4. Therefore, L <_ Ak if k >_ k . if TOr) = T{y) 9 x + y,
and & is such that k >_ k0 and x <_ kl - 1, y <_ kl - 1; then, since
4* = D(kl - l)/fc!,
Vtioofa:
We have L = Ak^
it follows that L <_ Ak <_ 1. If T is one-one, then it follows from the definition of the A1 s and the fs that Ak = 1 and e^ = 0 for all ?C, so L = 1.
It seems to be true, although possibly difficult to prove, that L < 1 if
each f(d, j) = /(<i) depends only on d and / satisfies the hypotheses of Theoren 2.1. It also seems to be the case that we should always have L > 0 under these hypotheses; this result again will be left to conjecture.
3.
j)
= 0(j!/j2 log2j)
The main drawback to Theorem 2.1 is the condition f(09 j) = 0. If we assume that f(d9 j) = 0(j!) uniformly in j for all "digits" d9 it seems to be
difficult to find a workable relationship between the quantities Ak9 but on
the other hand, it also seems to be difficult to find an example of an image
set <R which does not have density under this assumption. However, we do have
the following result.
JhzoJiQJtn 3.1:
of Q exists.
If f(d9
Vnoo{:
log 2 j) = 0(k\/k2
log2k).
j-i
k
9
then S(n) = ^ 0 ( j ! / j 2
j-i
1978]
log 2 (^ - 1))
fc-i
Z?(fe! - 1) = D(0) + ^Z?(q!, (q + 1)! - l) + 0(kl/k2
log2k).
(d + l)k\
- 1) = D(dkl,
log2k)
(d + l)kl
- 1) = D(dkl,
+ - + D(dkl + (fc - 1 ) ! ,
2
e r r o r s of s i z e 0(kl/k
k_
i.e.,
(3.4)
!-
(d+l)k\-l
2
log k),
+
Z?(d/c!, W + l)k\ - 1 ) = (dfc!, dfc!) + ^D{dk\
4! > ^ !
+ (q + 1)! - 1 ) % 0(fc!/fc l o g 2 k ) .
(q + 1)! - 1 ) ,
1 1 q <. k - 1
fc-i
0(dfe! , (d + 1) ! - 1) = 0(0) + X (<7! , (? + 1)! - 1)
+ 0(fc!/fc log2k)
D(dkl9
486
[Dec.
2k\ - 1) + D(2kl,
2
+ 0(k\/k
= 2D(k\ - 1) +D{2k\,
(k + 1)! - 1)
+ 0(k\/k
log k)
log2k)
log2k)
(k + 1)! - 1) + 0(k\/k
log2k).
- 1) + k + 1,
log k);
thus,
D((k + 1)1 - 1) = (k + l)(fc! - 1) + 0(ft + l)!/fe log 2 k).
(3.7)
+1
k +i = o + E
^ '
lo
g^')k
using the above equation, we may conclude that there exists a constant L such
that
(3.8)
Ak = L + 0(l/log fc).
Now l e t n = 2, dk.k\.
j
Then
= 1
+A - i ^ _ i
) + 0(fc!/*4 l o g 2 ^ ) .
By the same type of reasoning employed to get 3.4 and 3.7, we see that
D{dkmk\m - 1) = dKD(kl
Since dk
log2km)
+ D{dKk\m,
dknk\m+
Therefore,
- 1) + 0(kl/k
Z?() = ^ ^ ( L
/m-1
"1
'
\J=X
+ 0 ( l / l o g km))
'
'
+ 0{k\Jkm
log 2 /c m ) + D( Y.dk.kl)
= nL + 0(fc!/fcm l o g 2 ^ ) + o ( j ! / l o g
J =I
Hence,
>(n) = nL + Of/cl/log kj
so
Z?(n)/n = L + 0 ( l / l o g fcm) = L + 0 ( 1 ) ,
j).
1978]
487
RomoJtk 1: Theorem 3.1 has the drawback that the computability of the
density has been lost.
ROMCUik 2: If we assume that f(d9 j) = o(j\)
uniformly in j , then there
exists an image set < which does not have density. For example, let f(d9 j)
= 0 when j is even and f(d, j) = j! when j is odd. Then,
/
\zi
k-l
fe.i
if /c is odd, and
kl +
k-l
= kl +
X^ ! J
/
J=I
if fc is even. Therefore, the number of integers between k\ and 2kl that belong to ^ if fe is odd is at most 1 + (k - 2)! + (k - 4)! + + 1, and the
number of integers between k\ and 2k\ that belong to 9 if k is even is at
k\ - (k - 1)! - (k - 3)! - - 1!. Hence, if we let 6 and A denote the lower
and upper density of <, respectively, we see that
6 <. 0 + 0(1) and A >. 1 + 0(1),
so 6 = 0 and A = 1.
It is also interesting to note that, if we let f(d9 j) = o(jl)
uniformly
in j , there do exist image sets ^ of density 0. For example, if f(d, j) = 0
when d f 1 or j = 1 and f(d9 j) = 2j! if d =. 1 and j > 1, then no member of
(except 1) has the "digit" 1 anywhere in its factorial representation, and
the set
(3.8)
| n\n = J^ djJl>
dd + 19 1 <3
<k\
is easily seen to be the set of density 0.
Our next result is an immediate corollary of Theorem 3.1.
CotiolloAbj
3.9:
fid)
I f f(d9
= 0(j!/j
j)
= f(d)
log j)
488
[Dec.
1
2
3
4
5
6
7
8
9
10
0
0
0
2
6
8
14
17
26
39
fc
0
0
0
0.066667
0.041667
0.008929
0.002401
0.000375
0.000064
0.000009
REFERENCES
1.
2.
3.
G. Faber, "Uber die Adzahlberkei der Rationale)! Zahlen," Math, Ann., Vol.
60 (1905), pp. 196-203.
Rosalind Guaraldo, "On the Density of the Image Sets of Certain Arithmetic FunctionsI," The Fibonacci
Quarterly,
Vol. 16, No. 4 (Aug. 1978),
pp. 319-326.
Rosalind Guaraldo, "On the Density of the Image Sets of Certain Arithmetic FunctionsII," The Fibonacci
Quarterly,
Vol. 16, No. 5 (Oct. 1978),
pp. 428-434.
*****
H. W. GOULD
University,
Morgantown,
W. Va.
26506
ABSTRACT
INTRODUCTION
Neuman and Schonbach [9] have obtained a formula for the series of convolved powers
n
(1.1)
S(i,j;n)
= fcf(n - k)s'
fc = o
1978]
489
G(t;i,j)
Y,
A GENERATING FUNCTION
tnS(i,j;n).
n=0
Then
G(t;i,j)
k =0
n=*k
k=0
n=0
(2.2)
= kHk ^ n -
G(t;i,j)
k=0
n=0
2>v<
k
may be summed in a variety of ways. We shall use the methods of (i) Stirling
numbers of the second kind and (ii) Eulerian numbers. Our (2.2) is (3.4) in
Carlitz [3].
3.
(3.1)
(tDff(t) = X S(p9k)tkDkf(t),
k=0
(3.2)
k!S(p,k)
Expli-
= A V = E<-D*- J "0)fc P .
j-o
The formula dates back more than 150 years, but, for a recent example, see
Riordan [10, p. 45, ex. 18]. Riordan gives a full account of the properties
of Stirling numbers of both first and second kinds. Other historical remarks
and variant notations are discussed in [6], Applying the formula is easy because (tD)ptk
= kptk,
whence we have
490
(3.3)
^ V = V k\S(p,k)
[Dec.
^ .
This, too, is a very old formula. It converges for \t\ < 1, but we treat it
as a formal power series. Carlitz [2] gives a good discussion of formal power
series techniques.
Using (3.3) in (2.2), we find
i+j
(3.4)
, +2 Z ^ ! ( p -
=E "
G(t;iJ)
(1 - t)
r=0
k)lS(i,k)S(j,r-k).
k=0
(3.5)
5,(i,j) =2/c!(r-/c)!5(i,/c)5(j,2-/c).
k=o
i> = 0
n = 0
=Ei(i)-.'''
v = 0 n=r
n=0
r=0
In the next-to-last step here, the upper limit v = i, + j might as well have
been v = because of zero terms involved, since S(p,k)
= 0 when k > p. This
makes manipulation easier. Equating coefficients of tn and dropping some zero
terms, we find finally then our desired formula
(3.6)
S(iJ;n)
=J2[Hl)sr(i,j).
r =0
This simple expression may be compared with the bulky form of expression given
in [9] using Bernoulli numbers.
Having found our desired formula, we can next offer a much quicker proof.
Recall [10, p. 33] that
n
(3.7)
x" -
(X\v\S(n,r).
V =0
ky =Y(
?o
rlS(i,r)i2slS(j,s)(l)(n-k),
ZTo
\i /\ s /
S(i,j;n)
^riSfi.rl^lsy.s)!^!
r=Q
s=0
1978]
491
For the sake of completeness, we recall [10, p. 48] some of the values of
S(n,k):
0
0
1
2
3
4
5
6
7
1
1
1
1
1
1
1
1
3
7
15
31
63
4_
3_
1
15
140
1
21
7 -U
k_
1
1
6
25
90
301
1
10
65
350
n
Here, S(n,k)
= 0 when k > n and S(n,0)
= 0 for n >. 1.
For j = 0, formula (3.6) becomes the well known
1
(4.1)
S(i
>
'
(n+21),
5(2,0;) = ( n 2 1 )
+ 2
("31)'
- (n+21)
14(^)
Se^1)
2A(^1).
5(2,1;) = {
) + 2^
J=
12
>
3?r
5nd + 2n
60
^^^^crj^ftVKn )-^; )
2n6 - 5 ^ + 3n^
60
492
[Dec.
i =2
i =3
i
i
i
For j =
sa,3;n)-(n+31)+6(nl1)
6(n+51)t
(5.D
Afi
^ ( - D T ^ C J
k)K
k=o
These must not be confused with Euler numbers appearing in the power series
expansion of the secant function. The Eulerian numbers satisfy
A n>3 = An,n-j
and
An,
= jAn.1}j
+ i9
+ (n - j + D ^ - i , j - i ,
n
J = l
1_
3_
4_
5_
1
120
J_
1
1
1
1
1
1
1
1
1
4
1
1
11
11
26
66
26
1
57 302 302
57
120 1191 2416 1191
1978]
493
these numbers and tracing back to Euler. For our purposes, we need the wellknown expansion
(5.2)
x> K ^= a -
tr-lyEtkAntk.
k=0
k=0
This expansion is known to be valid for \t\ < 1, but again we treat all series
here as formal power series since we do not use the sums of any infinite series. We never assign t a value, but equate coefficients only.
Applying this to (2.2), we find
.
G(t;i,j)
s=0
r=0
i + 3
r=0N
n=0
A .
r=0
s=o
"
's=0
A .
i + j
.3)
w^-zfi^nE
A-
A
j,
r-8
Here we have again dropped some of the terms that are zero by noting that
An = 0 whenever j > n.
Formula (5.3) is (3.6) in Carlitz [3].
As with our previous Stirling number argument, we could obtain (5.3) by
another method. We recall that in fact
^n=i2(x+i~i)^,i
(5.4)
J=0
and form the product kz (n - k)J and sum from k = 0 to k = n to obtain a formula for (5.3) analogous to (3.8). We omit the details.
6.
(6.1)
S(i,0;n)
>l.
r =l
\
r= 0
i+1
)Ai,r,
By (5.3),
s ^0, v - s
s i n c e A0jr_s
= 0 f or r f
s,
k3k
iL \
)Ai,r
"+i
[Dec.
= 0 for i >_ 1,
r= 1
= Y^ Q + iV^-^+i' by P u t t i n S ^ - ^ + 1 for p,
r= 1
t
( '+i)Ai,r9
I = 1
5(1,0;n) = j ^ 1 ) ,
S(2,0;n)
=("^)
(^
),
W . O ^ - C ^ + l l ^ + ^ + ll^^J + f'J4),
etc.
For j = 1, we find
5(
2,i;)-(M;1)
s(3,i;n).(B;i)
and so on.
7.
+ (n:2)
*f,;2) + (n;3)
These again are a different way of saying what was found in [9],
ALTERNATIVE EXPRESSION OF THE STIRLING NUMBER EXPANSION
( n +1\
Formula (3.6) uses the values of I
-. ). We wish to show now that we can
transform this result easily into a formula using just f . i)> i-e? directly as a series of binomial coefficients in n rather than n + 1. We will need
to recall, see [10], the recurrence relation for Stirling numbers of the second kind
(7.1)
S(m,k) = kS(m-l,k)
S(m-l,k-l).
S(j+l,r-k)
= (r-k)S(j9r-k)
We get
+ S(j ,r - k - 1).
Now, by (3.6) and the usual recurrence for binomial coefficients, we have
1978]
i+
'
^95
D(")5,(i,;) + U
(Z)Sr-iV>3)
V=0
r=0
However, Sr(i,j)
Sr_1(iij)
V - l
= J^ kl (r -k)\S(i,k)S(j,r
k=o
-k)
k\ (r - 1 - k) \S(i,k)S(j,r
- 1 - k)
=o
fc
v- 1
(r - k)S(j,r
- k) + ^ fc! (p - 1 -fc)\S(i9k)SU,r
k =0
-1-k)
k =0
= X
feKr-l-WllW^W+l^-fe)
-5(i,/c)5'(j,3?-fe-l)},
by (7.2)
k=0
r- 1
+ J^kl(Y>-l
-k)\S(i,k)S(j,r-l
-k)
fe=0
v-1
= f c ! (r-fe-l)!5(t,W5(j+l,r-W + r\S(i,r)S(j
,0).
k=0
The extra term here may be dropped when we consider j >_ 1. . Therefore, we have
the new result that
i+j
(7.3)
S(i,j;n)
r=0
fe
Let J = 1 again.
ExampteA:
=0
We find
(5),
1
1
1
1
1
1
1
2_
1
3
7
15
31
63
3_
4_
5_
6_
2
12
6
50
60
24
180 390 360 120
602 2100 3360 3520
7 -
720
k36
(7.4)
S(i,0;n)
<-!)*(*. + ) E
fc = 0
[Dec.
( )) W + D* ,
j" = 0
(7.5)
S(i,0;n)
=E
n >. 0, i >. 1.
r= 0
This occurs, for example, as the solution to a problem [13] in the .American
Mathemat-ieal
Monthly.
Examples:
S(l,0;n)
= j^1),
5(2,0;n)-.("J1) + 2 ( " + 2 ) ,
^<3,0;) - ("J1) - 6( + 2 ) +
and so forth.
8.
D+ 24 ( W D'
FINAL REMARKS
It is interesting to note that the original sum (1.1) is a type of convolution. So also formulas (3.6), (5.3), and (7.3) involve convolutions of
the Stirling and Eulerian numbers. The formula found in [9] is not of this
type. This is so because of the way in which the binomial theorem was first
used. It would evidently be possible to obtain convolutions of the Bernoulli
numbers. To get such a formula using Bernoulli polynomials is easy. Let us
recall that
(8.D
~=Y,h~B(x)>
(8.2)
x" = ; r ^ T E r 1 H f e ) , n >
k-o
1978
h31
Form the product k1'(n - kY by using this formula to expand kl and (n - k)J'.
Sum both sides and we get
(8-3)
S(iJ;n} - ^t^^^jt^^tBAmAn
r=0
'
s=0
- k),
k = 0
(8.4)
{l)xn""B^
BAx) = X)
m = Q
it would be possible to secure a convolution of the Bernoulli numbers. However, the author has not reduced this to any interesting or useful formula
that appears to offer any advantages over those wev have derived here or those
in [9]. We leave this as a project for the reader.
It is also possible to obtain a mixed formula by proceeding first as in
[9] to get
'
k=0
apply one of our Stirling number expansions to the inner sum and get, e.g.,
J
i+r
(8.5)
+r k)
> >
fc=0
V=0
'
vv
(continues
Vol. 32
^n\
on page 560)
For the binomial coefficients in Pascal's triangle we write their smallest nonnegative residues modulo a base b. Then blocks of consecutive integers
within the rows may be interpreted as Z?-adic numbers. What Z?-adic numbers
can occur in the Pascal triangle modulo bl In this article we will give the
density of such numbers and determine the smallest positive integer h(b) ,
such that its b-adc representation does not occur (see [3] for b = 2).
We use the notation
m
t
Y1 a
i=0
(amam-l
'"
v \
>
' am $ >
numbers
498
L<imma 1:
gle modulo 2.
[Dec.
a)
for i = 0, 1, 2, gives n E k (mod 2 ) , & E 0 (mod 2 ) , and n E 1 (mod 2 ) , respectively, which is a contradiction.
Lomma 2:
(111)^ is not to be found within any row of Pascalfs triangle
modulo b with b > 2.
Psioofi:
We assume that
Together with (1), for i = 0 and = 1, we conclude that n E 2& + 1 (mod 2?),
and w E 2fc + 3 (mod b), respectively. However, both congruences are possible
only if b = 2.
We are now able to determine the density.
ThdOJiQjn 1: Almost all b-adc
Pascal's triangle modulo b.
VKOO^X As noted in [4], it is well known that the density of those -adic
integers not containing a given sequence of digits is 0 (see [2], p. 120).
Thus, the proof is given by Lemmas 1 and 2.
Jk2.0h.Qyn 2: Let h(b) be the smallest b-adc number not being found within any row of Pascal's triangle modulo b. Then, h(b) - b2 + b + 1 = (111)^
for b > 2, and 7z(2) = 11 = (1011) 2 .
We first prove two lemmas.
Lemma 3: Let b = bYb2 with (bl9 b2) = 1. Then (am ... a0)b occurs in the
Pascal triangle modulo b if and only if (aim ... a ^ ) ^ for i = 1, 2 occur in
the triangles modulo bi with a^j E a3- (mod b^) , j = 0, 1, . . . , m.
VJ100&: One direction of the proof is trivial.
, ,
In the following, we use the result of [1] and [6], that (^j (mod b) is
periodic for fixed k with the minimal period N being the product of all prime
powers p a + $ with pa from the canonical factorization of b and 3 from p$
^k
< p+ 1. Thus, N depends only on the prime factors of b and on k (see [5] for
further references). By reasons of symmetry, a corresponding periodicity of
I n _j_ \
length L holds for I 7, , p ) with fixed n and k.
From this and by the assumption, we are able to find n^ and k^ such that
for i = 1, 2,
with minimal periods L^ each being the lowest common multiple of m + 1 minimal
1978]
499
"t" X TJ I
has solutions o^ , # 2 .
N with
lni
fVo
'
vC lyLl p 5
+xiLi+yiNi\
x2J
_
(mod bi),
j = 0, 1,
. . . , m.
= n2 + x2L2
y2N2,
VK.00^',
(2)
(^/T)
( m d pa) f r pa3
fc ^ p a 3 .
rp a g - k
ppQ
y-Z
pi
i2-l
T. -
V%
t = a$ - a +1
(3)
Pj (r) -
TT
(r
+l
"y
ag
1, for v = 1, 2,
Let
i =1
K)= ^TT^-^)
\ ^
Ev
< m o d p)
j -1
pp
a$\
a - 1
( ae ) = TJ =Tl j ( )
a _ 1
^TT
j( )
J = l
TT J(P)
J-.l
^p
a_1+
~1
P T<7T=P 1J ( P )
<mod p a >-
500
Dec. 1978
, includ-
(5)
(,-- 2).
R. D. Fray, "Congruence Properties of Ordinary and ^-Binomial Coefficients," Duke Math, Journal,
Vol. 34 (1967), pp. 467-480.
G. H. Hardy & E. M. Wright, An Introduction
to the Theory of Numbers, 4th
ed. (Oxford: Oxford University Press, 1962).
H. Harborth, "Aufgabe P 424, Dualzahlen im Pascal-Dreieck," Praxis
der
Mathematik,
Vol. 13 (1971), pp. 76-77.
D. Singmaster, written communication to the author.
D. Singmaster, "Divisibility of Binomial and Multinomial Coefficients by
Primes and Prime Powers" (to appear).
W. F. Trench, "On Periodicities of Certain Sequences of Residues," Amer.
Math. Monthly,
Vol. 67 (I960), pp. 652-656.
E. Landau, Vorlesungen
uber Zahlentheorie
(New York: Chelsea, 1950).
###*#
San Jose
State
University,
San Jose,
California
95192
D i v i s i b i l i t y p r o p e r t i e s of t h e F i b o n a c c i sequence {Fn}
i n c l u d i n g t h e p r o p e r t y of g r e a t e s t common d i v i s o r s ,
\^m s ^n '
a r e w e l l known,
~ ^ (m, n) '
Here the derivation of the greatest common divisor of a sequence pair is extended to the Fibonacci polynomials, the Morgan-Voyce polynomials, the Chebyshev polynomials, and more general polynomials from a problem of Schechter [1] .
Moreover, all of these polynomials have coefficients which lie along rising
diagonals of Pascal's triangle, and all of these polynomials satisfy
(um(x),
un(x))
= U(m n x(x) with suitable adjustment of subscripts.
1.
INTRODUCTION
= 1, Bx(x)
Bn(x)
= bn.(x)
bn(x)
= xBn_1{x)
Bn(x)
= Bn.(x)
= x + 2; bQ(x)
= 1, b1(x)
= x + 1,
and
(1.1)
+ (1 +
+
+
x)Bn-!(x),
bn_(x),
bn(x).
= (x + 2)un+1(x)
un+2(x)
un(x),
=
=
=
=
=
=
=
=
1
1
x + 1
x + 2
x2 + 3x
x1 + kx
x3 + 5x2
x3 + 6x2
+ 1
+ 3
+ 6x + 1
+ 10^ 4- 4
Clearly, we see that the coefficients of this double sequence lie along the
rising diagonals of Pascal's triangle.
The Fibonacci polynomials are
(1.3)
f0(x)
= 0, f.ix)
= 1, fn+2{x)
= xfn+1(x)
= 1
= x
= x2 + 1
= x3 + 2x
501
fn(x),
502
= xh
= x5
= xe
= x7
f5(x)
fe(x)
f7(x)
f8(x)
3x2
kx2
5xh
6x5
+
+
+
+
[Dec.
+1
+ 3x
+ 6x2 + 1
+ lO^r3 + kx
Once again, we see that the coefficients lie along the rising diagonals of
Pascal's triangle.
It can be shown that [3], [4]
(1.4)
bn(x2)
= f2n
+1
xBn(x2)
= f2n
+2
(x)
(x),
and the fact that coefficients lie on the rising diagonals of Pascal's triangle follows from that property for the Fibonacci polynomials. The Fibonacci
polynomials obey
/+(*) = (x2 + 2)fn
(1.5)
+ 2(x)
fn(x),
For F i b o n a c c i p o l y n o m i a l s ,
(fm(x)>
Vtioofa:
fnW)
- xfm(x)
= d(x) .
fm + 1(x))
(#)-
= f(m,n)
fm_(x)9
Then, since d(x)\fm(x)
In turn, fm(x)
- xfm_1(x)
implies that
= 1. Therefore, d(x)
= 1,
/p +r(a0 = fP-i(x)fr(x)
+ / p (a0/ r+ 1 ( * ) ,
and
(1.8)
/ m ( # ) ! / (#)
i f
an
d only i f 777 | n .
fn(x)).
Since c = (rn,
ri) , by the Euclidean algorithm, there exist integers a and 2? such that c =
3??7 + bn. Since c <_ m, m9 n > 0, a 0 or b 0. Suppose a 0 and let & = -a.
Then bn = c + km applied to (1.7) gives
fhn (*) = fa
By ( 1 . 8 ) , fn(x)\fbn(x)
we have d(x) \fc (x)fkm
+ km(x)
= fc.1(x)fkm(x)
and fm (x)\fkm
+ 1(x)
+ fo(x)fkm
1(x).
implies
1978]
= 1, and d(x)\fa(x).
+ 1(x))
fa(x),
d(x)
503
+ Sk_2,
k even,
n)
S2(m,
ri)
S3(m,
= bQ(mn)
= m
= mB0(mn)
n) = mn + 1
ri)
= (mn)
ri)
= mBl(mn)
+ 3mn + 1
= m[(mn)
b1(mn)
b2(mn)
+ kmn + 3] = mB2(mn)
2k + 2 (jn, n) = mBk(mn) ,
2k + i (m>
bk(mn).
S2k^2(m2,
n2)
= f2k+2(mn)
= m2Bk(m2n2)
= % f2k+2
n2)
= m2B1(m2n2)
thus,
(mn).
= m2n2 + 2, and fh (mn) = (mn) 3. +
= ^(mn) (m2n2
+ 2)
Next, we state and prove a matrix theorem in order to derive further results for the polynomials Sk(m, ri).
Th2.0A.rn 2.1:
(AB)k
Let A = (j
/ bk(xy)
(l
o) *
Then
'
xBk_1(xy)\
=
\yBk_1(xy)
where bk(x)
o ) '5
and Bk(x)
bk_(xy)
504
{AB)
-\o
l
(AB) =/^
Assume that (AB)k
xybk(xy)
i
XyB.^xy)
, i
b_x(xy)
/ b Axy)
+ l *\ = f ^
+ bk(xy)
,
xBn(xy)\
^ M
+ xyBk_1(xy)
ybk(xy)
1/
[Dec,
x[(xy
+ yBk_1(xy)
Then,
+ l^.iGrz/)
xyBk_x(xy)
fcfc_i0c7/)r
+ fck-1(a:z/) ,
fcfc-iO*^) xBk{xy)
yBk(xy)
bk(xy)
by applying the mixed recurrences of (1.1), completing a proof by induction.
Now, returning to the matrices of Theorem 2.1, since the determinant of
AB is 1, it follows that
(2.4)
- xyB2k_1
bk(xy)bk_1(xy)
= 1.
045)* =
71
mS2k
$2k + l
"
777
^zk -1 .
2k-iS2k
+ i " w^2k
1B
(^2fc( w
S2k(m,
= nBk_1(mn)
n) -
mBk_1(mn).
Since the polynomials Sk(m, n), the Morgan-Voyce polynomials, and the
Fibonacci polynomials are interrelated by (1.4) and (2.3), which can be rewritten as
Sik + i(m> n) = f2k + 1Umn),
(2.7)
;
171
S2k(m9
n) = -7=/ 9 . (vmn)
1978]
505
and since the coefficients of the Fibonacci polynomials lie along the rising
diagonals of Pascal's triangle, we can write the following theorem.
ThdOtim 2.2: The coefficients of fk(x),
bk(x) , Bk (x) , and Sk(m, ri) are
all coefficients which lie along the rising diagonals of PascalTs triangle.
3.
(Bm(x),
Bn(x))
D
(ii)
and Bn(x)
(m + 1, n + l ) - l
(bm(x)9
bn(x))
= b{(Zm
(iii)
(Bm(x)9
bn(x))
= b ((2m + 2j
Vnoo{\
(i)
x(Bm(x2),
Bn(x2))
+ i, 2 n + l ) - l ) / 2 (x) ,
2n + 1)
_l)/2
= (f2m+2(x)9
= CCB(m +
f2n
(x).
+ 2(x))
n+l)-l(%Z)
lt
= (f2m
+ 1(ar),
/2n
+ 1(^))
Thus,
(2m + 2 ,
bn(x2))
= (xBm(x2)9
bn(x2))
(S m (x ), M * 2 ) ) = f2k
where
+ 1(x)
= bk(x2)
fe = ((2m + 2 , 2n + 1) - l ) / 2 .
2
ri))
VfiOO^i If i and j are both odd, (2,7) and Theorem 1.1 give the above result immediately. If i, and j are both even,
506
n ) , Sj(m,
n))
= (S2k(m,n),
S2h (m,
lK
\vmn
n))
Vmn
= Vmn
-p=(f2k2K (i/?)> f2h
(Sb)
zn
= Si{k,h)
1 so t h a t Vmn)(S2k
{SiQn,
n),
+1
Qn9 n),
S3-(m9 n))
and s i n c e f2k
= (S2k
= (S2k
(m, n) = S(2k,2h)
+ 1(m,
+1
()
(m, n) = ( i , j ) ( m ,
n).
a l s o ends i n 1 ,
S2h (m,
n),
/mnS2h(m,
m
~F=f2{kZKK>n) h^ff)
Vmn
n) always ends i n t h e c o n s t a n t
n),
[Dec.
n))
f2yl^m<^))
n))
=
(f2k
+ 1(i/mn)
>
f2h(Jrnn))
3.3:
Vnoo^\
If i\j,
S^im,
n)\Sj(m,
n) if and only if
i\j.
S^(m, n)\Sj(m,
n) with i )(j, then /^(^)|/. (#) where i^j,
a contradiction of
J
(1.8).
From all of this, we can also write divisibility properties for MorganVoyce polynomials.
TkzotiQM 3.4:
VKOO^:
Bm(x)\Bn(x)
bm(x)\bn(x)
bm(x)\Bn(x)
Bm(x)\Bn(x)
(Bm(x),
Bn(x))
Bn(x))
= Bm(x),
but
= B(
m+15rc+
1) - 1
or
2777 + 1 = (2777 + 1 , In + 2 ) .
1978]
507
= B ( n + 1, m + 1) -
Bm(x))
(B*(x),
and B*(x)
= Bn^1(x)
= Bm_1(x),
then
= Bfn>a)
B*(x))
= { 0 , 1, x -f 2, ...} obeys
and is in fact the Fibonacci polynomial, so to speak, for the auxiliary polynomial A2 - (x + 2)X + 1 = 0, since
Un _ 2
Tn + 1(x)
= 2xTn(x)
= Tm(Tn(x))
Tn{Tm(x))
Tn_(x),
=
Tmn(x).
The property (3.4) is easy to prove from the Binet form associated with the
auxiliary polynomial
(3.5)
A2 - 2x\ + 1 = 0,
Un + 1(x)
= 2xUn(x)
= 1, and
Un_(x).
Ax = Tn(x)
+ /(x2
xl
= Tn(x)
/(x2
- lJUn^iix)
- !)/_! ( a ) .
We prove only one part, since the second part is entirely similar.
/-i(aO = 0, and TQ(x)
= 1, A^ = Tn(x)
+ /(x
- l)Un_1(x)
for n = 0.
Since,
Assume
508
[Dec.
t h a t Ai = Tk(x) +
Then, by ( 3 . 5 ) ,
and
x\
- Tk(x))
/ ( x 2 - l)Uk
+1
= Tk + 1(x)
/(x2
+
/ \ x 2 - 1) U
- 1) (2xUk + 1(x)
(x).
Uk(x))
1(x),
T2(x)
- 1 = (ar2 - l)^_i(x).
Tn(x)
= (Ar+
= x + /x^ - 1.
Xn2)/2.
Replace x by Tm(x)
But X1XZ
+ 1 = 0.
- 2Tm(x)X1{Tm(x))
A2 (Tm (x) ) + 1
Tm(x) =
= 1, so
Tm(x)
2X1(Tm(x))
[A
i^-(x))
+1AI(^W)]/2.
= [ A ^ f r ) ) + A2(Tm(x))]/2.
and
A* = A 2 ( T m ( * ) ) .
Now,
^ n W
Xn2(Tm(x))]/2,
xl-
x\
where Ax and A2 are the roots of A2 - 2xX + 1 = 0. We now list the first few
polynomials and let
U*ix) = / n - i ( x )
U.l{x)
= 0
= U*{x)
U0(x)
= 1
= U*(x)
Ux(x)
= 2x
U2(x)
= hx
= U*(x)
- 1
U3(x)
= 8x
Uk(x)
= 16a: * - 12a;:2 + 1
- 4x = 4 x ( 2 x
= U*(x)
2
-- 1) = U*(x)
= U*(x)
1978]
UAx)
509
- 8x2 + 3) = U%(x)
= U*(x)
U*(x)9
U*(x)
U*mtn)(x).
= 2xU*(x)
{[/ n (x)} i s a s p e c i a l c a s e of
Hoggatt and Long [7] as
(3.11)
Un+2(x,
(3.12)
- C^.!^),
t h e p o l y n o m i a l sequence {Un (x,
y) = xUn + 1(x9
Since U*(x)
y) + yUn(x9
z/) = 1.
d e f i n e d by
y),
y)}
y)9
= U(m> n ) .
un)
k.
Define Sk(a9
b9 c9 d) by taking Sl
^ = aSk_1
+ bSk_l9
= 1, S2 = a,
k even,
(4.1)
fc = c5fc_! + dSk_29
k odd.
+ ^ * _ 2 , fe odd.
_ ^
Let
/#
&\
/<?
^\
( ac + b
li o j ' d o ) i
Now, {Kn}
^2fc + i
^2*
S*2k
dS2k.J
c
\
Xx - A 2 ?
ad\
_,
dp
then
bdKn_2.
/ %-k + i ~ d%k
\oKk
daKk
d(Kk -
bKk_x)
510
[Dec.
b,
c9
d)}.
S, = 1
52 = a
53 = ac + d
Sh = a2c + ad + ab
55 = a2c2
+ 2acd + abo + d2
56 = a3c2
+ 2a2cd + 2a2bo
7 = a3o3
+ 3a 2 c 2 J + 2a2bo2
+ a^2 + abd + a 2
+ 3aad2 + 2a&c^ + ab2o + ^ 3
We consider some special cases. If a - 0, then S2k + z = 0> and S2k + i ~ dk9
k >. 0. If fc= 0, ^2^ + 2= a(ac + d ) k and S2k + 1 = (ac + d)k, & >, 0. If c = 0,
then 5 2 k _! = c^-i and 2 k =a[(dk - bk)/ (d - b)]9 k > 1. If d = 0, then S2k =
a(ae + b)k~1 and /S^k + i = ao(aa + b)*"19
k >_1. The expansions of *(a, b9
c9 d) are not very interesting, since they are the same as those of Sk(a9
b9
c9 d) with the roles of a and c exchanged.
The special case of Sk(a9 b9 c9 d) where b = d proves fruitful. We list
the first few terms of {Sk(a9
b9 c)} below:
51 = 1
52 = a
53 = ac + b
S^ = a2c + lab
55 = a2c2
+ 3abc + b2
56 = a3c2
+ ha2bc
+ 3ab2 = a(ac
+ 3b)
2k + j
^2^ +i +j
j + lS2k
= $*k+i
bSjS2k-l>
Sj + iS 2k + i
+
+
bS^S2k'
that
a n d
aS
2k
cS
*k-
Sj(m + 1)
Sj + lSjm
bS^S
2k
= Sj
From (4.4),
since 5j = S^ for j odd. So, if Sj \Sj and 6^ l^-^, then Sj \Sj(m + i) for j odd.
Thus, for j and m both odd, we see that Sj\Sjm
for all odd m.
Next, suppose that j is odd and m is even; then, from (4.3),
S
2m'j
- Sf+1S2mlj
+ bS-S^.j^,
m =
2m'.
1978]
Now, i f SjlSj
and Sj\S2mlj9
Next, l e t j be e v e n ;
S
= s
2k+2j
t h e n Sj\S(2ml+1)j
*j' + iszk
+ bS2j'S2k-i
511
Sj(m+1).
and
Ik =
2j'm.
Since S-\S2-,
and Sj\S2{jtm
= S2k , we have Sj \S 2j m + zj' ~ -Sj'o + i) T h i s comp l e t e s t h e proof t h a t i f i\j9
t h e n S^Sj,
S i n c e , a l g e b r a i c a l l y , {Si} a r e of
i n c r e a s i n g d e g r e e i n t h e two v a r i a b l e s a and c c o l l e c t i v e l y , SJ-)(Si f o r i < j .
L a s t , u s i n g ( 4 . 3 ) and ( 4 . 4 ) , i t i s now s t r a i g h t f o r w a r d t o show
ThzoKQjm 4 A:
Sj(a9
b> o)\Si(a9
b9 c) i f and o n l y i f j \ i .
(Si(a9
b9 c),
Sj(a9
b9 c))
= S^t3-)(a9
b9
c).
VK.00^: Let P(x) be a monic polynomial of degree v + s with integral coefficients with two factors Q(x) and R(x) of degree v and s9 respectively.
Then,
br
+3
P(x/b)
P*(x9
brQ(x/b)bsR(x/b)
b) = Q*(x9 b)R*(x9
b).
In p a r t i c u l a r , i f P(x) i s of d e g r e e p , T(x)
(P(x)9
T(x)) = J/(ar), t h e n
{bpP(x/b)9
btrP{xlb))
bwW(x/b).
For a p p l i c a t i o n t o Theorem 4 . 2 :
(4.6)
c2S2m{a2,
(4-7).
S 2 r a + 1 ( a 2 , &2, c 2 ) =
Co6 7:
b2,
a2)
= aci2"-Y2BI(ac/Z>) ;
Both s u b s c r i p t s
2
(o SZm(a ,
fc2ffl/2m+1(ac/&).
even.
2
b,
o2S2n(a2,
a ),
b2,
= {acb2m-xfZm(,aclb),
a2))
aob2n-lfln{ae/b))
(2m, 2n) - 1 .
acb^m>^^f{2my2n)iac/b)
=
2
= c S(2mi
(a2, 2,
2n)
c2).
Therefore,
(S2m(a2,
Ca6 2 :
b2,
a2),
S2n(a2,
Both s u b s c r i p t s
2
b2,
S2n
2
+ 1(a ,
= (b2mf2m+1(ac/b),
Ca6 3 :
= S(2n>
b2,
a2))
odd.
2
( S 2 m + 1 ( a , Z? , a ),
c2))
b2nf2n+l(aa/b))
$ (2m + l, 2n + l) ($
s D , <3 ) .
(c S2m(a
c )9
S2n+1(a
h(2m,
2n + l)-1
(2m, Zn+1)
n
i.CC
b , c2))
b2nf2n
= (acb^-^Jac/b),
(rr^/h^
b
O )
even.
1(ao/b))
2B)
512
(c S2m(a ,
b,
c )9
= (S2m(a2,
b2,
c2)9
[Dec.
= 1,
S2n+1(a ,
b,
o2))
S2n+1(a2,
b2,
o2))
finishing the proof of Theorem 4.2 by replacing a2 with a, Z?2 with b9 and <?
with e.
Let f*(x) be a modified Fibonacci polynomial, with
fn to) = / to), n odd,
f2(x)
, n even.
f^W = 1
f*s(x)
= x2 + 1
f*(x)
= x2 + 2
/*(#) = ^ + 3jr2 + 1
f%(x)
+ kx2 + 3
= ^
/ * t o ) = ^ 6 + 5X4 + 6x 2 + 1
/*(a?) = x 6 + 6 ^ + 10;c2 + 4 .
Here,
=
fn + 2 <*>
f*+l
^)
fn^
'
even
'
= Bn(x),
v2n+1(x)
bn(x).
are
vi to) = l
y 2 to) = l
z;3to) = x + 1
b0(x)
B0(x)
bx(x)
7^ to) = a; + 3
= Bito)
v5(x)
= x2 + 3x +
y6to) = x
2
3
y 7 t o ) = a:
vQ(x)
= x
= 2? 2 to)
+ 4x + 3
+ 5x
+ 6x
= S2to)
+ 6^ + 1
b3(x)
+ lOx + 4 = to + 2) to + 4^ + 2) =
B3(x)
1978]
Since Vk(x)
this is {Sk(a,
513
satisfies
( vn(x)
= vn_1{x)
( vn(x)
= xvn_(x)
n even,
+ vn_2(x),
+ vn_2(x),
n odd,
vm(x))
= v
(m,n)
REFERENCES
1.
2.
3.
4.
5.
6.
7.
8.
9.
*****
University,
Philadelphia,
PA 19122
514
[Dec.
Kyrie
Kyrie
Kyrie
un 11
Christe
Christe
Christe
un 11
Kyrie
Kyrie
Kyrie
eleison.~\
e'leison.
eleison.
sect ion
eleison.
eleison.
eleison.
total chant
eleison.
eleison.
eleison.
Each chant falls into nine separate sections. The three repetitions of the
sections form three larger units which, in turn, make up the complete chant.
While there is considerable variety in the melodic treatment of the text, the
text itself had remained constant in the above form since ca. 900.
The actual nature of the rhythm of these chants is still open to question.
Because music is a time art, any analysis that does not account for the proportional movement of the pitches in time cannot pretend to be a statement
about the total nature of the music. In this sense, the following findings,
though factual, remain theoretical to the degree that while pitches in succession imply time, exact temporal proportions are not deducible from that
succession alone. In addition, the reader should be advised that there are
more than 200 "Kyrie" melodies known to exist.
In this light, the chants
analyzed for this study represent a sampling of the repertory.
METHOD OF ANALYSIS
Because different treatments of the same text are usually set to different pitches, 146 distinct musical sections are present in the 30 chants, the
remaining being exact repetitions of other sections. The pitches in each section were totaled, and <f> was determined for each section. A section was examined to determine if any significant musical event occurred at either the
major or minor mean. A significant event was defined as the beginning or ending of a musical phrase. The three statements of the "Kyrie," the "Christe,"
and the "Kyrie" tend to form larger units; these were analyzed according to
the same procedure. Finally, the pitches in the complete chant were totaled,
i.e., nine separate sections of text.
THE FINDINGS
Applying the analytical method described above revealed the presence of
the golden section in 105 of the 146 individual sections of the "Kyries" in
the Liber
Usualis.
These 105 sections make up .72 of the cases. The major
1978]
515
mean precedes the minor mean twice as often as the minor mean precedes the
major mean. Example 1 is a section of chant conforming to the M:m proportion,
Ky-n - e
|e_j_son
Example 1*
Example 2 shows the proportion in reverse.
e-
le-i-son
Example 2*
University
CA 94117
INTRODUCTION
Let F(n), L(n) denote the nth Fibonacci and Lucas numbers, respectively.
(This slightly unconventional notation is used to avoid the need for secondorder subscripts.) Consider the equation
(0)
F(m)
"Source:
Liber
GP,
Usualis
516
[Dec.
z(q).
If (x, y) = 1 and xy = zn, then a? = un and y = vn, where (u, v) = 1 and
wi> = z.
(1)
(2) F(2n) =
F(n)L(n).
/o\ /n/ \ r/\\ ( 2 if n = 0 (mod 3)
^0 (mod3)
(7)
2|F(w) -* 3|n.
(8)
(9)
(F(n), F(kn)/F(n))\k.
(10)
(11)
= 1.
Re.maA.kA: (1) through (8) and (11) through (14) are elementary and/or wellknown; for proof of (9), see [2], Lemma 16; (10) follows from (8) and (9).
THE MAIN THEOREMS
For a given prime, p, let m = m(p) >2 be the least integer such that, by
assumption, (0) has a nontrivial solution. By inspection,
777(2) = 1 2
a n d 7??(3) = 6.
Tk&Ofim 7: If m = 2n > 2 is the least integer such that F(m) = cv, where
p is prime, then either (i) m = 6, p = 3, or (ii) m = 12, p = 2.
Cctt>& J I f n 0 (mod 3), then by hypothesis, (1), (2), and (3), we have
1978]
517
F(n) = bp.
If b > 1, we have a contradiction, since n < m. If b = 1, then
hypothesis ->n = 2->m = 4-> F(m) = 3, a contradiction.
C<UZ 2 I f ft = 3 (mod 6 ) , then (4) -> F(n) = 2p, L(n) = 4s, with rs odd.
Now hypothesis and (2) -> F(/77) = 8PS = cp, so that (5) - p 13 -> p = 3. By [3] ,
we must have
c=2,n=3,m=6.
Co6e 3-If n = 0 (mod 6 ) , let n = n0 = 2J3fet, where j , fc > 1 and (6, ) =
1. Let m = 2~Vz0 for each i such that 1 <_ i <_ j . Let fr0 = rij = 3 k t, and let
hi = 3~^/z0 for each i such that 1 <_ i <_ k, so that t = hk.
By (6), we have
F(n) = 22+^r, L(n) = 2s, where rs is odd and (p, s) = 1. Now hypothesis, (1),
and (2) imply r = r%9 s = SQ, with r0s0 odd and (P 0 , s 0 ) = 1. Therefore, F(n)
= F(n0)
= 22+JTQ,
L{n)
= L(nQ)
= 2s p ,
^OSO
'
Since ni
peat our reasoning to obtain F(m) = F(rii + 1)L(ni + 1) = 2 Z + J'"2'p?, L(ft^) = 2s|
for i = 0, 1, 2, ..., j - 1. By (4) we have F(hQ) = F(rij) = 2r?, L(n3-) = 4sJ;
moreover, r^s^ = 2^ _ x is odd and (r^ , s^) = 1 for i = 1, 2, 3, ..., J. Now,
let rj = Mo, so that F(7z0) = 2up. We have F(hi-i)
= F{hi)^F{hi.{)/F{hi)
for
i = 1, 2, 3, . . . , k.
By (7), (10), and (1), if i < k, we have F(hi) = 2u\,
F{hi-i)lFQii)
= v\\ if i = k, we have F(t) = F(hk) = up, F(hk_x) /F(hk)
= 2^;
moreover, (ui, Vi) = 1 and u^Oi = u^_i is odd for i = 1, 2, 3, ..., fc.
But (11) + F(t) < F(6t) <. F(n) < F(m) = c p - * u k = l - * t = l .
If k >_ 2,
then F(hk-1)lF(hk_1)
= F(9)/F(3)
= 17 = v_ -> p = 1, a contradiction. Hence,
k = 1, ft0 = nJ = 3- I f C ^ 2> t h e n Li^o-2)
= (12) = 322 = 2s|_ 2 -* sf_ 2 =
161 -> p = 1, a contradiction. Therefore, k = j = 1, ft = 6, m = 12, p = 2.
Cofio&Lafiy: If (0) has a nontrivial solution for p _> 5, then it has a nontrivial solution such that m is odd.
VfLOO&i The proof follows directly from Theorem 1.
Tk&Oltm 2: If F(m) = cv > 1, where the prime p ^ 5, and m is odd, then
either (i) m = 1 (mod 12) and c = 1 (mod 8) or (ii) m = 5 (mod 12) and <? =
5 (mod 8; furthermore, if q is any prime factor of c, then 2/(q) >. 5, so that
q e {5, 13, 37, 73, 89, 97, 113, 149, 157, . . . } .
Vnooi'. If 2\c9 then 2 P | op -> 2 P |F(m) , so that by (6), 3*2p~2\m,
contradicting hypothesis. Now c is odd, so that F (m) is odd, and by (7), 3)(m. Therefore, m E 1 or 5 (mod 12). If q is any prime factor of o, then
(12) + #(4) |g(q) |TTZ.
Since ( 6 , m) = 1, we must have y (q) >_ 5 .
Co6 ?If w = 1 2 t 1, then (13) -> F(w?) = F ( 6 t ) 2
F(6t) = 0 (mod 8) and F(6t 1) i s odd, so F(6t l ) 2
cp E 0 + 1 El (mod 8 ) , and (14) i m p l i e s c = 1 (mod
Co6e 2 I f m = 1 2 t 5 , t h e n (13) > F(TT?) = F ( 6 t
Now, F ( 6 t 3) = 2 (mod 8) and F(6t 2) i s odd, so
T h e r e f o r e , c p E 4 + l E 5 (mod 8 ) , and (14) i m p l i e s
+ F ( 6 t l ^ 2 = c p . Now,
E l (mod 8 ) . T h e r e f o r e ,
8).
3 ) 2 + F(6t 2 ) 2 = cp.
F(6t 2 ) 2 = 1 (mod 8 ) .
c = 5 (mod 8 ) .
REFERENCES
1.
2.
3.
4.
Then,
^).
If ft + 1 is even, then (ft + l)/2 is an odd number _> 3, since ft >_ 3 in this
case; obviously, thereforey In 4-ftis an odd number >_ 5. Hence, m is a composite number. If ft + 1 is odd, then ft is even, and since ft >. 2, one must
have that In + ft is an even number >_ 4. Since m is an odd number, it follows
that {In + ft)/2 is an odd number > 2. The fact that ft + 1 >. 3 now shows that
m is a composite number.
VKOO^ (<=): Consider an arbitrary factorization m = kg,
(3 <_ft<_ g).
a = g
Hence, we have
E:
^e
m = / \ , (q + p - 1) ,
p = i
which i s a a - p a r t i t i o n ,
s i n c e a e IN and ft J>. 3 .
Then,
ft + l
, N
w = - (/ 20 n + 7ft).
If m = 2s , then, since ft >. 2, we must have that ft + 1 = 2t, _> 2, and In + ft
= 2", u >. 2. This is a contradiction, since ft 4- 1 = 2* would imply that ft is
518
Dec. 1978
519
CcU>& U) :
W r i t e fc = i; and g = 2 n .
m = kg = ~^{2g) = where
2(g - ^~Y^)
n >_ 1, y an
Then,
+k
- A= l(2a
+k
- 1} >
^ - 1
a
= g _
i$jm
Thus we have
k
]T (a + r - 1)
m =
r=1
Now,
k > ^-^k
<L
-^
^>
and we have
m = gk = f (2/0 = f [ 2 ^ - 2 - ^ A ) + ^ _ ij = | ( 2 a + ^ - 1)
where
Hence,
a = k - ^"2"^.
9
m = y^ (a + p - 1 ) ,
r =1
a a-partition.
Ik
2k
-Y9 = -y
where
a = ^ ^
( S L X^ - k) + 2k - l\ = *, ^ + r - 1)
[Dec.
520
x 3l
< x 3
^)
If g < 2k + 1, then k - ^^
Z
that
(2)
t[(*
~^)+r-l]
is a a-partition of 77?. Clearly, (2) coincides with the fixed partition (1)
if and only if k = g9 i.e. , the case where m is a square and is factored as
such,
If g _> 2k + 1, then clearly ^z
k E IN, and we obtain the a-partition
2k
(3)
17
E [(^HP - *) + * - ! ] f m-
The partitions (1) and (3), having different lengths, can never coincide.
Conversely, the indicated possible a-partitions corresponding to the particular type of factorization of m are the only possible ones m can have.
For, consider an arbitrary a-partition
r=1
V = 1
(2)
1978]
521
and
<}
Lv^ir1 ~ k) + r ~ X J
(if
a > 2k + i)
If g < 2k + 1, the two valid partitions (J) are (2) are different, except in
the case where g = k9 and if g >_ 2k + 1, the valid partitions (1) and (3) are
always different.
Now consider any two different factorizations (if they exist) of 777 into
two factors: m = kg (3 k <_ g) and m = kTgr (3 <_kr <_ gr) . Then, comparing
the lengths of the resulting a-partitions, one sees that every possible partition corresponding with m = kg differs from every possible one corresponding with 777 = kTgr.
Hence, if 777 admits t different factorizations of the form
m = kg (3 <_ k g) , then 7?? admits It different a-partitions, except in the
case where m is a square in which case the number is 2t - 1.
Finally, we consider the nature and number of a-partitions of even numbers other than powers of 2.
Tfieo^em 3: Let m be an even number other than a power of 2. Then, there
exists at least one factorization of the form m = kg (k < g) , where one of
the factors is an even number and the other an odd number J> 3. For each such
factorization, exactly one of the following three conditions holds:
(1) k is even, g is odd and g < 2k + 1;
(2) k is even, g is odd and g >_ 2k + 1;
(3) k is odd, g is even;
and only t h a t sum i n t h e
list
(V)
2k
r/
-1
(2<)
'
Finally,
PflOO&i First we observe that m may be written in the form m = 2nu, where
n >_ 1 and u is an odd number >_ 3. Now, consider an arbitrary factorization,
m = kg (k < g), into an even and an odd factor. If k is even, then the possibility (3) is ruled out and the ordering axiom ensures that exactly one of
(1) and/or (2) holds. If k is odd, then the first two possibilities are excluded and (3) obviously holds.
Concerning the next part of Theorem 3, we note that each of the indicated
sums is of the form
E',(a + v - 1),
V = 1
and that n >_ 3 and a e fft, providing that the condition to which the particular
sum corresponds holds. Moreover, in each of these cases the indicated summation results in the product kg. On the other hand, the validity of any given
condition (i) clearly results in a f IN in the sums (j f) 5 i 4 j . This concludes
this part of the proof.
Dec. 1978
522
F i n a l l y , we consider an a r b i t r a r y cr-partition
n
m = ^2 (a + v -" !)
j(2a
+ n - 1).
r=1
If n is even,
one must have
write k = n/2
Moreover, g >
the
Theory
of
Numbers,
G. L. COHEN
Wales Institute
of Technology,
Sydney,
Australia
(1)
- -<log2,
i = 1 r"z-
a result due to Suryanarayana [5], who also gave upper and lower bounds for
t
Z ^
< .673634,
the upper bound in [5] being .6 73770. We will also give a further improvement
in the upper bound when 5 is a factor and 3 is not; namely,
V- 1
> < .677637,
the upper bound in [6] being .678036. (These are six-decimal-place approximations to the bounds obtained.)
We assume henceforth that n is an odd perfect number.
An old result, due to Euler, states that we may write
i =l
523
524
(2)
\{pk
[Dec.
result
l)\n,
(3)
PX \ipk
+ 1).
We a l s o use t h e i n e q u a l i t y
1 + x + x2
(4)
0 < x <_ -k
To prove t h i s , n o t e t h a t
exp(x
+ ^x2^)
- (1 + x + x2)
= 1 + x + ~~ + - y L r + -)
_, x3 ^ xh , 1 /
+ - (1 + x + x2)
,x2V
x2\3
x2Y
Now,
f + f U n + Jiiir- 09
and
i / , * 2 V . i I , *2X*
a; + -7+
r uc + -7+
4
3!x2\
W
4!a; 2 \
<
^(*
1 /l3\
~ 18 \ 1 2 /
T )
3
36
23
<
T) + (*
T)2
'2'
Hence (4) is true. Other and better inequalities of this type can be established but the above is sufficient for our present purposes.
Now we prove (1). Since n is perfect,
t
In = o(n) = Y\ (1 + p. + p?
+ . . + p*1)
' - f t K - * )
By E u l e r ' s r e s u l t , ak >_ 1 and a^ >_ 2 (i- ^ k) , so
P k /i = i \
by ( 4 ) .
Pi
P?/
Hence,
log2>1 8 1+
( t) + 5'(t + ^)
Pfc/i.'i
VPi
4P;
1978]
t ,
525
>
~ti?i
ttx Pi
(Pk + D
4&P?
4Pl2
M P ,
2p 2
2p , 4 l P i
using ( 3 ) .
We end w i t h t h e
Th2.OA.2m: ( i )
I f 15 I n ,
then
x - l ^ l ^ l ^ l ^ l , ,
AFT
(ii)
I f 5 | n and 3 | n ,
l 3+ 6 l
+ l o g
= lr^
<
= a s a y
'
'
then
v - 1 ^ 1 ^ 1 ^ 1 ^ ,
Z-^7
2950753
2815321
31 + 61 +
l o g
293105
190861
= 6
'
Say
'
61 61 2
P;* 13,61
V^ 1 <, -,
o , 13 , 6
^ - < l o g 2 - l o g T - l o g - -
183
14076605
log ^
- log 1 3 8 4 5 8 4 1
i =i
+ 1 x 1 ,1 , 1
3
If ot/13N .> 4 , t h e n we s i m i l a r l y
13
61
'
obtain
ZjT
526
[Dec.
Then,
X ) ~ < lg 2 - logfl +f + ^
Log
(1+I
A+J_
< a
j. i
a. i
!
+
log T +
log 25 " !3 " j
" - !r
But p _>_ 13 (though we can easily demonstrate that in fact p >_ 17), so,
t
l ^
,
13
3 1 1 1 1 ^
< l o g 02 - l o g - nl o g + - + - + 333 < a.
i = l ^i
t^<iog2-iog(i
i)-io g (i
- i T + -i ?
-i,;
- i o g ( l + ^ + i + ^- + ^-= b.
If pk = 1 3 , t h e n , by ( 2 ) , p 2 = 7 .
A l s o , ot , j 3 .> 1 , so
L ^
a ( 7 2 ) = 57 = 3 - 1 9 , so a 2 >. 4 , s i n c e 3/fn.
logfl +-~
+~
31
+ ^j 3+ ~4
31 z
31 31
Since
1978]
If p
527
Then pk
M. Buxton & S. Elmore, "An Extension of Lower Bounds for Odd Perfect Numbers," Notices
Amer. Math. Soc., Vol. 23 (1976), p. A-55.
D. B. Gillies, "Three New Mersenne Primes and a Statistical Theory,"
Math. Comp. , Vol. 18 (1964), pp. 93-97.
Guiness Book of Records, 22nd ed., 1975, p. 81.
P. Hagis, Jr., "Every Odd Perfect Number Has at Least Eight Prime Factors,"
Notices
Amer. Math. Soc,
Vol. 22 (1975), p. A-60.
D. Suryanarayana, "On Odd Perfect Numbers II," Proc. Amer. Math.
Soc,
Vol. 14 (1963), pp. 896-904.
D. Suryanarayana & P. Hagis, Jr., "A Theorem Concerning Odd Perfect Numbers," The Fibonacci
Quarterly,
Vol. 8, No. 3 (1970), pp. 337-346, 374.
IRVING ADLER
Bennington,
VT 05257
1.
While working on some mathematical aspects of the botanical problem
of phyllotaxis, I came upon a property of simple continued fractions that is
simple, pretty, useful, and easy to prove, but seems to have been overlooked
in the literature. I present it here in the hope that it will be of interest
to people who have occasion to teach continued fractions. The property is
stated below as a theorem after some necessary terms are defined.
2. TeAmZnoZogy:
For any positive integer n, let n/0 represent . Let
us designate as a "fraction" any positive rational number, or 0, or , in the
form alb,
where a and b are nonnegative integers, and either a or b is not
zero. We say the fraction is in lowest terms if (a, b) = 1. Thus, 0 in lowest terms is 0/1, and in lowest terms is 1/0.,
If inequality of fractions is defined in the usual way, that is
alb
< eld
if ad < be,
528
[Dec.
3. Tfie M&cLiant: If a/b and c/d are fractions in lowest terms, and a/b
< c/d, the mediant between a/b and c/d is defined as (a + c) / (b + d) . Note
that a/b < {a + c)/(b + d) < c/d.
ExamptoJsThe mediant between 1/2 and 1/3 is 2/5. If n is a nonnegative
integer, the mediant between n and is n + 1. If n is a nonnegative integer
and m is a positive integer, the mediant between n and n + l/m is n +1/(777 + 1).
4. A Mediant N&>t: A mediant nest is a nest of closed intervals
..., T n , ... defined inductively as follows:
I0,Il9
^o = [0, " J .
For n >_ 0, if In = [r, s] , then .Tn+ 1 = either [r, m] or [m, s] , where 777 is the
mediant between v and s.
It is easily shown that: if at least one In for n >_ 1 has for form
[r,m],
then the length of Jn approaches 0 as n -* , so that such a mediant nest is
truly a nest of intervals, and it determines a unique number x that is contained in every interval of the nest. For the case where every In for n >_ 1
has the form [m, s] , let us say that the nest determines and "contains" the
number . Mediant nests are obviously related to Farey sequences.
5. Long Notation
^OH. a Mediant N&>t: A mediant nest and the number it
determines can be represented by a sequence of bits b1b2b3 ...b^
... , where,
for > 0, if Ii^i
= [P, s] and 77? is the mediant between r and s, bi = 0 if
Ii = [P, m] , and bi = 1 if _Z\- = [777, s] .
ExampZ&> 0 = 0; 1 = ; 10 = T, the golden section; where each of these
three examples is periodic, and the recurrent bits are indicated by the dots
above them.
6. kbbfizvtat&d
Notation
^OK. a Mediant N&6t: The sequence of bits representing a mediant nest is a sequence of clusters of ones and zeros,
b1b2b3
. . . bi
Tko.on.zmi
I. The nonterminating case, # = (a1, a 2 , ..., a^, . . . ) . Thus, x is irrational. Let p^/q^9
for >_ 1, be the principal convergents of a 2 + l/a2 +
... . Then a straightforward proof by induction establishes that for all
even _> 2,
1978]
529
x + --'+^ = [Pi/Rt>
Pi-i/Vi-i]Consequently, the nest determined by successive pairs of consecutive principal convergents of al + l/a2 + + l/an + defines the same number as the
mediant nest (ax, a2, . . . , a n , ... ).
II.
where
Pn + 1/qn
Since
+1
= [pn/qn,
Pn + 1/qn
= (pn_! + an + lPn)/(qn_1
l]
o r
[ P n + l / ^ n + 1* P / ? ] >
+ an + 1qn).
it follows that
x =
= p n /q n = a1 + l/a2 + + l/an.
III. The "conversely" in the theorem follows from the fact that the mapping of the set of mediant nests into the set of simple continued fractions
established in I and II is one-to-one and onto.
ExamptzThe mediant nest (0, 2, 3, ) and the continued fraction 0 + 1/2
+ 1/3 represent the same number. Verification:
a.
Interval
Interval
lo = [0/1, 1/0]
Ii = [0/1, 1/1]
(0 + 1)/(1 + 0) = 1/1
Ii
= [0/1, 1/2]
(0 + 1)/(1 + 2) = 1/3
I,
= [1/3, 1/2]
(1 + l)/(3 + 2) = 2/5
I,
= [2/5, 1/2]
(2 + l)/(5 + 2) = 3/7
I5
= [3/7, 1/2]
(3 + l)/(7 + 2) = 4/9
I,
= [3/7, 4/9]
(3 + 4)/(7 + 9) = 7/16
In = [3/7, mn-i]
(0 + 1)/(1 + 1) = 1/2
Since lim rnn_1 = 3/7, the number defined by this mediant nest is 3/7.
n --
b.
^ - = 3/7.
2+
West Virginia
University,
Morgantown,
WV 26506
Readers of this journal have long been interested in Bode's Rule, see,
e.g., [26] and [15]. Indeed attempts to solve it have ranged from those devoid of science but fairly accurate to those based on some physical principle (s) but rather inaccurate, such as Berlage's and 0. Schmidt's theories.
The problem of planetary motions was first tackled by Eudoxus, who proposed
rotating tilted concentric spheres, rather like a gyroscope, to explain each
planet. When Kepler solved the problem of their motions with the concept of
areal velocity, the area swept out in an invariable plane divided by the time
is a constant, the problem of a law for their spacing remained. Indeed, I
think the unit of angular momentun should be named after Kepler for his contribution of area as a vector. Bode's problem is of great value to the history and especially the philosophy of science. The qualities that distinguish
pure mathematics are succinctness, elegance, fertility, and relevance to the
unsolved problems. But to a scientist the first criterion is reproducibility.
The multifarious, variegated, and at times loquacious and mellifluous monographs on this aspect of cosmogony attest to man's persistent and insistent
attempts, at times based on specious assumptions, to find order in a theatre
of nature that may have no reason to be other than nearly random. Such is
one view. But while science corrects its mistakes (so far) it must be remembered that Boltzmann committed suicide because his contemporaries would not
accept his counting of molecules, Wagoner's 1911 theory of continental drift
was not believed until the 'sixties, and Newton's theory was not believed on
the continent until Clairaut's prediction of the return of Halley's comet in
1759 (P = 76.75 1.5 yr) came true. Such is the lag between prediction and
proof. The final answer to Bode's problem will be known within twenty years
when sophisticated computer simulations are finished. My work will probably
remain the most accurate, namely 1-percent with a few exceptions either way.
In any case, my work has led to some interesting mathematics, especially the
Self-Lucas property (see Section 2). May it be that Urania and Euterpe have
recessed a part of Nirvana to sequester all who have slaved over this vexing
problem. The impetus for this paper comes from Kowal's [17] recent discovery
of an object between Saturn and Uranus that prima faciae,
see Table 1, fits
my rule [6, 7] and does not fit any other rule published! Sequentially, I
present an overview of the history of Bode's Rule, Kowal's discovery and then
generalized folded sequences.
1.
530
Dec. 1978
Dermott
von Weizsacker
Greig
Quadranacci
Dermott
Cale
Pierucci & Dermott
Gaussin
Blagg
Dermott
Belot
von Weizsacker
Greig
Quadranacci
Titius
b
b
b
b
b
b
b
b
b
b
b
b
b
b
b
=
=
=
=
=
=
=
=
=
=
=
=
=
=
=
531
2
= 1.2599 (Saturnian moons)
1.370889 (10 eddies, inner planets),
1.378241 (Saturn's inner moons)
1.380 (see below)
3
= 1.4422 (Uranian moons)
0/J3 = 1.5115 (see [15])
41/3 = 1.5874 (Jovian moons)
1.72
1.73
61/3 = 1.817121
1.886
1.894427
1.8995476
1.905 (see below)
2.00
See Gould [15] for references I have not cited. Dermott (see [2] and [10] of
paper [7]) was forced to take Earth and Venus together to retain his period
factor of /6~. DermottTs arbitrary period factor of /2~ for Saturn's moons
misses both Rhea and Janus.
Indeed, I have myself happened upon some rather well-fitting arbitrary
rules. One such is a bisection of the Quadranacci recurrence Qn+i = Qn+Qn_3,
namely: 1, 1, 1, 1, 2, tf=3, 4, /=5, 7,4 = 10, 14, J = 19, . . . , = 95,
..., which is very good at representing reciprocal distances. Another one
for the distances begins: 4, 7, 10, 15, ... . The distance factor, b, is
(1.380278)2 = 1.905166. The most complex rule of which I am aware is Rothman's (see [9]), d = n(5. 5 + F2)/9 (1 +Fn ). It has seven parameters: n, 5.5,
2, 9, and 1, and two to determine the Fibonacci sequence, and since only 9
planets are fitted then only 2 degrees of freedom are left which is unscientific. All of the above rules are arbitrary, except von Weizsacker's and my
own. My own view is that any rule with more than two parameters violates
yvaetev
necessitatem."
Dermott {ibid.) proposed different period factors for each satellite system, whereas my theory is simpler since the same limiting ratio, 0 , applies
to all. He also ignored the outer Jovian and Saturnian moons. I chose to
emphasize them. This suggests the principle of Contrary Ignorability: whatever earlier researchers ignorethat is the path to pursue. My work indicates that outer Jovian moons should cluster well within 10 percent of 97,
257, 730, and 608 days (Table 2 of [6]). The announcement of Jupiter's XIII
moon [21] at 239 day and i = 27 [25] came after my initial work [4], [22]
and satisfies the above sequence. I also studied the relevance, if any, of
rotation periods and grazing periods of parent bodies in [4], using:
(35)
The period of a satellite just grazing the surface of Saturn, the Sun, Jupiter, and Uranus would be 0.167, 0.116, 0.12, and 0.11 day. The criterion
for coalescence against tidal forces may be written
(36)
pm > /M/distance3
or
where m, M are the satellite and parent masses, p m and pM their densities, Pm
the period of a satellite at this Roche limit, and / a factor between 2 and
10 [24, p. 18]. When (36) is not satisfied "rings" result. In [4] an inner
Uranian moon was suggested which would have a period of 1/1.3292 = 0.752 day
532
[Dec.
according to (17) of [6]. It has even been proposed that the separation, a,
between binary stars satisfies Bode's rule [27], but I am very skeptical of
that.
My own interest in Fibonacci numbers dates back at least to 1966 when I
obtained a copy of Vorobyev's book. I tried these numbers on the planets
with what seemed good accuracy and communicated this to Gould [13] , pointing
out the relevance of bisected Fibonacci sequences. After long arduous efforts, I thought I had put the problem to rest when news of Kowal's discovery
[17] of a planetoid between Saturn and Uranus, too big to be a comet nucleus
and too small to be a large planet, was announced. He calls the object Chiron
after one of the Greek half-man/half-horse animals. Is this discovery to
prognosticate that this Chinese year 4676 (see [18]), beginning 7 Feb. 1978,
should not be the year of Earth-Horse, but rather the year of the Centaur!?
One can see in Table 1 that Chiron fits very neatly into the bisected halfinteger sequence. I could have predicted this object three years ago [5, 6]
from the folded sequences I had discovered but it would have been considered
wildly delusionary at the time. The major body that should occur before Neptune in my sequence given in [5] and Table 2 of [6] is easily calculated to
have a reciprocal period corresponding to -1974 - 4558 = -6532. Its period
should then be 317816/6532 = 48.66 yrs. The agreement with ChironTs period
of 47 to 51 yrs is quite good. The ellipsis (...) in [5] and [6] clearly indicated that the sequence continued in both directions so that a body at Chiron's position was implied. In an earlier work [4] I had stated, " . . . one
should really ask why don't Jupiter, Uranus and Neptune have a plane of particulate matter [rings] inside Roche's limit since that is natural considering the pervasiveness of grains and cometesimals . . ." (p. 16). As we now
know, rings have been found around Uranus [16], [28]. The way to test my theory would be a computer simulation using reciprocal periods given by (23a) or
(23b) or [7] (or, perhaps, using part of a bisected odd-N folded sequence in
[6]) to see if the broad maxima in [10] can be sharpened.
The agreement of Chiron's period with Folded Fibonacci sequences is
reassuring but not perfect. I have been able to represent Neptune and outer
satellites in general more accurately than any other rule simply because I
worked with recursive sequences rather than naive power laws. Also, mine is
the only work to represent the several comet groups (see Table 1). So that,
in terms of completeness and goodness of fit, my hypothesis is the best. It
remains for a computer simulation to test whether my proposal gives maximum
stability. Such a simulation may solve the following question: Why are some
period ratios nearly but not quite small integers? Saturn:Jupiter is not 5:2
but is 6551:2638 to seven, significant digits. Accurate to only five digits
is 149:30. Neptune:Uranus is not 2:1 but is 51:26 to five digits. And Uranus : Saturn is 77:27 to nearly five digits. Similarly, Earth:Venus is not F7:
F6 but 1172:721 to eight-digit accuracy. These ratios suggest that low-order
commensurabilities (LOG) are avoided, except for the ratio 2 among the Galilean satellites. The Kirkwood gaps indicate that LOC are unstable if the
ratio >i 2, such as 11/5, 9/4, 7/3, 5/2, 8/3, and 3/1.
The problem is ancient. The Pythagoreans believed in orbits in arithmetic progression and added a Central Fire and a Counter-Earth [23] to obscure
that Fire so that the total number of moving bodies be the "magic" number
1 + 2 + 3 + 4 = 10. Yet Aristarchus placed the Sun in the center for reasons
of simplicity 18 centuries before Copernicus. Later, Ptolemy and others confounded the picture with equants and epicycles, until Kepler discovered that
blemished curvethe ellipse.
1978]
533
Period
(yrs)
-550 - 233
340 + 144
-210 - 89
0.999969
1.618
2.6183
130 + 55
-80 - 34
4.235
6.859
50 + 21
-30 - 13
20 + 8
-0 - 5
0+ 3
11.07
18.035
28.66
48.66
69.73
- 2
0+ 1
0- 1
20
30
50
80
130
210
340
550
890
1440
2330
3770
6100
9870
1
1
2
3
5
8
13
21
34
55
89
144
233
Solar System
161.00
123.0
521.0
99.5
83.55
45.42
29.42
17.85
11.11
6.849
4.237
2.6177
1.618
1.0
0.618
0.382
0.236
2.
null
Hungaria #434 (991 da = 2.71 yrs), etc. (Average
= 2.75 yrs)
null
Faye (7.35 yrs); Brooks II (6.72 yrs); d'Arrest
(6.67 yrs), Finlay (6.90 yrs), etc.
null
Neujmin (17.97 yrs)
null
Chiron (47 to 51 yrs), other Centaurs
Olbers (69.6 yrs); Brorson-Metcalf (69.1 yrs);
Pons-Brooks (71 yrs); Halley
Neptune (164.79 yrs); N + P (168.4 yrs)
Swift-Tuttle (119.6 yrs); Barnard II (128.3 yrs)
Planet X (464.? yrs)
null
Uranus (84.01 yrs)
null
Saturn (29.46 yrs)
null
(Jupiter 11.86 yrs) not meant to fit Jupiter.
null
Astrea (4.13 yrs); asteroids (0.23 yr )
null
(Mars)
(Earth)
(Venus, 0.615 yr)
null
(Mercury, 0.241 yr)
GENERALIZED FOLDED SEQUENCES
534
[Dec.
point differs from (1) of [6] wherein folded sequences were made infinite by
repeating the cycle ad infinitum.
It will help to display the Folded Pell
array {q>2} in Table .2:
TABLE 2
Folded Pell Sequences for N Odd
Sum
N
1
3
5
29
7
9
985
-1
-11
13
41
169
-69
31
-7
17
27
71
-407
171
-65
41
17
75
167
239
1393
409
11
5741
-2377
987
-403
181
-41
99
157
413
983
2379
13
-13859
5743
-2373
997
-379
239
99
437
973
2383
5739
-13855
5753
-2349
1055
-239
577
915
2407
5729
13865
5811
-2209
1393
577
2547
5671
13889
6+5/8
-2 -2/8
2+/8~
2 +0
6 +/8~
15
17
6 +/8 -2-0
-5/2
-7/2
-3/2
8119
13861
47321
-h
/2
3/z
5/2
7/2
+1
/{cPj,N)k
Define 0 = N modulo 4.
k = (N - l)/2 =
+1
Since r = (k - N/2)
then r = -h gives
[N/2].
= -1, from which one finds
= -z =
-Pk/P-k..
= i^^a
Then,
- g ) / a + z)
(38)
{S.-U/isA
= (a + D/(a - 1) = -{S$}/{$}-,,.
is suppressed.
1978]
535
(41)
(2 + d)P*+h
and Pn = (*n+h-
= 2*+ 1 + r
Qn.h)/jd,
= ( p + 5 ? ^ ) ,
(42)
Q r = (ar + i3 r )
and *
= (ah + ^3_/z) = (a
n+1
and analogous to Fn
(-l)
(45)
(5.
^ > Y>
'
^ 2
+ oT") =
F.n.
-i($}
) , where i,
-1,
we have
r+h
*PH
Now the initial values are determined by (37) or (38) and are
(46)
(2 + d) =
j,-fe
and
for all j .
J.?J
The bisections of the general Half-Integer sequence for N = 1 (mod 4) appear in (47) and (48), where t instead of j is used from now on for the parameter. Compare with the penultimate rows of Table 2. Also (49) is the
h
/5.
subscript r.
In Table 3 note that (2 + d) = Z,' = (a + a~n)/L and V
(47)
(th+dth+it2+3dt2+d+2)
(t2+dt2+d+2)
-(t3+dt3+3t+2dt)
(48)
2+d
(t2+2+d)
(t3+3t+dt)
-(t+dt)
(49)
TABLE 3
Parameters of Half-Integer -Fib Sequences
I-h/Ih
(ah + a-h)
(ah-a'h)
td
Jv - 2
/y + 2
1.272020
8.352410
2.014490
0.241187
1.000
1//I
3//2
1.414213
5.828427
2.030104
0.348311
1.500
/5
1.618034
4.236068
2.058171
0.485868
2.236
1.5
5/2
2.0
3.0
2.121320
0.707107
3.750
/8
2.414213
2.414213
2.197368
0.910180
5.657
/5
2.618034
2.236068
2.236068
1.0
6.708
8/3
10/3
3.0
2.0
2.309401
1.154700
8.888
/n
3.302776
1.868517
2.367605
1.267103
10.817
3.732050
1.732050
2.449490
1.414213
13.856
/20
4.236068
1.618034
2.544039
1.572303
17.888
/32"
5.828427
1.414214
2.828427
2.000000
33.941
2/3
536
[Dec.
Table 3 shows us that the Half-Integer Pell sequence is the only one in
which the ratio of the central pair, S_-h/Sh equals the characteristic root.
This is clear from (a + l)/(a - 1) = a whose solution is 1 + /l.
Sequences
given by (a^ + a~h) = a have a very simple Binet-like formula but the larger
root a = 2.1478990 is the solution to a cubic. Are there Half-Integer sequences with integer terms? Consider t-Fib sequences,
^t,
n+1
tPt
t n
we see that both bisections are integers (after dividing by common factors).
Furthermore, one bisection consists of every fourth Fibonacci number including F5 and the other consists of every fourth Lucas number including L5. Can
this be generalized? As a little algebra shows, yes,
(51)
di
= d\ - 2
or
d = t\
+ 2
or
= thdh
ti
where t 4 , dh refer to the sequence from which every fourth term is extracted
and di, ti refer to the chosen sequence. To illustrate this, the Root-32
Half-Integer sequence, see Table 3, has from (46)
Sh = t = k/2
and
S_h = 2 + d = 8.
..., 33461,
985, 29, 1,
5,
169,
5741, ...
(53)
which are every fourth Pell number as expected beginning with P5 = 29 and
P*/2 = 41. Proofs of statements above follow easily from (39) or (43). So
(23a, b) are every fourth term of the t = /(/5 - 2) sequences. The "F" sequence to 3 decimals is ..., 1.236, -0.486, 1, 0, 1, 0.485868, 1.236, 1.086,
1.764, 1.943, ... and clearly the ratios
1.236:1:1.764 = 3 + v:2 + v:3 + 2v = 20:0 + 1 : 0 + 3
show that every fourth term of "F" gives (23a) of [7] and equivalently a bisection of Table 1. The general recurrence of these bisections is
Bt,n
where P*
rence
+1
P't,>+Bj,n
P>t,n-1
= 34 for (52) and (53). A bisected t-Fib sequence has the recurPn+2
= (t2 + 2)P - P_ 2
or
&2Pn =
t2Pn.
Indeed the recurrence's middle term for ^-sectioning has a coefficient given
by the 777th rising diagonal of the Lucas triangle. The bi-bisection case is
Pn+1+ = (ah + ka2b + 2b2)Pn - bhPn.hi
and so on.
We come now to what I regard as the most important property of these sequences. The Self-Lucas property, (14) of [6.], remains unchanged in this
generalization, namely
(54)
(2 r + 1 + gr.j)/d
= (-l)r'\.P
and
(.i* +
1+S*.1)/d=(-lf
+h
S*1.,
1978]
537
where d = (a - 3) and (a + g) = j = t .
Now (54) may be proven from (43).
Note that terms with subscripts (r + 1) , (r - 1) and -r all belong to the
same bisection of t or . This is obvious since (r + 1) - (r - 1) = 2 and
(r - 1) - (-p) = 2k - (N + 1) which is also even. Taking ratios of (54) one
may form the triplet rule, for both t an<3 ?,
(55)
( r + Q>v + 2)/(Q>r
+ r _ 2 ) = S-r-i/S-r+i.
Now this is readily illustrated when J is an integer so consider (50). Obviously (10 + 65)/3 = 25 and (lit; + 76y)/3 = 29i> and both are members of (50).
Again from (52) one has (985 + 29)/6 = 169. From Table 1, I illustrate (55)
by ((80 - 2) + (30 - l))/((30 - 1) + (0 - 1)) = (0 + 5)/2 by using the trick
(0 + 2) = 0v.
But this last ratio, (0 + 5):2 is Chiron:Neptune.
I introduce a new operator capital lambda, A:
(56)
A E I + E,
where E and I are the forward shift and identity operators and, therefore,
(A - V) E (E + E~l).
Then the Self-Lucas property may be written
(57)
(A - V)$t,r = d(-l)"*St,-r,
where 2" is another notation for and *, respectively. We may also write
KBn - dB-n or -dB-n depending upon which bisection of Qt or * is being considered. The Self-Lucas property does not hold for F-like sequences: P t ) 1 =
1 = P t ) _ 1 ? or L-like sequences: P%t 1 = 1 = -P%t_1.
In this sense, the HalfInteger sequences are more important.
How are (54) and (57) to be interpreted geometrically?
Let alternate
terms of the bisection be made negative, then the recurrence is
(58)
6zNn = ~(t2
+ 4)tfn.
Further let the terms of Nn be reciprocal periods of planets and let a minus
sign mean retrograde (backward) motion. Then the Self-Lucas property may be
written
(59)
ANn =
-dN-n,
which in words says that the set of synodic (apparent) frequencies of a collection of alternately pro- and retrograde planets are simply proportional to
the negative of the sidereal (real) frequencies in reverse order.
3.
Why are the planetary frequencies not Folded or Half-Integer Pell sequences? The limiting distance ratio would be 3.2386766. One solution to this is
point (x) of [7], namely to bring the planets closer to each other, thereby
minimizing their potential energy which is negative; that is,
(60)
m^YJ{GMJmi/d2av),
i *J
where Mj is Jupiter's mass, m^ the mass of any planet except Jupiter, and dav
is a time-weighted distance from Jupiter. The Pell and, indeed, all t -Fib
sequences satisfy point (ix), the avoidance of low-order commensurabilities,
since gcd(Pt,n+i, Pt,n)
= 1 for all integers t and n. This can also be seen
by noting that the continued fraction of the roots of any t-Fib recurrence
consists of repeated (l/t) T s, so no one convergent is a great deal better
than another. The sequence, 11, 12, 16, 24, 38, ..., is an example who---
538
[Dec.
gcd(12, 16) ^ 1. Given that the total number of planets is a constant, then
minimization of the cumulative perturbation frequencies (synodic) occurs as
t becomes small (point xii). Of course, as t becomes small the average distance becomes smaller and the average perturbation force becomes very large.
k.
FINALE
The logic [7] of my rule suggests that other civilizations may be signaling us in binary code with 1/02 = 0.01100,00111,00100,01000,01100, .... But
let me assure you that my getting into Bode's Rule was not a matter of choice.
Its rewards, though, have been a large number of empyreal highs, some over
ideas I later rejected; but now I am glad to be through with this whirlpool.
Finally, we all know that the idea of the "music of the spheres" which dates
back to Eudoxus is poetic license, nonetheless I could not help noting that
though most of the "notes" in my scale are cacophpnous, the first note, 2+v,
corresponds to C-sharp two octaves high, since (2+v)
= 225 .
REFERENCES
1. Johannes Kepler, Mysterium Cosmographicum
(Tubingen, 1597).
2. Johann Elert Bode, Vom dem neuen3
zwischen
Mars und Jupiter,,
entdeckten
achten Haupt-Planeten
des Sonnensystems
(Berlin, 1802).
3. M. A. Blagg, Mon. Not. Roy. Soc. Astr.,
Vol. 73 (1913), pp. 414-422.
4. W. E. Greig, The Distance
Rule in Planetary
and Satellite
Systems
(1974),
55 pp. (private circulation).
5. W. E. Greig, "The Reciprocal Period Rule," Bull. Amer. Astron.
Soc. , Vol.
7 (1975), p. 449.
6. W. E. Greig, "Bodefs Rule and Folded Sequences," The Fibonacci
Quarterly,
Vol. 14 (1976), p. 129.
7. W. E. Greig, "The Reciprocal Period Law," The Fibonacci
Quarterly,
Vol.
15 (1977), p. 17.
8. W. E. Greig, "Abstract of Reference," Zentralblatt
fur Matematik,
1978 (in
press).
9. I. J. Good, "A Subjective Evaluation," J. Amer. Statistical
Assn. , Vol. 64
(1969), pp. 23-49.
10. J. G. Hills, Nature,
Vol. 225 (1970), p. 840.
11. Owen Gingerich, Icarus,
Vol. 11 (1969), pp. 111-113.
12. J.L.Brady, Astron.
Soc. Pacific,
Vol. 84 (1972), pp. 314-322.
13. W. E. Greig, letter to H. W. Gould, dated 30 Aug. 1973.
14. M. G. Monzingo, "Extending the Fibonacci Numbers," The Fibonacci
Quarterly,
Vol. 12 (1974), p. 292.
15. H.W.Gould, "Cale's Rule," Proc. W. Va. Acad. Sci. , Vol. 37 (1965), pp.
243-257.
16. Anon., "Rings of Uranus," Science News,Vol.
Ill (1976), pp. 180, 245,52.
17. C. T. Kowal, "Chiron the Centaur," Science News, Vol. 112 (1977), pp. 311,
388, 409, 419.
18. J. L. Friend, Sky and Telescope,
Vol. 32 (1963), p. 329.
19. D. E. Thomsen, "Toro! Toro!," Science
News, Vol. 103 (1972), pp. 186, 173.
20. W. E. Greig, "On Sums of Fibonacci-Type Reciprocals," The Fibonacci
Quarterly,
Vol. 15 (1977), pp. 356-358.
21. Anon., "Jupiter!s 13th Moon," Science
News, Vol. 109 (1974), pp. 195 and
367.
22. W. E. Greig, Bull. Amer. Astron.
Soc,
Vol. 7 (1974) , p. 337.
23. N. A. Pananides, Introductory
Astronomy
(Reading: Addison-Wesley, 1973).
1978]
24.
25.
26.
27.
28.
539
of California
at Irvine,
Irvine,
CA 92717
Sn
= Sn-1
+ Sn-2,
Si
= 0, S2
= 1.
^)
The probability of
Z2 "n^n>
and using the generating function (1 - x - x2) ~l for the Fibonacci numbers,
yields g{x) = x2/(1 - x - x2).
Hence, the expected number of trials is
J2nSn/2n
= (1/2)^(1/2) = 6.
n=l
G(x)
=Y,Fnxn
k
is given by G(x)
r
(3)
(px)
1 - x + -^
U-px)
^-(px)k
+1
is
=r
'
P\
5^0
[Dec.
q = 1 - p.
= qPn.l9
Pr(HT)Pr(On_2)
qpPn.2.
Continuing in this way, the last possibility to be considered is that a sequence in 0n begins with k - 1 #Ts followed by a T and then by a sequence in
(9n_k, the probability of which is qpk~1Pn_^.
Hence, the recursion:
(4)
Pn = qPn.
-- + ^ p ^ - 1 ? ^ ,
+ qpPn.2+
(Note that the probability of achieving k heads with k tosses is pk, while
with less than k tosses it is impossible.) The technique to find the generating function for the Fibonacci numbers applies to finding
00
G(X) = J^ Pn*"'
Consider
H{x) = ) Pn + 1xn;
n = k
then
xH(x) = Pn + 1xn
+1
n
= ^P]T
nx
Hence,
- Pkxk
(px)h
= G(x) -
H(x) = [G(x) -
(px) ]/x.
HM = p n
k
+ 1tf
and r e c a l l i n g t h a t Pj = 0 f o r j < k,
=
'
1978]
5^1
+ E(Z)
= E(Y)
+ tf(Z|y = l)Pr(Y
= E(Y)
+ J2E(Z\Y
= 1) + #(Z|j = 2)Pr(J = 2) + . . .
n)Pr(Y
= n) = E(Y)
+ Y^rLE{Xx)Pr{Y
n=l
= E(Y)
= n)
n=l
E{XX)E{Y).
But E(Y) = the expected number of tosses to observe a head = 1/p, and E(X{)
E7;,.!. Thus S7?, = 1/p + (l/p)Ek_,
which yields (3).
REFERENCE
1.
2.
L. E. Dickson, History
of the Theory of Numbers, Vol. I (1919; Chelsea reprint 1966).
W. Feller, Introduction
to Probability
Theory and Its Applications,
Vol.
I (New York: John Wiley & Sons, 1968).
University
IN
INITIAL TERM
47702
un
+k
= alun+k_l
+ .- - +
linear
...
akun9
where a19 a2, . . . , ak are rational integers and um\un whenever m\n.
Some divisibility sequences satisfy a stronger divisibility property, expressible in
terms of greatest common divisors as follows:
(um,
un)
= U(m^ n )
un+M
E un
(mod m)
n \J 9 I, ...
5^2
[Dec.
U ;
+ U
ir+
U V
K(t)
(1 - Xit)a
- x2t)
-. (1 -
Lmma 1:
Vh.00^:
xkt)'
Let q = x x
for w = 1, 2, ... .
Since um\uQm
q= p ^ p s z '
Then, since um\qmuQ
?""" , m = 0, 1, 2, .... .
az(j)
(l if K | = 1
a(2) =
\a%U
a
(j)
, if | cf | > 1
~ 1) if \z\ l i - 1
:1 <. r < j -
l\
(if K l > j - 1.
Thus, either the sequence a (l), a (2), a (3), ... is strictly increasing and
unbounded, or else it is strictly increasing up to some point and constant
thereafter, or else it is the constant sequence 1, 1, ... .
Lemma 2:
1978]
Lemma, 4:
Vtiooi:
(uW9 ua) = U(Wta)9
so mn{iZtW,
by Lemma 3, so iZi(<Wya)
= ^i,w
Lmma. 5: Suppose 1 <_ <_'z; and j >_ 2.
positive integer. Then
(ti,b
a +
i,
t [ , k + 2j
Since az(r)
ig,a}
Suppose a = az(j)
^A.fca + a-l)
(^Z,ly
\a9
= a, then
Now
= ii,(w,a)
5^3
iitWiita,
>_ 2 and H s a
^ 1, 2 > s ^ , a - l )
Vh.00^1
+ l ) ^l,bL
+ 2 J
^ , 2 L - l ) = (^Jl, 1 ^ , 2 J J
^Jt,2>L-l)
(^,2L+1 s J ^ , 3L-1 )
(^.(fc-DL + l ,
---J ^,L-l)5
1 , 2 , ..., fc.
For
Pl
'p2
...
pt'
. . . UQet
5hk
[Dec.
Finally, u0 = up,
and un\u0
for n = 0, 1, ... .
VK.00^'. By Lemma 7 and the fact that {un} is a kth order recurrent sequence, the sequence {un) is periodic with period M. Letting p be the fundamental period, we now show that the denominator of the generating function
H(t)/K(t)
must be of the form 1 - tp:
Y^r
+ + Up.^"1
= u0 + uxt
= u0(l
+ t
= (u0 + uYt
+ t
2p
+ uQtp
+ u1tQ
p
+ ...) + uxt(X
+ t
p 1
+1
+ t
+ t
2p
2p
+
+ ...) +
+ )
.
A. "0
= tfW 8l U
fll
. . . U et
un\up.
akup.
REFERENCES
1.
2.
W. HERGET
Braunschweig,
Universitat,
Fed.
Rep.
Germany
(1)
*-d^2( n :>
(m>o)
i =0
(2)
(-iy
(r.)^i^
= oiodp",
with u: = pe~1{p
- 1), where r >_, e >_ 1, m > re, p prime such that p - 1^/7
With r = 1 we get, in particular
1978]
(3)
mod p e ,
=
m+y6-1^
5^5
- 1)
p - i n t e g r a l =>p -
l\li.
Since B0 = 1, B1 = - 1 / 2 and B 2 i + 1 = 0 f o r i e N, we g e t
(4)
be the smallest integer greater than 1 with the following prop3 77?0 V/c, m >_ m0:
(5)
^ Bm ?z-integral f\Bk
E 5 m mod n ) .
and
J:
of
We will prove
F(2e) = F(3e) = 2
(c)
L(pe) = pe(p
(d)
V(p )
<.e + 1.
(6)
+ p e - 1 ( p - 1)
ni
1
1 j _ moc | pe. hence,
m + tpe (p - 1)
k = m + spe(p - l)Ap --\mhm
> e =^Bk ~ Bm mod p e .
546
[Dec.
Then, Bk = Bs l
k: = s + (p - l ) p * = (2 + (<7p2 + 3 ) ( p - D ) p * = 2 + t(p - 1) ,
where
t : = 2 r _ ^ + (qp* + 3 ) p i e ff;
P
t h e n by (3) w i t h e = 1 and m = 2 ,
52 _
B2 + ( p - l )
2 + (p - T )
_ Bk
_
=
" = x
m o d p
'
L(pe)
= pe(p - 1)
by (6).
Now we may improve this last inequality as follows:
Th&QJim 2:
1. V(p) = 2 for p prime.
2. Let p be a prime, p > 3 and e e (2, 4, 6, . . . } . Then,
(a) Be f 0 mod p A p - l f e =>V(pe) = e + 1.
(b) fc maximal such t h a t
V0<.i<.k:
3.
= e - 1 - 2k.
V O i <_ k: {Be_l_2i
=>F(pe) = e - 2 - 2k.
VKOOl: By Theorem 1(d), we have V(p) <. 2. But V(p) < 2 is impossible
since B1 - -1/2 ^ 0 mod p and S 1+ L ( P ) = 0, thus V(p) = 2.
For the proof of the other assertions we note that [4, p. 321, Cor.]:
Y,(-iy(l)Bm + iv(l
i =0
fc
1978]
5^7
F(p ) e
- 3 <=> (p\eAp
^ e
l|e
l|e
- l\e
= 2 p mod p(p
- 2 ) V ( p - l\eAp3\e
- 1) V e = 2 - 2p
- 2)
3
mod p 3 ( p - 1)
If 2 | e
>. e - 5 .
then
F ( p ) = e <=^p2\e
- l A p - \\e
- 1
F ( p ) e - 2 < = p | e - l V p - l\e
e
- 1
<=> e = 2p
h
= -2p
- 3) V (p - l\e
2
+ 1 mod p ( p -
- hph\e
- 3)
l)Ve
+ 3 mod p ^ ( p - 1)
V(pe)
V(pe)
>_ e -
6.
S 2 i E 0 mod p e =^>pe|2i.
The assertions of Corollary 1 with "<==" are also valid for any irregular
prime p.
By Corollary 1, you may see that only for greater integers pe, the value
e
V(p ) differs from e and e + 1, respectively. We get
CoKolLcUty 2: For prime p , p > 3 , l e t e
2 p 2 + l , e5 = 252, e6 = 1253. Then we have
(a)
F(peO < e
- i ,
= p - l ,
e2 = p ,
e3 = 2p,
e^ =
{l, . . . , 4}.
548
2.
[Dec.
that e < et and V(pe) < . - , where i e {-2, 3, 4}? (By the computational results in [5] we may conclude that this does not happen when
p < 30,000.)
Is there a power n = pe of some irregular prime such that
V(pe)
<e-5?
VAJIGJL RojncUik: Professor L. Carlitz and Jack Levine in [3] asked similar
questions about Euler numbers and polynomials. Analogous results about the
periodicity of the sequence of the Bernoulli polynomials reduced modulo n and
the polynomial functions over Z generated by the Bernoulli polynomials will
be derived in a later paper.
REFERENCES
1.
2.
3.
4.
5.
HANSRAJ GUPTA
University,
Chandigarh,
1.
India
INTRODUCTION
T
T
e
o
)
)
c)
<)
1978]
5^9
is four deep or is of depth four. It is clear that a graph of depth k cannot have less than k nodes.
Denote by x the number of nodes on the horizontal, and by yi the number
of those on the vertical section of the ith. right angle, starting with the
outermost right angle as the first. Then, the partition can be very conveniently represented by the 2 x k matrix:
^2
^3
ys
hi
or simply by
.yi
Evidently, we must have
(1.1)
xi
>_ xi
+1,
+1
_> yi
+1
+ 1,
i <_ k - 1.
ykl
y^
yi<.
= I** -
yih>
which we shall call the rank-vector both of the graph and of the partition it
represents.
The number of nodes in the graph is given by
^
(1.3)
+ yi
(xi
- 1).
2"
.9
1.
1].
yi
2.
Xi >_ max(xi+ 1 + 1, xi
+ ri
+1
+1
+ 1).
>_ max(rk
+ 1 , 1).
- ri
Y^Xi = Un + k
i=i
+Y,rA.
i=i
+1
+ 1.
Hence, for
550
[Dec.
+ E rt
i= 1
Moreover, i
There are an infinite number of Ferrars graphs which have the same rankvector. All such graphs have the same depth but not the same number of nodes
necessarily.
ThJL0h&m\
Among the graphs with the same rank-vector,
with the least number of nodes.
VJWOfa: Using the equality sign in place of the sign _>_ in (2.2) and ( 2 . 1 ) ,
we obtain the least value of each of the x 's, i k.
(1.3) and (1.4) then
give n0 that is the least n for which a graph with the given rank-vector e x ists. This proves the theorem.
Incidentally, w e also get the unique partition with the given rank-vector
and the least number of nodes. W e call this unique partition the basis of
the given rank-vector.
Example:
Let us find the basis of the rank-vector [-2 3 0 1 ] . With
the equality sign in place of the of the inequality sign, (2.2) gives xh = 2.
With the equality sign in place of _> , (2.1) n o w gives, in succession,
= 3, x2 - 7, and x1 = 8.
x3
yh
= 10.
2"
.10
1_
If
+ K
y<i
lyi
according as h is positive or negative.
(b)
If hi >_ h2 _> ... >_ hk >_ 0 are integers, then the graphs of
XJS
+ 7z,-l
and
.Mi
Mi
have the same rank-vector.
7*1
1978]
k.
Let
551
in the basis. For our n to have any partitions with the given rank-vector,
it is necessary that n has the same parity as n0 and n>_nQ. Assume that this
Write
m = j(n
- n0).
m = hl + h2 + + hk,
with hi >_ h2 >_
L b>k O J t> e a
y s u c n p a r t i t i o n of m.
+ hi
(4.1)
2/t
+ /z,
" 7 6 4 1 "
_10
2.
so that n 0 = 33 and m = 5.
The partitions of 5 into at most 4 parts are:
5; 4 + 1 , 3 + 2 ; 3 + 1 + 1 ,
2+2+1;
2+1+1+1.
[12
Ll5
l2 > 1"11
J-4
"10
-13
2_
10
2,
13
l"
2.
12
J-2
552
[Dec.
We show that these ranks are bounded both above and below. Since x1 <. a,
the number of rows a partition of n can occupy is not less than u, where
u - 1 < n/a <_ u.
Hence, none of the ranks can exceed (a - u).
Similarly, none of the ranks can fall short of (v - b), where
v - 1 < n/b <_ v.
Of course, for n to have a partition of said type, it is necessary to have
n <_ ab,
REFERENCE
1.
J.
University,
Candigarh,
India
(n;p) for [J .
The two relations
(1.1)
(1.2)
S(n,r) = Y,(n;j),
J
where j runs over all nonnegative integers which are E p (mod 5 ) .
As a consequence of this definition and the relations (1.1) and (1.2) we
have
(1.4)
S(n,r)
(1.5)
S(n9r
2.
=S(n,n
- r), and
- 1) + S(n,r)
= S(n + l,r).
(2.1)
Fl = 1 = F29 and
(2.2)
1978]
553
n =H
(-l)'(n
- 1 - 5j)/2]);
- l;[(n
(2.4)
(-DJ'(n;[(n
Fn = Y,
- 1 - 5j)/2]);
or
2m + 1 (mod 10)
Fn = S(n - l9m)
(2.6)
Fn = S(n9m)
- S(n - l,m
- 2);
- 1).
- S(n9m
We first assert that (2.5) and (2.6) are equivalent and prove the assertion as follows:
For any integer j , we have
(n;m + 5j) - (n - ;m + 5j) = (n - l;m + 5j - 1).
Also
(n;m - 1 + 5j) - (n - l;m - 2 + 5j)
= (n;n
- m + 1 - 5j) - (n - l;n
- m + 1 -
5j)
- m - 5j)
= (n - l;n
= (n - ;m + 5j - 1).
Hence, letting J vary suitably, we get
S(n,m)
- Sin
- l,m)
- 1) - S(n - l9m - 2 ) ,
= S(n,m
Ft = S(t,m)
(3.3)
- 1)
- S(t9m
= S(t,m)
= S(t9m
- S(t9m
- 2) for t odd,
+ 1) - S(t9m
- 1) f o r t e v e n .
S(t9m)
= S(t9t
- m) = S(t910k
+ m + 1) = S(t9m
+ 1).
I f t i s e v e n , l e t t = 10/c + 2TT7 + 2 ; t h e n
(3.5)
S(t9m
- 1) = S(t9t
-777 + 1) = S(t910k
+ m + 3) = S(t9m
- 2);
so t h a t
(3.6)
Ft + 1 = S(t9m
+ 1) - S(t9m
- 2) f o r t odd as w e l l as t e v e n .
554
[Dec.
From ( 3 . 1 ) and ( 3 . 6 ) , we g e t
Ft + Ft + 1 = {S(t,m)
= S(t
+ S(t9m
+ 1)} - {S(t,m
+ 1,777 + 1) - S(t
- 1) + S(t,m
2)}
+ 1,777 - 1 ) .
Thus,
F
t + 2 = S(t
+ 1,777 + 1) - S(t
+ 1,777 - 1 ) .
and so also
[p2/2] = 0 (mod 5 ) .
Hence,
= S(p 2 , 0) - S(p 2 , 4 ) ,
Fp2
V +i
and t h e r e f o r e ,
5(
P2'
0)
5(
"
P2'
3);
F p , . ! = 5 ( p 2 , 4) - 5 ( p 2 , 3 ) .
(ii)
and so a l s o
is
[ p 2 / 2 ] E 4 (mod 5 ) .
Hence
F p 2 = 5 ( p 2 , 4) - 5 ( p 2 , 3 ) ,
*Wi
and
= (p 2 > 4) - ^ ( p 2 ,
);
therefore,
= 5(
^2'
3)
Lmmcii
For l ^ p
"
5(
2)
P2'
All that we need now for our purpose is the
*V-i
-l,
(p ; h) E ( - l ^ ' V A (mod p 2 ) .
P/LOOfJ:
We have
2.
*-P '
(p
fe)
n>
j . P2 - 1 . P2 - 2
h
1
2
Since f o r 1 <. r ft - 1 ,
^-f^t h e lemma f o l l o w s
= -1 (mod p 2 )
immediately.
p 2 - ft + 1
h - 1
1978]
555
(4.1)
E 0 (mod p 2 ) ;
otherwise,
(p2;h)/p
(4.2)
5 (-D^p/fc (mod p ) .
We have, of course,
(p2;p2).
(p2;0) = 1 -
(4.3)
(4.4)
S(p\m)
On the right of the sigma in (4.4), we need consider only those nonnegative
values of j for which
m + 5j <_ p 2 and 77? + 5j E 0 (mod p) ;
(ii)
when 77? = 0 , we h a v e ,
S(p2,0)
(4.5)
so t h a t
(4
*(p2,0) -
.6)
(p2;5j)
- 1 =
(mod p 2 ) ,
~1
E^
(-l^-VSj
(mod p ) ,
Thus
^121 ~ 1
1
1
1
1
11
^ - ^ + t - i E
(mod 1 1 ) .
Therefore,
Fl21
E 89 (mod 1 2 1 ) .
REFERENCE
1.
G. E . A n d r e w s , The Fibonacci
Quarterly,
West Virginia
University,
W. GOULD
Morgantown,
W.Va.
26506
Operational formulas can play a fascinating role in finding transformations and sums of series. For instance, by using the differential operator
D (=d/dx) we can transform
(1)
2 , * *
k = 0
into
YJkxk~1
ri
(1 - x)2
556
[Dec.
The operator 0 = xD is even more interesting. It has the basic property that
0pxk = kvxk, so that (1) can be transformed into
<2>
S* v - 6 '{5-M-
(4)
k\S(p,
k) = A*0* = (-D""J'(5)jp5
X J /
j-o
then series (2) can be found in closed form for it is trivial to find the
higher derivatives needed in (3). The result is a very old and well-known
formula. In [7] is given an extension of (3) applied to generalized Hermite
polynomials. There are numerous similar generalized expansions involving the
D operator. Here we propose to examine some rather unusual variations that
are not too well known, and which have applications to Fibonacci numbers
among other things.
We shall need several, other well-known operational formulas whose proofs
involve some calculus and/or mathematical induction, and we tabulate these
below:
(5)
9 = Dz, where x = ez 9
(6)
(7)
e" = D?,
xD-nl(Dz),
e D = 1 + A = E,
where
A/(a:) = fix + 1) - f{x)
More generally
(9)
etD* = f(x + t) =
Efix).
x ,h
The ^-operator
(10)
fiqx)
= Qf(x),
where Q = q* .
This was used, e.g., in [10], and is very convenient when working with basic
hypergeometric series.
In the references at the end are several papers, viz. [1], [2], [4], [5],
from the older literature where properties of a great number of familiar and
unfamiliar operators were developed. The master calculator was almost certainly George Boole. The English literature for the period from about 1830
to 1890 is especially rich in papers on unusual operators.
In [1], Boole gave the pair of very remarkable operational expansions
(11)
f(x + 0riD))uix)
e^D)fix)e-*{D)uix),
f(D + 0rix))uix)
e-4(x)f(D)e0(xy>uix),
and
(12)
1978]
557
*L-t,WLK
(13)
k=o
'
(14)
= ZE^E(-i)- (^ k /(x)
n =0 j =0
k =0
where Cj are Stirling numbers of the first kind, i.e., coefficients in the
expansion of a binomial coefficient:
j =0
DDf={(D-l)+l}Df=Y,(Dn)(D-l)nf
(16)
n =0
=Z
= ESt(-D"-k(?)o^).
n =0
' k =0
n=0
" k=0
k=0
k =0
= eax.
Then
558
= zCCaa =
SERIES
[Dec.
aae<ZX
>
n=0
so that w e have the attractive formula
DDeax
(17)
aaeax,
(18)
"
n- l
Y,kkek*tk
n-l
= ^tkDDekx
n-l
= DD ^(te*)k
= DD {*"*"* " H .
In particular,
(19)
y>*** = W ^ \\
All that would be necessary to sum (19) would be to find a different method
of attaching a meaning to the right-hand member.
For a Fibonacci-Lucas application, recall the general Lucas function
Ln = Ln(a,
Then
(20)
b) = an + bn.
>**%* = W g p ,+ V V
fc=o
* ap - 1
and for the generalized Fibonacci function
Fn = Fn(a,
bpex
then
(21)
- 1
aP"gn* _
(a - b)^kkek*Fpk
= W -
p x
ae
- 1
2)Pn5nX - 1
bpex
(22)
J2kkek*u%lLk,
k=0
(23)
LLnn = e-xDDeLnX
F*" = e-xDDeF"x
(24)
!>*"" =
rc = 0
so that we
e-'DD^ttneL-x,
n 0
(25)
S(t,
u) = ] tnuLn ,
1978]
559
and
(26)
T(t9 u) = ]T tnuFn .
n= 0
These a r e o f f e r e d a s r e s e a r c h p r o j e c t s ; t h e a u t h o r would b e i n t e r e s t e d i n
h e a r i n g of any s u c c e s s by o t h e r s . DUS\U=1 and DUT\U=1 a r e known.
The o p e r a t o r 0 e may be c o n s i d e r e d f i n a l l y .
We f i n d
"
-tfi*- ( - i ) B ~*(*) e ^
n=0
k=0
so we have
(27)
=Y,^D^(-l)n-k(n)&Kf(x).
Bf(x)
n
=0
"
k=0
^'
Let / ( x ) = x , t h e n
)exp
QBxp =
ppxp.
(29)
**** = eeJ>* = e 9 ! ^ } .
We wish to remark that even stranger formulas have been published. Cayley [4], [5] expressed the Lagrange series inversion formula in the most
curious operational form
(30)
F(x) =
(DJ^'^F'^e^^)},
By differentiation,
Cayley says these are well known, and goes on to write similar formulas for
functions of several variables.
Bronwin [2] writes
(32)
f(a + x) =
D!{f(a)ex)
560
[Dec.
from page
497)
1978]
561
EVERY INTEGER IS A
State
HARVEY J. HIND IN
College,
Stony Brook,
NY 11790
In this note we prove the following: Every positive integer n can be expressed in an infinite number of ways as a quadratic function for each of the
infinite number of figurate number types.
The nth figurate r-sided number pr is given by
(1)
p* = n((r
- 2)n - r + 4)/2,
pvsn
= sn((r
- 2)sn
- r + 4)/2.
However, (2) is a quadratic in n. Solving for n and taking the positive root
yields
(p - 4) + /(r - 4 ) 2 + 8(r - 2)psrn
which allows us to express n as stated above. A special case of (3) for pentagonal numbers (r = 5) was obtained by Hansen [1].
REFERENCE
1.
Quarterly,
University
New Mexico
87131
Proposed
by Herta
T. Freitag,
Roanoke,
VA.
+ T
+ i
Show that
l + 3 + 5 2 + + (2n - l ) 2
s + '" + 2n-i
Proposed
by Gregory
Wulczyn,
Bucknell
University,
Lewisburg,
PA.
Find the complete solution, with two arbitrary constants, of the difference equation
(n2 + 3n + 3)Un + 2 - 2(n2 + n + l)Z7w+ i + (n2 - n + l)Un = 0.
B-390
Proposed
by V. E. Hoggatt,
State
University,
San Jose,
CA
B~391
Proposed
by M. Wachtel,
Zurich,
Switzerland.
B-392
Proposed
by Phil
Mana, Albuquerque,
Proposed
by V. E. Hoggatt,
Let Tn = ( n 2 X ) '
hFn + kLn.
Jr., San Jose
o = X> Pn = TiT2
"Tn
f o r i n t e g e r s k and n w i t h 0 . k n .
fnl
lk\
NM.
State
^r
Show t h a t
1
(n\(n + \
n - k + l U / U + 1 /
562
University,
San
n > 0 , and [ j ] =
Jose,CA.
Pn/PkPn.k
Dec. 1978
563
SOLUTIONS
1NCONTIGUOUS ZERO DIGITS
B-364
Proposed
by George Berzsenyi,
Lamar University,
Beaumont,
TX.
Find and prove a formula for the number R(n) of positive integers less
than 2 whose base 2 representations contain no consecutive O's. (Here n is
a positive integer.)
Solution
by C. B.A.
Peck,
State
College,
PA.
Let Sn be the number of integers 777 with 2n~ <_ m < 2n and having a binary
representation 5(777) with no consecutive pair of 0Ts. Clearly Sn = Rn - i?n_i
for n > 1 and S1 = R. Also,
Sn = n_i + 5 n _ 2 for w > 2,
since Sn_1 counts the desired 77? for which B(m) starts with 11 and Sn_2 counts
the desired 77? for which 5(777) starts with 101. It follows inductively that
S
= Fn+1>
and
then
Rn
= S,
+ S2
+ + Sn
= F2 + F3 + . . + Fn + 1 = Fn
+3
- 2.
Prielipp,
proposer.
CONGRUENT TO A G.P.
B-365
Proposed
by Phil
Mana, Albuquerque,
NM
Show that there is a unique integer m > 1 for which integers a and r exist with Ln E avn (mod 777) for all integers n >_ 0. Also, show that no such 777
exists for the Fibonacci numbers.
Solution
by Graham Lord,
Universite
Laval,
Quebec.
2 5
Since 7 = LiiLl = a r
E LZL3 = 12 (mod 777), then 77? divides 5, hence 777 = 5.
Furthermore, a = ar E L0 = 2 (mod 5). And finally, or2 =L2 = L1 + LQ E av +
a (mod 5) together with a E 2 (mod 5) implies v1 E v + 1 (mod 5), i.e., v E 3
(mod 5). In all, 777 = 5, and a and v can be taken equal to 2 and 3, respectively. Note for any n >_ 1, Ln + 1 = Ln + Ln_1 E ap n + a p n _ 1 E ar n+ 1 (mod 5).
For the Fibonacci numbers, if 777 were to exist, then
3 = FxFh
Bob Prielipp,
A. G.
Shannon,
LUCAS CONGRUENCE
B-366
Proposed
Slippery
E LhLk
of Tennessee,
Rock, PA.
Knoxville,
TN and
(mod 5) when i + j - h + k.
Solution
by Paul S. Bruckman, Concord,
lege,
Clarion,
PA
(independently)
Using the result of B-365,
Clarion
State
Col-
564
L.L-
= 2 3i
- LhLk
+j
[Dec.
- 2 3h + k = 0 (mod 5 ) ,
since i + j = h + k.
Also solved by George Berzsenyi, Herta T. Freitag,
Graham Lord, T.
Bob Priellpp,
A. G. Shannon, Gregory Wulczyn, and the proposer.
Ponnudurai,
ROUNDING DOWN
B-367
Dakota State
University,
by George Berzsenyi,
T
In view of Binet s
Lamar University,
_ 2)2n~1
7
a - b
Prove
Beaumont, TX.
2n
a2n-i
l
and n >. 1.
formula,
ccFr>n_ T F 0y. = a
zn
Brookings,
Zn
b2n
7
a - b
-
Similarly,
2p
^ '
Since - 1 < b =
follow.
j?
"
+1
<
~2&
~ O
a - b
"
t h a t 0 < -b2n~x
implies
a - b
< 1,
the
" ~b
2n-l
desired
results
Priellpp,
and h(n)
Roanoke, VA.
such that
^2^FiLn-i=
g(n)Fn + h(n)Ln
i-l
by Sahib Singh,
Clarion
State
College,
Clarion,
PA.
Let An = 2_\^i^n-i'
i=i
and
h(n)
= n/10
so t h a t
(5n2 + l\)Fn + nLn = 0 (mod 10).
This also gives us nLn - Fn (mod 5 ) .
Also solved
poser.
by Paul S. Bruckman,
and the
pro-
1978]
565
B-369
Proposed
by George Berzsenyi,
Lamar University,
Beaumont,
TX.
-^2n+3J
\^Zn+l9
-^2n+5/
by Graham Lord,
Universite
Laval,
Quebec.
2
+ 5 = L2 n + 2 5
and
L
2n + l L 2n+5 +
2r7 + 3 '
lL6/c + i a n d 2 lL6fe+5
4= L L
3 l 6fc+3
4|L2k
Of these, (l) is somewhat well known and the latter three are consequences of the results in "A Note on Fibonacci Numbers," The Fibonacci
Quarterly,
Vol. 2, No. 1 (February 1964), pp. 15-28, by L. Carlitz.
By (T) and \2) there is exactly one even number, L6fc+3 > in the set Sn,
n >. 0. So if {2} U 5 n has the desired property, then zLsk+3 + 5 will be an
odd square and thus congruent to 1 modulo 8. This implies that 2, if it exists, is odd.
Now the other two members of Sn are either:
(a)
L6k_i,
(b) L
Lsk+1;
s k + 5
L6k+7;
o r (c) L6k+l9
LBk+5-
Each of these is odd by (l), and hence the sum of 5 and any one of them
multiplied by z will equal an even square. Thus, in case (a) [and similarly
in case (b)]:
zLsk_1
+ 5 = 0
(mod 4 ) ;
i. e. ,
zL&k
= z(L6k
+1
= (zLsk+5
+ 5) - (zL6k+1
+ 5) = 0 (mod 4 ) ,
T. Ponnudurai,
Bob
Prielipp,
PA 17745
Proposed
by Gregory
Wulczyn,
Bucknell
University,
Lewisburg,
PA.
Show that:
(a)
(b)
H-291
FkFLsr
Proposed
~ ?U,r
^2r +l^k+8r
(-Dk + 1F2Zr(Fk
by George Berzsenyi,
+ 8p
+h ~ ^k+kr+2'
>
+ 2F k + l f r ) .
Lamar University,
Beaumont,
TX
Prove that there are infinitely many squares which are differences of
consecutive cubes.
H-292
Proposed by F. S. Cater
land, OR.
and J. Daily,
Portland
State
University,
Port-
Find all real numbers v (0,1) for which there exists a one-to-one function fr mapping (0,1) onto (0,1) such that
(1) fr and f~ are infinitely many times differentiable on (0,1), and
(2) the sequence of functions fr, frofP9
frofrofr,
fv^fvfvfv
> '
converges pointwise to v on (0,1).
H-293
Proposed
by Leonard
Carlitz,
Duke University,
Durham, NC.
defined by
'
- z) (k = 0, 1, 2, . . . ) .
'
<x>
Z-4+fc^Jr = E ^ ^ J r 4 ( x - z) (k = 0, 1, 2, . . . ) ,
n =0
where f0(x)
satisfy
n=0
566
Dec. 1978
H-294
Proposed
by Gregory
Wulczyn,
Bucknell
University,
567
Lewisburg,
PA.
Evaluate
^2r+l
^6r+3
-^10r+5
^lhr+7
^ I 8r +9
* hr +2 ^ 1 2 r + 6
^ 2 OP + 1 0 ^ 2 8 r + 1 h^ Z Sv + 1 8
^6r+3
^ 3 6 r + 1 5 ^ 4 2 ^ + 2 1^5 4 r + 2 7
^18r+9
^ 7 8 r + L f - ^ 7 2 i + P + 1 2 ^ 7 i + 0 P + 2 0 ^ 5 6r + 2 8 ^ 7 2 r + 3 6
^ 1 Or + 5 ^ 2 Or + 1 5 ^ 5 Or +2 5 ^ 7 Or + 3 6 ^ 5 Or+1+5
SOLUTIONS
SYMMETRIC SUM
H-272 ( C o r r e c t e d )
Proposed
by Leonard
Carlitz,
Duke University,
Durham,
NC.
Show t h a t
EmcAX^rv^A)
is symmetric in p, q,
Solution
C m (p, (7, r)
v.
by Paul Bruckman,
Concord,
CA.
Define
a)
Clearly,
therefore
/?? - j in
volving p
j =0
".<p-"-')-i;(.^)(5)(j)/(;)fc^))
k=0
j-o
m
m -J
E0t(A)UKi)(.?oc;')/(j)-
j=0 k=0
7 = 0fe= 0
Replacing k by m - k in the
last expression yields
te
last
expression yields:
m
m
'.<*'>-E(.-*)(5)(.?*)(0(*)/(;)
J
= 0fc= j
(-,)/(?)-O/rTO-O^'/fT1)Therefore,
777
*.<*.. w C)E(?)(.f*)?:w)'(J)(5)/(-r1)-
568
<
[Dec.
Identities
D-u'UK^/cr)/-
fc-0
1 \/q+r-m
k
+ k\ I IT -m + k\
) /\
k
)
k=0
<p.,.rt=E(?)(^-)L-J-)r'r+j')/()
3
= Cm(q, p, p ) .
Also solved
by the
Q.E.D.
proposer.
A RAY OF LUCAS
Rock State
College,
Slippery
Rock, PA.
Consider, after Hoggatt and H-257, the array D, indicated below, In which
^2n+i (n = 0, 1, 2, ...) is written in staggered columns:
i.
ii.
iii.
1
4
1
11
4
1
29
11 4
1
76
29 11 4 1
Show that the row sums are L2n + 2 ~ 2;
Show that the rising diagonal sums are F2n+3 ~ 1 where L2n+1
is the
largest element in the sum.
Show that if the columns are multiplied by 1, 2, 3, ... sequentially to
the right then the row sums are L2n+3
- {In + 3).
Solution
by A. G. Shannon,
of Technology,
Australia.
2Ll ^2n-2j +1
^2n+2
~ 2;
J-0
[n/2]
i1,
i;
J-0
n
iii.
(j + l)L2n_2j+1
= L2n + 3 - {In + 3 ) .
j-o
n
n
(^
Z^L2^-2J + 1
J=0
I-** \
j - 0
J2n+2 " ^ O J a
2n-2(j-
1)
s
"
2n-2j)
required.
/ ^ L2n - 2 (j - 1)
j=0
n+1
"
,L2n-2(,i-l)
j=l
1978]
(ii)
2 ^
[n/2]
L
2n-Hj
[n/2]
+ 1 = 2-J {F2n-Hj
J-o
+3 "
2n - 4 j - 1 ) =
j=o
if n\m9
if n\m.
[n/2] + l
F
2n-kj
2-j\^2n-2i+
=0
2n-kj+3
n + 1) = F 2 n + 3 - 1
j-i
+3 ~ X ,
^.Q
'
= F2n+3
in which
569
*Yj 2in+i)-n+i - Ll
)
[from (i)]
3 ~ -^2n-2i+l ~ ^j
- (L2n + 2 - 2) - 2(n + 1)
i =Q
= L2n+,
- 2 - 1 - L 2 n + 2 + 2 - 2(n + 1)
= L2n+3
- (In + 3 ) , as required.
H. Freitag,
B. Prielipp,
Dinh
The'Hung,
H-225
Proposed
by G. A. R. Guillotte,
Quebec,
Canada.
VOLUME INDEX
ADLER, IRVING. "A Simple Continued Fraction Represents a Mediant Nest of Intervals," Vol. 16, No. 6, pp. 527-529.
ANDERSON, 0. D. "Some More Patterns from PascalTs Triangle," Vol. 16, No. 4,
pp. 296-301.
ANDERSON, PETER G. "On the Formula IT = 2Z arcot /2fc+i>"
118.
Vo1
16
>
No
> P-
AUSTIN, RICHARD. "Binary Sequences without Isolated Ones," Vol. 16, No. 1,
pp. 84-86 (co-author, Richard Guy).
BEITER, MARION. "Coefficients of the Cyclotomic Polynomials F3qr(x)9"
16, No. 4, pp. 302-306.
Vol.
BERGUM, G. E. "A Combinatorial Problem Involving Recursive Sequences and Tridiagonal Matrices," Vol. 16, No. 2, pp. 113-118 (co-author, V. E. Hoggatt,
Jr.); "A Family of Tridiagonal Matrices," Vol. 16, No. 3, pp. 285-288 (coauthor, V. E. Hoggatt, Jr). Problem proposed: B-367, Vol. 16, No. 6, p.
564. Problem solved: B-367, Vol. 16, No. 6, p. 564BERNSTEIN, LEON. "An Invariant for Combinatorial Identities," Vol. 16, No. 4,
pp. 354-369.
BERZSENYI, GEORGE. Problems proposed: B-351, Vol. 16, No. 1, p. 91; H-266,
Vol. 16, No. 1, p. 94; B-378, Vol. 16, No. 2, p. 184; B-387, Vol.16, No. 5,
p. 473, H-269, Vol. 16, No. 5, p. 478; B-364, B-369, Vol. 16, No. 6, pp.
563, 565, H-291, Vol. 16, No. 6, p. 566. Problems solved: B-346, B-347,
B-248, B-349, B-350, B-351, Vol. 16, No. 1, pp. 89-91, H-266, Vol. 16, No.
1, p. 94; B-352, B-356, Vol. 16, No. 2, pp. 185, 186; B-358, B-359, B-362,
Vol. 16, No. 5, pp. 474, 475, H-269, Vol. 16, No. 5, p. 478; B-364, B-365,
B-366, B-367, B-369, Vol. 16, No. 6, pp. 563-565.
BEVERAGE, D. Problems Solved: H-264, H-265, Vol. 16, No. 1, p. 93.
Proposed: H-283, Vol. 16, No. 2, p. 188.
Problem
BICKNELL-JOHNSON, MARJORIE. "A Golden Double Crostic," Vol. 16, No. 1, pp.
67-69; "A Golden Double Crostic Solution," Vol. 16, No. 1, p. 83; "A Primer
for the Fibonacci Numbers XVII: Generalized Fibonacci Numbers Satisfying
un+iUn-i
- Un = 1," Vol. 16, No. 2, pp. 130-137 (co-author, V. E. Hoggatt,
Jr.); "Properties of Generating Functions of a Convolution Array," Vol. 16,
No. 4, pp. 289-295 (co-author, V. E. Hoggatt, Jr.); "Convolution Arryas for
Jacobsthal and Fibonacci Polynomials," Vol. 16, No. 5, pp. 385-402 (coauthor, V. E. Hoggatt, Jr.); "Fibonacci Chromotology or How To Paint Your
Rabbit," Vol. 16, No. 5, pp. 426-428; "Divisibility Properties of Polynomials in Pascal's Triangle," Vol. 16, No. 6, pp. 501-513 (co-author, V. E.
Hoggatt, Jr.).
BRADY, WRAY G. "More on Benford's Law," Vol. 16, No. 1, pp. 51-52. Problems
proposed: B-366, Vol. 16, No. 6, p. 563, H-273, Vol. 16, No. 6, p. 568.
Problems solved: B-347, B-348, B-349, Vol. 16, No. 1, pp. 89-91; B-352, B354, B-355, Vol. 16, No. 2, pp. 185-186; H-270, Vol. 16, No. 5, p. 479; B366, Vol. 16, No. 6, p. 563, H-273, Vol. 16, No. 6, p. 568.
BRIDGER, CLYDE A. Problem solved:
570
Dec. 1978
VOLUME INDEX
571
(co-proposer,
COHEN, G. L. "On Odd Perfect Numbers," Vol. 16, No. 6, pp. 523-527.
CULL, PAUL. "Knight's Tour Revisited," Vol. 16, No. 3, pp. 276-284 (co-author,
Jeffery De Curtins).
DAILY, J. Problem Proposed:
Cater).
DE CURTINS, JEFFERY. "Knight's Tour Revisited," Vol. 16, No, 3, pp. 276-284
(co-author, Paul Cull).
572
VOLUME INDEX
[Dec.
Problem Solved:
GLADWIN, A. S. "Expansion of the Fibonacci Numbers Fnm in nth Powers of Fibonacci or Lucas Numbers," Vol. 16, No. 3, pp. 213-215.
GOULD, H. W. "Remark on Problem H-123," Vol. 16, No. 2, p. 189; "Evaluation
of Sums of Convolved Powers Using Stirling and Eulerian Numbers," Vol. 16,
No. 6, pp. 488-497, 560-561; "Operational Formulas for Unusual Fibonacci
Series," Vol. 16, No. 6, pp. 555-560. Problem Proposed: H-282, Vol. 16,
No. 2, p. 188 (co-proposer, W. E. Greig).
GRASSL, RICHARD M. Problems Proposed: B-349, B-350, Vol. 16, No. 1, pp. 9091. Problems Solved: B-349-B-350, Vol. 16, No. 1, pp. 90-91.
1978]
VOLUME INDEX
573
GREGORY, M. B. "Fibonacci Sine Sequences," Vol. 16, No. 2, pp. 119-120 (coauthor, J. M. Metzger).
GREIG, W. E. "On Generalized G^k Numbers," Vol. 16, No. 2, pp. 166-170;
"Folded Sequences and Bode's Problem," Vol. 16, No. 6, pp. 530-539. Problem Proposed: H-282, Vol. 16, No. 2, p. 188 (co-proposer, H. W. Gould).
GUARALDO, ROSALIND. "On the Density of the Image Sets of Certain Arithmetic
FunctionsI," Vol. 16, No. 4, pp. 318-326; "On the Density of the Image
Sets of Certain Arithmetic FunctionsII," Vol. 16, No. 5, pp. 428-434; "On
the Density of the Image Sets of Certain Arithmetic FunctionsIII," Vol.
16, No. 6, pp. 481-488.
GUERIN, E. E. "Matrices and Convolutions of Arithmetic Functions," Vol. 16,
No. 4, pp. 327-334.
GUILLOTTE, G. A. R. Problem Proposed: H-225, Vol. 16, No. 6, p. 569 . Problems Solved: H-272(corrected), H-273, Vol. 16, No. 6, pp. 567-568.
GUPTA, HANSRAJ. "The Rank-Vector of a Partition," Vol. 16, No. 6, pp. 548552; "The Andrews Formula for Fibonacci Numbers," Vol. 16, No. 6, pp. 552555.
GUY, RICHARD. "Binary Sequences without Isolated Ones," Vol. 16, No. 1, pp.
84-86 (co-author, Richard Austin).
HANSEN, RODNEY T. "General Identities for Linear Fibonacci and Lucas Summations," Vol. 16, No. 2, pp. 121-128.
HARBORTH, HEIKO. "&-Adic Numbers in Pascal's Triangle Modulo &," Vol. 16, No.
6, pp. 497-500.
HEED, JOSEPH L. "Entry Points of the Fibonacci Sequence and the Euler 0
Function," Vol. 16, No. 1, pp. 47-50 (co-author, Lucille Kelly).
HENDY, M. D. "Stolarsky's Distribution of Positive Integers," Vol. 16, No. 1,
pp. 70-80.
HENSLEY, DOUGLAS. "Fibonacci Tiling and Hyperbolas," Vol. 16, No. 1, pp. 3740.
HERGET, W. "Minimum Periods Modulo n for Bernoulli Numbers," Vol. 16, No. 6,
pp. 544-548.
HICKERS0N, DEAN R. "Identities Relating the Number of Partitions into an Even
and Odd Number of Parts," Vol. 16, No. 1, pp. 5-6; "An Identity Relating
Compositions and Partitions," Vol. 16, No. 1, pp. 23-26.
HIGGINS, FRANK. Problem Proposed: B-357, Vol. 16, No. 2, p. 186.
Solved: B-357, Vol. 16, No. 2, p. 186.
Problem
HILLMAN, A. P. "A Property of Wythoff Pairs," Vol. 16, No. 5, p. 472 (coauthor, V. E. Hoggatt, Jr.). Editor: "Elementary Problems and Solutions,"
Vol. 16, No. 1, pp. 88-91, 96; Vol. 16, No. 2, pp. 184-187; Vol. 16, No. 5,
pp. 473-476; Vol. 16, No. 6, pp. 562-565.
HOGGATT, V. E., JR. "A Primer for the Fibonacci Numbers, Part XVI, The Central Column Sequence," Vol. 16, No. 1, pp. 41-46 (co-author, J. L. Brown,
Jr.); "A Combinatorial Problem Involving Recursive Sequences and Tridiagonal Matrices," Vol. 16, No. 2, pp. 113-118 (co-author, G. E. Bergum); "A
Primer for the Fibonacci Numbers XVII: Generalized Fibonacci Numbers
574
VOLUME INDEX
[Dec.
Satisfying un+1un_1
- u\
= 1," Vol. 16, No. 2, pp. 130-137 (co-author,
Marjorie Bicknell-Johnson); "Generalized Eulerian Numbers and Polynomials,"
Vol. 16, No. 2, pp. 138-146, 151 (co-author Leonard Carlitz); "A Family
of Tridiagonal Matrices," Vol. 16, No. 3, pp. 285-288 (co-author, G. E. Bergum); "Properties of Generating Functions of a Convolution Array," Vol. 16,
No. 4, pp. 289-295 (co-author Marjorie Bicknell-Johnson); "Convolution Arrays for Jaeobsthal and Fibonacci Polynomials," Vol. 16, No. 5, pp. 385402 (co-author, Marjorie Bicknell-Johnson); "A Property of Wythoff Pairs,"
Vol. 16, No. 5, p. 472 (co-author, A. P. Hillman); "Divisibility Properties
of Polynomials in Pascal!s Triangle," Vol. 16, No. 6, pp. 501-513 (co-author,
Marjorie Bicknell-Johnson). Problems Proposed: B-375, Vol. 16, No. 1, p.
88, B-346, B-347, Vol. 16, No. 1, p. 89; H-278, H-265, Vol. 16, No. 1, pp.
91, 94; B-381, B-352, B-353, Vol. 16, No. 2, pp. 184-185, H-281, Vol. 16,
No. 2, p. 188, H-267(corrected), Vol. 16, No. 2, p. 190; H-285, Vol. 16,
No. 5, p. 477; B-390, B-393, Vol. 16, No. 6, p. 562. Problems Solved: B346, Vol. 16, No. 1, p. 89; B-352, B-353, Vol. 16, No. 2, p. 185; H-267,
Vol. 16, No. 2, p. 190.
H0RADAM, A. F. "Diagonal Functions," Vol. 16, No. 1, pp. 33-36; "Wythoff
Pairs," Vol. 16, No. 2, pp. 147-151.
H0RADAM, E. M. "Solved, Semi-Solved, and Unsolved Problems in Generalized
Integers: A Survey," Vol. 16, No. 4, pp. 370-381.
HOWARD, F.
HUNG, DINH THE'. Problems Solved: B-348, B-351, Vol. 16, No. 1, pp. 90, 91;
B-352, B-353, B-354, Vol. 16, No. 2, pp. 185-186; H-273, Vol. 16, No. 6,
p. 568.
HUNTER, J. A. H.
407-411.
IV'IE, JOHN.
Problems Solved:
JONES, BURTON W. "A Second Variation on a Problem of Diophantus and Davenport," Vol. 16, No. 2, pp. 155-165.
KELLY, LUCILLE. "Entry Points of the Fibonacci Sequence and the Euler 0 Function," Vol. 16, No. 1, pp. 47-50 (co-author, Joseph J. Heed).
KIMBERLING, CLARK. "Strong Divisibility Sequences with Nonzero Initial Term,"
Vol. 16, No. 6, pp. 541-544.
KLAMIN, M.
Problem Solved:
Problem
1978]
VOLUME INDEX
575
LARSON, PAUL. "The Golden Section in the Earliest Notated Western Music,"
Vol. 16, No. 6, pp. 513-515.
LORD, GRAHAM. Problems Solved: B-346, B-347, B-349, B-350, B-351, Vol. 16,
No. 1, pp. 89-91; B-352, B-353, B-354, B-355, B-356, Vol. 16, No. 2, pp.
185-186; B-358, B-359, B-361, B-362, B-363, Vol. 16, No. 5, pp. 474-476;
B-365, B-366, B-367, B-368, B-369, Vol. 16, No. 6, pp. 563-565.
MANA, P. L. Problems Proposed: B-373, Vol. 16, No. 1, p. 88; B-354, Vol. 16,
No. 2, p. 185; B-358, Vol. 16, No. 5, p. 474; B-392, B-365, Vol. 16, No. 6,
pp. 562, 563.
METZGER, J. M. "Fibonacci Sine Sequences," Vol. 16, No. 2, pp. 119-120 (coauthor, M. B. Gregory) .
MILSOM, JOHN W.
pp. 185-186.
MOHANTY, S. P.
384.
MULLEN, GARY L. Problem Proposed: B-376, Vol. 16, No. 2, p. 184 (co-proposer,
Frank Kocher).
MULLIN, A.
Problem Proposed:
0'CALLAHAN, T. Problem Proposed: B-360, Vol. 16, No. 5, pp. 474-475. Problem Solved: B-360, Vol. 16, No. 5, pp. 474-475.
PATERSON, MUCHAEL S. "Identities from Partition Involutions," Vol. 16, No. 3,
pp. 198-212 (co-author, Donald E. Knuth).
PECK, C. B. A. Problems Solved: B-346, B-348, Vol. 16, No. 1, pp. 89-90; B352, B-353, B-354, B-355, Vol. 16, No. 2, pp. 185-186; B-364, Vol. 16, No.
6, p. 563.
PIGNO, LOUIS. "Interpolation of Fourier Transforms on Sums of Fibonacci Numbers," Vol. 16, No. 3, pp. 193-194 (co-author, Robert E. Dressier); "Topological, Measure Theoretic and Analytic Properties of the Fibonacci Numbers," Vol. 16, No. 3, pp. 195-197 (co-author, Robert E. Dressier).
PONNUDURAI, T. Problems Solved:
PR I EL IPP, A. G. Problems Solved: B-346, B-347, B-351, Vol. 16, No. 1, pp.
89-91; B-352, B-353, B-354, B-355, B-356, Vol. 16, No. 2, pp. 185-186;
B-364, B-365, B-366, B-367, B-369, Vol. 16, No. 6, pp. 563-565. H-273,
Vol. 16, No. 6, p. 568.
ROBBINS, NEVILLE.
pp. 515-517.
576
VOLUME INDEX
Dec. 1978
SINGH, SAHIB. Problems Solved: B-358, B-360, B-361, B-362, B-363, Vol. 16,
No. 5, pp. 474-476 (co-solver, Roger Engle); B-359, Vol. 16, No. 5, p. 474
(co-solvers, Benjamin Freed & Roger Engle).
SINNAMON, GORDON.
Problems Solved:
Problem Solved:
Problem
STEINER, RAY. "On Nth Powers in the Lucas and Fibonacci Series," Vol. 16,
No. 5, pp. 451-469; "On fcth Numerical Centers," Vol. 16, No. 5, pp. 470471.
STERN, FREDRICK.
Problem Proposed:
SUN, HUGO S. "Embedding a Group in the pth Powers," Vol. 16, No. 1, p. 4.
SUTTENFIELD, JAMES M., JR.
Problem Solved:
"The Fibonacci Sequence Modulo /I/," Vol. 16, No. 5, pp. 403-
Problem Proposed:
Problem Proposed:
SUSTAINING MEMBERS
*H. L. Alder
*J. Arkin
Do A* Baker
Murray Berg
Gerald Bergum
Je Berkeley
George Berzsenyi
C. A. Bridger
John L. Brown, Jr.
Paul Bruckman
Paul F. Byrd
Co R0 Burton
Lo Carlitz
Go Do Chakerian
P. J, Cocuzza
Mo Jo DeLeon
Harvey Diehl
Jo Lo Ercolano
^Charter Members
Do R* Farmer
*James Maxwell
Harvey Fox
R. Ko McConnell, Jr.
Eo To Frankel
*Sister M. DeSales McNabb
Ro Mo Giuli
Lo P. Meissner
*Ho Wo Gould
Leslie Miller
Nicholas Grant
F. J. Ossiander
William Greig
Fo Go Rothwell
Vo Co Harris
Co Eo Serkland
Ao Po Hillman
A. G o Shannon
*A F9 Horadam
Jo A. Sehumaker
*Verner E Hoggatt, Jr 8 D. Singmaster
Virginia Kelemen
Co Co Styles
Ro Po Kelisky
L. Taylor
Co Ho Kimberling
Ho Lo Umansky
Jo Krabacker
*Lo Ao Walker
George Ledin, Jr e
Marcellus Waddill
*C. To Long
Paul Willis
Jo Ro Ledbetter
Co Fo Winans
Do Po Mamuscia
Eo L o Yang
Typed by
~*J0 ANN VTSEh
Campbell, California