シルベスターの数列の具体例・性質・漸化式について

はるか
はるか
シルベスターの数列って知ってる?
ふゅか
ふゅか
漸化式がいろんな形になるから面白いよね!

1. シルベスターの数列とは

シルベスターの数列 {an}\lbrace{a_n}\rbrace は、次のような漸化式になります。a1=2 a_1 = 2 an+1=an2an+1 a_{n+1} = a_n^2 – a_n + 1

この漸化式によって、シルベスターの数列(Sylvester’s sequence) ana_n が以下のように計算されます。

1. a1=2 a_1 = 2
2. a2=222+1=42+1=3 a_2 = 2^2 – 2 + 1 = 4 – 2 + 1 = 3
3. a3=323+1=93+1=7 a_3 = 3^2 – 3 + 1 = 9 – 3 + 1 = 7
4. a4=727+1=497+1=43 a_4 = 7^2 – 7 + 1 = 49 – 7 + 1 = 43
5. a5=43243+1=184943+1=1807 a_5 = 43^2 – 43 + 1 = 1849 – 43 + 1 = 1807

このように、次の項を計算することができます。

2. シルベスターの数列の性質

2.1. 漸化式が積の形になる

シルベスターの数列{an}\lbrace{a_n}\rbrace の漸化式はan+1=k=1nak+1 a_{n+1} = \prod_{k=1}^{n} a_k + 1 となる。
はるか
はるか
実は、この数列の漸化式、積の形にもできるんだ。
ふゅか
ふゅか
うん、確かに!an+1=k=1nak+1a_{n+1} = \prod_{k=1}^{n} a_k + 1 っていう形だよね。

数学的帰納法を利用して証明します。

n=1n=1のとき、次の項 a2 a_2 は次のようになります。

a2=a12a1+1=222+1=42+1=3 a_2 = a_1^2 – a_1 + 1 = 2^2 – 2 + 1 = 4 – 2 + 1 = 3

また、漸化式より、a2a_2は次のようになります。

a2=222+1=42+1=3a_2 = 2^2 – 2 + 1 = 4 – 2 + 1 = 3

n=k n = k のときに ak+1=i=1kai+1 a_{k+1} = \prod_{i=1}^{k} a_i+ 1 であると仮定します。

次に、n=k+1 n = k+1 のときの式を考えます。ak+2 a_{k+2} は以下のように表されます。

ak+2=ak+12ak+1+1 a_{k+2} = a_{k+1}^2 – a_{k+1} + 1

ここで、仮定により ak+1=i=1kai+1 a_{k+1} = \prod_{i=1}^{k} a_i + 1 ですので、ak+2 a_{k+2} を次のようになります。

ak+2=(i=1kai+1)2(i=1kai+1)+1 a_{k+2} = \left(\prod_{i=1}^{k} a_i + 1\right)^2 – \left(\prod_{i=1}^{k} a_i + 1\right) + 1

=(i=1kai)2+2i=1kai+1i=1kai1+1= \left(\prod_{i=1}^{k} a_i\right)^2 + 2\prod_{i=1}^{k} a_i + 1 – \prod_{i=1}^{k} a_i – 1 + 1

=i=1kai×(i=1kai+1)+1= \prod_{i=1}^{k} a_i \times \left(\prod_{i=1}^{k} a_i + 1\right) + 1

となります。ak+1=i=1kai+1 a_{k+1} = \prod_{i=1}^{k} a_i+ 1 より、

ak+2=i=1kai×ak+1+1 a_{k+2} = \prod_{i=1}^{k} a_i \times a_{k+1} + 1

ak+2=i=1k+1ai+1\therefore a_{k+2} = \prod_{i=1}^{k+1} a_i + 1

となります。

よって、n=k+1 n = k+1 の場合にも、与えられた漸化式 が成り立ちます。したがって、数学的帰納法により、与えられた数列が次の形で表されることが証明されました。

an+1=k=1nak+1 a_{n+1} = \prod_{k=1}^{n} a_k + 1

2.2. 逆数の和

シルベスター数列 an a_n について、以下の等式を示します。

1a1+1a2++1an=11k=1nak \frac{1}{a_1} + \frac{1}{a_2} + \dots + \frac{1}{a_n} = 1 – \frac{1}{\displaystyle\prod_{k=1}^{n} a_k}

ここで、k=1nak\prod_{k=1}^{n} a_ka1×a2××an a_1 \times a_2 \times \dots \times a_n を表す積の記号です。

数学的帰納法を利用して証明します。

まず、初項 n=1 n = 1 の場合を考えます。この場合、a1=2a_1=2より、式は次のようになります。

1a1=11a1 \frac{1}{a_1} = 1 – \frac{1}{a_1}

この等式は明らかに成り立ちます。

次に、n=k n=k に対して次の等式が成り立つと仮定します。

1a1+1a2++1ak=11i=1kai \frac{1}{a_1} + \frac{1}{a_2} + \dots + \frac{1}{a_k} = 1 – \frac{1}{\prod_{i=1}^{k} a_i}

n=k+1 n=k+1 のとき、式は次のようになります。

1a1+1a2++1ak+1ak+1 \frac{1}{a_1} + \frac{1}{a_2} + \dots + \frac{1}{a_k} + \frac{1}{a_{k+1}}

仮定より、次のように変形します。

=(11i=1kai)+1ak+1 = \left(1 – \frac{1}{\prod_{i=1}^{k} a_i}\right) + \frac{1}{a_{k+1}}

次に、先ほど示したシルベスター数列の性質 ak+1=i=1kai+1 a_{k+1} = \prod_{i=1}^{k} a_i + 1 を利用します。これにより、

=11i=1kai+1i=1kai+1 = 1 – \frac{1}{\displaystyle\prod_{i=1}^{k} a_i} + \frac{1}{\displaystyle\prod_{i=1}^{k} a_i + 1}

=1+i=1kai1+i=1kaii=1k+1ai = 1 + \frac{-\displaystyle\prod_{i=1}^{k} a_i – 1 + \displaystyle\prod_{i=1}^{k} a_i }{\displaystyle\prod_{i=1}^{k+1} a_i}

=11i=1k+1ai = 1 – \frac{1}{\displaystyle\prod_{i=1}^{k+1} a_i}

これにより、n=k+1 n=k+1 でも等式が成り立つことが示されました。

したがって、数学的帰納法により次の等式が成り立つことが証明されました。

1a1+1a2++1an=11k=1nak \frac{1}{a_1} + \frac{1}{a_2} + \dots + \frac{1}{a_n} = 1 – \frac{1}{\displaystyle\prod_{k=1}^{n} a_k}

はるか
はるか
逆数の和が 11k=1nak1 – \frac{1}{\prod_{k=1}^{n} a_k} にもなる。

2.3. エジプト分数

 

自然数の組 (a1,a2,,an)(a_1, a_2, \dots, a_n) が、k=1n1ak<1\displaystyle\sum_{k=1}^{n} \frac{1}{a_k} < 1 の条件を満たしつつ、その和を最大にする。このとき、自然数の組はシルベスターの数列になる。
はるか
はるか
さらに、この数列はエジプト分数の組み合わせにも使われる。
ふゅか
ふゅか
そうそう!自然数の組み合わせで、和が1未満で最大の値を取るとき、それがシルベスターの数列になるんだよね!例えば、a1=2a_1 = 2, a2=3a_2 = 3のときとか。

以下に、具体的な例を挙げて解説します。

n=2 n = 2 の場合
1a1+1a2<1 \frac{1}{a_1} + \frac{1}{a_2} < 1
この条件を満たす自然数 a1<a2a_1 < a_2 を求めると、最大値は a1=2a_1 = 2, a2=3a_2 = 3 のときです。このときの和は
12+13=56 \frac{1}{2} + \frac{1}{3} = \frac{5}{6}
となります。

n=3 n = 3 の場合
1a1+1a2+1a3<1 \frac{1}{a_1} + \frac{1}{a_2} + \frac{1}{a_3} < 1
この条件を満たす自然数 a1<a2<a3a_1 < a_2 < a_3 を求めると、最大値は a1=2a_1 = 2, a2=3a_2 = 3, a3=7a_3 = 7 のときです。このときの和は
12+13+17=4142 \frac{1}{2} + \frac{1}{3} + \frac{1}{7} = \frac{41}{42}
となります。