베르누이 수·거듭제곱 합 계산기
베르누이 수 Bₙ을 유리수로 정확히 내고, 그 값으로 1^k + 2^k + … + n^k의 닫힌 식(파울하버 공식)을 만들어 줍니다. 만든 식을 실제 합과 대조하는 검산표도 함께 보여 줍니다.
1^3 + 2^3 + … + n^3
1/4·n^4 + 1/2·n^3 + 1/4·n^2
k = 3이므로 n에 대한 4차 다항식입니다. 최고차항 계수는 언제나 1/(k+1) = 1/4입니다. 아래 검산표에서 n = 1..12에 대해 실제로 하나씩 더한 값과 같은지 확인할 수 있습니다.
닫힌 식과 실제 합이 같은가
| n | 닫힌 식 | 하나씩 더한 값 | 일치 |
|---|---|---|---|
| 1 | 1 | 1 | ● |
| 2 | 9 | 9 | ● |
| 3 | 36 | 36 | ● |
| 4 | 100 | 100 | ● |
| 5 | 225 | 225 | ● |
| 6 | 441 | 441 | ● |
| 7 | 784 | 784 | ● |
| 8 | 1296 | 1296 | ● |
| 9 | 2025 | 2025 | ● |
| 10 | 3025 | 3025 | ● |
| 11 | 4356 | 4356 | ● |
| 12 | 6084 | 6084 | ● |
왼쪽은 위에서 만든 식에 n을 넣은 값, 오른쪽은 13부터 n3까지 실제로 하나씩 더한 값입니다. 유리수로 계산하므로 반올림 오차 때문에 어긋날 일이 없습니다.
베르누이 수 B₀ ~ B12
| n | Bₙ | 소수로 |
|---|---|---|
| 0 | 1 | 1.000000000 |
| 1 | -1/2 | -0.5000000000 |
| 2 | 1/6 | 0.1666666667 |
| 3 | 0 | 0 |
| 4 | -1/30 | -0.03333333333 |
| 5 | 0 | 0 |
| 6 | 1/42 | 0.02380952381 |
| 7 | 0 | 0 |
| 8 | -1/30 | -0.03333333333 |
| 9 | 0 | 0 |
| 10 | 5/66 | 0.07575757576 |
| 11 | 0 | 0 |
| 12 | -691/2730 | -0.2531135531 |
B₁을 뺀 홀수 번째는 모두 0입니다. 생성함수 x/(eˣ−1)에 x/2를 더하면 짝함수가 되고, 짝함수의 전개에는 홀수 차수 항이 없기 때문입니다.
계산 방법
- 1보고 싶은 베르누이 수의 번호를 넣으면 B₀부터 그 번호까지 유리수로 나옵니다.
- 2거듭제곱 차수 k를 넣으면 1^k + 2^k + … + n^k의 닫힌 식이 나옵니다.
- 3검산표에서 만든 식에 n을 넣은 값과 실제로 하나씩 더한 값이 같은지 확인합니다.
- 4B₁₂ = −691/2730처럼 갑자기 큰 소수가 나타나는 자리를 눈으로 확인해 봅니다.
자주 묻는 질문
1^k + 2^k + … + n^k을 n에 대한 닫힌 식으로 쓸 때 그 계수를 만들어 내는 유리수 수열입니다. k가 무엇이든 이 합은 n에 대한 (k+1)차 다항식이 되는데, 그 계수가 베르누이 수로 표현됩니다. B₀ = 1, B₁ = −1/2, B₂ = 1/6, B₄ = −1/30 식으로 이어집니다.
재귀식 Σ_{j=0}^{m} C(m+1, j)·B_j = 0에서 B_m을 풀면 됩니다. 정리하면 B_m = −1/(m+1) × Σ_{j<m} C(m+1, j)·B_j이고, B₀ = 1에서 출발해 차례로 채웁니다. 이때 반드시 유리수로 정확히 계산해야 합니다 — 부동소수점으로 하면 B₂₀ 근처에서 이미 값이 무너집니다.
둘 다 쓰이며 어느 쪽도 틀린 것이 아닙니다. 표에 실리는 「알려진 값」은 B₁ = −1/2이고 이 계산기의 표도 그쪽입니다. 다만 그 값을 그대로 파울하버 공식에 넣으면 Σ_{i=1}^{n}이 아니라 Σ_{i=0}^{n−1}이 나옵니다. 그래서 거듭제곱 합을 낼 때만 j = 1 항의 부호를 뒤집습니다. 나머지 Bₙ은 두 관례에서 완전히 같고 오직 B₁의 부호만 다릅니다.
B₁을 뺀 모든 홀수 번째가 0입니다. 생성함수 x/(eˣ−1)에 x/2를 더하면 짝함수가 되기 때문입니다. 짝함수의 테일러 전개에는 홀수 차수 항이 없으므로 그에 대응하는 베르누이 수도 사라집니다.
691은 분모가 아니라 분자에 있습니다 — B₁₂ = −691/2730입니다. 분모 2730 = 2·3·5·7·13은 폰 슈타우트–클라우젠 정리로 설명됩니다. B₂ₙ의 분모는 (p−1)이 2n을 나누는 모든 소수 p의 곱이고, 12를 나누는 1·2·4·6·12에 1을 더한 2·3·5·7·13이 정확히 그것입니다. 분자에 갑자기 나타나는 691은 리만 제타 함수와 모듈러 형식 이야기에서 되풀이해 등장하는 유명한 소수입니다.
검산표에서 n = 1부터 12까지 두 값을 나란히 보여 줍니다. 왼쪽은 닫힌 식에 n을 넣은 값이고 오른쪽은 1^k부터 n^k까지 실제로 하나씩 더한 값입니다. 두 값이 모두 같아야 하며, 유리수 연산으로 계산하므로 반올림 오차 때문에 어긋날 일이 없습니다.
전송되지 않습니다. 모든 계산은 브라우저 안에서 이뤄지고, 입력값은 이 기기에만 남습니다.
알아두면 좋은 점
- 검증은 두 갈래로 했습니다. B₀부터 B₁₂까지가 알려진 값(1, −1/2, 1/6, 0, −1/30, 0, 1/42, 0, −1/30, 0, 5/66, 0, −691/2730)과 정확히 같은지, 그리고 B₁₄ = 7/6·B₁₆ = −3617/510·B₃₀ = 8615841276005/14322까지 맞는지 고정했습니다.
- 폰 슈타우트–클라우젠 정리를 따로 검사합니다. B₂ₙ의 분모가 「(p−1)이 2n을 나누는 소수들의 곱」과 같은지 n = 1..15에서 소수를 실제로 훑어 대조합니다.
- 거듭제곱 합의 닫힌 식은 k = 0부터 20까지, n = 0부터 40까지 하나씩 더한 값과 정확히 같은지 전수로 확인했습니다. 유리수 연산이라 결과의 분모가 언제나 1이 되는 것도 함께 검사합니다.
- 검증에서 잡은 것이 있습니다. B₁ = −1/2를 그대로 파울하버 공식에 넣었더니 1+2+…+10 = 55가 아니라 45(= 0+1+…+9)가 나왔습니다. 두 관례가 정확히 이 자리에서 갈립니다 — B₁ = −1/2는 Σ_{i=0}^{n−1}을, +1/2는 Σ_{i=1}^{n}을 줍니다. 표는 −1/2로 보이되 합 공식에서만 부호를 뒤집도록 고쳤습니다.
- k = 3의 합이 (1+2+…+n)²과 같은 것, 최고차항 계수가 언제나 1/(k+1)인 것, 상수항이 없는 것(n = 0에서 합이 0이어야 하므로)을 따로 고정했습니다.
- 베르누이 수는 40번까지, 거듭제곱 차수는 20까지 다룹니다. 그보다 크면 분자·분모의 자릿수가 화면에 담기지 않습니다.
함께 보면 좋은 도구
마지막 검증: 2026년 9월 1일 · 결과는 참고용 추정치입니다.