베르누이 수

testwiki
둘러보기로 이동 검색으로 이동

틀:위키데이터 속성 추적 수론에서, 베르누이 수(Bernoulli數, 틀:Llang)는 거듭제곱수의 합이나 탄젠트멱급수 전개 따위의 다양한 공식에 등장하는 유리수 수열이다. 야코프 베르누이세키 다카카즈가 비슷한 시기에 독립적으로 발견하였다.

정의

베르누이 수열 Bn은 다음과 같은 생성함수로 정의할 수 있다.

texp(t)1=n=0Bntnn!

보다 일반적으로, 다음이 성립한다.

textet1=ntnn!k=0n(nk)Bnkxk

따라서, 베르누이 수열 B0,B1,B2,은 다음과 같다.

1, −1/2, 1/6, 0, −1/30, 0, 1/42, 0, −1/30, …

여기서 분자는 틀:OEIS이며, 분모는 틀:OEIS이다.

일부 저자들은 첫 베르누이 수를 1/2 대신 +1/2로 사용하기도 한다. 여기서는 이를 B~로 표기하자.

B~i=(1)iBi={+1/2i=1Bni1
t1exp(t)=n=0B~ntnn!

베르누이 다항식

베르누이 다항식. 점차 사인 또는 코사인에 수렴하는 것을 알 수 있다.

베르누이 수 Bn가 주어졌을 때, 이로부터 다음과 같은 베르누이 다항식(틀:Llang) Bn(x)[x]을 정의할 수 있다.

Bn(x)=k=0n(nk)Bnkxk

이는 아펠 다항식열을 이룬다.

그렇다면, 베르누이 수는 베르누이 다항식으로부터 다음과 같이 얻을 수 있다.

Bn=Bn(0)
B~n=Bn(1)

성질

B1=1/2을 제외하고, 홀수차 베르누이 수는 모두 0이다. 4의 배수차 베르누이 수는 음수이며, 4의 배수가 아닌 짝수차 베르누이 수는 양수이다.

sgnBn={1n=101n1,3(mod4)1n0(mod4)+1n2(mod4)

생성 함수

베르누이 다항식의 생성 함수는 다음과 같다.

n=0Bn(x)tnn!=texp(xt)exp(t)1

따라서, 베르누이 수의 생성 함수는 다음과 같다.

n=0Bntnn!=texpt1
n=0B~ntnn!=texptexpt1

대칭

베르누이 다항식은 다음 항등식들을 만족시킨다.

Bn(1x)=(1)nBn(x)n0
(1)nBn(x)=Bn(x)+nxn1

삼각 함수의 멱급수

베르누이 수는 탄젠트코탄젠트매클로린 급수에 등장한다.

tanx=n=1(1)n122n(22n1)B2nx2n1(2n)!(|x|<π/2)
cotx=n=0(1)n22nB2nx2n1(2n)!

마찬가지로, 쌍곡 탄젠트의 매클로린 급수는 다음과 같다.

tanhx=n=122n(22n1)B2n(2n)!x2n1(|x|<π/2)

푸리에 급수

베르누이 다항식의 푸리에 급수는 다음과 같다.

Bn(x)=n!(2πi)nk=0e2πikxkn=2n!k=1cos(2kπxnπ2)(2kπ)n

따라서, 베르누이 다항식은 다음과 같이 사인 또는 코사인으로 수렴하는 것을 알 수 있다.

(2π)nn!Bn(x)=2cos(2πxnπ2)+O(2n)

제타 함수와의 관계

xζ(1x)의 그래프

베르누이 수는 리만 제타 함수의 특수한 값이다.

B~n=()nBn=nζ(1n)=2(2π)nsin(π(1n)/2)n!ζ(n)

따라서, 리만 제타 함수는 베르누이 수의 복소수 n에 대한 해석적 연속으로 볼 수 있다.

n=0일 경우 리만 제타 함수는 을 갖지만, 이 경우 다음과 같이 극한을 취할 수 있다.

B~0=B0=1=limn0nζ(1n)
B~1=B1=2limn1sin(π(1n)/2)(2π)nn!ζ(n)=1/2

따라서, 스털링 공식에 따라 짝수차 베르누이 수의 절댓값에 대하여 다음과 같은 점근 공식이 성립한다.

|B2n|=2(2n)!(2π)2nζ(2n)4πn(nπe)2n(n1)

보다 일반적으로, 베르누이 다항식은 후르비츠 제타 함수 ζ(s,q)의 특수한 값이다.

Bn(x)=nζ(1n,x)
B~n=Bn(1)=nζ(1n,1)=nζ(1n)

따라서, 후르비츠 제타 함수는 베르누이 다항식의 복소수 n에 대한 해석적 연속으로 볼 수 있다.

스털링 수와의 관계

베르누이 다항식과 베르누이 수는 제2종 스털링 수하강 포흐하머 기호로 다음과 같이 나타낼 수 있다.

Bn+1(x)=Bn+1+k=0nn+1k+1{nk}xk+1_
Bn=k=0k(1)kk!k+1{nk}

이 경우 B1=1/2이 된다.

반대로, 하강 포흐하머 기호를 베르누이 다항식으로 다음과 같이 나타낼 수 있다.

xn+1_=k=0nn+1k+1[nk](Bk+1(x)Bk+1)

여기서 [nk]제1종 스털링 수이다.

거듭제곱수의 합

임의의 자연수 mn에 대하여, 처음 n개의 m제곱수들의 합은 다음과 같다.

k=1nkm=1m+1k=0m(1)k(m+1k)Bknm+1k=1m+1(Bm+1(n+1)Bm+1(0))

이를 베르누이 공식(Bernoulli公式, 틀:Llang) 또는 파울하버 공식(틀:Llang)이라고 한다. 예를 들어, m=1일 경우, 1부터 n까지의 자연수의 합인 삼각수 공식을 얻을 수 있다.

1+2+3++n=12(B0n22B1n1)=12(n2+n)

m=2,3일 경우, 제곱수·세제곱수의 합 공식을 얻을 수 있다.

12+22+32++n2=13(B0n33B1n2+3B2n1)=13(n3+32n2+12n)
13+23+33++n3=14(B0n44B1n3+6B2n24B3n1)=14(n4+2n3+n2)

베르누이 공식은 음계산법을 통하여 간단하게 적을 수 있다. 우선, 음변수 𝖻에 대하여 다음과 같은 선형 범함수를 정의하자.

L:𝖻nBn

그렇다면

L:(a+𝖻)nBn(a)

가 된다. 따라서,

k=1nkm=1m+1L((n+1+𝖻)m+1𝖻m+1)=L𝖻𝖻+n+1xmdx

이다.

알고리즘

베르누이 수를 계산하는 효율적인 알고리즘들이 알려져 있으며,[1][2][3][4] Bn를 계산하는 가장 빠른 알려진 알고리즘의 시간 복잡도는

O(n2log2+ϵn)

이다.[4]

수론적 성질

폰 슈타우트-클라우젠 정리(틀:Llang)에 따르면, 모든 양의 정수 n에 대하여, 만약 Bn0이라면 다음이 성립한다.

Bn+(p1)n1p

여기서 (p1)np1n의 약수가 되는 모든 소수 p에 대한 합이다. 즉, Bn0의 분모는 (p1)np이다.

소수 p에 대하여 다음 두 조건이 서로 동치이며, 이를 만족시키는 소수를 정규 소수이라고 한다.

  • pB2n2n{2,4,6,,p3}
  • ph(ζp)). 여기서 hKK유수이며, (ζp)원분체이다.

이를 정규 소수의 쿠머 조건(틀:Llang)이라고 한다.

낮은 차수의 베르누이 수는 다음과 같다. 틀:OEIS, 틀:OEIS

n Bn
0 1
1 −1/2
2 1/6
4 −1/30
6 1/42
8 −1/30
10 5/66
12 −691/2730
14 7/6
16 −3617/510
18 43867/798

낮은 차수의 베르누이 다항식은 다음과 같다.

B0(x)=1
B1(x)=x12
B2(x)=x2x+16
B3(x)=x332x2+12x
B4(x)=x42x3+x2130
B5(x)=x552x4+53x316x
B6(x)=x63x5+52x412x2+142

역사

세키의 《괄요산법》
베르누이의 《추측술》 97쪽 (현대적으로 재조판)
러브레이스의 베르누이 수 계산 프로그램. 이는 세계 최초의 컴퓨터 프로그램이다.

1631년에 요한 파울하버(틀:Llang, 1580~1635)는 거듭제곱수의 합을 계산하는 알고리즘을 출판하였다.[5] 파울하버의 알고리즘은 효율적이지만, 파울하버는 이를 베르누이 수를 통하여 나타내지 않았다. 이에 대하여 도널드 커누스는 다음과 같이 적었다. 틀:인용문2

세키 다카카즈(1642~1708)는 1712년 사후에 출판된 《괄요산법》(틀:Llang틀:Lang)[6]에 거듭제곱수의 합에 대한 일반 공식 및 베르누이 수를 제시하였다. 《괄요산법》에는 산가지로 표기된 파스칼 삼각형 밑에 다음과 같이 처음 12개의 베르누이 수가 수록돼 있다.

一級
二級 取㆓二分之一㆒爲㆑加
三級 取㆓六分之一㆒爲㆑加
四級
五級 取㆓三十分之一㆒爲㆑減
六級
七級 取㆓四十二分之一㆒爲㆑加
八級
九級 取㆓三十分之一㆒爲㆑減
十級
十一級 取㆓六十六分之五㆒爲㆑加
十二級

즉, 분수 k/nn틀:Langk로 표시돼 있고, 부호는 양수일 경우 틀:Lang, 음수일 경우 틀:Lang로 표시돼 있다. (㆑, ㆒, ㆓와 같은 부호는 간분의 가에리텐(틀:Llang)이다.)

세키와 거의 동시에, 야코프 베르누이는 1713년 사후에 출판된 저서 《추측술》(틀:Llang)[7] 2부 3장 97쪽에 거듭제곱수의 합의 일반 공식을 제시하였지만, 증명하지 않았다. 이 책에서 베르누이는 다음과 같이 적었다. 틀:인용문2

1830년대에 레온하르트 오일러콜린 매클로린오일러-매클로린 공식을 발견하면서 베르누이 수를 재발견하였다. 아브라암 드무아브르레온하르트 오일러는 "베르누이 수"라는 표현을 최초로 사용하였다.[8] 1834년에 카를 구스타프 야코프 야코비는 베르누이 공식을 엄밀하게 증명하였다.[9]

에이다 러브레이스는 1843년에 찰스 배비지해석기관에 대한 책[10]의 주석 G(틀:Llang)에 베르누이 수를 계산하는 알고리즘을 기술하였다. 이 주석 G는 세계 최초의 컴퓨터 프로그램으로 여겨진다.

폰 슈타우트-클라우젠 정리는 카를 게오르크 크리스티안 폰 슈타우트(틀:Llang)[11]와 토마스 클라우젠(틀:Llang)[12] 이 1840년에 독자적으로 발견하였다. 정규 소수의 쿠머 조건은 에른스트 쿠머가 1850년에 증명하였다.[13]

같이 보기

참고 문헌

틀:각주

외부 링크

틀:전거 통제