シルベスター数列

数論において、シルベスター数列 (英語: Sylvester's sequence) とは、各項がそれまでの項の総積に1を足したものであるような数列である。最初のいくつかの項は次のようになる： 2, 3, 7, 43, 1807, 3263443, 10650056950807, ... (オンライン整数列大辞典の数列 A000058) 「シルベスター数列」という名称は、1880年にこの数列を最初に調査したジェームス・ジョセフ・シルベスターに由来している。数列の各項の値はシルベスター数 (英: Sylvester number) と呼ばれることもある。 シルベスター数列の値は二重指数関数的に増加し、その逆数の和は他の単位分数の級数よりも速く1に収束する。シルベスター数列を定義する漸化式は、各項の数の素因数分解をより容易にさせるが、数列の増加速度が急速であるために、完全な素因数分解はいくつかの項に対してしか知られていない。 シルベスター数は1のエジプト式分数表現、佐々木・アインシュタイン多様体、オンラインアルゴリズムの実装などに応用されている。 シルベスター数列の第 n 項は次の式で定義される： s n = 1 + ∏ i = 0 n − 1 s i . {\displaystyle s_{n}=1+\prod _{i=0}^{n-1}s_{i}.} 0個の項の積 (空積) は 1 であるため、s0 = 2 である。 あるいは、次のような漸化式で定義してもよい： s i = s i − 1 ( s i − 1 − 1 ) + 1 , ( s 0 = 2 ) {\textstyle \displaystyle s_{i}=s_{i-1}(s_{i-1}-1)+1,\quad (s_{0}=2)} シルベスター数列は n の関数として、二重指数関数的に増加する。具体的には s n = ⌊ E 2 n + 1 + 1 2 ⌋ {\displaystyle s_{n}=\left\lfloor E^{2^{n+1}}+{\frac {1}{2}}\right\rfloor } という形式で書くことができる。ここで E はおよそ 1.2640847353...である (オンライン整数列大辞典の数列 A076393)。 シルベスター数列が二重指数関数的増加を示すことは、フェルマー数列 Fn と比較すると驚くべきことではない。フェルマー数は通常、二重指数関数による表式 F n = 2 2 n + 1 {\textstyle F_{n}=2^{2^{n}}+1} によって定義されるが、これはシルベスター数列のような総積を使った漸化式によっても定義が可能である： F n = 2 + ∏ i = 0 n − 1 F i . {\displaystyle F_{n}=2+\prod _{i=0}^{n-1}F_{i}.} シルベスター数列の逆数和による無限級数 ∑ i = 0 ∞ 1 s i = 1 2 + 1 3 + 1 7 + 1 43 + 1 1807 + ⋯ . {\displaystyle \sum _{i=0}^{\infty }{\frac {1}{s_{i}}}={\frac {1}{2}}+{\frac {1}{3}}+{\frac {1}{7}}+{\frac {1}{43}}+{\frac {1}{1807}}+\cdots .} について考える。この級数の部分和は次のような単純な形で書ける： ∑ i = 0 j − 1 1 s i = 1 − 1 s j − 1 = s j − 2 s j − 1 . {\displaystyle \sum _{i=0}^{j-1}{\frac {1}{s_{i}}}=1-{\frac {1}{s_{j}-1}}={\frac {s_{j}-2}{s_{j}-1}}.} これは、帰納法によって、またはより直接的に、任意の i に対する次の等式 1 s i − 1 − 1 s i + 1 − 1 = 1 s i {\displaystyle {\frac {1}{s_{i}-1}}-{\frac {1}{s_{i+1}-1}}={\frac {1}{s_{i}}}} から、級数の畳み込みを行うことによって示される： ∑ i = 0 j − 1 1 s i = ∑ i = 0 j − 1 ( 1 s i − 1 − 1 s i + 1 − 1 ) = 1 s 0 − 1 − 1 s j − 1 = 1 − 1 s j − 1 . {\displaystyle \sum _{i=0}^{j-1}{\frac {1}{s_{i}}}=\sum _{i=0}^{j-1}\left({\frac {1}{s_{i}-1}}-{\frac {1}{s_{i+1}-1}}\right)={\frac {1}{s_{0}-1}}-{\frac {1}{s_{j}-1}}=1-{\frac {1}{s_{j}-1}}.} 部分和が j → ∞ で1に収束することから、級数全体が1に収束することがわかる。従って私たちは、これによって無限の長さを持つ1のエジプト式分数表示を得たことになる。 1 = 1 2 + 1 3 + 1 7 + 1 43 + 1 1807 + ⋯ {\displaystyle 1={\frac {1}{2}}+{\frac {1}{3}}+{\frac {1}{7}}+{\frac {1}{43}}+{\frac {1}{1807}}+\cdots } このとき、級数を適当な長さで切って、最後の分母を1小さいものに置き換えることで、任意の長さを持つ1のエジプト式分数表示を得ることができる。 1 = 1 2 + 1 3 + 1 6 , 1 = 1 2 + 1 3 + 1 7 + 1 42 , 1 = 1 2 + 1 3 + 1 7 + 1 43 + 1 1806 , … . {\displaystyle 1={\tfrac {1}{2}}+{\tfrac {1}{3}}+{\tfrac {1}{6}},\quad 1={\tfrac {1}{2}}+{\tfrac {1}{3}}+{\tfrac {1}{7}}+{\tfrac {1}{42}},\quad 1={\tfrac {1}{2}}+{\tfrac {1}{3}}+{\tfrac {1}{7}}+{\tfrac {1}{43}}+{\tfrac {1}{1806}},\quad \dots .} また、無限級数の最初の k 項の和は、k 項からなるエジプト式分数表示のうち1未満で最も大きいものを提供する。例えば、最初の4項の和は 1805/1806 となるため、開区間 (1805/1806, 1) に含まれる数のエジプト式分数表示は少なくとも5つの項が必要になる。 シルベスター数列は、各ステップごとにその部分和が1以下となるような最小の分母を選ぶ貪欲法の結果として見ることもできる。あるいは、(初項である1/2を除外して、) 1/2の奇数貪欲法によるエジプト式分数展開 (Odd greedy expansion) だと考えることもできる。 シルベスター数の逆数和が有理数である1に収束することから、数列が二重指数関数的に増加することは、その逆数和による級数が無理数となること、さらに数列が irrationality sequence[定訳なし] となることの十分条件ではないことがわかる。 一方で Badea (1993) の結果から、シルベスター数列 (の漸化式) は、二重指数関数的な増加を持ち逆数和が有理数に収束するような数列に対する、ある種の唯一性を持っている。つまり、数列 {an} が不等式 a n ≥ a n − 1 2 − a n − 1 + 1 , {\displaystyle a_{n}\geq a_{n-1}^{2}-a_{n-1}+1,} を満たし、逆数和が有理数に収束するならば、十分に大きなすべての n に対して、その漸化式はシルベスター数列と同じ a n = a n − 1 2 − a n − 1 + 1 {\displaystyle a_{n}=a_{n-1}^{2}-a_{n-1}+1} となる。 エルデシュとグラハムは、数列 {an} が lim n → ∞ a n a n − 1 2 = 1. {\displaystyle \lim _{n\rightarrow \infty }{\frac {a_{n}}{a_{n-1}^{2}}}=1.} を満たすとき、逆数和が有理数となるならば、十分大きな全ての n に対して a n = a n − 1 2 − a n − 1 + 1 {\textstyle a_{n}=a_{n-1}^{2}-a_{n-1}+1} が成り立つと予想した。Badea (1995) ではこの予想についていくつかの結果が纏められている。 i < j のとき、定義から sj ≡ 1 (mod si) が成り立つ。従って、任意の2つのシルベスター数は互いに素である。任意の素数に対して、それで割り切れるシルベスター数が高々1つとなるため、これを用いて素数が無限に存在することを証明できる。より強い事実として、シルベスター数の素因数に6を法として5と合同である数 (p = 6k + 5 と表せるような素数) は存在しないこと、12を法として7と合同である素数が無限に存在することがシルベスター数列を用いて証明できる。 シルベスター数の素因数分解については、多くのことが未解決のままである。例えば、全ての数が平方因子をもたない整数であるかどうか不明である (既知の数は全て平方因子を持たないことがわかっている)。 Vardi (1991)が説明しているように、与えられた素数 p で割り切れるシルベスター数がどれかを決定することは簡単である：p を法として0に合同となる数か、周期的な列のいずれかが見つかるまで、法 p の下でのシルベスター数列を計算すればよい。この方法を使って、Vardiは最初の300万個の素数のうち1166個がシルベスター数の素因数であること、それらの数の平方がシルベスター数の因数にならないことを突き止めた。シルベスター数の因数として現れる素数の集合は、全ての素数の集合に対して密度0 (Natural density の意味で) である。実際、x より小さいそのような数は O ( x / ( log ⁡ x ⋅ log ⁡ log ⁡ log ⁡ x ) ) {\textstyle O(x/(\log x\cdot \log \log \log x))} のオーダーとなる。 次の表は、シルベスター数の既知の因数分解を示している。ただし最初の4つは全て素数であるため除いている。また、一部の数は大きすぎるため桁数のみを表しており、特に明記されていなければそれらの数は素数とは限らない (unfactoredである)。 Boyer, Galicki & Kollár (2005) ではシルベスター数列の性質を使って、奇数次元の球面またはエキゾチック球面（英語版）に対する非同値な佐々木・アインシュタイン多様体の族を与えている。彼らは論文中で、2m－3 次元の球面またはエキゾチック球面上で少なくとも 1 3 ( s m − 1 − 1 ) {\textstyle {\tfrac {1}{3}}(s_{m-1}-1)} 個の非同値な佐々木・アインシュタイン多様体の族を与える方法を示し、従ってそれらの数は m に対して二重指数関数的な増加を見せる。 Galambos & Woeginger (1995)が説明しているように、 Brown (1979)とLiang (1980)はオンラインのビンパッキングアルゴリズムについて、アルゴリズムの評価の下界を与えるためにシルベスター数列の値を利用した。Seiden & Woeginger (2005) でも同様に、2次元カッティングストック問題に対して、3-stageギロチンカット解の下界を与えるためにシルベスター数列が用いられている。 Známの問題（英語版）は、集合の各要素がそれぞれ他の数の総積に1を足した数の真の約数であるような集合についての問題である。「真の」約数であるという条件を除くと、シルベスター数列の値はこの問題を解決する。本来の条件の下では、シルベスター数列を定めるものと同様の再帰的な手続きによって、問題の解を得ることができる[要出典]。Známの問題の解は表面特異点の分類や、unary n-cyclic 正規言語を通して非決定性有限オートマトンの理論への応用が存在する。 Curtiss (1922) は単位分数の k 項和で1に最も近い近似が、完全数の約数の数に下界を与えることが記されている。Miller (1919) では同じ性質が k 個の部分群を持つ群の位数の上界を与えるために使われている。 カエン定数（英語: Cahen's constant） Primary pseudoperfect number - オンライン整数列大辞典の数列 A054377 レオナルド数（英語: Leonardo number）- オンライン整数列大辞典の数列 A001595 ^ 数列の増加速度と級数の無理性については、例えば数列 {an} が十分に速く増加するとき、 ∑ k 1 a k a k + 1 {\textstyle \sum _{k}{\tfrac {1}{a_{k}a_{k+1}}}} が無理数となることが知られている (Erdős & Graham 1980, p. 64)。 ^ Andersenはこの区間で1167の素因数を見つけたため、おそらくこれは誤記である。 ^ p < 5 × 107 かつ n ≦ 200 を満たす範囲において、全てのシルベスター数 sn の素因数 p はVardiによってリストされている。Ken Takusagawa は s9 までの素因数分解と s10 の素因数分解をブログに記している。それ以外の因数分解については、Jens Kruse Andersen によるリストを出典としている。 ^ 佐々木多様体（英語版）でもあるアインシュタイン多様体 ^ 論文中でSeidenとWoegingerは、シルベスター数列を Salzer (1947) の仕事にちなんで「Salzer's sequence」という名前で言及している。 ^ Finch (2003, p. 444) ^ Golomb (1963), Aho & Sloane (1973, §2.5) ^ Curtiss (1922)、Miller (1919) あるいは Soundararajan (2005) を参照。 ^ Erdős & Graham (1980) ^ Guy & Nowakowski (1975). ^ a b Andersen (2007–20) ^ Jones (2006). ^ Odoni (1985). ^ Takusagawa (2006a) ^ Takusagawa (2006b) ^ Boyer, Galicki & Kollár (2005, §7)、特に Prop. 44–Prop. 48 ^ Brenton & Hill (1988) ^ Domaratzki et al. (2005). Aho, A. V.; Sloane, N. J. A. (11 1973). “Some Doubly Exponential Seqiences”. The Fibonacci Quarterly 11 (4): 429–437. https://www.fq.math.ca/11-4.html.  Badea, Catalin (1993). “A theorem on irrationality of infinite series and applications”. Acta Arithmetica 63 (4): 313–323. doi:10.4064/aa-63-4-313-323. MR1218459.  Badea, Catalin (1995). On some criteria of irrationality for series of positive rationals : a survey (PDF). Actes des Rencontres Arithmétiques de Caen. 2008年9月11日時点のオリジナル (pdf)よりアーカイブ。2006年9月27日閲覧。 Boyer, Charles P.; Galicki, Krzysztof; Kollár, János (2005). “Einstein metrics on spheres”. Annals of Mathematics 162 (1): 557–580. arXiv:math.DG/0309408. doi:10.4007/annals.2005.162.557. MR2178969.  Brenton, Lawrence; Hill, Richard (1988). “On the Diophantine equation 1=Σ1/ni + 1/Πni and a class of homologically trivial complex surface singularities”. Pacific Journal of Mathematics 133 (1): 41–67. doi:10.2140/pjm.1988.133.41. MR0936356. http://projecteuclid.org/euclid.pjm/1102689567.  Brown, D. J. (1979). A lower bound for on-line one-dimensional bin packing algorithms (Technical report). Coordinated Science Laboratory Report. Vol. R-864. Coordinated Science Lab., Univ. of Illinois, Urbana-Champaign. hdl:2142/74218。 Curtiss, D. R. (1922). “On Kellogg's diophantine problem”. American Mathematical Monthly 29 (10): 380–387. doi:10.2307/2299023. JSTOR 2299023.  Domaratzki, Michael; Ellul, Keith; Shallit, Jeffrey; Wang, Ming-Wei (2005). “Non-uniqueness and radius of cyclic unary NFAs”. International Journal of Foundations of Computer Science 16 (5): 883–896. doi:10.1142/S0129054105003352. MR2174328. http://www.cs.umanitoba.ca/~mdomarat/pubs/DESW_dcfs.ps.  Erdős, Paul、Graham, Ronald L.『Old and new problems and results in combinatorial number theory』Monographies de L'Enseignement Mathématique, No. 28, Univ. de Genève、1980年。MR0592420。  Finch, Steven R.『Mathematical constants』Cambridge University Press、New York〈Encyclopedia of Mathematics and its Applications〉、2003年。ISBN 978-1-107-26691-9。ISSN 0953-4806。OCLC 847526740。https://www.worldcat.org/oclc/847526740。  Galambos, Gábor; Woeginger, Gerhard J. (1995). “On-line bin packing — A restricted survey”. Mathematical Methods of Operations Research 42 (1): 25. doi:10.1007/BF01415672. MR1346486.  Golomb, Solomon W. (1963). “On certain nonlinear recurring sequences”. American Mathematical Monthly 70 (4): 403–405. doi:10.2307/2311857. JSTOR 2311857. MR0148605.  Graham, R.、Knuth, D. E.、Patashnik, O.『Concrete Mathematics』（2nd）Addison-Wesley、1989年、Exercise 4.37頁。ISBN 0-201-55802-5。  Guy, Richard K.「E24 Irrationality sequences」『Unsolved Problems in Number Theory』（3rd）Springer-Verlag、2004年、346頁。ISBN 0-387-20860-7。Zbl 1058.11001。https://books.google.com/books?id=1AP2CEGxTkgC&pg=PA346。  Guy, Richard; Nowakowski, Richard (1975). “Discovering primes with Euclid”. Delta (Waukesha) 5 (2): 49–63. MR0384675.  Jones, Rafe (2006). “The density of prime divisors in the arithmetic dynamics of quadratic polynomials”. Journal of the London Mathematical Society 78 (2): 523–544. arXiv:math.NT/0612415. Bibcode: 2006math.....12415J. doi:10.1112/jlms/jdn034.  Liang, Frank M. (1980). “A lower bound for on-line bin packing”. Information Processing Letters 10 (2): 76–79. doi:10.1016/S0020-0190(80)90077-0. MR0564503.  Miller, G. A. (1919). “Groups possessing a small number of sets of conjugate operators”. Transactions of the American Mathematical Society 20 (3): 260–270. doi:10.2307/1988867. JSTOR 1988867.  Odoni, R. W. K. (1985). “On the prime divisors of the sequence wn+1 =1+w1⋯wn”. Journal of the London Mathematical Society. Series II 32: 1–11. doi:10.1112/jlms/s2-32.1.1. Zbl 0574.10020.  Rosenman, Martin; Underwood, F. (1933). “Problem 3536”. American Mathematical Monthly 40 (3): 180–181. doi:10.2307/2301036. JSTOR 2301036.  Salzer, H. E. (1947). “The approximation of numbers as sums of reciprocals”. American Mathematical Monthly 54 (3): 135–142. doi:10.2307/2305906. JSTOR 2305906. MR0020339.  Seiden, Steven S.; Woeginger, Gerhard J. (2005). “The two-dimensional cutting stock problem revisited”. Mathematical Programming 102 (3): 519–530. doi:10.1007/s10107-004-0548-1. MR2136225.  Soundararajan, K. (11 February 2005). "Approximating 1 from below using n Egyptian fractions". arXiv:math/0502247。 Sylvester, J. J. (1880). “On a point in the theory of vulgar fractions”. American Journal of Mathematics 3 (4): 332–335. doi:10.2307/2369261. JSTOR 2369261.  Vardi, Ilan『Computational Recreations in Mathematica』Addison-Wesley、1991年、82–89頁。ISBN 0-201-52989-0。  Andersen, Jens Kruse (2007–20). "Factorization of Sylvester's sequence". 2022年4月7日時点のオリジナルよりアーカイブ。2022年9月4日閲覧。 Takusagawa, Ken (19 January 2006a). "factoring Sylvester's sequence". Ken's blog. 2022年8月18日時点のオリジナルよりアーカイブ。2022年9月4日閲覧。 Takusagawa, Ken (2 April 2006b). "Sylvester 10th factored". Ken's blog. 2022年8月18日時点のオリジナルよりアーカイブ。2022年9月4日閲覧。