数列の和
・1 から n までの自然数の和 : 1+2+・・・+n =n(n+1)/2
・1 から n までの自然数の平方の和 : 12+22+・・・+n2 =n(n+1)(2n+1)/6
・1 から n までの自然数の立方の和 : 13+23+・・・+n3 =n2(n+1)2/4
上記の和の公式は、高校で学ぶ、よく知られた公式である。しかし、その証明は、恒等式
を使うなど、あまりに代数的で味気ない。このページでは、もっと直感的に分かりやすい、味
のある証明についてまとめたいと思う。まずは、教科書の記述を復習してみよう。
1 から n までの自然数の和
教科書では、恒等式 (k+1)2−k2 =2k+1 を利用して、k に1,2,3,・・・,n を
代入して求めるのが普通である。
n 個の式を辺々加えて、 (n+1)2−1 =2( 1+2+3+・・・+n )+n
よって、 2( 1+2+3+・・・+n )=(n+1)2−(n+1)=n(n+1)
この式より、上記の公式を得る。
(追記) 令和元年6月17日付け
恒等式 (k+1)2−k2 =2k+1 を辺々加えずに、次のような導き方も知られている。
Σk=1n {(k+1)2−k2} =Σk=1n (k+1)2−Σk=1n k2=(n+1)2−1=n2+2n
よって、 2Σk=1n k+n=n2+2n より、 Σk=1n k=n(n+1)/2
1 から n までの自然数の平方の和
教科書では、恒等式 (k+1)3−k3 =3k2+3k+1 を利用して、k に1,2,3,・・・,n
を代入して求めるのが普通である。
n 個の式を辺々加えて、
(n+1)3−1 =3(12+22+32+・・・+n2)+3( 1+2+3+・・・+n )+n
よって、
3(12+22+32+・・・+n2)
=(n+1)3−3×n(n+1)/2−(n+1)=n(n+1)(2n+1)/2
この式より、上記の公式を得る。
(追記) 令和元年6月17日付け
恒等式 (k+1)3−k3 =3k2+3k+1 を辺々加えずに、次のような導き方も知られてい
る。
Σk=1n {(k+1)3−k3} =Σk=1n (k+1)3−Σk=1n k3=(n+1)3−1
よって、 3Σk=1n k2+3n(n+1)/2+n=(n+1)3−1 より、
3Σk=1n k2
=(n+1)3−3n(n+1)/2−n−1
={(n+1)/2}{2(n+1)2−3n−2}
={(n+1)/2}{2n2+n}
=n(n+1)(2n+1)/2
よって、 Σk=1n k2=n(n+1)(2n+1)/6
(追記) 令和3年10月13日付け
上記の証明が一番よく知られていると思われるが、2項係数の「総和の公式」を用いた証
明も斬新で興味深い。
(追記) 令和3年9月7日付け
公式:12+22+・・・+n2 =n(n+1)(2n+1)/6 の覚え方について、
12+22+・・・+n2 =n(n+0.5)(n+1)/3
とした方がいいという意見もある。nから0.5刻みに増えていって、3で割るという解釈なの
だが、皆さんはいかがだろか?
例 12+22+・・・+202 =20×20.5×21/3=410×7=2870
(コメント) この発想が活かせる場面はどんなだろう?
1 から n までの自然数の立方の和
教科書では、恒等式 (k+1)4−k4 =4k3+6k2+4k+1 を利用して、k に
1,2,3,・・・,n を代入して求めるのが普通である。
その計算やまとめ方のコツは、平方数の場合と全く同じなので、ここでは、その記述を省
略することにする。
この1 から n までの自然数の立法の和が、1 から n までの自然数の和 N を用いて、
13+23+33+・・・+n3 =N2
と書けるということは、私自身高校時代に不思議に思った事実である。
このことについて、当HPがいつもお世話になっているS(H)さんからの情報(平成21年6月
29日付け)によると、
Ulm 生まれの Johann Faulhaber (1580〜1635)が、その著
Academia Algebra (1631)
において、 15+25+35+・・・+n5 =(4N3−N2)/3
等の公式(117+217+317+・・・+n17 までの奇数冪について)を述べているという。
一般の場合については、1834年に、Jacobi が証明したらしい。
15+25+35+・・・+n5 を求めてみよう。
まず、恒等式 (k+1)5−k5 =5k4+10k3+10k2+5k+1 を利用して、
(n+1)5−15 =5Σk4+10Σk3+10Σk2+5Σk+Σ1
これより、 Σk4=n(n+1)(2n+1)(3n2+3n−1)/30
また、恒等式 (k+1)6−k6 =6k5+15k4+20k3+15k2+6k+1 を利用して、
(n+1)6−16 =6Σk5+15Σk4+20Σk3+15Σk2+6Σk+Σ1
これより、 Σk5=n2(n+1)2(2n2+2n−1)/12
よって、 n(n+1)=2N を代入すれば、
15+25+35+・・・+n5=4N2(4N−1)/12=(4N3−N2)/3
であることが示された。
例 15+25+35+・・・+105 を計算せよ。
(解) N=55 なので、
(4N3−N2)/3=N2(4N−1)/3=3025・219/3=3025・73=220825
より、 15+25+35+・・・+105=220825 である。 (終)
(コメント) 比較的簡単に求められますね!
連続する整数の積の形で、k乗和の公式を求める方法もある。(→ 参考)
このような数列の和の公式に対する別解として世間一般では、2S法とか3S法、4S法
などの方法が紹介されている。(この方法にも限界があることが自ずから分かるであろう。)
2S法・・・1 から n までの自然数の和
![]() |
左図より明らかに、2S=(n+1)がn個=n(n+1) が成り |
3S法・・・1 から n までの自然数の平方の和
上記の計算において、各三角形の上からk段目、左からm列目の数を求めてみると、そ
れぞれ k 、n−k+m 、n−m+1 なので、その和は、2n+1 となる。
したがって、上図より明らかに、3S=(2n+1が1+2+・・・+n=n(n+1)/2 個) が成
り立ち、直ちに公式を得ることができる。
#自然数の平方の和の公式は、すでに、アルキメデス(紀元前287?〜212)が知ってい
たと言われている。
4S法・・・1 から n までの自然数の立方の和
4S法は、2S法、3S法と同様な発想であるが、上記の方法と比べて煩雑な計算が避けら
れず、エレガントさに欠けるので、紹介は、あらすじだけにとどめたい。
基本的には、S’=1+2(1+2)+3(1+2+3)+・・・+n(1+2+・・・+n) を計算する
わけであるが、次のような立体を4個用いる。用い方は、3S法と同様である。
一つの点に注目すると、和はすべて 3n+1 であることに注意する。
1 から n までの自然数の立方の和
立方の和の計算について、実は美しい解法が知られている。私が知っている和の計算の
中でも、その美しさは群を抜いている。それは、群数列による解法である。
![]() |
左図のように、自然数列 1、2、・・・、n を、2倍、3倍、・・・ したものを、下の行に書き足していき、青線で区切り、群数列 を作る。左図にあらわれる自然数の総和は、明らかに、 (1+2+3+・・・+n)2=n2(n+1)2/4 そこで、第k群に入っている自然数の和を求めると、 (1+2+・・・+k)k+(1+2+・・・+k)k−k2=k3 である。 |
従って、第1群から第n群までの和を求めることにより、公式は示される。
(別解) 1 から n までの自然数の立方の和について、次のような群数列を用いる解法も
知られている。
(参考文献:金子昌信 著 べき乗和からベルヌーイ数,ゼータ函数へ
(数学セミナー ’03年10月号) (日本評論社))
![]() |
1 から始まる奇数列の初項から第 n 項までの和は、n2 で与えられる。 第 k 段には、k 個の奇数が含まれる ので、第1段から第 n 段までに含まれる 奇数の総数は、 1+2+3+・・・+n =n(n+1)/2 |
したがって、第1段から第 n 段までに含まれる奇数の総和は、 n2(n+1)2/4
ところで、第 k 段は、初項が k2−k+1、末項が k2+k−1、項数が k の等差数列であ
るから、その和は、簡単な計算から、k3 となる。
したがって、 13+23+33+・・・+n3 =n2(n+1)2/4 (終)
1 から n までの自然数の平方の和(ライプニッツの方法)
平方の和の計算について、階差数列を利用する方法が知られている。この方法は、ライ
プニッツが発見したと言われている。
|
左図において、 bn=an−an-1、cn=bn−bn-1、・・・ である。 このとき、たとえば、 |
a3=a2+b3=(a1+b2)+(b2+c3)=a1+2b2+c3 ・・・・・(*)
=(a0+b1)+2(b1+c2)+(c2+d3)=a0+3b1+3c2+d3 ・・・・・(**)
このことから、a3 の値は、対角線の成分から求められることが分かる。さらに、その各項
の係数は、その項から、a3 へ行く最短経路の数に等しい。
たとえば、b1 から a3 へ行くには、横に2マス、縦に1マス進めばよいから、その方法の
数は、順列・組合せの理論により、3C1 =3(通り)である。
上の計算で、(*)(**)のところを是非確かめて欲しいと思う。
定理(ライプニッツ)
an=nC0a0+nC1b1+nC2c2+nC3d3+・・・
実際に、自然数の平方の和の計算に応用してみよう。
|
左図において、第1行は、 an=12+22+32+・・・+n2 で与えられる数列であり、第2行以下は、その階差数列で ある。ただし、a0=0 とする。 |
このとき、第5行以下は、全て0であることに注目する。 定理により、
an=nC0×0+nC1×1+nC2×3+nC3×2
となり、これを計算することにより、公式を得ることができる。
(参考文献:遠山 啓 著 数学入門(岩波新書)
久保季夫 著 数列と級数(科学新興社))
( → 数列の和の求め方で重心の概念を利用した方法については、「ギブスの三角形」を参照)
(追記)1 から n までの自然数の4乗の和
4乗の和の計算について、立方の和と同様に、群数列を用いる方法がある。
![]() |
左図のように、数列 1、4、9、・・・、n2 を、2倍、3倍、・・・ したものを、下の行に書き足していき、青線で区切り、群数列 を作る。 左図にあらわれる数列の総和は、明らかに、 (1+2+3+・・・+n)(12+22+32+・・・+n2) =n2(n+1)2(2n+1)/12 そこで、第k群に入っている自然数の和を求めると、 |
(1+2+・・・+k)k2+(12+22+・・・+k2)k−k3=(5/6)k4+(1/6)k2
したがって、第1群から第n群までの和を求めることにより、求める和を S とすると、
(5/6)S+n(n+1)(2n+1)/36=n2(n+1)2(2n+1)/12
よって、 S=n(n+1)(2n+1)(3n2+3n−1)/30
(参考文献:中村義作・阿邊恵一 著 代数を図形で解く(講談社ブルーバックス))
(追記) 上記では、 Sk(n)=1k+2k+3k+・・・+nk の和の公式のいろいろな求め方
をみてきた。
S1(n)=1+2+3+・・・+n=n(n+1)/2=(1/2)n2+(1/2)n
S2(n)=12+22+32+・・・+n2=n(n+1)(2n+1)/6=(1/3)n3+(1/2)n2+(1/6)n
S3(n)=13+23+33+・・・+n3=n2(n+1)2/4=(1/4)n4+(1/2)n3+(1/4)n2
・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・
このような式をぼんやりと眺めていると、
S3’(n)=n3+(3/2)n2+(1/2)n=3・S2(n) (形式的に、n で微分!)
という関係式が成り立つことに気づく。このような関係は、
半径 r の円の面積 S=πr2 を、r で微分すると、S’=2πr は、半径 r の円周の長さを
表す
半径 r の球の体積 V=(4/3)πr3 を、r で微分すると、V’=4πr2 は、半径 r の球の表
面積を表す
に似ている。(これらの関係を高校時代に気づいたときには、とても感動した!)
実は、このような関係は、一般的に成立するらしい。浅井哲也さんが、2004年9月号の
数学セミナー(日本評論社)の中の、「エレガントな解答をもとむ」に出題されている。
(解答は、12月号を参照)
奇数 k に対して、 Sk’(n)=k・Sk-1(n) (形式的に、n で微分!)
「k は奇数」という限定条件が入るが、とても美しい関係である!
(補足) 平成19年8月13日付け
上記の公式を証明させる問題が、平成18年度東京大学 後期試験 数学に出題され
た。(一部表現等修正)
東京大学後期(2006)
p、n を自然数とする。次の問いに答えよ。
(1) p+1 次多項式 Sp(x)があって、
と表されることを示せ。
(2) q を自然数とする。(1)の多項式 S1(x)、S3(x)、・・・、S2q-1(x)に対して、
が恒等式となるような定数 a1、・・・、aq を q を用いて表せ。
(3) q を2以上の自然数とする。(1)の多項式 S2(x)、S4(x)、・・・、S2q-2(x)に対して、
が恒等式となるような定数 c と b1、・・・、bq-1 を q を用いて表せ。
(4) p を3以上の奇数とする。このとき、
を示せ。
(略解) (1)は、数学的帰納法を用いる。教科書でも示されている手法
(k+1)p+1−kp+1 の展開式を k=1、2、・・・、n について和をとる
ことにより、容易に示される。
(2)も、 F(x)=xq(x+1)q に対して、F(k)−F(k−1) の展開式を k=1、2、・・・、n
について和をとることにより示される。
(q+1)/2≦j≦q のとき、 aj=2qC2q-2j+1
上記以外のとき、 aj=0
(3)が少し大変だが、用いる手法は(1)(2)と同じである。 S2j(n)−S2j(n−1)=n2j
に注意して計算することにより、 c=2q であることが分かる。このとき、
G(x)=xq-1(x+1)q-1(2qx+q) に対して、G(k)−G(k−1) の展開式を
k=1、2、・・・、n について和をとることにより示される。
(q−1)/2≦j≦q−1 のとき、 bj=2(2j+1)qC2q-2j-1
上記以外のとき、 bj=0
(4)は、数学的帰納法を用いる。前述で示したように S3’(n)=3・S2(n) が成り立つので、
p=3 のとき、与えられた命題は成り立つ。2q−1 未満の奇数 p について、
が成り立つものとする。p=2q−1 のとき、 p≧3 なので、 q≧2
(2)で求められた恒等式 a1S1(x)+a2S3(x)+・・・+aqS2q-1(x)=xq(x+1)q
において、 q≧2 なので、 a1=0 である。a2S3(x)+・・・+aqS2q-1(x)=xq(x+1)q
の両辺を微分して、
a2・3S2(x)+・・・+aq-1・(2q−3)S2q-4(x)+aqS2q-1’(x)=qxq-1(x+1)q-1(2x+1)
(2)(3)の aj 、 bj の関係から、
b1S2(x)+・・・+bq-2S2q-4(x)+aqS2q-1’(x)=xq-1(x+1)q-1(2qx+q)
が成り立つ。(3)の恒等式と比較して、 aqS2q-1’(x)=bq-1S2q-2(x)
ここで、 aq=2qC1=2q 、 bq-1=2(2q−1)qC1=2q(2q−1) なので、
S2q-1’(x)=(2q−1)S2q-2(x) が成り立つ。
よって、p=2q−1 のときも命題は成り立つ。
したがって、3以上の奇数 p に対して、
が成り立つ。(終)
(コメント) 平成16年9月に出版されて、平成18年度入試に出題ということを考えると、多
分出題者は数学セミナーの問題を意識していたのかな?それとも独立に考えら
れた問題なのかな?
何れにしても、この問題を時間内に解くことはほぼ困難だということは言える。3題150分
の試験だが、受験生には辛い出題だ!
(追記) 平成17年7月25日付けで、神奈川県の S さんという方から、
kが奇数の時、 Sk’(n)=k・Sk-1(n) が成り立つが、これを一般化して、Sk(n)
から Sk+1(n) を求める方法を考えました
というメールを頂戴した。S さんが、この話題について、ある研究会誌に掲載された原稿の
PDFファイルも参考にさせていただいた。S さんに大変感謝いたします。
S さんが考えられた方法は次のようなものである。
(1) Sk−1(n) に対して、両辺を k 倍する。(最高次の係数が1になるようにする)
(2) (1)で得られたものを、 n に関して形式的に積分する。積分定数は、0
とする。
(3) 各項の係数の和が 1 になるように、適宜 1 次の項を付け加える。
この操作により、Sk(n) が得られる。
例 S1(n) = (1/2)n2+(1/2)n 、S2(n) = (1/3)n3+(1/2)n2+(1/6)n であるが、
S2(n) は、S1(n) から次のようにして得られる。
まず、(1) から、 2・S1(n) = n2+n で、 (2)から、(1/3)n3+(1/2)n2
1−{(1/3)+(1/2)}=1/6 なので、 (3)から、(1/6)n を付け加えて、
S2(n) = (1/3)n3+(1/2)n2+(1/6)n ( ← ぴったしカンカン!!)
例 S2(n) = (1/3)n3+(1/2)n2+(1/6)n 、S3(n) = (1/4)n4+(1/2)n3+(1/4)n2
であるが、S3(n) は、S2(n) から次のようにして得られる。
まず、(1) から、 3・S2(n) = n3+(3/2)n2+(1/2)n で、(2)から、
(1/4)n4+(1/2)n3+(1/4)n2
1−{(1/4)+(1/2)+(1/4)}=0 なので、 (3)で、付け加える項はない。
よって、 S3(n) = (1/4)n4+(1/2)n3+(1/4)n2 ( ← ぴったしカンカン!!)
(コメント) k が奇数の時、Sk’(n)=k・Sk-1(n) が成り立つので、たとえば、k=3 のとき、
S3’(n)=3・S2(n) が成り立つ。 S2(n) = (1/3)n3+(1/2)n2+(1/6)n なので、
S3’(n) = n3+(3/2)n2+(1/2)n となるが、両辺を積分して、
S3(n) = (1/4)n4+(1/2)n3+(1/4)n2+C (C は積分定数)
ここで、S3(1) = 1 なので、 C=0 でなければならない。
よって、 S3(n) = (1/4)n4+(1/2)n3+(1/4)n2
このように、k が奇数の時は、不定積分と初期条件(?) Sk(1) = 1 という事実から、
Sk(n) を決定することができるが、k が偶数の時は、上記のような漸化式が存在しないの
で、積分で Sk(n) を決定することは無理だろうと諦めていた。今回、S さんの鮮やかな方
法(特に、(3)!)に大きな感動を覚えました。
ここで、S さんの方法を若干検証しておこうと思う。
Sk(n)=1k+2k+3k+ ・・・ +nk において、
Sk(n)−Sk(n−1) = nk 、 Sk(1) = 1
という関係式が本質的である。
上式は、 n の恒等式なので、上式により、 n の多項式 Sk(n) は一意的に決定される。
ところで、n の多項式 Sk(n) は特徴ある性質を持っている。
上式において、 n=1 とおいてみると、 Sk(1)−Sk(0) = 1 で、Sk(1)=1 から、
Sk(0)=0 が成り立つとするのが自然だろう。
これは、Sk(n) が定数項を持たないことを意味する。
すなわち、Sk(n) は、n で割り切れる。
さらに、n=0 とおいてみると、 Sk(0)−Sk(−1) = 0 で、Sk(0)=0 から、
Sk(−1)=0 が成り立つとするのが自然だろう。
これは、Sk(n) が、n+1 で割り切れることを意味する。
したがって、任意の k に対して、Sk(n) は、n(n+1) で必ず割り切れる。
ところで、上記の関係式を形式的に n で微分すると、
Sk’(n)−Sk’(n−1) = k・nk−1 より、 Sk’(n)/k−Sk’(n−1)/k = nk−1
このとき、 Sk’(n)/k が、Sk−1(n) と何らかの関係があることは上述のことから予想
される。
ここで、n=1 のとき、Sk’(1)/k =1 が成り立つとは一般に言えないので、適当な定
数 a を定めることにより、
Sk’(1)/k + a = 1
となるようにできる。 n=1 のとき、上式が成り立つということは、n=0 のとき、0にな
らなければいけないので、定数 a は、Sk’(n)/k の定数項を打ち消すようにおけばよい。
このとき、 Sk−1(n) =Sk’(n)/k + a である。
これより、 k・Sk−1(n) =Sk’(n) + ak ・・・ (1)の手順
k・∫Sk−1(n)dn=Sk(n) + akn ・・・ (2)の手順
よって、 Sk(n)=k・∫Sk−1(n)dn + akn において、
Sk(1)=1 となるように、 a の値を定める。 ・・・ (3)の手順
以上の手順から、確かに、Sk(n) が求められるわけである。
例 S4(n) = (1/5)n5+(1/2)n4+(1/3)n3−(1/30)n から、S5(n) を求めてみるのも
楽しい!
5S4(n)=n5+(5/2)n4+(5/3)n3−(1/6)n を積分して、
(1/6)n6+(1/2)n5+(5/12)n4−(1/12)n2
1−((1/6)+(1/2)+(5/12)−(1/12))=0 より、
S5(n)=(1/6)n6+(1/2)n5+(5/12)n4−(1/12)n2
(コメント) S さんからいただいた資料を参考に話を展開したが、上記のような理解でいい
のかどうか少し不安である。是非 S さんの検証をお願いしたいと思う。
このように、単に、1+2+3+・・・+n=n(n+1)/2 という認識しかないと、つまらない
公式も、1+2+3+・・・+n=(1/2)n2+(1/2)n という見方もあることに気がつくと、俄然
公式が光り輝いて見える。因数分解されたものは美しいが、ただそれだけである。
(追記) 平成19年11月22日付け
上記の計算をさらに実行して、
S6(n)=(1/7)n7+(1/2)n6+(1/2)n5−(1/6)n3+(1/42)n
S7(n)=(1/8)n8+(1/2)n7+(7/12)n6−(7/24)n4+(1/12)n2
・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・
このとき、 S5(n)+S7(n)=2{S1(n)}4
という美しい関係式が成り立つ。(ひたすら計算すると確かめられる!)
この関係式は、ヤコビ(1804〜1851)により発見されたものだという。
(追記) 平成21年6月4日付け
上記の関係式について、HN「ポルテ」さんより、階差を利用した証明法をご教示いただい
た。ポルテさんに感謝します。
(k+1)4−(k−1)4=8k3+8k であることを利用して、
k4・{(k+1)4−(k−1)4}=k4・(8k3+8k)=8(k5+k7)
すなわち、 {k(k+1)}4−{k(k−1)}4=8(k5+k7) が成り立つ。
そこで、 an={n(n+1)}4 とおくと、上式より、k5+k7=(ak−ak-1)/8
よって、 S5(n)+S7(n)=(an−a0)/8={n(n+1)}4/8=2{S1(n)}4
(コメント) ひたすら計算すると思っていたら、上記のような簡明な解法があるのですね。感
動しました。ポルテさんによれば、階差を利用してΣ計算が綺麗に解ける問題は
とても爽快だとのこと。私も同感です!
ポルテさんの手法を用いると、1 から n までの自然数の立方の和の公式
13+23+33+・・・+n3 =n2(n+1)2/4
の証明が簡単にできることを、HN「FN」さんよりご教示いただいた。
(平成22年4月15日付け...泊まり込みの出張で直ぐレス出来ませんでしたf(^^;)。)
FNさんによれば、この公式の証明の中で、最も簡単な証明の一つだという。その簡明さ
に私も感動しました。FNさんに感謝します。
(証明) (k+1)2−(k−1)2=4k であることを利用して、
k2・{(k+1)2−(k−1)2}=4k3 即ち、{k(k+1)}2−{k(k−1)}2=4k3
が成り立つ。 k に1,2,・・・,n を代入し、 n 個の式を辺々加えて示される。(証終)
さらに、当HPがいつもお世話になっているS(H)さんから類題を提起いただいた。S(H)さ
んに感謝します。
S5(n)+10S7(n)+5S9(n)=16{S1(n)}5
これも、ポルテさんよりご教示いただいた階差の手法を用いれば容易だろう。
(k+1)5−(k−1)5=2(1+10k2+5k4) であることを利用して、
k5・{(k+1)5−(k−1)5}=2k5・(1+10k2+5k4)=2(k5+10k7+5k9)
即ち、 {k(k+1)}5−{k(k−1)}5=2(k5+10k7+5k9) が成り立つ。
そこで、 an={n(n+1)}5 とおくと、上式より、k5+10k7+5k9=(ak−ak-1)/2
よって、S5(n)+10S7(n)+5S9(n)=(an−a0)/2={n(n+1)}5/2=16{S1(n)}5
が成り立つ。
上記の2例から察するに、左辺が与えられて、それを計算した結果、右辺となるのではなく、
逆に、右辺から左辺が自然に導かれるようだ。
(追記) 階差数列の一般項の計算で、次の数列の和の公式
12+22+32+・・・+n2 =n(n+1)(2n+1)/6
よりも、むしろ、次の和の公式:12+22+32+・・・+(n−1)2 =n(n−1)(2n−1)/6
の方が多用される。このとき、 n(n+1)(2n+1)−n(n−1)(2n−1)=6n2
という等式の存在に気づかされる。このことから、次の積分の等式
の成り立つことが分かる。(両辺を、k に関して微分してみるとよい。)
したがって、「1 から n までの自然数の立方の和」について、積分を用いた新しい証明
法を得ることができる。すなわち、
において、
より、
したがって、 13+23+33+・・・+n3 =n2(n+1)2/4 が成り立つ。
(追記) 平成21年6月2日付け
「自然数の和」や「平方数の和」を求める場合に定積分を用いる手法は新鮮である。
たとえば、「自然数の和」を求める場合は、 F1(x)=x+1/2 について、
となることから、
ところで、F1(x)=x+1/2 は次のようにして求められる。
1次式 F1(x)=ax+b について、
が成り立つものとして、両辺を k で微分すると、 ak+b−{a(k−1)+b}=1
これより、 a=1 である。このとき、定積分の計算から、
k−1/2+b=k より、 b=1/2
よって、条件式を満たす1次式 F1(x) は、 F1(x)=x+1/2 と確定する。
同様にして、2次式 F2(x)=ax2+bx+c について、
が成り立つものとして、両辺を k で微分すると、
ak2+bk+c−{ a(k−1)2+b(k−1)+c}=2k
よって、 2ak−a+b=2k より、 a=b=1
このとき、定積分の計算から、 k2+c−1/6=k2 より、 c=1/6
よって、条件式を満たす2次式 F2(x) は、 F2(x)=x2+x+1/6 と確定する。
この F2(x)=x2+x+1/6 を用いて、
と、平方の和が求められる。
(追記) 当HP読者のS.Yさんから、上記の話題についてメールにてご投稿いただいた。
(平成29年10月23日付け)
道具としてxの関数 X[n](x) を以下のように定義します。
X[n](x)=x^n - (x-1)^n (nは自然数)
これを形式的に、X^n = x^n -(x-1)^n と表し、形式的なXのn次式と呼ぶことにします。
(実際は x のn-1次式です)
以下、∫は、x=0 からx までの積分、Σは、x=1 から n までの和とします。
このとき、次の式が成り立ちます。
(1) X=x-(x-1)=1
(2) ∫ X^n dx = 1/(n+1){X^(n+1)+(-1)^(n+1)}
(3) Σ X^m = n^m
(2)の補足 (xとXに注意してください)
∫ X^n dx
=1/(n+1)[x^(n+1)-(x-1)^(n+1)] [0,x]
=1/(n+1){x^(n+1)-(x-1)^(n+1)-0+(-1)^(n+1)}
=1/(n+1){X^(n+1)+(-1)^(n+1)}
(3)の補足
Σ X^m =Σ[1,n]{x^m -(x-1)^(m-1)}= n^m
x^m をXの形式的なm+1次式で表すことができれば、(3)より和の式が求まります。
例えば、x=1/2 X^2 + 1/2 X ⇔ Σx =1/2 n^2 + 1/2 n などです。
順に、和の式を求めていきます。
x =∫ 1 dx =∫ X dx [(1)より]=1/2 (X^2 +1) [(2)より]=1/2 X^2 + 1/2
=1/2 X^2 + 1/2X [(1)より] …(4)
∴ Σx = 1/2 n^2 +1/2 n [(3)より]
さらに、
x^2=∫ 2x dx =∫ (X^2+X) dx [(4)より]=1/3 (X^3-1)+1/2 (X^2 +1) [(2)より]
=1/3 X^3 +1/2 X^2 + 1/6=1/3 X^3 +1/2 X^2 + 1/6X [(1)より] …(5)
∴ Σx^2=1/3 n^3+ 1/2 n^2 +1/6 n [(3)より]
続けて、
x ^3= ∫ 3x^2 dx =∫ (X^3 +3/2X^2 +1/2X) dx [(5)より]
=1/4 (X^4+1)+1/2 (X^3-1)+1/4 (X^2+1) [(2)より]=1/4X^4 +1/2 X^3 +1/4 X^2
…(6)
∴ Σx^3=1/4n^4 +1/2 n^3 +1/4 n^2 [(3)より]
4乗の和を求めると、(説明略)
x ^4= ∫ 4x^3 dx =∫ (X^4 +2X^3 +X^2) dx =1/5 (X^5-1)+1/2 (X^4+1)+1/3 (X^3-1)
=1/5X^5 +1/2 X^4 +1/3 X^3-1/30=1/5X^5 +1/2 X^4 +1/3 X^3-1/30X
∴ Σx^4=1/5n^5 +1/2 n^4 +1/3 n^3-1/30n
一般に、m≧0 の整数について、x^m = S[m](X) (Xの形式的なm+1 次式) とおくと、
(イ) Σx^m=S[m](n) が成り立つ。
(ロ) x^m = ∫ mx^(m-1) dx より、S[m](X)= m ∫ S[m-1](X) dx を満たす。
(ハ) S[0](X)=1=X から、x^m =S[m](X)=m ! ∫∫… X (dx)^m と表すことができる。
これらの積分は(1)、(2)、(3)を用いて行う。
(1)の X=1 が ポイントだったと思いました。
(補足1) x^m =m ! ∫∫… (dx)^m とも表せる。
(補足2) x^m =S[m](X) の左辺にx = 1 を代入すると1 になる。
また、右辺の各項 の X^k=X[k](x)=x^k -(x-1)^kにx=1を代入すると X[k](1)=1
X=1 より S[m](X)= a[m+1] X^(m+1)+… +CX の形になり、Cは mS[m-1](X) を形式
的な不定積分をして(最後の項は …CX) 、x =1 (すなわちX^k=1)を代入して決定して
もよい。(Sさんの解法の根拠になります。)
(コメント) 解法が奇麗ですね!S.Yさんに感謝します。
(追記) 平成19年3月7日付け
上記で、 n(n+1)(2n+1)−n(n−1)(2n−1)=6n2 という等式の存在が
脚光を浴びたが、この等式を背景とする問題が平成19年度筑波大学前期入試で出題され
た。筑波大学は毎年質の高い入試問題を出題するという点でいつも注目している!
筑波大学前期(2007)
(1) 一般項 an が、 an3+bn2+cn で表される数列 { an }において、
n2=an+1−an (n=1、2、3、・・・)
が成り立つように、定数 a 、 b 、 c を定めよ。
(2) (1)の結果を用いて、
となることを示せ。 ((3)は略)
上記の等式によれば、(1)は、 a=1/3 、 b=−1/2 、 c=1/6 と即答だろう!
(ただし、入試問題の答案としては、恒等式の性質を用いてきっちり論述しないと点数には
ならないかもしれない。)
(2)も、問題の背景に、 n≧1 に対して
ということが分かっていれば、
(a2−a1)+(a3−a2)+・・・(an+1−an)=an+1−a1=an+1
から、明らかに等式が成り立つ訳だが、やはり、ここもグッと堪えて、細々とした計算を答
案に示してあげないと合格点はつかないだろう。
(追記) 平成22年3月14日付け
当HPがいつもお世話になっているHN「GAI」さんが、「ただひたすらの方法」と題して
「1 から 10 までの自然数の10乗の和」について、ひたすら紙と鉛筆だけを使用し、
手間暇かけて求めていく方法について、平成22年3月10日付けで当HPの掲示板「出
会いの泉」に書き込まれたものを紹介したいと思う。
計算手法のコツ
次の2つの公式もどきを覚えておく。
<公式1> 1k+2k+3k+・・・+nk={(n+1+B)k+1−Bk+1}/(k+1)
<公式2> Bn は (B+1)k=Bk の展開式 (k>1)から漸次決定していく。
実際の計算は次のように行う。
(1) まず、公式2 において、k=2 とし、(B+1)2=B2 を展開して、B1=−1/2
同様に、 (B+1)3=B3 を展開して、B2=1/6
(B+1)4=B4 を展開して、B3=0
(B+1)5=B5 を展開して、B4=−1/30
・・・・・・・・・・・・・・・・・・・・・・・・・
これを地道に続けていくと、ベルヌーイ数が出現する。
B5=0 、B6=1/42 、B7=0 、B8=−1/30 、B9=0 、B10=5/66 、・・・
そこで、試しに、公式1を用いて、自然数の平方の和を求めてみよう。
12+22+32+・・・+n2
={(n+1+B)3−B3}/3
={(n+1)3+3(n+1)2B+3(n+1)B2}/3
=(n+1){2(n+1)2−3(n+1)+1}/6=n(n+1)(2n+1)/6
同様の計算を、
110+210+310+・・・+1010={(11+B)11−B11}/11
について行えばよい。
{(11+B)11−B11}/11
=1110+1110B+5・119B2+15・118B3+30・117B4+42・116B5
+42・115B6+30・114B7+15・113B8+5・112B9+11B10
=1110−1110/2+5・119/6−117+115−113/2+5/6
=1110/2+5・119/6−117+115−113/2+5/6
ここで、 1110=25937424601 、119=2357947691 、117=19487171
115=161051 、113=1331 なので、代入して計算すると、
14,914,341,925 となる。
(コメント) 上記の計算の検証に、Excel さんのお世話になりました...f(^^;)
また、上記の計算に次のホームページを参考にさせていただきました。
数学的思考と科学的思考 (1998年2学期 木曜1限)
ベルヌーイ数のはなし −整数のべき乗の和と交代順列−
(兵庫教育大学の大学院生向けの講義録) 兵庫教育大学自然系数学 松山 廣 教授
また、GAI さんの解法とは違って、Sさんの手法を活用して、順次
S1(n) = (1/2)n2+(1/2)n
S2(n) = (1/3)n3+(1/2)n2+(1/6)n
S3(n) = (1/4)n4+(1/2)n3+(1/4)n2
S4(n) = (1/5)n5+(1/2)n4+(1/3)n3−(1/30)n
S5(n)=(1/6)n6+(1/2)n5+(5/12)n4−(1/12)n2
S6(n)=(1/7)n7+(1/2)n6+(1/2)n5−(1/6)n3+(1/42)n
S7(n)=(1/8)n8+(1/2)n7+(7/12)n6−(7/24)n4+(1/12)n2
S8(n)=(1/9)n9+(1/2)n8+(2/3)n7−(7/15)n5+(2/9)n3−(1/30)n
S9(n)=(1/10)n10+(1/2)n9+(3/4)n8−(7/10)n6+(1/2)n4−(3/20)n2
S10(n)=(1/11)n11+(1/2)n10+(5/6)n9−n7+n5−(1/2)n3+(5/66)n
と、ほぼ機械的に求められるので、こちらを利用してもいいだろう。
このとき、 S10(10)=14,914,341,925 であることが分かる。
自然数のべき乗和の式について、GAI さんから問題提起です。(平成24年1月29日付け)
S(k)=1k+2k+3k+・・・+nk を表す式は、それぞれ公式として記述されているが、これ
らには、次の性質が成り立つことを知りました。
k:奇数の場合 → S(k)は、S(1)の多項式で表される。
k:偶数の場合 → S(k)は、S(2)で割れ、その商は、また、S(1)の多項式で表される。
これを、具体的に式にしてみました。ただし、表記が便利なように、A=S(1)、B=S(2)で表
します。
S(3)=A2
S(4)=(6A−1)B/5
S(5)=(4A−1)A2/3
S(6)=(12A2−6A+1)B/7
S(7)=(6A2−4A+1)A2/3
S(8)=(40A3−40A2+18A−3)B/15
S(9)=(2A−1)(8A2−6A+3)A2/5
S(10)=(2A−1)(24A3−28A2+20A−5)B/11
・・・・・・・・・・・・・・・
これらを一般的に表せる式として作れないものでしょうか?ご存知の方はお知らせ下さい。
よおすけさんからのコメントです。(平成24年1月29日付け)
上記に、ある程度のことは書かれていると思います。または、「お気に入り集」から、らす
かるさんのサイト「k乗和の公式」に行くか…。
空舟さんからのコメントです。(平成24年1月30日付け)
k:奇数の場合 → S(k)は、S(1)の多項式で表される。
k:偶数の場合 → S(k)は、S(2)で割れ、その商は、また、S(1)の多項式で表される。
という性質は、S(k)=f(n) が、 f(n+1)−f(n)=(n+1)k 、f(0)=0 を満たすこと等を考え
ると証明できました。それほど難しくないと思いますのでぜひ考察を...。
これらを一般的に表せる式というのは難しそうですが、検索してみると、ベルヌーイ数とい
う数を使って、一般的に、(nの多項式としてですが)表す方法があるそうです。
S(H)さんからのコメントです。(平成24年1月30日付け)
あの「Donald E. Knuth」先生が、関連して;
「Johann Faulhaber and Sums of Powers」、「Johann Faulhaber」、「数列の和」
(追伸) Σ計算をぐっと睨んで、個別に、例えば、169+269+・・・+n69 からS(70)
また、S(15) 、S(17) 、・・・、S(169) をぐっと睨んで。
当HPの読者のKSさんより、平成24年7月28日付けで累乗和の公式について、メールを
頂いた。
(1) ベルヌーイの公式を利用
1k+2k+3k+・・・+(n−1)k
={(nk+1+k+1C1B1nk+・・・+k+1CkBkn}/(k+1)
={(n+B)k+1−Bk+1}/(k+1) (← 形式的表記)
ただし、 B0=1 、Σk=0〜n n+1CkBk=0 (n=1、2、3、・・・)
上記で計算したように、順次 (B+1)n=Bn を展開することにより、ベルヌーイ数が得
られる。
B0=1 、B1=−1/2 、B2=1/6 、B3=0 、B4=−1/30 、B5=0 、B6=1/42
B7=0 、B8=−1/30 、B9=0 、B10=5/66 、・・・
ベルヌーイ数を用いて、
xex/(ex−1)=Σn=0〜∞ Bnxn/n!
Σm=0〜n mk=Σm=0〜k kCmBm・(n+1)k+1-m/(k+1−m)
例 k=3 のとき、
Σm=0〜n m3
=Σm=0〜3 3CmBm・(n+1)4-m/(4−m)
=B0・(n+1)4/4+3C1B1・(n+1)3/3+3C2B2・(n+1)2/2+3C3B3・(n+1)
=(n+1)4/4−(n+1)3/2+(n+1)2/4
(2) 微分の利用
Sk’(n)=k・Sk-1(n)+Sk’(0) ただし、 Sk’(0)=Bk
下位の次数から積分し、積分定数を Sk(1)=1 により確定する。また、ベルヌーイ数
Bk=0 (kは3以上の奇数) を利用する。
(コメント) 類似の計算が上記にあるので参考にしてください。
(3) 階乗関数
係数が、第一スターリング数から、第二スターリング数にかわる昇階乗のときは、符号が
交互になる。
[X]k=X(X+1)・・・(X+k−1) 、 [X]k=X(X−1)・・・(X−k+1)
Σx=1〜n[X]k=[n]k+1/(k+1) 、 Σx=1〜n[X]k=[n]k+1/(k+1)
[X]n=Σk=0〜n nCkXk 、 [X]n=Σk=0〜n (−1)n-knCkXk
Xn=Σk=0〜n (−1)n-knCk[X]k 、 Xn=Σk=0〜n nCk[X]k
Σx=1〜nXk=Σk=0〜n (−1)n-knCk[n]k/(k+1)
Σx=1〜nXk=Σk=0〜n nCk[n+1]k/(k+1)
(4) 連立方程式の利用(展開式)
Σm=1〜n mk=a・nk+1+b・nk+・・・+c とおいて、n=1、2、・・・、k+1を代入し、a、
b、・・・、cの連立方程式を解く。
(5) 教科書にある方法
(n+1)k−nk=k・nk-1+・・・+1
nk−(n−1)k=k・(n−1)k-1+・・・+1
・・・・・・・・・・・・・・・・・・・・・・
2k−1k=k・1+・・・+1
辺々加えて、 (n+1)k−1=kΣm=1〜n mk-1+・・・+n
後は、k−1乗以下の和の公式を利用する。
(6) ファウルハーバーの定理
χ=n(n+1)/2 とおくと、
Σr≧1 k2r+1=(aχr+1+・・・+bχ2の形)
Σr≧1 k2r=(2n+1)(aχr+・・・+bχの形)
(7) 高階差分法を利用(組合せの数で表現する)
Σk=1〜n k3=nC1+7・nC2+12・nC3+6・nC4
ただし、係数は、
nDr=n-1Dr-1+r・n-1Dr 、nD0=1 、nDn=n! 、nDr=0 (n<r)
(8) 組合せの数(C)で対称的に表現する。
オイラーの係数で与えられる。
n2=nC1+2・nC2=n+1C2+nC2 より、 Σk=1〜n k2=n+2C3+n+1C3
n3=nC1+6・nC2+6・nC3=nC3+4・n+1C3+6・n+2C3 より、
Σk=1〜n k3=n+1C4+4・n+2C4+n+3C4
ここで、 mDr=(m+1−r)m-1Dr-1+r・m-1Dr 、mD1=mDm=1 (オイラーの数)
このとき、 Σm=1〜n mk=kD1・n+1Ck+1+kD2・n+2Ck+1+・・・+kDk・n+kCk+1
(9) 一番目と二番目次数はすぐにわかるので、他の係数は以下の比較で求めることがで
きる
例 Σk=1〜n k4=n5/5+n4/2+an3+bn2+cn
=(n+1)5/5−(n+1)4/2+a(n+1)3+b(n+1)2+c(n+1)
を解いて、a=1/3 、b=0 、c=−1/30
以上から、
Σk=1〜n k3=(n+1)4/4−(n+1)3/2+(n+1)2/4 (ベルヌーイ)
(スターリング)
Σk=1〜n k3=nC1+7・nC2+12・nC3+6・nC4 (スターリング)
Σk=1〜n k3=n+1C4+4・n+2C4+n+3C4 (オイラーの対称性)
Σk=1〜n k3=n4/4+n3/2+n2/4 (未定係数)
Σk=1〜n k3=χ2 (ファウルハーバー)
(コメント) K.S.さん、いつも投稿ありがとうございます。「ファウルハーバー」という名前は
初見でした!
(補足) 令和5年7月24日付け
上記では、差分の形で数列の和を求めることを種々行った。このことを一般化して、次の
ようなマニュアルがあることを最近知ることができた。
連続するN個の整数の積を、GN(k)=k(k+1)(k+2)・・・(k+N−1)で表すとする。
このとき、 G1(k)=k、G2(k)=k(k+1)、・・・ である。
「整数値多項式」で示した等式 : Fn+1(X+1)−Fn+1(X)=(n+1)Fn(X+1) から
GN+1(k)−GN+1(k−1)=(N+1)GN(k)
が成り立つ。
例 N=1 のとき、 k(k+1)−(k−1)k=2k
N=2 のとき、 k(k+1)(k+2)−(k−1)k(k+1)=3k(k+1)
・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・
この例からも分かるように、Σk、Σk2、・・・ の公式が順次得られる。
また、GN+1(k)−GN+1(k−1)=(N+1)GN(k) において、
k=1、2、・・・、n を順次代入し、辺々加えることにより、
GN+1(n)−GN+1(0)=(N+1)(GN(1)+GN(2)+・・・+GN(n))
もちろん、GN+1(0)=0 なので、
Σk=1〜n GN(k)=GN+1(n)/(N+1)
が成り立つ。(何となく、不定積分の公式に似ていますね!)
この公式を用いると、
Σk=1〜n k=Σk=1〜n G1(k)=G2(n)/2=n(n+1)/2
Σk=1〜n k(k+1)=Σk=1〜n G2(k)=G3(n)/3=n(n+1)(n+2)/3
Σk=1〜n k(k+1)(k+2)=Σk=1〜n G3(k)=G4(n)/4=n(n+1)(n+2)(n+3)/4
Σk=1〜n k(k+1)(k+2)(k+3)
=Σk=1〜n G4(k)=G5(n)/5=n(n+1)(n+2)(n+3)(n+4)/5
Σk=1〜n k(k+1)(k+2)(k+3)(k+4)
=Σk=1〜n G5(k)=G6(n)/6=n(n+1)(n+2)(n+3)(n+4)(n+5)/6
などが得られる。
よって、
Σk=1〜n k2=Σk=1〜n k(k+1)−Σk=1〜n k
=n(n+1)(n+2)/3−n(n+1)/2=n(n+1)(2n+4−3)/6=n(n+1)(2n+1)/6
Σk=1〜n k3=Σk=1〜n (k−1)k(k+1)+Σk=1〜n k
=Σk=1〜n-1 k(k+1)(k+2)+Σk=1〜n k
=(n−1)n(n+1)(n+2)/4+n(n+1)/2=n2(n+1)2/4
Σk=1〜n k4=Σk=1〜n (k4−k2)+Σk=1〜n k2
=Σk=1〜n (k−1)k(k+1)(k+2)−2Σk=1〜n (k−1)k(k+1)+Σk=1〜n k2
=Σk=1〜n-1 k(k+1)(k+2)(k+3)−2Σk=1〜n-1 k(k+1)(k+2)+Σk=1〜n k2
=(n−1)n(n+1)(n+2)(n+3)/5−(n−1)n(n+1)(n+2)/2+n(n+1)(2n+1)/6
=(n−1)n(n+1)(n+2)(2n+1)/10+n(n+1)(2n+1)/6
=n(n+1)(2n+1)(3n2+3n−1)/30
Σk=1〜n k5=Σk=1〜n (k5−k3)+Σk=1〜n k3
=Σk=1〜n k(k+1)(k−1)(k2−4+4)+Σk=1〜n k3
=Σk=1〜n (k−2)(k−1)k(k+1)(k+2)+4Σk=1〜n (k−1)k(k+1)+Σk=1〜n k3
=Σk=1〜n-2 k(k+1)(k+2)(k+3)(k+4)+4Σk=1〜n-1 k(k+1)(k+2)+Σk=1〜n k3
=(n−2)(n−1)n(n+1)(n+2)(n+3)/6+(n−1)n(n+1)(n+2)+n2(n+1)2/4
=n2(n+1)2(2n2+2n−1)/12
(追記) 平成25年4月9日付け
Sk(n)=1k+2k+3k+・・・+nk について、名古屋大学で次のような問題が出題された。
名古屋大学 理系(2013)
k、m、nは整数とし、n≧1とする。nCkを二項係数として、Sk(n)、Tm(n)を以下のように
定める。
Sk(n)=1k+2k+3k+・・・+nk 、Sk(1)=1 (k≧0)
Tm(n)=mC1S1(n)+mC2S2(n)+mC3S3(n)+・・・+mCm-1Sm-1(n)
=Σk=1〜m-1 mCkSk(n) (m≧2)
(1) Tm(1)とTm(2)を求めよ。
(2) 一般のnに対して、Tm(n)を求めよ。
(3) pが3以上の素数のとき、Sk(p−1) (k=1、2、3、・・・、p−2) は、pの倍数である
ことを示せ。
#和の公式を導くときに出てくる式に関する問題で一見難しそうだが、頑張ってみる価値は
ありそう...ということで、詳解に挑戦してみた。
(解)(1) Tm(1)=mC1S1(1)+mC2S2(1)+mC3S3(1)+・・・+mCm-1Sm-1(1)
=mC1+mC2+mC3+・・・+mCm-1
=mC0+mC1+mC2+mC3+・・・+mCm-1+mCm−(mC0+mCm)
=2m−(1+1)=2m−2
Tm(2)=mC1S1(2)+mC2S2(2)+mC3S3(2)+・・・+mCm-1Sm-1(2)
ここで、 Sk(2)=1k+2k=1+2k なので、
Tm(2)=(mC1+mC2+mC3+・・・+mCm-1)
+mC1・2+mC2・22+mC3・23+・・・+mCm-1・2m-1
=(2m−2)+mC1・2+mC2・22+mC3・23+・・・+mCm-1・2m-1
X=mC1・2+mC2・22+mC3・23+・・・+mCm-1・2m-1 とおくと、
X=(1+2)m−mC0−mCm・2m=3m−1−2m
よって、 Tm(2)=(2m−2)+3m−1−2m=3m−3
(2) (1)より、Tm(n)=(n+1)m−(n+1) と類推される。この類推が正しいことを数学
的帰納法により示す。
n=1 のとき、(1)より明らかに成り立つ。
n=k(k≧1)のとき成り立つと仮定する。すなわち、Tm(k)=(k+1)m−(k+1)
このとき、
Tm(k+1)=Σ k=1〜m-1 mCkSk(k+1)
=Σ k=1〜m-1 mCk{Sk(k)+(k+1)k}
=Tm(k)+Σ k=1〜m-1 mCk(k+1)k
=Tm(k)+{1+(k+1)}m−1−(k+1)m
=(k+1)m−(k+1)+(k+2)m−1−(k+1)m=(k+2)m−(k+2)
よって、n=k+1のときも成り立つ。
以上から、すべての自然数nについて、 Tm(n)=(n+1)m−(n+1) が成り立つ。
(3) Sk(p−1) (k=1、2、3、・・・、p−2) は、pの倍数であるであることを、kについ
ての数学的帰納法で示す。
k=1 のとき、 S1(p−1)=1+2+・・・+(p−1)=p(p−1)/2
ここで、pは3以上の素数すなわち奇数であるので、(p−1)/2は整数。
よって、S1(p−1)はpの倍数である。
k=1、2、3、・・・、q (q≦p−3)のとき、Sk(p−1)はpの倍数であると仮定する。
Tm(n)=mC1S1(n)+mC2S2(n)+mC3S3(n)+・・・+mCm-1Sm-1(n)
=Σ k=1〜m-1 mCkSk(n) (m≧2)
において、m=q+2、n=p−1 とおくと、
Tq+2(p−1)=q+2C1S1(p−1)+q+2C2S2(p−1)+・・・+q+2Cq+1Sq+1(p−1)
ここで、(2)より、Tq+2(p−1)=pq+2−p なので、
q+2Cq+1Sq+1(p−1)=pq+2−p−(q+2C1S1(p−1)+・・・+q+2CqSq(p−1))
帰納法の仮定により、右辺は、pの倍数。すなわち、(q+2)Sq+1(p−1)は、pの倍数。
ここで、q+2≦p−1 で、p は素数なので、q+2 は、p で割り切れない。
よって、Sq+1(p−1)は、pの倍数でなければならない。
これは、k=q+1のときも命題が成り立つことを示す。
以上から、pが3以上の素数のとき、
Sk(p−1) (k=1、2、3、・・・、p−2) は、pの倍数である
ことが示された。 (終)
自然数の冪乗和計算について、GAI さんからの続報です。(平成27年6月18日付け)
自然数の冪乗和公式を眺めていたら、
奇数冪では、因数 n(n+1) 、偶数冪では、因数 n(n+1)(2n+1)
を共通に含むことから、
奇数冪では、s[1]=1+2+3+・・・+n 、偶数冪では、s[2]=12+22+32+・・・+n2
だけを頼りに求まらないものかと挑戦してみました。(冪数が大きくなると公式も結構複雑)
1+2+3+・・・+n=n(n+1)/2の値をxとしたら、自然数の奇数冪の和
s[r]=1r+2r+3r+・・・+nr (rは奇数)
がxの値を元に次の計算で求まる。
s[3]=x2 、s[5]=(1/3)x2(-1+4x) 、s[7]=(1/3)x2(1-4x+6x2) 、s[9]=(1/5)x2(-1+2x)(3-6x+8x2)
s[11]=(1/3)x2(5-20x+34x2-32x3+16x4)
s[13]=(1/105)x2(-691+2764x-4720x2+4592x3-2800x4+960x5)
s[15]=(1/3)x2(105-420x+718x2-704x3+448x4-192x5+48x6)
s[17]=(1/45)x2(-10851+43404x-74220x2+72912x3-46880x4+21120x5-6720x6+1280x7)
s[19]=(1/105)x2(219335-877340x+1500334x2-1474592x3
+950544x4-433536x5+145824x6-35840x7+5376x8)
・・・・・・・・・
これも複雑といえば複雑だな〜。これに関連させて、
12+22+32+・・・+n2=n(n+1)(2n+1)/6 の値をyとしたとき、自然数の偶数冪の和
s[s]=1s+2s+3s+・・・+ns (sは偶数)
がyの値を元に計算できませんかね?いろいろ試みたんですが、頓挫しました。
S(H)さんからのコメントです。(平成27年6月18日付け)
x = n (1 + n)/2、y = n (1 + n) (1 + 2 n)/6 として、Q(x,y)の元として、有理式で表現可能です!
例えば、s[4]=(-7y + 972y3)/(5 (7 + 42x + 144x2))
s[6]=(5y - 22xy + 108y3)/(7 (5 + 8x))
s[8]=(-21y + 86xy - 360y3 + 600xy3)/(5 (21 + 40x))
s[10]=(2125y - 8550xy + 30420y3 - 41040xy3 + 55728y5 -31104xy5)/(11 (425 + 840x - 576y2))
s[10]= (1/44)(20y - 120xy + 315x2y - 387y3 + 216xy3)∈Q[x,y] とも書ける。
GAI さんからのコメントです。(平成27年6月18日付け)
xも用いていいなら、
s[4]=(6x-1)y/5 、s[6]=(12x2-6x+1)y/7 、s[8]=(40x3-40x2+18x-3)y/15 、
s[10]=(2x-1)(24x3-28x2+20x-5)y/11
の方が遥かにスッキリです。
GAI さんからのコメントです。(平成27年6月19日付け)
確か以前、こんな事を考えたことがあったなと思いましたが、どうして出したか、今思い出
そうとしてもできませんでした。確かいろいろ苦労しながら一個ずつ求めていった記憶だけ
があります。あえて言えば、"アイデアル"を用いたです。
ところで、S(H)さんが求められた式は複雑ですが、コンピュータで確認しますとどれもちゃ
んと s[4]、s[6]、・・・、s[16] を満たしました。これらを求めるのに"イデアル"を利用されたと
のことですが、当方全くこの概念に疎くチンプンカンプンです。しかし、その威力は強烈であ
ることは感じられます。S(H)さんがこの結果を出されたのはとても手作業では有り得ず、何
らかのソフトを利用されたと推測します。もし、Mathematica、もしくは、PARI/GP のソフトで
この計算をさせるためには、どのようなコマンドでやったら可能なのかを教えてください。こ
れが自由自在に手に入れることができたら、活用範囲が一気に広がりそうです。
S(H)さんからのコメントです。(平成27年6月18日付け)
失礼いたしました。
s[10]=(1/44) (20y - 120xy + 315x2y - 387y3 + 216xy3)
s[12]=(-11056y + 66336xy - 176205x2y + 226485y3 - 158760xy3 + 60480x2y3)/7280
s[14]=(1/320)(2240y - 13440xy + 35805x2y - 46533y3 + 34344xy3 -17280x2y3 + 5184y5)
s[16]=(1/21760)(-925952y + 5555712xy - 14811615x2y
+ 19302615y3 - 14429880xy3 + 7767360x2y3 - 3058560y5 + 622080xy5)
S(H)さんからのコメントです。(平成27年6月19日付け)
s:偶数の場合、例えば、s[2*3] は、 s[6]=(1 + 2n)(1/21)(x - 6 x2 + 12x3)
らすかるさんからのコメントです。(平成27年6月19日付け)
x=[(9y2)(1/3)/2] ([ ]はガウス記号)なので、
s[4]=(6[(9y2)(1/3)/2]-1)y/5
s[6]=(12[(9y2)(1/3)/2]^2-6[(9y2)(1/3)/2]+1)y/7
s[8]=(40[(9y2)(1/3)/2]^3-40[(9y2)(1/3)/2]^2+18[(9y2)(1/3)/2]-3)y/15
s[10]=(2[(9y2)(1/3)/2]-1)(24[(9y2)(1/3)/2]^3-28[(9y2)(1/3)/2]^2+20[(9y2)(1/3)/2]-5)y/11
のようには書けます。ガウス記号を使わずに、xをyで厳密に表すと、
x={{7776y2-1+72y√(11664y2-3)}(1/3)+{7776y2-1-72y√(11664y2-3)}(1/3)-1}/24
ですから、これを使えばガウス記号を使わずにも書けます。
よって、例えば、s[4]をyで表すと
s[4]={{7776y2-1+72y√(11664y2-3)}(1/3)+{7776y2-1-72y√(11664y2-3)}(1/3)-5}y/20
のようになりますので、おそらく簡単な形にはならないでしょうね。
GAI さんからのコメントです。(平成27年6月19日付け)
x、yの間に知るよしもないが、まず次の関係式があることに驚きました。
(自分も2つの間の関係式を探っていたんですが、わけの分からない3次方程式が出現していました。)
x=[(9y2)(1/3)/2] ([ ]はガウス記号)
確かに、コンピュータで数値確認していくと、ピタリと一致するではありませんか。
次の、ガウス記号を使わない場合も、一体どこからやってきたか検討もつきませんが、平
方根があるので、nを実数の世界で動かしてやはりコンピュータで超複雑計算をさせてみた
ら、中には、「29.99999」などど返すこともありますが、xの列と一致するものが出現しました。
全く想像もしていない関係式で繋がれるものだと二度驚きました。
また、s[4]の式を見て、これでは覚えておく気になりませんから、もうこれ以上これに挑戦す
ることは時間の無駄であることがハッキリ認識できました。いつもモヤモヤをスッキリと払い
のけて下さるので有り難いです。
らすかるさんからのコメントです。(平成27年6月19日付け)
xとyの関係式は、y=n(n+1)(2n+1)/6=x(2n+1)/3 から、
9y2=x2(4n2+4n+1)=x2(8x+1)=8x3+x2
となります。 8x3<8x3+x2<8(x+1)3 から、8x3<9y2<8(x+1)3 より、x3<9y2/8<(x+1)3
すなわち、 x<(9y2)(1/3)/2<x+1 なので、 x=[(9y2)(1/3)/2] となり、また、
8x3+x2-9y2=0
という三次方程式をゴリゴリ解いた結果が、
x={{7776y2-1+72y√(11664y2-3)}(1/3)+{7776y2-1-72y√(11664y2-3)}(1/3)-1}/24
です。
GAI さんからのコメントです。(平成27年6月19日付け)
この部分見直したら、 y=x(2n+1)/3 から、 n=3y/x-1 とやってしまっていた。
(n=(3y/x-1)/2 だった!)
これを、x=n(n+1)/2 へ代入したので、8x3+x2-3y2=0 になっていた。これを一応
3次方程式で解いたが(ここはコンピュータソフトのお世話になった。)余りにも複雑怪奇
な式が出現したので、ここから怖じ気づいた。でも計算ミスをしなくても、あのガウス記号
でxを表現する技は思ってもできない。
GAI さんからのコメントです。(平成27年6月19日付け)
聞きたかったのは偶数冪の時のやり方の方でした。ちなみに、奇数冪の場合は次の様な
方法で求めました。
恒等式 (k+1)2-(k-1)2=4k より、両辺にk2を掛けて、 {(k+1)k}2-{k(k-1)}2=4k3
この式に、k=1、2、3、・・・、n としたものを全て加えて、 {2s[1]}2=4s[3]
よって s[3]=s[1]2=x2・・・・・@
同様に、 (k+1)3-(k-1)3=2(1+3k2) の両辺にk3を掛けて、{(k+1)k}3-{k(k-1)}3=2(k3+3k5)
これから同様な手順で、 s[3]+3s[5]=4s[1]3=4x3・・・・・A
以下、同じ手順を冪を一つずつ上げながら進むと、
s[5]+s[7]=2s[1]4=2x4・・・・・B (ヤコビの関係式というらしい。)
s[5]+10s[7]+5s[9]=16x5・・・・・C
3s[7]+10s[9]+3s[11]=16x6・・・・・D
s[7]+21s[9]+35s[11]+7s[13]=64x7・・・・・E
s[9]+7s[11]+7s[13]+s[15]=16x8・・・・・F
s[9]+36s[11]+126s[13]+84s[15]+9s[17]=256x9・・・・・G
5s[11]+60s[13]+126s[15]+60s[17]+5s[19]=256x10・・・・・H
以上から@→A→B→・・・→Hを解くことで奇数冪の場合の結果を、xを用いて出しました。
これを、教わったプログラムで作業しますと、2時間ぐらいかけてやっていた作業が2分くら
いで片づいてしまいました。
所で偶数冪の場合がどうしたものか分からなくなったので、これをプログラムで出したいの
ですが偶数バージョンの方を教えて下さい。
DD++さんからのコメントです。(平成27年6月19日付け)
実際に試してないので、どの程度楽になるかわかりませんが、n=(m-1)/2 と変換すると見
通しが立ちやすくなるのでは、と思います。
つまり、x = n(n+1)/2 = (m2-1)/8 より、 m2 = 8x+1 で、y = n(n+1)(2n+1)/6 = m(m2-1)/24
として、奇数次項は y×偶数次 という形にして、残りは、m2 に上を代入して終わり。
次数ごとの元々の式があればという前提ですが。
GAI さんからのコメントです。(平成27年6月20日付け)
指示された手順を実行してみました。ただし、最後に出た結果に3/Aを掛けることになるよ
うですが、s[4]〜s[10]まで出していた結果が次々に現れてきました。コンピュータで自動化さ
せる力がないので、結構手間がかかりましたが、上記のテクニックであのnだけの複雑きわ
まりない式が、x、y だけの式に置き換わっていけることに感動します。
一応、s[20]までの作業を終えましたので掲載してみます。(ミスしていないかな?)なお、x、
y では紛らわしかったので、s[1]=n(n+1)/2=A、s[2]=n(n+1)(2n+1)/6=B にしています。
s[4]= (-1+6 A) B/5
s[6]= (1-6 A+12 A^2) B/7
s[8]= (-3+18 A-40 A^2+40 A^3) B/15
s[10]= (-1+2 A) (-5+20 A-28 A^2+24 A^3) B/11
s[12]= (-691+4146 A-9440 A^2+11480 A^3-8400 A^4+3360 A^5) B/455
s[14]= (105-630 A+1436 A^2-1760 A^3+1344 A^4-672 A^5+192 A^6) B/15
s[16]= (-3617+21702 A-49480 A^2+60760 A^3-46880 A^4+24640 A^5-8960 A^6+1920 A^7) B/85
s[18]= (219335-1316010 A+3000668 A^2-3686480 A^3+2851632 A^4
-1517376 A^5+583296 A^6-161280 A^7+26880 A^8) B/665
s[20]= (-3666831+22000986 A-50165840 A^2+61638680 A^3-47710320 A^4
+25465440 A^5-9926400 A^6+2914560 A^7-633600 A^8+84480 A^9)
B/1155
GAI さんが実験されました。(平成27年6月20日付け)
今、s[5]=15+25+35+・・・+1005 をどうしても求めねばならない場面に遭遇したと想像してみ
よう。しかし、歳のせいもあり、s[5]の公式があったことは見たことがあるも、その正確な式は
思い出せない。手元に計算機もコンピュータもなく、あるのは紙と鉛筆のみ。
勿論、S(H)さんの方法(参考1、参考2)で、A[0]=0とあらかじめ分かっていたとしても6変
数の連立を解くことになる。実際上記の方針でコツコツ進んでみた。
まず、15、15+25、15+25+35、15+25+35+45、15+25+35+45+65 の計算と、出来上がる6
個の連立を解く作業は、実際やってみると結構面倒ですよ。
(コンピュータがあれば行列でチョチョイと行けますが・・・)
もし行列でやろうとしても、手作業で逆行列を出すのは面倒(上手い方法があるのかな?)。
これに対し奇数冪の和は、S[1]の値のみでOKという認識さえあれば、S[3]=S[1]2 はお馴
染みだし、(k+1)3-(k-1)3=2+6k もチョロイ。後は、これから、s[3]+3s[5]=4s[1]3 も容易く導ける。
正確に記憶していなくても、この2つを組み合わせれば、 s[5]=(4s[1]-1)s[1]2/3
(s[5]=n6/6+n5/2+5n4/12-n2/12=n2(n+1)2(2n2+2n-1) に比べ遥かに簡潔)
は容易に入手可能。
当然、s[1](n=1〜100)=5050 なので、(4*5050-1)*50502/3=171708332500
(ここも1002*1012*(20000+200-1)の計算より楽?)
で、めでたしめでたし。
私は断然元のs[5]をnの式で覚えるより(忘れても復元するのにも骨が折れる。)、トイレに
s[5]=(4x-1)*x2/3 を貼り付けて万一に備えます。
「張り紙無用」と題して、GAI さんからのコメントです。(平成27年6月21日付け)
パスカルの三角形に対し、重み付き(○印数)三角形を
1@
1@ 1A
1@ 3A 2B (← 先頭の「1」は、上の段の1×@、真ん中の「3」は、上の段の
1×@+1×A、右の「2」は、上の段の1×Aから)
1@ 7A 12B 6C
1@ 15A 50B 60C 24D
1@ 31A 180B 390C 360D 120E
・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・
と構成していく。
(上の「180」は、上の段の15*A+50*B=30+150から計算、「390」は、50*B+60*C=150+240
から計算)
この上から6段目に並んだ数字を用いてこれを二項係数nCkに貼り付けていく。
s[5]=15+25+35+・・・・・+n5
=1*nC1+31*nC2+180*nC3+390*nC4+360*nC5+120*nC6
=n+31*n(n-1)/2+180*n(n-1)(n-2)/6+390*n(n-1)(n-2)(n-3)/24
+360*n(n-1)(n-2)(n-3)(n-4)/120+120*n(n-1)(n-2)(n-3)(n-4)(n-5)/720
=n6/6+n5/2+5n4/12-n2/12
これなら、s[19]、s[20]でも張り紙無しで対応できそうだ。この重み付き三角形の計算が、
Mathematicaコマンドで、
Table[Table[Sum[(-1)^(i+k)*(i+1)^n*Binomial[k,i],{i,0,k}],{k,0,n}],{n,0,10}]//ColumnForm
で構成できました。{n,0,10}なら11段までですので21段必要なら{n,0,20}で実行。
更に面白かったのが、この三角形の数字と重みを組み合わせると、ベルヌーイ数Bnが続
々と出現するところ。
B0=1/@=1
B1=1/@-1/A=1/2 (B1=-1/2を第1ベルヌーイ数、B1=1/2を第2ベルヌーイ数というらしい。)
B2=1/@-3/A+2/B=1/6
B3=1/@-7/A+12/B-6/C=0
B4=1/@-15/A+50/B-60/C+24/D=-1/30
B5=1/@-31/A+180/B-390/C+360/D-120/E=0
B6=1/@-63/A+602/B-2100/C+3360/D-2520/E+720/F=1/42
・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・・
自然数の冪乗和にはベルヌーイ数がセットで出現していましたが、こんな風に隠れていた
とは思いませんでした。先人たちの仕事に感謝します。これで万が一の場合に備えるための
トイレの張り紙はしなくて済みました。
(追記) 「自然数の冪乗和計算からの進化」と題して、GAI さんからの投稿です。
(平成27年10月17日付け)
自然数の和の公式
Σk=1〜n k=n(n+1)/2
Σk=1〜n k2=n(n+1)(2n+1)/6
Σk=1〜n k3=n2(n+1)2/4
Σk=1〜n k4=n(n+1)(2n+1)(3n2+3n-1)/30
Σk=1〜n k5=n2(n+1)2(2n2+2n-1)/12
・・・・・・・・・・・・・・・・・・・・・・・
などの式を用いることで、
Σk=1〜10 k=55
Σk=1〜10 k2=385
Σk=1〜10 k3=3025
Σk=1〜10 k4=25333
Σk=1〜10 k5=220825
・・・・・・・・・・・・・・・・・・・・・・
が計算できる。
ここで、 f(x)=(x-1)(x-2)(x-3)(x-4)(x-5)(x-6)(x-7)(x-8)(x-9)(x-10) とおいて、x6・f’(x) を
f(x) で割った商を計算すると(なお、f’(x)は導関数)、
10x5+55x4+385x3+3025x2+25333x+220825
と、上記の一連の和の値が商の各係数として出現する。この原理を利用すると、
素数 2、3、5、7、11、13、17、19 に対する冪乗和でも
f(x)=(x-2)(x-3)(x-5)(x-7)(x-11)(x-13)(x-17)(x-19)
=x8-77x7+2451x6-41817x5+414849x4-2429223x3+8130689x2-14117683x+9699690
を利用することにより、x11・f’(x)をf(x)で割った商を計算することで、
8x10+77x9+1027x8+15803x7+260167x6+4448507x5+77915887x4
+1387350083x3+24995572327x2+454280348267x+8311148375647
が起こるので、
2+3+5+・・・+17+19=77
2^2+3^2+5^2+・・・+17^2+19^2=1027
2^3+3^3+5^3+・・・+17^3+19^3=15803
2^4+3^4+5^4+・・・+17^4+19^4=260167
2^5+3^5+5^5+・・・+17^5+19^5=4448507
・・・・・・・・・・・・・・・・・・・・・・・・・・
2^10+3^10+5^10+・・・+17^10+19^10=8311148375647
が一気に求められる。
(コメント) 面白い計算ですね!
f(x)=(x-1)(x-2)(x-3)(x-4)=x4-10x3+35x2-50x+24 において、 x2・f’(x) を f(x) で割った商
を計算すると、 4x+10
この係数を見て、 Σk=1〜4 1=4 、Σk=1〜4 k=10 とするのかな?
f(x)=(x-1)(x-2)(x-3)(x-4)(x-5)(x-6)=x6-21x5+175x4-735x3+1624x2-1236x+720 において、
x4・f’(x)=6x9-105x8+700x7-2205x6+3248x5-1236x4 を f(x) で割った商を計算すると、
6x3+21x2+91x+441
この係数を見て、
Σk=1〜6 1=6 、Σk=1〜6 k=21 、Σk=1〜6 k2=91 、Σk=1〜6 k3=441
とするのかな?
(追記) 「組合せだけの冪乗和への挑戦」と題して、当HPがいつもお世話になっているHN
「GAI」さんからの投稿です。(平成28年8月7日付け)
組合せ数を計算する道具として、nCr(=binomial(n,r)) なる記号がある。
よく n を固定して r を変化させて、
n-1Cr-1+n-1Cr=nCr 、nC0+nC1+nC2+・・・+nCn=2n
などの変形式をよく目にする。
ここでは、逆に r を固定して n を変化させることを考える。
例えば、r=5 で固定して、n を 0〜10 で変化させ和を取れば、
Σ n=0〜10 nC5= 1 + 6 + 21 + 56 + 126 + 252 = 462 = 11C6
また、r=3 で固定して、n を 0〜8 で変化させると、
Σ n=0〜8 nC3= 1 +4 +10 +20 +35 +56 = 126 =
9C4
一般に、 Σ k=0〜n kCr= n+1Cr+1 ・・・・・・・・・(*)
と、いくつかのkCrの和が一つのn+1Cr+1の値に対応してくれる。
そこで、一般に、k2 を kC2、kC1で作ってみると、 k2=2kC2+kC1 より
Σ k=0〜nk2=2Σ k=0〜n kC2+Σ k=0〜n kC1=2n+1C3+n+1C2
よって、 左辺=2(n+1)n(n-1)/6+(n+1)n/2=n(n+1)(2n+1)/6
同様に、 k3=6kC3+6kC2+kC1 より、
Σ k=0〜nk3=6Σ k=0〜n kC3+6Σ k=0〜n kC2+Σ k=0〜n kC1
=6n+1C4+6n+1C3+n+1C2=n2(n+1)2/4
以下、
k4=24kC4+36kC3+14kC2+kC1
k5=120kC5+240kC4+150kC3+30kC2+kC1
k6=720kC6+1800kC5+1560kC4+540kC3+62kC2+kC1
k7=5040kC7+15120kC6+16800kC5+8400kC4+1806kC3+126kC2+kC1
・・・・・・・
(これらの係数は恒等式とみて,k=1,2,3,4,・・・を両辺に代入することですぐ手に入る。)
を利用することで、これから自然数の冪乗和を組合せだけの道具立てで構成していける。
(コメント) 通常は、恒等式 (k+1)3−k3=3k2+3k+1 を使うところを新しい恒等式
k3=6kC3+6kC2+kC1
を使うわけですね!係数については、k3=xkC3+ykC2+zkC1 において、
k=1 として、 1=z
k=2 として、 8=y+2z より、 y=6
k=3 として、 27=x+3y+3z より、 x=6
が得られ、k3=6kC3+6kC2+kC1 であることが分かる。
GAI さんからのコメントです。(平成28年11月8日付け)
ある人のホームページを見ていたら何気に、
Σ k=1〜nk5=n+1C6+26n+2C6+66n+3C6+26n+4C6+n+5C6
を片隅に掲示してあり、最初はデザインかなと思っていたんですが、計算機で確認してみる
と、何とすべての自然数nでバッチリ公式が成立しているではないか!
こんなにきれいに自然数の冪乗和と組合せ関数が融合できることに驚き、何もこれは5乗
に限る等式では無いはずだと他の冪乗で探ってみました。
すると、
Σ k=1〜nk4=n+1C5+11n+2C5+11n+3C5+n+4C5
Σ k=1〜nk3=n+1C4+4n+2C4+n+3C4
Σ k=1〜nk2=n+1C3+n+2C3
Σ k=1〜nk=n+1C2 (正に、=(n+1)*n/2))
が成立してきた。そこで、上昇へ向かって、
Σ k=1〜nkr (r=6,7,8,9,10,・・・)
での関係式を見つけてみて下さい。
Seiichi Manyama さんからのコメントです。(平成28年11月8日付け)
昼休み中なので、ちらっと見た程度ですが、オイラー多項式と関係があるのではないでしょ
うか?数列にすれば、「A008292」じゃないでしょうか?
DD++ さんからのコメントです。(平成28年11月8日付け)
n=5 までを下から順に構成してみます。まずは、当たり前の等式から出発。
x = x
2x = (x-1)+(x+1) を両辺にかけて、2x^2 = (x-1)x + x(x+1)
3x = (x-2)+2(x+1) = 2(x-1)+(x+2) を各項にかけて、
6x^3 = (x-1)x{(x-2)+2(x+1)} + x(x+1){2(x-1)+(x+2)} = (x-2)(x-1)x + 4(x-1)x(x+1)
+ x(x+1)(x+2)
4x = (x-3)+3(x+1) = 2(x-2)+2(x+2) = 3(x-1)+(x+3) を各項にかけて、
24x^4 = (x-2)(x-1)x{(x-3)+3(x+1)} + 4(x-1)x(x+1){2(x-2)+2(x+2)} + x(x+1)(x+2){3(x-1)+(x+3)}
= (x-3)(x-2)(x-1)x + 11(x-2)(x-1)x(x+1)+ 11(x-1)x(x+1)(x+2) + x(x+1)(x+2)(x+3)
5x = 略 を各項にかけて、
(5!)x^5 = (x-4)(x-3)(x-2)(x-1)x + 26(x-3)(x-2)(x-1)x(x+1) + 66(x-2)(x-1)x(x+1)(x+2)
+ 26(x-1)x(x+1)(x+2)(x+3) + x(x+1)(x+2)(x+3)(x+4)
さて、ここで、各項にさらに、
6 = (x+1)-(x-5) = (x+2)-(x-4) = (x+3)-(x-3) = (x+4)-(x-2) = (x+5)-(x-1) をかけると、
(6!)x^5
= (x-4)(x-3)(x-2)(x-1)x(x+1) - (x-5)(x-4)(x-3)(x-2)(x-1)x + 26(x-3)(x-2)(x-1)x(x+1)(x+2)
- 26(x-4)(x-3)(x-2)(x-1)x(x+1) + 66(x-2)(x-1)x(x+1)(x+2)(x+3) - 66(x-3)(x-2)(x-1)x(x+1)(x+2)
+ 26(x-1)x(x+1)(x+2)(x+3)(x+4) - 26(x-2)(x-1)x(x+1)(x+2)(x+3) + x(x+1)(x+2)(x+3)(x+4)(x+5)
- (x-1)x(x+1)(x+2)(x+3)(x+4)
これに、x=1,2,3,……,n を代入して合計し、6! で割ると、
Σ k=1〜nk5=n+1C6+26n+2C6+66n+3C6+26n+4C6+n+5C6
もっと高次の式を出したければ、6x=略 をかけ、7x=略 をかけ、……、とやってから最終段
階に入ればいいはずです。もっとも、最終段階をやる前に係数は全部見えていますが。
さらに面倒ならば、
A[i,1] = A[i,i] = 1, A[i,j] = (i-j+1)A[i-1,j-1] + jA[i-1,j]
という漸化式で係数だけ計算できます。
GAI さんからのコメントです。(平成28年11月9日付け)
1 1
2;2
1 4 1
3;2 2;3
1 11 11 1
4;2 3;3 2;4
1 26 66 26 1
5;2 4;3 3;4 2;5
1 57 302 302 57 1
・・・・・・・・・・・・・・・・・・・・
57=1*5+26*2 、302=26*4+66*3 などとやれる訳ですね。
この係数を、Eulerian numbers と呼ぶんですか。(Eluer数とはまた違う)ほう、1700年
代でこの関係式を見つけているとはオイラーは面白くてしかたなかったでしょうね。失明する
まで研究に没頭する気持ちも分からなくはないです。しかし、それでも13人もの子供を作る
なんて、爪の垢でも煎じたいです。ひょっとして、DD++さんも子沢山では?
Σk=1〜100 k5 =101C6+26・102C6+66・103C6+26・104C6+105C6=171708332500
でスイスイ計算できるんだ!例の複雑になる5乗での公式を憶えなくてもこれさえあれば十
分です。
DD++ さんからのコメントです。(平成28年11月9日付け)
1^5+2^5+3^5+・・・+100^5=98*99*100*101*102*103/6+5*(100*101/2)^2-4*100*101/2
=171708332500
の方が、この場合は早いかも。
(追記) 平成30年4月2日付け
数学セミナー’18 3月号(日本評論社)のNOTEで、長岡京市のクスコ氏(91歳)が立方
数の和の公式の簡単な証明を発表されている。とても簡便な方法で感動した。
Sn=1+2+3+・・・+n とおくと、 Sn−Sn-1=n で、さらに、
Sn+Sn-1=(1+(n−1))+(2+(n−2))+・・・+((n−1)+1)+n=n2
よって、 n3=(Sn−Sn-1)(Sn+Sn-1)=Sn2−Sn-12
したがって、 Σk=1〜n k3=Sn2=(Σk=1〜n k)2 が成り立つ。
(コメント) 教科書だと、(k+1)4−k4 の恒等式を用いて示されますが、煩雑さがなくエレ
ガントですね!
(追記) 令和2年7月21日付け
Σk(k+1)=n(n+1)(n+2)/3 を示すのに、次のように階差型に変形すると、証明が
容易である。
Σk(k+1)=(1/3)Σ(k(k+1)(k+2)−(k−1)k(k+1))=n(n+1)(n+2)/3
(追記) 令和4年5月27日付け
これまで、
1 から n までの自然数の平方の和 : 12+22+・・・+n2 =n(n+1)(2n+1)/6
1 から n までの自然数の立方の和 : 13+23+・・・+n3 =n2(n+1)2/4
がいろいろな方法で示されている。最近、面積の考えを用いた次のような証明があることを
知った。
まず、自然数の平方の和から考えよう。
上図のn個の正方形の面積の総和が、 12+22+・・・+n2 である。これを横長の長方
形の面積の総和と考える。面積の総和をSとすると、
S=(1+2+・・・+n)+(2+3+・・・+n)+(3+4+・・・+n)+・・・+n
=n(n+1)/2+(n−1)(n+2)/2+(n−2)(n+3)/2+・・・+(n−(n−1))(n+n)/2
=(n2+n)/2+(n2+n−2)/2+・・・+(n2+n2−n(n−1)−n(n−1))/2
=(n2+n)/2+(n2+n−2)/2+・・・+(n2+n−n(n−1))/2
=(n2+n−1・0)/2+(n2+n−2・1)/2+・・・+(n2+n−n(n−1))/2
=(n2+n−(12−1))/2+(n2+n−(22−2))/2+・・・+(n2+n−(n2−n))/2
=n(n2+n)/2−S/2+n(n+1)/4
よって、 (3/2)S=n(n2+n)/2+n(n+1)/4=n(n+1)(2n+1)/4 より、
S=n(n+1)(2n+1)/6
次に、自然数の立方の和を考えよう。
上図の水色部分の面積は、
((1+2+・・・+n)+(1+2+・・・+(n-1)))n=(n(n+1)/2+n(n-1)/2)n=n3
なので、正方形全体の面積は、 13+23+・・・+n3 である。
ところで、正方形の1辺の長さは、 1+2+・・・+n=n(n+1)/2 なので、
13+23+・・・+n3=(n(n+1)/2)2
が成り立つ。
(追記) 令和5年3月7日付け
令和5年度大学入学共通テスト追試験の数学U・Bで、自然数の平方数の和に関する問
題が出題された。これまでとは異なる趣で、ハッとさせられた。
第2問[2] 12+22+・・・+102 をある関数の定積分で表すことを考える。次の問いに
答えよ。
(1) すべての実数 t に対して、 ∫tt+1 F(x)dx=t2 となる2次関数F(x)を求めよ。
(2) 12+22+・・・+102 をある関数の定積分で表せ。
(解)(1) F(x)=ax2+bx+c (a≠0) とおくと、
∫tt+1 F(x)dx=(a/3)(3t2+3t+1)+(b/2)(2t+1)+c=t2 が恒等式なので、
a=1 、a+b=0 、a/3+b/2+c=0
これを解いて、 a=1 、b=−1 、c=1/6 より、 F(x)=x2−x+1/6
(2) 12+22+・・・+102
=∫12 F(x)dx+∫23 F(x)dx+・・・+∫1011 F(x)dx
=∫111 (x2−x+1/6)dx と表される。 (終)
(コメント) 計算すると、 ∫111 (x2−x+1/6)dx=385 である。
この問題設定がどのような発展をもたらすのか、興味あるところである。
(追記) 令和5年6月29日付け
いくつかの数列の和の問題を考えてみた。
問題 数列の和 1+11+111+・・・+111111111 を求めよ。
計算結果は、当然のごとく「123456789」になるわけであるが、もっともらしく計算して
みよう。
(解) 第n項は、 1+10+・・・+10n-1=(10n−1)/9 なので、
与式=Σk=19 (10k−1)/9=10(109−1)/81−1=(1010−91)/81 (終)
上記の問題に若干手を入れて、次の問題を設定した。
問題 数列の和 1+22+333+・・・+999999999 を求めよ。
(解) 与式=(1+2+・・・+9)+10(2+3+・・・+9)+102(3+4+・・・+9)+
103(4+5+・・・+9)+104(5+6+・・・+9)+105(6+7+・・・+9)+
106(7+8+9)+107(8+9)+108・9
=900000000+
170000000+
24000000+
3000000+
350000+
39000+
4200+
440+
45
=1097393685 (終)
(追記) 令和6年8月25日付け
数列の和 Sn=Σk=1n ak と Tn=Σk=1n ak2 について、次の東北大学 理系(1979)の
問題は面白い。
第2問 各項が正の等差数列 a1,a2,・・・,an,・・・ において、
Sn=Σk=1n ak 、Tn=Σk=1n ak2 とする。任意の自然数nに対して、
Sn≦n2+n−1、Tn≧(4n3−n)/3 が成り立つとき、この等差数列の初項と公差を求めよ。
(解) a1=S1≦1 で、a12 =T1≧1 すなわち、 a1≧1 から、 a1=1
よって、初項 a は、1 となる。交差をdとおく。
Sn=n(2+(n−1)d)/2≦n2+n−1 から、
(d−2)n2−dn+2≦0 すなわち、 ((d−2)n−2)(n−1)≦0
n≧1 より、 (d−2)n−2≦0 なので、 d≦2+2/n
n → ∞ のとき、 d≦2
また、 1+(1+d)2≧10 より、 1+d≧3 すなわち、 d≧2
よって、d=2 でなければならない。
逆に、a=1、d=2 のとき、
Sn=n2 で、明らかに、 Sn≦n2+n−1 が成り立つ。
Tn=Σk=1n (2k−1)2 =(2/3)n(n+1)(2n+1)−2n(n+1)+n=(4n3−n)/3
より、 Tn≧(4n3−n)/3 が成り立ち、条件を満たす。
以上から、求める等差数列の初項は1で、公差は2である。 (終)
kuiperbelt さんからのコメントです。(令和6年8月26日付け)
数列の和
Σ_{k=1〜n}k=n(n+1)/2=C(n+1,2)
Σ_{k=1〜n}k(k+1)=n(n+1)(n+2)/3=2C(n+2,3)
Σ_{k=1〜n}k(k+1)(k+2)=(1/4)Σ_{k=1〜n}k(k+1)(k+2)(k+3)-(k-1)k(k+1)(k+2)=n(n+1)(n+2)(n+3)/4=3!C(n+3,4)
Σ_{k=1〜n}k(k+1)(k+2)(k+3)(k+4)=(1/5)Σ_{k=1〜n}k(k+1)(k+2)(k+3)(k+4)-(k-1)k(k+1)(k+2)(k+3)
=n(n+1)(n+2)(n+3)(n+4)/5=4!C(n+4,5)
・・・・・・
Σ_{k=1〜n}k(k+1)(k+2)(k+3)(k+4)...(k+m-1)=(m-1)!C(n+m-1,m)
と、
C(n,k)=C(n-1,k)+C(n-1,k-1) より、C(n-1,k-1)=C(n,k)-C(n-1,k) から、
C(n,k)=C(n+1,k+1)-C(n,k+1)=C(n+2,k+2)-2C(n+1,k+2)+C(n,k+2)
=C(n+3,k+3)-3C(n+2,k+3)+3C(n+1,k+3)-C(n,k+3)
=Σ_{j=0〜m}(-1)^j*C(m,j)*C(n+m-j,k+m)
より、
Σ_{k=1〜n}k^2=Σ_{k=1〜n}k(k+1)-Σ_{k=1〜n}k
=2C(n+2,3)-C(n+1,2)=2C(n+2,3)-(C(n+2,3)-C(n+1,3))=C(n+2,3)+C(n+1,3)
Σ_{k=1〜n}k^3=Σ_{k=1〜n}k(k+1)(k+2)-3Σ_{k=1〜n}k(k+1)+Σ_{k=1〜n}k
=3!C(n+3,4)-3*2C(n+2,3)+C(n+1,2)
=6C(n+3,4)-6C(n+2,3)+C(n+1,2)
=6C(n+3,4)-6(C(n+3,4)-C(n+2,4))+C(n+3,4)-2C(n+2,4)+C(n+1,4)
=C(n+3,4)+4C(n+2,4)+C(n+1,4)
Σ_{k=1〜n}k^4
=Σ_{k=1〜n}k(k+1)(k+2)(k+3)-6Σ_{k=1〜n}k(k+1)(k+2)+7Σ_{k=1〜n}k(k+1)-Σ_{k=1〜n}k
=4!C(n+4,5)-6*3!C(n+3,4)+7*2C(n+2,3)-C(n+1,2)
=24C(n+4,5)-36C(n+3,4)+14C(n+2,3)-C(n+1,2)
=24C(n+4,5)-36(C(n+4,5)-C(n+3,5))+14(C(n+4,5)-2C(n+3,5)+C(n+2,5))-(C(n+4,5)-3C(n+3,5)+3C(n+2,5)-C(n+1,5))
=C(n+4,5)+11C(n+3,5)+11C(n+2,5)+C(n+1,5)
Σ_{k=1〜n}k^5
=Σ_{k=1〜n}k(k+1)(k+2)(k+3)(k+4)-10Σ_{k=1〜n}k(k+1)(k+2)(k+3)+25Σ_{k=1〜n}k(k+1)(k+2)
-15Σ_{k=1〜n}k(k+1)+Σ_{k=1〜n}k
=5!C(n+5,6)-10*4!C(n+4,5)+25*3!C(n+3,4)-15*2C(n+2,3)+C(n+1,2)
=120C(n+5,6)-240C(n+4,5)+150C(n+3,4)-30C(n+2,3)+C(n+1,2)
=120C(n+5,6)-240(C(n+5,6)-C(n+4,6))+150(C(n+5,6)-2C(n+4,6)+C(n+3,6))
-30(C(n+5,6)-3C(n+4,6)+3C(n+3,6)-C(n+2,6))+(C(n+5,6)-4C(n+4,6)+6C(n+3,6)-4C(n+2,6)+C(n+1,6))
=C(n+5,6)+26C(n+4,6)+66C(n+3,6)+26C(n+2,6)+C(n+1,6)
Σ_{k=1〜n}k^6
=Σ_{k=1〜n}k(k+1)(k+2)(k+3)(k+4)(k+5)-15Σ_{k=1〜n}k(k+1)(k+2)(k+3)(k+4)
+65Σ_{k=1〜n}k(k+1)(k+2)(k+3)-90*Σ_{k=1〜n}k(k+1)(k+2)
+31*Σ_{k=1〜n}k(k+1)-Σ_{k=1〜n}k
=6!C(n+6,7)-15*5!C(n+5,6)+65*4!C(n+4,5)-90*3!C(n+3,4)+31*2C(n+2,3)-C(n+1,2)
=720C(n+6,7)-1800C(n+5,6)+1560C(n+4,5)-540C(n+3,4)+62C(n+2,3)-C(n+1,2)
=C(n+6,7)+57C(n+5,7)+302C(n+4,7)+302C(n+3,7)+57C(n+2,7)+C(n+1,7)
Dengan kesaktian Indukmu さんからのコメントです。(令和6年8月27日付け)
比べてみたらとても興味深かったです。(→ 数研通信)
(コメント) 東北大学 文系(1979) で、次の問題が出題されている。kuiperbelt さんの投
稿を拝見して、この問題の重要性がヒシヒシと感じられました。
問題 Sn=Σk=1n k (n=1、2、・・・) とするとき、次の和を求めよ。
(1) Tn=Σk=1n Sk
(2) Un=Σk=1n Tk
(解) Sn=n(n+1)/2 から、
(1) Tn=(1/2)Σk=1n k(k+1)=(1/6)n(n+1)(n+2)
(2) Un=(1/6)Σk=1n k(k+1)(k+2)=(1/24)n(n+1)(n+2)(n+3) (終)
(追記) 令和6年10月26日付け
次の東北大学 文系(1983)の問題も興味深い。
問題 (1) F(t)は t の4次式で、任意の実数 x に対して ∫x-1x F(t)dt=x4 を満たす
という。F(t)を求めよ。
(2) (1)を利用して Σk=1n k4 を求めよ。
(解)(1) F(t)=at4+bt3+ct2+dt+e とおくと、条件より、
(a/5)(x5−(x−1)5)+(b/4)(x4−(x−1)4)+(c/3)(x3−(x−1)3)+(d/2)(x2−(x−1)2)+e
=(a/5)(5x4−10x3+10x2−5x+1)+(b/4)(4x3−6x2+4x−1)
+(c/3)(3x2−3x+1)+(d/2)(2x−1)+e
=ax4+(−2a+b)x3+(2a−3b/2+c)x2+(−a+b−c+d)x+a/5−b/4+c/3−d/2+e
=x4 が x の恒等式なので、
a=1 、−2a+b=0 より、b=2 、2a−3b/2+c=0 より、c=1
−a+b−c+d=0 より、d=0 、a/5−b/4+c/3−d/2+e=0 より、e=−1/30
以上から、 F(t)=t4+2t3+t2−1/30 である。
(2) Σk=1n k4=∫0n (t4+2t3+t2−1/30)dt=n5/5+n4/2+n3/3−n/30 (終)
(コメント) 因数分解まで頑張れば、
n5/5+n4/2+n3/3−n/30
=n(6n4+15n3+10n2−1)/30
=n(n+1)(6n3+9n2+n−1)/30
=n(n+1)(2n+1)(3n2+3n−1)/30
となりますね!でも、ここまで求めよ、とも指示されていないので...。
(追記) 令和6年10月29日付け
次の東北大学 文理共通(1984)の問題は、不等式の評価に関する問題である。
問題2 数列{an}において、an≧0 (n=1、2、・・・)とし、sn=Σk=1n an とおく。このと
き、次の不等式を数学的帰納法を用いて証明せよ。
(1+a1)(1+a2)・・・ (1+an)≦1+sn+sn2/2!+・・・+snn/n!
(解) n=1 のとき、左辺=1+a1 、右辺=1+s1=1+a1 なので、左辺≦右辺
よって、n=1 のとき成り立つ。
n=k(k≧1) のとき、成り立つと仮定する。すなわち、
(1+a1)(1+a2)・・・ (1+ak)≦1+sk+sk2/2!+・・・+skk/k!
このとき、
(1+a1)(1+a2)・・・ (1+ak) (1+ak+1)≦1+sk+1+sk+12/2!+・・・+sk+1k+1/(k+1)!
が成り立つことを示す。
1+sk+1+sk+12/2!+・・・+sk+1k+1/(k+1)! において、
p=1、2、・・・、k+1 に対して
sk+1p/p!=(sk+ak+1)p/p!≧(skp+pskp-1ak+1)/p!
なので、
1+sk+1+sk+12/2!+・・・+sk+1k+1/(k+1)!
≧(1+sk+sk2/2!+・・・+skk/k!)+ak+1(1+sk+sk2/2!+・・・+skk/k!)
≧(1+a1)(1+a2)・・・ (1+ak)+ak+1(1+a1)(1+a2)・・・ (1+ak)
=(1+a1)(1+a2)・・・ (1+ak) (1+ak+1)
よって、n=k+1 のときも成り立つ。
したがって、すべての自然数nに対して、不等式は成り立つ。 (終)
DD++ さんからのコメントです。(令和6年11月3日付け)
こんな解き方が思いついた場合、どうすればいいんでしょうね。
(別解) 任意の自然数 k に対し 1+a[k]>0 なので、相加平均と相乗平均の大小関係により、
{(1+a1)+(1+a2)+・・・+ (1+an)}/n≧{(1+a1)(1+a2)・・・ (1+ak)}^(1/n)
すなわち、 Π[k=1→n] (1+a[k]) ≦ (1+s[n]/n)^n = Σ[k=0→n] nCk/n^k*s[n]^k
ここで、nCk/n^k * k! = nPk/n^k ≦ 1 なので、nCk/n^k ≦ 1/k!
したがって、Π[k=1→n] (1+a[k]) ≦ Σ[k=0→n] s[n]^k/k! (終)
#数学的帰納法を使うよりも本質を捉えた証明ではないかと思いますが、帰納法を使ってな
いので 0 点になるんでしょうか?
証明法の探究の場合はさておき、入試で証明方法に制限をつけるのは、ナンセンスと感じ
ますね。
(コメント) 見事な証明で、感動しました!DD++ さんに感謝します。ただ、証明方法に制限
がついているので、減点は免れないかも...。
以下、工事中