도구스학업·수학

베르누이 수·거듭제곱 합 계산기

베르누이 수 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에 대해 실제로 하나씩 더한 값과 같은지 확인할 수 있습니다.

B12-691/2730
소수로-0.2531135531
거듭제곱 합의 차수4차 다항식
검산 (n = 1..12)모두 일치

닫힌 식과 실제 합이 같은가

n닫힌 식하나씩 더한 값일치
111
299
33636
4100100
5225225
6441441
7784784
812961296
920252025
1030253025
1143564356
1260846084

왼쪽은 위에서 만든 식에 n을 넣은 값, 오른쪽은 13부터 n3까지 실제로 하나씩 더한 값입니다. 유리수로 계산하므로 반올림 오차 때문에 어긋날 일이 없습니다.

베르누이 수 B₀ ~ B12

nBₙ소수로
011.000000000
1-1/2-0.5000000000
21/60.1666666667
300
4-1/30-0.03333333333
500
61/420.02380952381
700
8-1/30-0.03333333333
900
105/660.07575757576
1100
12-691/2730-0.2531135531

B₁을 뺀 홀수 번째는 모두 0입니다. 생성함수 x/(eˣ−1)에 x/2를 더하면 짝함수가 되고, 짝함수의 전개에는 홀수 차수 항이 없기 때문입니다.

B₁₂ = −691/2730. 앞의 값들이 1/6, −1/30, 1/42처럼 얌전하다가 갑자기 691이라는 큰 소수가 분자에 나타납니다. 분모 2730 = 2·3·5·7·13은 폰 슈타우트–클라우젠 정리로 설명됩니다 — B₂ₙ의 분모는 (p−1)이 2n을 나누는 모든 소수 p의 곱이고, 12를 나누는 1·2·4·6·12에 1을 더한 것이 정확히 그 소수들입니다. 691은 리만 제타 함수와 모듈러 형식 이야기에서 되풀이해 등장하는 유명한 소수입니다.
B₁이 −1/2인가 +1/2인가는 관례가 갈리는 자리이고, 둘 다 틀린 것이 아닙니다. 위 표는 널리 실리는 B₁ = −1/2로 냈습니다. 다만 그 값을 그대로 파울하버 공식에 넣으면 Σi=1..n이 아니라 Σi=0..n−1이 나옵니다 — k = 1에서 n(n+1)/2 대신 n(n−1)/2가 됩니다. 그래서 이 계산기는 표는 −1/2로 보이되 합 공식에서만 그 항의 부호를 뒤집습니다. 나머지 Bₙ은 두 관례에서 완전히 같습니다.
베르누이 수는 반드시 유리수로 정확히 계산해야 합니다. 부동소수점으로 하면 B₂₀ 근처에서 이미 값이 무너집니다 — 분자·분모가 커지는 속도가 배정밀도의 자릿수를 금세 넘기 때문입니다. 이 계산기는 분자·분모를 BigInt로 들고 매 연산마다 약분합니다.

계산 방법

  1. 1보고 싶은 베르누이 수의 번호를 넣으면 B₀부터 그 번호까지 유리수로 나옵니다.
  2. 2거듭제곱 차수 k를 넣으면 1^k + 2^k + … + n^k의 닫힌 식이 나옵니다.
  3. 3검산표에서 만든 식에 n을 넣은 값과 실제로 하나씩 더한 값이 같은지 확인합니다.
  4. 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일 · 결과는 참고용 추정치입니다.