베르누이 수열

역사 raw
대문 랜덤 문서 최근 토론

1. 개요2. 상세3. 역사4. 정의
4.1. 일반항
5. 성질6. 이용7. 오일러 수열과의 관계
7.1. 오일러 수열을 이용한 베르누이 수 표현7.2. 베르누이 수열을 이용한 오일러 수열 표현

1. 개요

Bernoulli Numbers

음이 아닌 정수 nn에 대해 BnB_n으로 나타내어지는 수열이다.

2. 상세

이 수열 자체에 대해서는 그다지 많이 알려지지 않았지만 tanx\tan x, cotx\cot x, tanhx\tanh x, cothx\coth x 등 다수의 삼각함수, 쌍곡선 함수의 테일러 급수에 B2nB_{2n}이라는 형태로 자리를 차지하고 있어 미친 존재감을 자랑한다. BnB_n이 아닌 B2nB_{2n}을 쓰는 이유는 B2n+1=0 (n1)B_{2n+1} = 0~ (n\ge1), 즉 3\bf3항 이상의 홀수항이 모조리 0\bf0이라는 독특한 성질이 있기 때문이다.[1]11항에 대해서도 간혹 B1=12B_1 = \dfrac12라 나타내는 문헌이 존재하는데 이것은 베르누이 수열이 (1)nBn(-1)^nB_n로 정의되었기 때문이다.[2] 혼동을 피하기 위해 일반적인 베르누이 수열을 BnB^-_n로, (1)n(-1)^n을 곱한 베르누이 수열을 Bn+B^+_n로 표기하기도 한다. 대략 제1818항까지의 값은 다음과 같다.
nn
00
11
22
33
44
55
66
77
88
99
1010
1111
1212
1313
1414
1515
1616
1717
1818
BnB_n
11
12-\dfrac12
16\dfrac16
00
130-\dfrac1{30}
00
142\dfrac1{42}
00
130-\dfrac1{30}
00
566\dfrac5{66}
00
6912730-\dfrac{691}{2730}
00
76\dfrac76
00
3617510-\dfrac{3617}{510}
00
43867798\dfrac{43867}{798}

3. 역사

역사적으로는 다음과 같은 거듭제곱 합의 계수에 대한 연구에서 시작되었다.
k=1nk=12n2+12nk=1nk2=13n3+12n2+16nk=1nk3=14n4+12n3+14n2k=1nk4=15n5+12n4+13n3130n\displaystyle\begin{aligned} \sum_{k=1}^n k &= \frac12n^2 + \frac12n \\ \sum_{k=1}^n k^2 &= \frac13n^3 + \frac12n^2 + \frac16n \\ \sum_{k=1}^n k^3 &= \frac14n^4 + \frac12n^3 + \frac14n^2 \\ \sum_{k=1}^n k^4 &= \frac15n^5 + \frac12n^4 + \frac13n^3 - \frac1{30}n \end{aligned}
훗날 야코프 베르누이nn의 거듭제곱에 붙은 계수들에 대해 일반항을 제시하기 전까지, 이 공식에 대해 열심히 연구하던 당대 수학자들[3] 중, 요한 파울하버(Johann Faulhaber)가 무려 k17k^{17}에 대한 합의 공식까지 제시하여 빼어난 기록을 남겼기에 오늘날에도 이 거듭제곱 합의 공식은 파울하버의 공식으로 알려져 있으나, 일반식을 제시한 건 베르누이이기 때문에 종종 베르누이의 공식, 또는 단순히 거듭제곱 합의 공식이라고도 불린다. 자세한 것은 해당 문서 참조.

이와는 별개로 일본의 세키 다카카즈가 그의 저서 《괄요산법》(括要算法, 1712)에서 n=12n=12까지에 대해 구체적인 값을 제시하였으나 일반식을 제시한 건 아니기에 수열 이름에 포함될 정도의 업적으로 보지는 않는 듯 하다.[4]

4. 정의

다음 생성함수를 이용하여 정의된다.
xex1=x2(cothx21)=n=0Bnn!xn\displaystyle\frac x{e^x -1} = \frac x2 \left( {\rm coth}\,\frac x2 - 1 \right) = \sum_{n=0}^\infty \frac{B_n}{n!}x^n
B1=12B_1 = \dfrac12Bn+B^+_n의 경우, 위의 테일러 급수에서 xx만큼을 더한 급수의 계수에 해당하므로
xex1+x=xexex1=x1ex=x2(cothx2+1)=n=0Bn+n!xn\displaystyle\frac x{e^x -1} + x = \frac{xe^x}{e^x - 1} = \frac x{1 - e^{-x}} = \frac x2 \left( \coth\,\frac x2 + 1 \right) = \sum_{n=0}^\infty \frac{B^+_n}{n!}x^n
으로 정의된다.

물론 위 식들을 직접 미분하고 x=0x=0을 대입하는 미친짓(……)으로 값을 계산하진 않고, 각 식의 역수들이 테일러 급수식으로 용이하게 나타낼 수 있다는 점을 이용, 점화식을 유도하여 계산하는 것이 일반적이다. 예를 들어 BnB_n의 경우 양변에
ex1x=1x(n=0xnn!1)=1xn=1(xnn!)=n=1(xn1n!)=n=0{xn(n+1)!}=1+x2!+x23!+x34!+\displaystyle\begin{aligned} \frac{e^x -1}x &= \frac1x \left(\sum_{n=0}^\infty \frac{x^n}{n!} -1 \right) = \frac1x \sum_{n=1}^\infty \left( \frac{x^n}{n!} \right) = \sum_{n=1}^\infty \left( \frac{x^{n-1}}{n!} \right) = \sum_{n=0}^\infty \left\{ \frac{x^n}{(n+1)!} \right\} \\ &= 1 + \frac x{2!} + \frac{x^2}{3!} + \frac{x^3}{4!} + \cdots\cdots \end{aligned}
를 곱하면 우변이 11이 되므로 좌변의 급수식을 적절하게 변형해주면 점화식이 얻어진다.
xex1ex1x=(B0+B11!x+B22!x2+B33!x3+)(1+x2!+x23!+x34!+)=n=0r=0nBrxrr!xnr(nr+1)!=n=0r=0nBrxnr!(nr+1)!=n=0r=0n(n+1)!r!(nr+1)!Brxn(n+1)!=n=01(n+1)!r=0n(n+1r)Brxn=B0+12!r=01(2r)Brx+13!r=02(3r)Brx2+14!r=03(4r)Brx3+=1\displaystyle\begin{aligned} \frac x{e^x -1}\frac{e^x -1}x &= \left( {B_0} + \frac{B_1}{1!}x + \frac{B_2}{2!}x^2 + \frac{B_3}{3!}x^3 + \cdots\cdots \right) \left( 1 + \frac x{2!} + \frac{x^2}{3!} + \frac{x^3}{4!} + \cdots\cdots \right) \\ &= \sum_{n=0}^\infty \sum_{r=0}^n \frac{B_rx^r}{r!} \frac{x^{n-r}}{(n-r+1)!} = \sum_{n=0}^\infty \sum_{r=0}^n \frac{B_r x^n}{r!(n-r+1)!} = \sum_{n=0}^\infty \sum_{r=0}^n \frac{(n+1)!}{r!(n-r+1)!} \frac{B_r x^n}{(n+1)!} \\ &= \sum_{n=0}^\infty \frac1{(n+1)!} \sum_{r=0}^n \binom{n+1}r B_rx^n \\ &= B_0 + \frac1{2!} \sum_{r=0}^1 \binom2rB_rx + \frac1{3!} \sum_{r=0}^2 \binom3rB_rx^2 + \frac1{4!} \sum_{r=0}^3 \binom4rB_rx^3 + \cdots\cdots \\ &= 1 \end{aligned}

항등식이므로 r=0n(n+1r)Br=δ0,n\displaystyle\sum_{r=0}^n \binom{n+1}rB_r = \delta_{0,\,n}이며(단, δ0,n\delta_{0,\,n}크로네커 델타) 이 식으로부터 점화식이 얻어진다.
r=0n(n+1r)Br=(n+1)Bn+r=0n1(n+1r)Br=δ0,nBn=δ0,n1n+1r=0n1(n+1r)Br\displaystyle\sum_{r=0}^n \binom{n+1}rB_r = (n+1)B_n + \sum_{r=0}^{n-1} \binom{n+1}rB_r = \delta_{0,\,n} \\ \therefore B_n = \delta_{0,\,n} - \frac1{n+1} \sum_{r=0}^{n-1} \binom{n+1}rB_r
제1항이 δ0,nn+1\dfrac{\delta_{0,\,n}}{n+1}이 아닌 이유는, n=0n=0이면 11이고 n1n\ge1이면 00이므로 사실상 값이 같기 때문이다. 보통은 n1n\ge1이라는 조건을 붙이지만 공합(empty sum)[5]00으로 약속하는 일반적인 정의에 따르면 위 식은 음이 아닌 정수에 대해 성립한다.

한편, cothx\coth x는 정의에 따라 다음과 같이 나타내어지는데, 바로 위의 생성함수를 이용하여 표현할 수 있다.
cothx=coshxsinhx=ex+ex2exex2=ex+exexex=e2x+1e2x1=1+2e2x1=1+1x2xe2x1=1+1xn=0Bnn!(2x)n=1+n=02nBnn!xn1\displaystyle\begin{aligned} \coth x &= \frac{\cosh x}{\sinh x} = \frac{\dfrac{e^x + e^{-x}}2}{\dfrac{e^x - e^{-x}}2} = \frac{e^x + e^{-x}}{e^x - e^{-x}} = \frac{e^{2x} + 1}{e^{2x} - 1} = 1 + \frac2{e^{2x} - 1} = 1 + \frac1x \frac{2x}{e^{2x} - 1} \\ &= 1 + \frac1x \sum_{n=0}^\infty \frac{B_n}{n!}(2x)^n \\ &= 1 + \sum_{n=0}^\infty \frac{2^n B_n}{n!}x^{n-1} \end{aligned}
쌍곡선 함수를 복소평면으로 확장시키면 coshix=cosx\cosh ix = \cos x, sinhix=isinx\sinh ix = i \sin x의 관계가 있음을 알 수 있고, 이로부터 cothix=icotx\coth ix = -i\cot x임을 알 수 있으므로 위의 테일러 전개식에 ixix를 대입하면
cothix=1+n=02nBnn!(ix)n1=1+2n=0(2i)n1Bnn!xn1=1+2{n=0(2i)2n1B2n(2n)!x2n1+n=0(2i)2nB2n+1(2n+1)!x2n}={1+2n=0(4)nB2n+1(2n+1)!x2n}i{n=0(4)nB2n(2n)!x2n1}=icotx\displaystyle\begin{aligned} \coth ix &= 1 + \sum_{n=0}^\infty \frac{2^nB_n}{n!}(ix)^{n-1} = 1 + 2\sum_{n=0}^\infty \frac{(2i)^{n-1}B_n}{n!}x^{n-1} \\ &= 1 + 2 \left\{ \sum_{n=0}^\infty \frac{(2i)^{2n-1}B_{2n}}{(2n)!}x^{2n-1} + \sum_{n=0}^\infty \frac{(2i)^{2n}B_{2n+1}}{(2n+1)!}x^{2n} \right\} \\ &= \left\{ 1 + 2\sum_{n=0}^\infty \frac{(-4)^nB_{2n+1}}{(2n+1)!}x^{2n} \right\} - i\left\{ \sum_{n=0}^\infty \frac{(-4)^nB_{2n}}{(2n)!}x^{2n-1} \right\} \\ &= -i\cot x \end{aligned}
실수부가 00이 되어야하며 허수부의 급수가 곧 cotx\cot x의 테일러 급수가 된다. 위 식의 실수부에 대해
1+2n=0(4)nB2n+1(2n+1)!x2n=1+2B1+n=1(4)nB2n+1(2n+1)!x2n=0\displaystyle1 + 2\sum_{n=0}^\infty \frac{(-4)^nB_{2n+1}}{(2n+1)!}x^{2n} = 1 + 2B_1 + \sum_{n=1}^\infty \frac{(-4)^nB_{2n+1}}{(2n+1)!}x^{2n} = 0
이므로 B2n+1=12δ0,nB_{2n+1} = -\dfrac12\delta_{0,\,n}이 얻어지며 이 식으로부터 33 이상의 홀수항은 00이 된다는 것을 알 수 있다.
이 사실을 이용하면 전술했던 베르누이 수열의 점화식도 다음과 같이 축약시킬 수 있게 된다.
B2n=δ0,n12n+1r=02n1(2n+1r)Br=δ0,n+12(1δ0,n)12n+1r=0n1(2n+12r)B2r=1+δ0,n212n+1r=0n1(2n+12r)B2r\displaystyle\begin{aligned} B_{2n} &= \delta_{0,\,n} - \frac1{2n+1} \sum_{r=0}^{2n-1} \binom{2n+1}rB_r = \delta_{0,\,n} + \frac12(1 - \delta_{0,\,n}) - \frac1{2n+1} \sum_{r=0}^{n-1} \binom{2n+1}{2r}B_{2r} \\ &= \frac{1 + \delta_{0,\,n}}2 - \frac1{2n+1} \sum_{r=0}^{n-1} \binom{2n+1}{2r}B_{2r} \end{aligned}
12(1δ0,n)\dfrac12(1-\delta_{0,\,n})은 합의 기호 부분에서 r=1r=1일 때, 즉 B1B_1이 곱해진 항을 계산하여 빼낸 부분인데, n=0n=0이면 r=1r=1인 항이 존재하지 않으므로 해당 항이 00이 되면서 n1n\ge1이면 12\dfrac12로 남아있도록 변형한 것이다. 이를 정리하면
Bn{B2n=1+δ0,n212n+1r=0n1(2n+12r)B2rB2n+1=12δ0,n\displaystyle\therefore B_n \begin{cases} \begin{aligned} B_{2n} &= \frac{1 + \delta_{0,\,n}}2 - \frac1{2n+1} \sum_{r=0}^{n-1} \binom{2n+1}{2r}B_{2r} \\ B_{2n+1} &= -\frac12\delta_{0,\,n} \end{aligned} \end{cases}

4.1. 일반항

Bn=k=0n1k+1r=0k(kr)(1)rrn\displaystyle B_n = \sum_{k=0}^n \frac1{k+1} \sum_{r=0}^k \binom kr(-1)^rr^n
식이 복잡해 보이지만 두 번째 합의 식은 제2종 스털링 수의 일반항에서 유래했다. 즉, 제2종 스털링 수 표기를 이용해서 나타내면 다음과 같다.
Bn=k=0nk!(1)kk+1S(n,k)\displaystyle B_n = \sum_{k=0}^n \frac{k!(-1)^k}{k+1} S(n,\,k)
베르누이 수열의 일반항은 조금 특이한 과정을 거쳐서 구해진다. 우선 각 생성함수를 다음과 같은 조건 하에 치환을 거쳐 식을 변형해준다.
i){\rm i}) x1ex=n=0Bn+xnn!=n=0(1)nBnxnn!\displaystyle\frac x{1 - e^{-x}} = \sum_{n=0}^\infty B_n^+\frac{x^n}{n!} = \sum_{n=0}^\infty (-1)^nB_n\frac{x^n}{n!}에서 1ex=t1 - e^{-x} = t로 치환하면 x=ln(1t)x = -\ln(1-t)가 되는데 x>0x>0일 때, 0<t<10<t<1이므로 해당 식에 대해 매클로린 급수를 적용할 수 있다. 따라서
x1ex=ln(1t)t=1tdt1t=1tk=0tkdt=1tk=0tk+1k+1=k=0tkk+1=k=0(1ex)kk+1=k=0(ex1)kk!k!(1)kk+1\displaystyle\begin{aligned} \frac x{1 - e^{-x}} &= \frac{-\ln(1 - t)}t = \frac1t \int \frac{{\rm d}t}{1-t} = \frac 1t \int \sum_{k=0}^\infty t^k\,{\rm d}t = \frac1t \sum_{k=0}^\infty \frac{t^{k+1}}{k+1} = \sum_{k=0}^\infty \frac{t^k}{k+1} \\ &= \sum_{k=0}^\infty \frac{(1 - e^{-x})^k}{k+1} = \sum_{k=0}^\infty \frac{(e^{-x} - 1)^k}{k!} \frac{k!(-1)^k}{k+1} \end{aligned}
(ex1)kk!\dfrac{(e^{-x} - 1)^k}{k!}제2종 스털링 수의 생성함수이므로 생성함수 식으로 바꾼 뒤 일반항을 대입한다.
x1ex=n=0(1)nBnxnn!=k=0(ex1)kk!k!(1)kk+1=k=0n=0{nk}(x)nn!k!(1)kk+1=n=0k=0n{nk}(x)nn!k!(1)kk+1=n=0(1)n{k=0n{nk}k!(1)kk+1}xnn!=n=0(1)n{k=0n1k!r=0k(kr)(1)krrnk!(1)kk+1}xnn!=n=0(1)n{k=0n1k+1r=0k(kr)(1)rrn}xnn!Bn=k=0n1k+1r=0k(kr)(1)rrn\displaystyle \begin{aligned} \frac x{1 - e^{-x}} &= {\color{blue}\sum_{n=0}^\infty (-1)^n}{\color{red}B_n}{\color{blue}\frac{x^n}{n!}} \\ &= \sum_{k=0}^\infty \frac{(e^{-x} - 1)^k}{k!} \frac{k!(-1)^k}{k+1} = \sum_{k=0}^\infty \sum_{n=0}^\infty \begin{Bmatrix} n \\ k \end{Bmatrix} \frac{(-x)^n}{n!} \frac{k!(-1)^k}{k+1} = \sum_{n=0}^\infty \sum_{k=0}^n \begin{Bmatrix} n \\ k \end{Bmatrix} \frac{(-x)^n}{n!} \frac{k!(-1)^k}{k+1} \\ &= \sum_{n=0}^\infty (-1)^n \left\{ \sum_{k=0}^n \begin{Bmatrix} n \\ k \end{Bmatrix} \frac{ k!(-1)^k}{k+1} \right\} \frac{x^n}{n!} \\ &= \sum_{n=0}^\infty (-1)^n \left\{ \sum_{k=0}^n \frac 1{\cancel{k!}} \sum_{r=0}^k \binom kr(-1)^{k-r} r^n \frac{\cancel{k!}(-1)^k}{k+1} \right\} \frac{x^n}{n!} \\ &= {\color{blue}\sum_{n=0}^\infty (-1)^n}{\color{red}\left\{ \sum_{k=0}^n \frac1{k+1} \sum_{r=0}^k \binom kr(-1)^rr^n \right\}}{\color{blue}\frac{x^n}{n!}} \end{aligned} \\ \therefore B_n = \sum_{k=0}^n \frac 1{k+1} \sum_{r=0}^k \binom kr(-1)^rr^n
ii){\rm ii}) xex1=n=0Bnxnn!\displaystyle\frac x{e^x - 1} = \sum_{n=0}^\infty B_n\frac{x^n}{n!}에서 ex1=te^x - 1 = t로 치환하면 x=ln(1+t)x = \ln(1+t)가 되는데 x<0x<0일 때, 1<t<0-1<t<0이므로 마찬가지로 매클로린 급수를 적용한다.
xex1=n=0Bnxnn!=ln(1+t)t=1tdt1+t=1tk=0(t)kdt=1tk=0(1)ktk+1k+1=k=0tk(1)kk+1=k=0(ex1)k(1)kk+1=k=0(ex1)kk!k!(1)kk+1=k=0n=0{nk}xnn!k!(1)kk+1=n=0k=0n{nk}xnn!k!(1)kk+1=n=0{k=0n{nk}k!(1)kk+1}xnn!=n=0{k=0n1k!r=0k(kr)(1)krrnk!(1)kk+1}xnn!=n=0{k=0n1k+1r=0k(kr)(1)rrn}xnn!Bn=k=0n1k+1r=0k(kr)(1)rrn\displaystyle\begin{aligned} \frac x{e^x - 1} &= {\color{blue}\sum_{n=0}^\infty}{\color{red}B_n}{\color{blue}\frac{x^n}{n!}} \\ &= \frac{\ln(1+t)}t = \frac1t \int \frac{{\rm d}t}{1+t} = \frac1t \int \sum_{k=0}^\infty (-t)^k\,{\rm d}t = \frac1t \sum_{k=0}^\infty \frac{(-1)^kt^{k+1}}{k+1} = \sum_{k=0}^\infty t^k\frac{(-1)^k}{k+1} \\ &= \sum_{k=0}^\infty (e^x - 1)^k\frac{(-1)^k}{k+1} = \sum_{k=0}^\infty \frac{(e^x - 1)^k}{k!} \frac{k!(-1)^k}{k+1} \\ &= \sum_{k=0}^\infty \sum_{n=0}^\infty \begin{Bmatrix} n \\ k \end{Bmatrix} \frac{x^n}{n!} \frac{k!(-1)^k}{k+1} = \sum_{n=0}^\infty \sum_{k=0}^n \begin{Bmatrix} n \\ k \end{Bmatrix} \frac{x^n}{n!} \frac{k!(-1)^k}{k+1} \\ &= \sum_{n=0}^\infty \left\{ \sum_{k=0}^n \begin{Bmatrix} n \\ k \end{Bmatrix} \frac{k!(-1)^k}{k+1} \right\} \frac{x^n}{n!} \\ &= \sum_{n=0}^\infty \left\{ \sum_{k=0}^n \frac1{\cancel{k!}} \sum_{r=0}^k \binom kr(-1)^{k-r}r^n \frac{\cancel{k!}(-1)^k}{k+1} \right\} \frac{x^n}{n!} \\ &= {\color{blue}\sum_{n=0}^\infty}{\color{red}\left\{\sum_{k=0}^n \frac1{k+1} \sum_{r=0}^k \binom kr(-1)^rr^n \right\}}{\color{blue}\frac{x^n}{n!}} \end{aligned} \\ \therefore B_n = \sum_{k=0}^n \frac1{k+1} \sum_{r=0}^k \binom kr(-1)^rr^n

Bn+B_n^+의 경우 생성함수x1ex=n=0Bn+xnn!\displaystyle \frac x{1 - e^{-x}} = \sum_{n=0}^\infty B^+_n \frac{x^n}{n!}에서 좌변의 식은 xex1ex\dfrac x{e^x -1}e^x와 같다. 즉 같은 방식으로 식을 전개해나가면 제2종 스털링 수의 생성함수 식이 ex(ex1)kk!\dfrac{e^x(e^x - 1)^k}{k!}로 주어지고 ex(ex1)kk!=n=0{n+1k+1}xnn!\displaystyle \frac{e^x(e^x - 1)^k}{k!} = \sum_{n=0}^\infty \begin{Bmatrix} n+1 \\ k+1 \end{Bmatrix} \frac{x^n}{n!}이므로
Bn+=k=0n{n+1k+1}k!(1)kk+1=k=0n1k+1r=0k(kr)(1)r(r+1)n\displaystyle\begin{aligned} B^+_n &= \sum_{k=0}^n \begin{Bmatrix} n+1 \\ k+1 \end{Bmatrix} \frac{k!(-1)^k}{k+1} \\ &=\sum_{k=0}^n \frac1{k+1} \sum_{r=0}^k \binom kr(-1)^r(r+1)^n \end{aligned}

5. 성질

  • k=0n(n+1k)Bk=δn,0\displaystyle \sum_{k=0}^n \binom{n+1}k B_k = \delta_{n,\,0}
상기 점화식을 구하는 과정에서 유도된 것이다. δn,0\delta_{n,\,0}크로네커 델타δn,m={1 (n=m)0 (nm)\delta_{n,\,m} = \begin{cases} 1~(n=m) \\ 0~(n \ne m) \end{cases}를 만족하는 함수이다.
  • k=0n(n+1k)Bk+=n+1\displaystyle \sum_{k=0}^n \binom{n+1}k B^+_k = n+1
파울하버의 공식 k=1nkc=n=0c(1)kc+1(c+1k)Bknn+1k\displaystyle \sum_{k=1}^n k^c = \sum_{n=0}^c \frac{(-1)^k}{c+1} \binom{c+1}kB_kn^{n+1-k}에서 n=1n=1을 대입하고 Bn+=(1)nBnB^+_n = (-1)^nB_n를 이용하면 된다.
두 식을 더하면 베르누이 수열의 짝수 항만 남고 좌변이 2배가 되므로
k=0n2(n+12k)B2k=n+1+δn,02\displaystyle\sum_{k=0}^{\left\lfloor{\frac n2}\right\rfloor} \binom{n+1}{2k}B_{2k} = \frac{n+1 + \delta_{n,\,0}}2
로 간략화할 수 있다. \lfloor \cdot \rfloor바닥 함수이다.

6. 이용

주로 테일러 급수에서 많이 쓰이고, 전술한대로 거듭제곱 합의 공식에도 쓰인다. 아래 목록에 없는 secx\sec xsechx{\rm sech}\, x오일러 수열을 이용해서 표현한다. 베르누이 수열이 오일러 수열과 서로 합연산[6] 관계에 있기는 하나(후술) 이걸 이용해서 두 급수를 표현하려면 식이 엄청 복잡해진다. 자세한 것은 항목 참조.
  • k=1nkc=k=0c(1)kc+1(c+1k)Bknc+1k\displaystyle\sum_{k=1}^n k^c = \sum_{k=0}^c \frac{(-1)^k}{c+1}\binom{c+1}kB_kn^{c+1-k}
파울하버의 공식이라고 한다. 식의 유도 과정은 해당 문서 참조.
  • cotx=n=0(4)nB2n(2n)!x2n1=1x13x145x32945x5\displaystyle\cot x = \sum_{n=0}^\infty \frac{(-4)^nB_{2n}}{(2n)!}x^{2n-1} = \frac1x - \frac13x - \frac1{45}x^3 - \frac2{945}x^5 - \cdots\cdots
  • tanx=n=1{(4)n(16)n}B2n(2n)!x2n1=x+13x3+215x5+17315x7+\displaystyle\tan x = \sum_{n=1}^\infty \frac{\{(-4)^n - (-16)^n\}B_{2n}}{(2n)!}x^{2n-1} = x + \frac13x^3 + \frac2{15}x^5 + \frac{17}{315}x^7 + \cdots\cdots
cotxtanx=cosxsinxsinxcosx=cos2xsin2xsinxcosx=cos2x12sin2x=2cot2x\cot x - \tan x = \dfrac{\cos x}{\sin x} - \dfrac{\sin x}{\cos x} = \dfrac{\cos^2x - \sin^2x}{\sin x\cos x} = \dfrac{\cos2x}{\dfrac12\sin2x} = 2\cot2x에서 tanx=cotx2cot2x\tan x = \cot x - 2\cot 2x라는 관계를 유도할 수 있어 위의 식이 자연스럽게 얻어진다.
  • cscx=n=0{2(1)n(4)n}B2n(2n)!x2n1=1x+16x+7360x3+3115120x5+\displaystyle\csc x = \sum_{n=0}^\infty \frac{\{2(-1)^n - (-4)^n\}B_{2n}}{(2n)!}x^{2n-1} = \frac1x + \frac16x + \frac7{360}x^3 + \frac{31}{15120}x^5 + \cdots\cdots
12(tanx+cotx)=12(cosxsinx+sinxcosx)=cos2x+sin2x2sinxcosx=1sin2x=csc2x\dfrac12(\tan x + \cot x) = \dfrac12\left(\dfrac{\cos x}{\sin x} + \dfrac{\sin x}{\cos x}\right) = \dfrac{\cos^2x + \sin^2x}{2\sin x\cos x} = \dfrac1{\sin2x} = \csc 2x에서 cscx=12(tanx2+cotx2)\csc x = \dfrac12\left(\tan\dfrac x2 + \cot\dfrac x2\right)를 이용하면 된다.
  • cothx=n=04nB2n(2n)!x2n1=1x+13x145x3+2945x5\displaystyle\coth x = \sum_{n=0}^\infty \frac{4^nB_{2n}}{(2n)!}x^{2n-1} = \frac1x + \frac13x - \frac1{45}x^3 + \frac2{945}x^5 - \cdots\cdots
위에서 유도한 식의 형태와 조금 다른데, 베르누이 수열에서 33 이상의 홀수 항이 00이 된다는 점을 적용해서 간략화시킨 형태이기 때문이다. cothx=icotix\coth x = i\cot ix를 이용해서도 유도할 수 있다.
  • tanhx=n=1(16n4n)B2n(2n)!x2n1=x13x3+215x517315x7+\displaystyle\tanh x = \sum_{n=1}^\infty \frac{(16^n - 4^n)B_{2n}}{(2n)!}x^{2n-1} = x - \frac13x^3 + \frac2{15}x^5 - \frac{17}{315}x^7 + \cdots\cdots
tanhx=itanix\tanh x = -i \tan ix를 이용해서 유도할 수 있다.
  • cschx=n=0(24n)B2n(2n)!x2n1=1x16x+7360x33115120x5+\displaystyle{\rm csch}\,x = \sum_{n=0}^\infty \frac{(2 - 4^n)B_{2n}}{(2n)!}x^{2n-1} = \frac1x - \frac16x + \frac7{360}x^3 - \frac{31}{15120}x^5 + \cdots\cdots
cschx=icscix{\rm csch}\,x = i \csc ix를 이용해서 유도할 수 있다.

7. 오일러 수열과의 관계

삼각함수 및 쌍곡선 함수가 각종 사칙연산을 통해 서로 연관되어있기 때문에, 베르누이 수열과 오일러 수열 역시 서로 무관하지는 않다. 다만, 아무래도 각 함수의 곱(즉, 테일러 급수끼리의 곱)이 반드시 포함되어 있기에 서로 합연산의 관계에 있어서 손계산이 그렇게 간단한 형태로 나오지는 않는다. 차라리 서로 점화식의 관계에 있다고 이해하는 편이 빠를 것이다.

7.1. 오일러 수열을 이용한 베르누이 수 표현

sechxsinhx=tanhx{\rm sech}\,x\sinh x = \tanh x이므로
{n=0E2n(2n)!x2n}{n=0x2n+1(2n+1)!}=n=1(16n4n)B2n(2n)!x2n1n=0r=0nE2rx2r(2r)!x2n2r+1(2n2r+1)!=n=0r=0n1(2n+1)!(2n+1)!E2r(2r)!(2n2r+1)!x2n+1=n=0r=0n1(2n+1)!(2n+12r)E2rx2n+1=n=1r=0n11(2n1)!(2n12r)E2rx2n1(16n4n)B2n(2n)!=r=0n11(2n1)!(2n12r)E2rB2n=2n16n4nr=0n1(2n12r)E2r\displaystyle\begin{aligned} \left\{ \sum_{n=0}^\infty \frac{E_{2n}}{(2n)!}x^{2n} \right\}\left\{ \sum_{n=0}^\infty \frac{x^{2n+1}}{(2n+1)!} \right\} &= {\color{blue}\sum_{n=1}^\infty}{\color{red}\frac{(16^n - 4^n)B_{2n}}{(2n)!}}{\color{blue}x^{2n-1}} \\ \sum_{n=0}^\infty \sum_{r=0}^n \frac{E_{2r}x^{2r}}{(2r)!} \frac{x^{2n-2r+1}}{(2n-2r+1)!} &= \sum_{n=0}^\infty \sum_{r=0}^n \frac 1{(2n+1)!} \frac{(2n+1)!E_{2r}}{(2r)!(2n-2r+1)!}x^{2n+1} = \sum_{n=0}^\infty \sum_{r=0}^n \frac 1{(2n+1)!} \binom{2n+1}{2r}E_{2r}x^{2n+1} \\ &= {\color{blue}\sum_{n=1}^\infty}{\color{red}\sum_{r=0}^{n-1} \frac 1{(2n-1)!} \binom{2n-1}{2r}E_{2r}}{\color{blue}x^{2n-1}} \end{aligned} \\ \frac{(16^n - 4^n)B_{2n}}{(2n)!} = \sum_{r=0}^{n-1} \frac 1{(2n-1)!} \binom{2n-1}{2r}E_{2r} \\ \therefore B_{2n} = \frac{2n}{16^n - 4^n}\sum_{r=0}^{n-1} \binom{2n-1}{2r}E_{2r}
오일러 수열이 정수 수열이고 조합도 자연수이기 때문에 결과적으로 연산 자체는 정수의 사칙연산이 된다. 분수끼리 더하고 빼야하는 베르누이 수열의 점화식 계산보다는 훨씬 수월할 것이다.

7.2. 베르누이 수열을 이용한 오일러 수열 표현

coshxsinhxtanhx=sechx\cosh x - \sinh x\tanh x = {\rm sech}\,x이므로, sinhxtanhx\sinh x\tanh x부분에 대해
{n=0x2n+1(2n+1)!}{n=1(16n4n)B2n(2n)!x2n1}=n=1r=1n(16r4r)B2rx2r1(2r)!x2n2r+1(2n2r+1)!=n=1r=1n16r4r(2n+1)!(2n+1)!B2r(2r)!(2n2r+1)!x2n=n=11(2n+1)!r=1n(16r4r)(2n+12r)B2rx2n\displaystyle\begin{aligned}&\left\{ \sum_{n=0}^\infty \frac{x^{2n+1}}{(2n+1)!} \right\} \left\{\sum_{n=1}^\infty \frac{(16^n - 4^n)B_{2n}}{(2n)!}x^{2n-1} \right\} \\ &= \sum_{n=1}^\infty \sum_{r=1}^n \frac{(16^r - 4^r)B_{2r}x^{2r-1}}{(2r)!} \frac{x^{2n-2r+1}}{(2n-2r+1)!} = \sum_{n=1}^\infty \sum_{r=1}^n \frac{16^r - 4^r}{(2n+1)!} \frac{(2n+1)!B_{2r}}{(2r)!(2n-2r+1)!}x^{2n} \\ &= \sum_{n=1}^\infty \frac1{(2n+1)!}\sum_{r=1}^n (16^r - 4^r)\binom{2n+1}{2r}B_{2r}x^{2n}\end{aligned}
따라서 sechx{\rm sech}\, x에 관한 등식은 다음과 같이 되며
coshxsinhxtanhx=sechx=n=0x2n(2n)!{n=11(2n+1)!r=1n(16r4r)(2n+12r)B2rx2n}=n=0E2n(2n)!x2n=1+n=1E2n(2n)!x2n=1+n=11(2n)!x2n{n=11(2n+1)!r=1n(16r4r)(2n+12r)B2r}x2n=1+n=1{1(2n)!1(2n+1)!r=1n(16r4r)(2n+12r)B2r}x2n1(2n)!1(2n+1)!r=1n(16r4r)(2n+12r)B2r=E2n(2n)!E2n=112n+1r=1n(16r4r)(2n+12r)B2r\displaystyle\begin{aligned}&\cosh x - \sinh x\tanh x = {\rm sech}\,x \\ &= \sum_{n=0}^\infty \frac{x^{2n}}{(2n)!} - \left\{ \sum_{n=1}^\infty \frac1{(2n+1)!}\sum_{r=1}^n (16^r - 4^r)\binom{2n+1}{2r}B_{2r}x^{2n}\right\} = \sum_{n=0}^\infty \frac{E_{2n}}{(2n)!}x^{2n} = {\color{blue}1 + \sum_{n=1}^\infty}{\color{red}\frac{E_{2n}}{(2n)!}}{\color{blue}x^{2n}} \\ &= 1 + \sum_{n=1}^\infty \frac1{(2n)!}x^{2n} - \left\{ \sum_{n=1}^\infty \frac1{(2n+1)!}\sum_{r=1}^n (16^r - 4^r)\binom{2n+1}{2r}B_{2r}\right\}x^{2n} \\ &= {\color{blue}1 + \sum_{n=1}^\infty}{\color{red}\left\{ \frac1{(2n)!} - \frac1{(2n+1)!} \sum_{r=1}^n (16^r - 4^r)\binom{2n+1}{2r}B_{2r} \right\}}{\color{blue}x^{2n}}\end{aligned} \\ \frac1{(2n)!} - \frac1{(2n+1)!} \sum_{r=1}^n (16^r - 4^r)\binom{2n+1}{2r}B_{2r} = \frac{E_{2n}}{(2n)!} \\ \therefore E_{2n} = 1 - \frac1{2n+1} \sum_{r=1}^n (16^r - 4^r)\binom{2n+1}{2r}B_{2r}
r=0r=0이면 (16r4r)(2n+12r)B2r=0(16^r - 4^r)\dbinom{2n+1}{2r}B_{2r} = 0이므로 합의 기호 부분은 r=0r=0부터 더해주는 것으로 바꿔도 무관하다. 즉
E2n=1+12n+1r=0n(4r16r)(2n+12r)B2r\displaystyle E_{2n} = 1 + \dfrac1{2n+1} \sum_{r=0}^n (4^r - 16^r)\binom{2n+1}{2r}B_{2r}
한편 12n+1(2n+12r)=1(2n+1)(2n+1)!(2r)!(2n2r+1)!=(2n)!(2r)!(2n2r+1)(2n2r)!=12n2r+1(2n2r)\dfrac1{2n+1}\dbinom{2n+1}{2r} = \dfrac1{(2n+1)}\dfrac{(2n+1)!}{(2r)!(2n-2r+1)!} = \dfrac{(2n)!}{(2r)!(2n-2r+1)(2n-2r)!} = \dfrac1{2n-2r+1}\dbinom{2n}{2r}이므로
E2n=1+r=0n4r16r2n2r+1(2n2r)B2r\displaystyle E_{2n} = 1 + \sum_{r=0}^n \frac{4^r - 16^r}{2n-2r+1}\binom{2n}{2r}B_{2r}
로도 나타낼 수 있다. 어느 식이든 베르누이 수열이 유리수 수열이기 때문에 오일러 수열로 나타낸 베르누이 수열과는 달리 이쪽은 오히려 계산이 복잡해진다.

[1] 이것과 비슷한 성질의 수열로서 오일러 수열이 있는데 이 수열은 모든 홀수항이 0\bf0이다.[2] 당초 이 수열의 발견자인 야콥 베르누이 본인이 B1=12B_1 = \dfrac12인 수열을 BnB_n이라 정의했었는데, 훗날 연구를 통해 생성함수로 더 엄밀하게 정의될 수 있다는 것이 알려진 뒤 베르누이가 최초로 정의한 BnB_n은 사실 Bn+B^+_n임이 밝혀졌다.[3] 페르마도 이를 연구했었다! 사실 그는 구적법 때문에 거듭제곱 합의 중요성에 대해 인지하고 있었고, 그 일반식을 얻었으며 증명까지 해냈다고 했으나, 그 내용에 대해 자세한 기록을 남기지는 않았다(……)페르마가 또[4] 일본에서 출판되는 일부 교양 수학서들 중 세키 - 베르누이 수열이라는 명칭을 쓰는 게 있긴 하다.[5] 더해지는 수열 ana_n의 종류에 관계없이 α<β\alpha<\beta에 대해 합의 범위가 i=βαan\displaystyle\sum_{i=\beta}^\alpha a_n으로 주어지는 것.[6] 점화식이라고 생각하는 게 차라리 낫다.