1. 개요
Faulhaber's formula자연수 [math(k)]에 대한 거듭제곱 [math(k^c)]의 합
| [math(\displaystyle \sum_{k=1}^n k^c)] |
2. 상세
일반식은 다음과 같이 주어진다.| [math(\begin{aligned} \sum_{k=1}^n k^c &= \sum_{k=0}^c \frac{(-1)^k}{c+1} \binom{c+1}k B_k n^{c+1-k} \\ &= \frac1{c+1}B_0n^{c+1} - B_1n^c + \frac c2B_2n^{c-1} - \frac{c (c-1)}6B_3n^{c-2} + \cdots\cdots + \frac{(-1)^{c-1}c}2B_{c-1}n^2 + (-1)^cB_cn \\ &= \frac1{c+1}n^{c+1} + \frac12n^c + \frac c{12}n^{c-1} + \cdots\cdots + \frac{(-1)^{c-1}c}2B_{c-1}n^2 + (-1)^cB_cn \end{aligned})] |
식 자체는 복잡해보이지만 [math(B_k = b^k)]로 치환하면 위의 식은 이항정리를 풀어서 쓴 형태와 유사하다는 걸 알 수 있다.
| [math(\displaystyle \sum_{k=0}^c \frac{(-1)^k}{c+1} \binom{c+1}k B_k n^{c+1-k} = \sum_{k=0}^c \frac 1{c+1} \binom{c+1}k (-b)^k n^{c+1-k} = \frac{(n-b)^{c+1} - ( -b)^{c+1}}{c+1})] |
| [math(\displaystyle B_c(n) = \sum_{k=0}^c \binom ck B_k n^{c-k})] |
| [math(\displaystyle \sum_{k=1}^nk^c = \frac{(-1)^{c+1}}{c+1}\{B_{c+1}(-n) - B_{c+1}(0)\})] |
파울하버는 베르누이가 공식을 발견하기 전에 [math(c)]가 홀수일 경우에 대한 규칙성을 발견하고 [math(c=17)]까지의 식을 제시한 인물로, 공식 자체를 증명한 사람은 아니지만 이와 관련이 있는 '파울하버 다항식'을 먼저 발견한 업적이 있어서인지 파울하버의 이름이 붙은 쪽이 더 유명하다.[3]
2.1. 간략화
베르누이 수열 문서에도 나와있듯이 [math(k \ne 1)]인 홀수이면 [math(B_k = 0)]이라는 성질이 있고 [math(k=1)]일 때의 항은 수식 구조상 반드시 [math(\dfrac12n^c)]이 나오기 때문에 [math(c)]가 홀수일 경우와 짝수일 경우에 따라 다음과 같이 간략화할 수 있다. 자연수 [math(m)]에 대해- [math(c = 2m-1)]이면 [math(B_{c+2} = 0)]이고 항은 [math(k=c-1=2m-2)]까지 존재하므로
[math(\displaystyle \sum_{k=1}^n k^{2m-1} = \frac12n^{2m-1} + \sum_{k=0}^{m-1} \frac1{2m} \binom{2m}{2k} B_{2k} n^{2(m-k)})] - [math(c = 2m-2)]이면 [math(B_c \ne 0)]이므로
[math(\displaystyle \sum_{k=1}^n k^{2m-2} = \frac{1-\delta_{1,\,m}}2 n^{2m-2} + \sum_{k=0}^m \frac1{2m-1} \binom{2m-1}{2k} B_{2k} n^{2(m-k)-1})]
3. 역사
당초 베르누이 자신의 표기법은 다음과 같았는데, 재미있는 것은 그가 죽고 난 뒤 출판된 《추측술》(Ars Conjectandi, 1713)이란 저서에 공식만 덩그러니 놓여있었을 뿐 증명이 같이 실려있지 않았다는 점이다. 엄밀한 증명은 후대에 야코비에 의해 이루어졌다.| [math(\displaystyle \sum_{k=1}^n k^c = \frac{n^{c+1}}{c+1} + \frac12n^c + \sum_{k=2}^c \frac{B_k}{k!} c^{\underline{k-1}} n^{c+1-k})] |
| [math(\displaystyle \sum_{k=1}^n k^c = \sum_{k=0}^c \frac{B_k}{k!} c^{\underline{k-1}} n^{c+1-k})] |
| [math(\begin{aligned} \sum_{k=1}^n k^c &= \sum_{k=0}^c \frac{B_k^+}{k!} c^{\underline{k-1}} n^{c+1-k} = \sum_{k=0}^c \frac{(-1)^k B_k}{k!} c^{\underline{k-1}} n^{c+1-k} = \sum_{k=0}^c (-1)^k B_k \frac{c!}{k! (c-k+1)!} n^{c+1-k} \\ &= \sum_{k=0}^c \frac{(-1)^k}{c+1} \binom{c+1}k B_k n^{c+1-k} \end{aligned})] |
4. 유도
| [math(\begin{aligned} \sum_{k=1}^n \left( e^x \right)^k &= \frac {e^x \left(e^{nx}-1 \right)}{e^x-1} = \frac{e^{nx}-1}{1 - e^{-x}} = \frac{e^{nx}-1}x \frac x{1 - e^{-x}} = \frac 1x \left\{ \sum_{c=0}^\infty \frac{(nx)^c}{c!} - 1 \right\} \sum_{k=0}^\infty \frac{B_k^+ x^k}{k!} = \sum_{c=1}^\infty \frac{n^c x^{c-1}}{c!} \sum_{k=0}^\infty \frac{B_k^+ x^k}{k!} \\ &= \sum_{c=0}^\infty \frac{n^{c+1} x^c}{(c+1)!} \sum_{k=0}^\infty \frac{B_k^+ x^k}{k!} = \sum_{c=0}^\infty \sum_{k=0}^c \frac{B_k^+ x^k}{k!} \frac{n^{c-k+1} x^{c-k}}{(c-k+1)!} = \sum_{c=0}^\infty \sum_{k=0}^c \frac{(c+1)!}{k! (c-k+1)!} \frac{B_k^+ n^{c+1-k} x^c}{(c+1)!} \\ &= \sum_{c=0}^\infty \sum_{k=0}^c \binom{c+1}k \frac{B_k^+ n^{c+1-k}}{c+1} \frac{x^c}{c!} = \sum_{c=0}^\infty \left\{\sum_{k=0}^c \frac{(-1)^k}{c+1} \binom{c+1}k B_k n^{c+1-k} \right\} \frac{x^c}{c!} \end{aligned})] |
| [math(\displaystyle \sum_{k=1}^n \left( e^x \right)^k = \sum_{k=1}^n e^{kx} = \sum_{k=1}^n \sum_{c=0}^\infty \frac{(kx)^c}{c!} = \sum_{c=0}^\infty \sum_{k=1}^n k^c \frac{x^c}{c!} = \sum_{c=0}^\infty \left(\sum_{k=1}^n k^c \right) \frac{x^c}{c!})] |
5. 예시
| [math(\begin{aligned} \sum_{k=1}^nk^0 &= \frac11B_0n^1 \\ &= n \\ \sum_{k=1}^n k^1 &= \frac12 \left( B_0 n^2 - 2 B_1 n \right) = \frac12 n^2 + \frac12n \\ &= \frac{n(n+1)}2 \\ \sum_{k=1}^n k^2 &= \frac13 \left( B_0 n^3 - 3 B_1 n^2 + 3 B_2n \right) = \frac13n^3 + \frac12n^2 + \frac16n = \frac{2n^3 + 3n^2 + n}6 \\ &= \frac{n(n+1)(2n+1)}6 = \frac{2n+1}3\frac{n(n+1)}2 \\ \sum_{k=1}^n k^3 &= \frac14 \left( B_0 n^4 - 4 B_1 n^3 + 6 B_2 n^2 \right) = \frac14n^4 + \frac12n^3 + \frac14n^2 \\ &= \left\{ \frac{n(n+1)}2 \right\}^2 \\ \sum_{k=1}^n k^4 &= \frac15 \left( B_0 n^5 - 5 B_1 n^4 + 10 B_2 n^3 + 5 B_4 n \right) = \frac15n^5 + \frac12n^4 + \frac13n^3 - \frac1{30}n = \frac{6n^5 + 15n^4 + 10n^3 - n}{30} \\ &= \frac{n(n+1)(2n+1) \left(3n^2+3n-1 \right)}{30} = \frac{2n+1}5{\left[2{\left\{\frac{n(n+1)}2\right\}}^2 - \frac13\frac{n(n+1)}2\right]}\\ \sum_{k=1}^nk^5 &= \frac16 \left( B_0n^6 - 6B_1n^5 + 15B_2n^4 + 15B_4n^2 \right) = \frac16n^6 + \frac12n^5 + \frac5{12}n^4 - \frac1{12}n^2 = \frac{2n^6 + 6n^5 + 5n^4 - n^2}{12} \\ &= {\left\{\frac{n(n+1)}2\right\}}^2\frac{2n^2+2n-1}3 = \frac43{\left[{\left\{\frac{n(n+1)}2\right\}}^3 -\frac14{\left\{\frac{n(n+1)}2\right\}}^2\right]} \\ \sum_{k=1}^nk^6 &= \frac17 \left( B_0n^7 - 7B_1n^6 + 21B_2n^5 + 35B_4n^3 + 7B_6n \right) = \frac17n^7 + \frac12n^6 + \frac12n^5 - \frac16n^3 + \frac1{42}n = \frac{6n^7 + 21n^6 + 21n^5 - 7n^3 + n}{42} \\ &= \frac{n(n+1)(2n+1)(3n^4 + 6n^3 - 3n + 1)}{42} = \frac{2n+1}7{\left[{4\left\{\frac{n(n+1)}2\right\}}^3-2{\left\{\frac{n(n+1)}2\right\}}^2+\frac13\frac{n(n+1)}2\right]} \\ \sum_{k=1}^nk^7 &= \frac18 \left( B_0n^8 - 8B_1n^7 + 28B_2n^6 + 70B_4n^4 + 28B_6n^2 \right) = \frac18n^8 + \frac12n^7 + \frac7{12}n^6 - \frac7{24}n^4 + \frac1{12}n^2 = \frac{3n^8 + 12n^7 + 14n^6 - 7n^4 + 2n^2}{24} \\ &= \frac{(n^4+2n^3+n^2)(3n^4+6n^3-n^2-4n+2)}{24} = 2{\left\{\frac{n(n+1)}2\right\}}^4 - \frac43{\left\{\frac{n(n+1)}2\right\}}^3 + \frac13{\left\{\frac{n(n+1)}2\right\}}^2 \end{aligned})] |
6. 대한민국 교육과정
대한민국 고등학교 교육과정에서 수열의 합을 다룰 때 [math(c\ge4)] 범위는 안 다룬다. [math(c = 1,\,2,\,3)]일 때 공식을 유도하는 원리조차 다르며[5] 애초에 이 공식을 외운다 하더라도 베르누이 수열을 같이 외워두지 않으면 아무런 쓸모가 없으므로 자연수 거듭제곱의 합에 대해 이런 공식이 있다는 것 정도로만 알아두면 좋을 것이다.7. 제타 함수와의 관계
거듭제곱 합의 공식에서 [math(c = -s)]를 대입하면| [math(\displaystyle \sum_{k=1}^n k^{-s} = \sum_{k=1}^n \frac1{k^s})] |
앞서 베르누이 다항식으로 표현했던 식을 곱씹어보자.
| [math(\begin{aligned} \sum_{k=1}^nk^c &= \frac{(-1)^{c+1}}{c+1}\{B_{c+1}(-n) - B_{c+1}(0)\} = \frac{(-1)^{c+1}}{c+1}B_{c+1}(-n) - \frac{(-1)^{c+1}}{c+1}B_{c+1}\\ &= \frac{(-1)^{c+1}}{c+1}\sum_{i=0}^{c+1}\binom{c+1}iB_i{\cdot}(-n)^{c+1-i} + \frac{(-1)^c}{c+1}B_{c+1} \\ &= \frac1{c+1}\sum_{i=0}^{c+1}\binom{c+1}i(-1)^iB_in^{c+1-i} + \frac{(-1)^c}{c+1}B_{c+1} \\ &= \frac1{c+1}\sum_{i=0}^{c+1}\binom{c+1}iB_i^+n^{c+1-i} + \frac{(-1)^c}{c+1}B_{c+1}\end{aligned})] |
| [math(\begin{aligned} \sum_{k=1}^nk^{-s} &= \sum_{k=1}^n\frac1{k^s} \\ &= \frac1{1-s}\sum_{i=0}^{1-s}\binom{1-s}iB_i^+n^{1-s-i} + \frac{(-1)^s}{1-s}B_{1-s} \\ &= \frac1{1-s}\left\{B_0n^{1-s} + \frac{(1-s)}{1!}B_1^+n^{-s} + \frac{(1-s)(-s)}{2!}B_2n^{-s-1} + \cdots\right\} + \frac{(-1)^s}{1-s}B_{1-s}\end{aligned})] |
| [math(\begin{aligned} \sum_{k=1}^n\frac1{k^s} &= \frac1{1-s}\sum_{i=0}^{1-s}\binom{1-s}i\frac{B_i^+}{n^{s+i-1}} + \frac{(-1)^s}{1-s}B_{1-s} \end{aligned})] |
| [math(\displaystyle \zeta(s) = \sum_{k=1}^\infty\frac1{k^s} = \frac{(-1)^s}{1-s}B_{1-s})] |
| [math(\zeta(-n) = \dfrac{(-1)^n}{n+1}B_{n+1})] |
| [math(B_n^+ = -n\zeta(1-n))] |
| [math(\zeta(-1) = -\dfrac12B_2 = -\dfrac12{\cdot}\dfrac16 = -\dfrac1{12})] |
[1] 후술하겠지만 증명하지는 않았다.[2] 사실 베르누이 수열을 발견한 야콥 베르누이 본인도 후자의 수열을 [math(B_k)]라 정의했다. 오늘날에는 생성함수를 통해 더 엄밀하게 정의할 수 있기 때문에 [math(B_1 = -\dfrac12)]이 되는 베르누이 수열을 [math(B_k)]라고 정의한다.[3] 그래서 파울하버 다항식에 등장하는 각 항의 계수는 베르누이 수열을 이용해서 표현할 수 있으나 당시 파울하버 본인은 이러한 연관성을 발견하지 못했다.[4] 고등학교 수학에서 '순열'이라는 이름으로 등장하는 연산이다.[5] 보통은 항등식 [math(\displaystyle (k+1)^c - k^c = \sum_{i=0}^{c-1} \binom cik^i = ck^{c-1} + \sum_{i=0}^{c-2}\binom cik^i)]을 [math(k=1)]부터 [math(k=n)]까지 모조리 더해서 [math(\displaystyle (n+1)^c -1 = c\sum_{k=1}^n k^{c-1} + \sum_{k=1}^n\left(\sum_{i=0}^{c-2} \binom cik^i\right))]로 만들고, 이항하면 [math(\displaystyle \sum_{k=1}^n k^{c-1} = \dfrac1c\left\{(n+1)^c -1 -\sum_{k=1}^n\left(\sum_{i=0}^{c-2} \binom cik^i\right)\right\})]이 되므로 우변의 [math(\displaystyle\sum_{k=1}^n\left(\sum_{i=0}^{c-2} \binom cik^i\right))]에 대해 이미 알고있는 합의 공식들을 모두 대입해서 유도하는 방식을 쓴다.그야말로 개노가다 [math(k^{c-1} = k^2)] 또는 [math(k^3)] 즉 [math(c = 3,\,4)]라면 이 방법으로 충분히 구할 수 있으나 [math(c\ge5)]가 되면 그야말로 개노가다의 산물이 된다. 한편, 본 문서에서 소개한 과정을 이해하려면 최소한 테일러 급수는 알고 있어야 하는데 이미 여기서부터 고등학교 범위를 한참 벗어났다.