6
エンタメ解説
文献あり

中央二項係数っていつもそうですよね?私の事どう思っているんですか?

510
0

中央ニ項係数って素敵♡

んちゃ!
今回は中央二項係数で作って遊ぶよ!
この記事は書きかけだよ!仕事しながらだから許して...

ベータ関数とに項係数を含む和

B(m,n)=(m1)!(n1)!)(m+n1)!=01xm1(1x)n1dx(m,nN)

(nk)1=(n+1)01xk(1x)nkdx

Lehmerの仕事

{n=0an(2nn)n=0an(2nn)

y=(arcsinx)2y=2arcsinx1x2y=21x2(xarcsinx1x2+1)
ゆえに以下の微分方程式を得る。
(1x2)yxy=2

解析的解

y=n=0anxnと置きこれを先の微分方程式に代入すると次の様に計算できる。
2a2+n=1{(n+2)(n+1)an+2n2an}=2a2=1an+2=n2(n+2)(n+1)an
a2n=(2n2)22n(2n1)a2(n1)=22n1{(n1)!}2(2n)!a2=22n1(n!)2n2(2n)!=22n1n2(2nn)
a0=0a1=a3==0
(arcsinx)2=12n=122nn2(2nn)x2n

1(2nn)=nB(n,n+1)

n=0xn(2nn)=01n=0nxntn1(1t)ndt=01x(1t)n=0n(xt(1t))n1dt=01x(1t)(1xt(1t))2

n=0xn(2nn)=4x(4x)32arctan(x4x)+x4x+1

n=1xnn(2nn)=4x4xarctanx4x

n=1xnn(2nn)=01n=1xntn1(1t)ndt=01x(1t)n=1(xt(1t))n1dt=01x(1t)1xt(1t)dt
xt2xt+1=x(t12)2+4x4
01x(1t)1xt(1t)dt=011t(t12)2+4x4xdt(t12=4x4xtanθ)=4x4xarctanx4xarctanx4x(124x4xtanθ)dθ=4x4x0arctanx4xdθ=4x4xarctanx4x=2x4xarctanx4x

ddxn=1xnn(2nn)=n=1xn1(2nn)=4(4x)32xarctanx4x+24x
n=0xn(2nn)=4x(4x)32arctanx4x+x4x+1

x=1の場合
{n=11n(2nn)=π33n=01(2nn)=2327π+43

k=0xk(2k+1k)=24x+8(4x)4xx2arctanx4xx2

k=0xk(2k+1k)=01k=0(2k+2)xktk(1t)k+1dt=2xt01k=0k(xt(1t))kdt=2xt01xt(1t)(1xt(1t))2dt=2011t(1xt(1t))2dt(1t=τ)=201τ(1xτ(1τ))2dτ(τ12=4x4xtanθ)=2x216x2(4x)24x4xarctanx4xarctanx4x(4x4xtanθ+12)cos2θdθ=16x(4x)320arctanx4xcos2θdθ=8x(4x)32arctanx4x+24x=8(4x)4xx2arctanx4xx2+24x(0<x<4)=4(4x)4x+x2log14x+x21+4x+x2(4<x<0)
最後の変形は次の事実を用いた。
{θ=arctanλsinθ=λ1+λ2cosθ=11+λ21iarctanθi=12log1θ1+θ

k=0xk(2k+1k)=D

m=0xm(2m+1m)(2m+3)=(4x1)D4x

実際に計算して係数比較するだけ。
4(k+1)(k+2)(2k+3)(2k+2)=2(k+1)

f(x)=n=0anxnとすると次の事が成り立つ。
n=1an1xnn(2nn)=01x(1t)f(xt(1t))dt

S(k;x)=n=0nkxnとすると次の様に書ける。
n=0nk(2nn)=201S(k+1;x(1x))dx+01S(k;x(1x))dx

n=0nk(2nn)=n=0nkn!n!(2n)!=n=0nkΓ(n+1)Γ(n+1)Γ(2n+1)=n=0nk(2n+1)B(n+1,n+1)=201n=0nk+1xn(1x)ndx+01n=0nkxn(1x)ndx=201S(k+1;x(1x))dx+01S(k;x(1x))dx

S(k;x)=n=0nkxnxddxS(k;x)=S(k+1;x)

xddxS(k;x)=xddxn=0nkxn=xn=0nk+1xn1=n=0nk+1xn=S(k+1;x)

S(k;x)=(xddx)kS(0;x)=(xddx)k11x

(xddx)kdidxi(i=1,2,...,k)の係数は次の様に書ける。
ただし、第k行、第i列は第k1行の第i1,i列の1,i倍して足し合わせたもの。

1
11
131
1761
11525101

k=0とすると次の事が成り立つ。
n=01(2nn)=201x(1x)(1x(1x))2dx+0111x(1x)dx=2011(1x(1x))2dx0111x(1x)dx=8327π239π=2327π+43

k=0nakbnk(nk)=n+1(a+b)(1a+1b)nk=1n+1(ak+bk)(1a+1b)k1k

n=0xnk=0nakbnk(nk)=n=0k=0n(n+1)xnakbnk01tk(1t)nkdt=ddx(x01k=0n=k(axt)k(bx(1t))nkdt)=ddx(x0111axt11bx(1t))dt=ddx(1a+babx01(ax1axt+bx1bx+bxt)dt)=ddx(log(1ax)log(1bx)a+babx)=ddx[1a+bn=0anbnxn(a+b)nk=1ak+bkkxk]=ddx[x1a+bn=0anbnxn(a+b)nk=1ak+bkkxk1]=ddx[x1a+bn=0anbnxn(a+b)nk=0ak+1+bk+1k+1xk]=ddx[x1a+bλ=0xλk=0λakbk(a+b)kaλk+1+bλk+1λk+1]=ddx[1a+bλ=0xλ+1k=0λaλkbλk(a+b)λkak+1+bk+1k+1]=ddx[1a+bn=0anbnxn+1(a+b)nk=0nak+1+bk+1k+1(a+b)kakbk]=1a+bn=0(n+1)anbnxn(a+b)nk=1n+1ak+bkk(a+b)k1ak1bk1
xの各次数の係数比較を行うと次に結論を得る。
k=0nakbnk(nk)=n+1(a+b)(1a+1b)nk=1n+1(ak+bk)(1a+1b)k1k

a=b=1とすると次式が成り立つ。
k=0n1(nk)=n+12n+1k=1n+12kk=(n+1)k=0n1(n+1k)2k

a=x,b=1とおくと次式を得る。
k=0nxk(nk)=n+1(1+z)(z1+z)nk=1n+11+zkk(1+zz)k

逆べき生成変換

次の変換を逆べき生成変換と呼ぶことにする。
I(f|x)=0xf(t)f(0)tdt
また、nNに対して次の様な変換を定義する。これをn階逆べき生成変換と呼ぶことにする。
In(f|x)=0xdx10x1dx20xn1dxnf(xn)f(0)xn

f(x)=n=0anxnk回逆べき生成変換を行った結果は次の様に書ける。
Ik(f|x)=n=1annkxn

q級数とに項係数を含む総和

q級数のEuler展開

Πk=1(1qk)=1+k=1(1)k(qk(3k1)2+qk(3k+1)2)

[1]f(q;x)=1k=1qkxk+1Πl=1k1(1qlx)と置く。するとこれは次の式を満たす。
f(q;x)=1qx2q2x3f(q;qx)
これは次の様にして証明できる。
f(q;x)=1qx2k=2qkxk+1Πl=1k1(1qlx)=1qx2k=1qk+1xk+2Πl=1k(1qlx)=1qx2k=1qk+1xk+2Πl=2k(1qlx)+k=1qk+2xk+3Πl=2k(1qlx)=1qx2q2x3(1+k=1qkxk(Πl=2k+1(1qlx)Πl=2k(1qlx)))=1qx2q2x3(1k=1qk(qx)k+1Πl=2k(1ql1(qx)))=1qx2q2x3(1k=1qk(qx)k+1Πl=1k1(1ql(qx)))=1qx2q2x3f(q;qx)
[2]
f(q;x)=1qx2q2x3f(q;qx)=1qx2q2x3(1q3x2q5x3f(q;q2x))=1qx2q2x3+q5x5+q7x6f(q;q2x)=1qx2q2x3+q5x5+q7x6(1q5x2q8x3f(q;q3x))=1qx2q2x3+q5x5+q7x6q12x8q15x9f(q;q3x)
ここで少し考察をしてみる。とりあえず符号は抜きにして、qx2,q5x5,q12x8,...;q2x3,q7x6,q15x9,...のべきをテーブルとして並べてみる。

(1,2)(2,3)
(5,5)(7,6)
(12,8)(15,9)
このテーブルを見るとわかると思うが、前者は(n(3n1)2,3n1)後者(n(3n+1)2,3n)なる規則で並んでいることが分かる。よって、これを示す。
[3]
これを示すには、上記方法をN回使用した際に出てくるf(q;qNx)にかけられているqμxνμ,νを求めればよい。
f(q;qNx)だけに着目し符号を無視すると、q2x3q2(qx)3q2(qN1x)3=qN2(3N+1)x3N
ゆえに、q2x3,q7x6,q15x9,...のグループのべきの組は(N2(3N+1),3N)また、qx2,q5x5,q12x8,...のグループの場合はf(q,x)を始まりとして先の結果にq(qN1x)2=q2N1x2をかけたものなので、N12(3N2)+2N1=N(3N1)2
[4]
よって結論
f(q;x)=1+n=1(1)n(qn(3n1)2x3n1+qn(3n+1)2x3n)
[5]
帰納法を使うと次の式を簡単に証明できる。
1n=1NΠk=1n1(1qk)qn=Πn=1N(1qn)
[6]
[5]よりx=1とすれば次の式が成り立つ。
f(q;1)=Πn=1(1qn)=n=(1)nqn(3n+1)2

(α)n=Πk=0n1(1qk1a)

(α)n=(α)(qnα)

q-超幾何級数

2ψ1(α,β;γ|x)=n=0(α)n(β)n(γ)n(q)nxn

(az)(z)=m=0(α)m(c)mzm

Jackson積分

0cf(t)dqt=c(1q)n=0f(cqn)qn

k=0(xk)(ynk)=(x+yn)

x+y個の異なるものをn個選ぶ方法の数はx個のものA、y個の異なるものBに分け、Aからk個、Bからnk個選択する方法の数で分け、k=0,1,...,nすべてに渡り総和を取ったものになるので証明完了。

k=0n(nk)(x+ynxk)=(x+yx)

k=0n(nk)(1)kkm={0if  m<n(1)nn!if  m=n

f(x)=(1x)nとする。するとこの式は次の様に書ける。
f(x)=k=0n(1)k(nk)xk
(xddx)mf(x)=k=0n(1)k(nk)kmxk
定理15より(xddx)m=1++xmdmdxmとかけるので、m<nの場合は(ddx)mf(x)1xで割り切れるのでx=1を代入すれば0。m=nの場合は、1xで割り切れない項は(1)nn!のみ。
以上より証明は完了した。

k=0n(nk)(1)k(xk+y)m={0if  m<n(1)nxnn!if  m=n

の中のkの最大べきにのみ着目すればいい。kの最大べきがn未満の場合はすべて0
m=nの場合はkn部分だけが0でない値を持つ。ゆえにその値は(1)nxnn!

f(t)=c0+c1t++cmtmの場合
k=0n(nk)(1)kf(k)m={0if  m<n(1)ncnn!if  m=n

超幾何関数と二項係数

x=0近傍の超幾何関数は下記の様に定義される。
y=2F(α,β;γ|x)=n=0(α)n(β)n(γ)nn!xn(|x|<1)

=ddx,D=xとすると下記の式が成り立つ。
((D+a)(D+b)D(D+c1)x)y=0

D2=x22+xより
(D+a)(D+b)=D2+(a+b)D+ab=x22+(a+b+1)x+ab
D(D+c1)=x22+cx
より
((D+a)(D+b)D(D+c1)x)2F1(a,b;c|x)=(x(x1)2+((a+b+1)xc)+ab)n=0(a)n(b)n(c)nn!xn=n=0[((n+1)n+(n+1)c)(a+n)(b+n)(c+n)(n+1)+(n(n1)+(a+b+1)n+ab)](a)n(b)n(c)nn!xn=0

超幾何微分方程式

(x(1x)2+(c(a+b+1)x)ab)y=0

両辺をx(1x)で割ると
(2+(cx(1x)a+b+11x)ab)y=(2+(cx1+a+bc1x)abx(1x))y=0
よりx=1の近傍での解はe=1+a+bcとおけば2F1(a,b;e|1x)
これを示したいなら1xxなる変数変換を行えばいい。

超幾何関数y=2F1(a,b;c|x)は次の微分方程式を満たす。
ψ(x)=xc(1x)e(e=1+a+bc)とすると次式が成り立つ。
(ψ(x)abψ(x)x(1x))y=0

f(x)=f(x)+f(x)
ψ(x)=(xc+cxc1)(1x)e=xc((1x)ee(1x)e1)+cxc1(1x)e=ψ(x)(e1x+cx)
これを代入すると
(ψ(x)abψ(x)x(1x))y=ψ(x)(2+(cxe1x)abx(1x))y=0

D=f(x)g(x)とし、Dy=0の解をF(x)と置く。

F(z(x))の満たす微分方程式は下記の様に書ける。
(f(z(x))z(x)z(x)g(z(x)))F(z(x))=0

z=1z(x)
より
(f(z(x))z(x)z(x)g(z(x)))F(z(x))=z(x)(zf(z)zg(z))F(z)=0

ψ(x)=xc(1x)e(e=1+a+bc)
(ψ(x)abψ(x)x(1x))y=0(y(0)=1,y(0)=abc)
に対してz(x)=xsとすると次式が成り立つ。
(xcss+1(1xs)es2abxcs1(1xs)e1)y(xs)=0

z(x)=1x1+(r1)xとすると以下の計算が成り立つ。
1z=rx1+(r1)xz(x)=r(1+(r1)x)2
(xe(1x)c)re1(1+(r1)x)c+e2+αβxe1(1x)c1re(1+(r1)x)c+e=re1((xe(1x)c)1(1+(r1)x)c+e2αβrxe1(1x)c11(1+(r1)x)c+e)
これより、以下の式が成り立つ。
ρ(x)=11+(r1)x(xe(1x)cρc+e2(x)rabxe1(1x)c1ρc+e(x))y=0

z(x)=(1x1+(r1)x)s=(1xρ(x))sとすると、z(x)=rs(1x)s1ρs+1(x)ゆえに次の微分方程式を得る。
(r,s)=(2,2),(3,3),(4,2)
{(1+x)2(1x)2=4x4x(1x)2(1+2x)3(1x)3=9x1x31x9x(1x)(1x)(1+2x)2(1+3x)2(1x)2=8x1x21x8x(1x2)(1x)(1+3x)2

Di=fi(x)gi(x)(i=1,2)に対する解をF1,F2とする。このとき、ある適当な関数hが存在して次の事が成り立つとする。
h(x)F2(x)=F1(x)
このとき、hD1h=D2が成り立つとすると
{f2(x)=f1(x)h2(x)g2(x)=g1(x)h2(x)(f1(x)h(x))h(x)

hh=hを用いる。
hf1h=(hh)f1(h+h)=f1h2+f1hhf1hhf1(h)2=f1h2+(f1h)h
ゆえに
hD1h=f1h2+(f1h)hg1h2=f1h2(g1h2(f1h)h)

上記hが存在するとすると下記の式が成り立つ。
D1h(x)F2(x)=1h(x)D2F2(x)=0

一次変換

{(1x)a+bc2F1(a,b;c|x)=2F1(ca,cb;c|x)(1x)a2F1(a,b;c|x)=2F1(a,cb;c|xx1)

[1]次の様に記号を決める。
{h=(1x)a+bcF1=2F1(ca,cb;c|x)F2=2F1(a,b;c|x)
する
{D1=xc(1x)cab+1(ca)(cb)xc1(1x)cabD2=xc(1x)a+bc+1abxc1(1x)a+bc
f1h2=xc(1x)cab+1(1x)2a+2b2c=xc(1x)a+bc+1
g1h2(f1h)h=(ca)(cb)xc1(1x)a+bc+(a+bc)cxc1(1x)a+bc=abxc1(1x)a+bc=g2
初期条件検査
{F1(0)=(hF2)(0)=1F1(0)=(hF2)(0)=(ca)(cb)c
[2]次の様に記号を定める。
{h=(1x)aF1=2F1(a,cb;c|xx1)F2=2F1(a,b;c|x)
{D1=xc(1x)a+bc+1+a(cb)xc1(1x)a+bc1D2=xc(1x)a+bc+1abxc1(1x)a+bc
以上より、
f1h2=xc(1x)a+bc+1(1x)2a=xc(1x)a+bc+1=f2
g1h2(f1h)h=a(cb)xc1(1x)a+bc1axc1(1x)a+bc1(bxc)=abxc1(1x)a+bc=g2
初期条件検査
{F1(0)=(hF2)(0)=1F1(0)=(hF2)(0)=a(bc)c

{(1+x)a2F1(a2,ab+12;b+12|x2)=2F1(a2,b2;b|1(1x1+x)2)(1+2x)a2F1(a3,a+13;a+56|x3)=2F1(a3,a+13;a+12|1(1x1+2x)3)(1+3x)a22F1(a4,a+24;a+56|x2)=2F1(a4,a+24;a+23|1(1x1+3x)2)

A.(r,s)=(2,2),ρ=1+x
{z=(1x)2ρ21z=4xρ2z=41xρ3
{h=(1+x)a=ρaF1=2F1(a2,b2;b|1(1x1+x)2)F2=(1+x)a2F(a2,ab+12;b+12|x2)
{D1=xb(1x)ab+1ρab+1abxb1(1x)ab+1ρab1 =xb(1x2)ab+1ρ2aabxb1(1x2)ab+1ρ2a2D2=xb(1x2)ab+1a(ab+1)xb(1x2)ab
f1h2=xb(1x2)ab+1=f2
(f1h)=axb(1x2)ab+1)(1+x)a1(bx2(ab+1)x1x2a+11+x)
(f1h)h=axb1(1x2)abρ1((ab+1)x2(a+1)x+b)
g1h2(f1h)h=axb1(1x2)ab+1ρ2+axb1(1x2)abρ1((ab+1)x2+(a+1)xb)=axb1(1x2)abρ2(ab+1)x(1+x)2=a(ab+1)xb(1x2)ab=g2
初期条件検査
{F1(0)=F2(0)=1F1(0)=(hF2)(0)=a
B.(r,s)=(3,3),ρ=1+2xとする。
{h=(1+2x)a=ρaF1=2F1(a3,a+13;a+12|1(1x1+2x)3)F2=2F1(a3,a+13;a+56|x3)
{D1=xa+12(1x3)a+12ρ2aa(a+1)xa12(1x)2(1x3)a12ρ2a2D2=xa+12(1x3)a+12a(a+1)xa+32(1x3)a12
f1h2=xa+12(1x3)a+12=f2
(f1h)=a(a+1)xa+12(1x3)a+12(1+2x)a+1(12x3x22(1x3)2(1+2x))
(f1h)h=a(a+1)xa12(1x3)a12ρ2(4x44x32x+1)
g1h2(f1h)h=a(a+1)xa12(1x)2(1x3)a12ρ2+a(a+1)xa12(1x3)a12ρ2(4x4+4x3+2x1)=a(a+1)xa+32(1x3)a12=g2
初期条件検査
{F1(0)=hF2(0)=1F1(0)=(hF2)(0)=2a
C.(r,s)=(4,2),ρ=1+3xとする。
{h=(1+3x)a2=ρa2F1=2F1(a4,a+24;a+23|1(1x1+3x)2)F2=2F1(a4,a+24;a+56|x2)
{D1=xa+23(1x2)a+23ρaa(a+2)2xa13(1x)(1x2)a13ρa2D2=xa+23(1x2)a+23a(a+2)4xa+23(1x2)a13
f1h2=xa+23(1x2)a+23=f2
f1h=3a2xa+23(1x2)a+23(1+3x)a+22
(f1h)=3a(a+2)2xa+23(1x2)a+23ρa+22(13x2x3(1x2)32(1+3x))
(f1h)h=a(a+2)4xa13(1x2)a13ρ2(9x36x23x+2)
g1h2(f1h)h=a(a+2)4xa+23(1x2)a13=g2
初期条件検査
{F1(0)=hF2(0)=1F1(0)=(hF2)(0)=3a2

K(x)=0π2dθ1x2cos2θ=π22F1(12,12;1|x2)

{11x=n=0(2n)!22n(n!)2xn0π2cos2nθdθ=(2n)!22nn!π2(12)n=(2n)!22nn!
を用いればよい。
(12)n(12)n(1)n=((2n)!)224n(n!)3

K(x)=π22F1(12,12;1|x2)=π2(1+x)2F1(12,12;1|4x(1x)2)

差分と加速

数列{an}nに関して次の演算子Δabを一般化差分と呼ぶことにする。
Δaban=aan+ban+1

{Δab0an=anΔabkan=Δab(Δabk1an)(kN)
の様に定義すると、次の式が成り立つ。
Δabkan=l=0kkClaklblan+l

k=1の場合は明らか。
1,2,...,k1場合まで成り立つと仮定する。次にkの場合を考える。
Δabkan=Δ(l=0k1k1Clalbk1lan+l) =k1C0akb0an+k1C1ak1b1an+1+k1C2ak2b2an+2++k1Ck1a1bk1an+k2  +k1C0ak1b1an+1+k1C1ak2b2an+2++k1Ck2a1bk1an+k2+k1Ck1a0bkan+k1=akb0+l=1k1kClak1blan+l+a0bkan+k=l=0kkClaklblan+l

{f(x)=n=0Anxnf(l)(x)=n=0n!An(nl)!xnl=l!n=lnClAnxnl
とすると次式が成り立つ。
f(λΔab)a0=k=0AkλkΔabka0=k=0Akλkl=0kaklblkClal=l=0(λb)lall!l!k=lkClAk(λa)kl=l=0(λb)lall!f(l)(λa)

f(x)=11xとするとf(l)(x)=l!(1x)l+1なので
11λΔaba0=k=0(λΔab)ka0=l=0(λb)lal(1λa)l+1
a=b=1,λ=12とすると下記の式を得る。
k=0(Δ2)ka0=2l=0(1)lal
すなわち、下記の式を得る。
l=0(1)lal=12k=0(Δ2)ka0

f(x)=log(1x)=n=1xnnとするとf(n)(x)=(n1)!(1x)n
log(1λΔab)a0=n=1(λΔab)nna0=l=1all(λb(1a))l

数列変換

数列{sn}nNから数列{tn}nNへの写像T:{sn}{tn}を数列変換と呼ぶ。

Richarson補外

数列変換T:{sn}{tn}を次の様に定める。
まず次のような漸化式を数列{sn}を満たすものとする。
sn+1sλ(sns)(λ1)
λを既知とし、stnとする。すると次式が成り立つ。
tn=sn+1λsn1λ
この数列変換T:{sn}{tn}をRichardson補外という。

先のRichardson補外の場合でλが未知とした場合次の式が成り立つ。
λ=sn+1ssns=sn+2ssn+1s
この式でsについて解き、stnの様に置きなおす。
tn=snsn+2sn+12sn2sn+1+sn+2=sn(sn+1sn)2sn2sn+1+sn+2=sn(Δsn)2Δ2sn
この数列変換T:{sn}{tn}AitkenΔ2法という。
ちなみに、簡単な計算からtnは次の様に書ける事も分かる。
tn=sn+1(sn+1sn)(sn+2sn+1)sn2sn+1+sn+2=sn+2(sn+2sn+1)2sn2sn+1+sn+2=sn+2(sn+2)22sn+2
ただしan+1=an+1anとした。

Euler-Maclaurin展開

k=mnf(k)=mnf(x)dx+12(f(m)+f(n))+k=2N(1)kbk(f(k1)(n)f(k1)(m))+RN

k=mn101ψ(x)f(x+k)dx=k=mn1([ψ(x)f(x+k)]0101ψ(x)f(x+k)dx)=ψ(1)f(m)+ψ(0)f(n)+k=mn(ψ(1)ψ(0))f(k)k=mn101ψ(x)f(x+k)dx
ここで、ψ(x)=x+aとおきψ(0)=ψ(1)=12とするとa=12
k=mn101ψ(x)f(x+k)dx=12(f(m)+f(n))+k=mnf(k)k=mn101f(x+k)dx=12(f(m)+f(n))+k=mnf(k)mnf(x)dx
よって、下記の様に書ける。
k=mnf(k)=mnf(x)dx+12(f(m)+f(n))+k=mn101ψ(x)f(x+k)dx
ここで次のような記号を定める。
{R1=k=mn101ψ(x)f(x+k)dxRN=(1)N1k=mn101BN(x)f(N)(x+k)dxB1(x)=ψ(x)BN(0)=BN(1)=bNBN(x)=BN1(x)
すると下記の式を得る。
RN=(1)N1bN(fN1(n)fN1(m))+RN1
ゆえに、以下の式を得る。
R1=b2(f(1)(n)f(1)(m))+R2=b2(f(1)(n)f(1)(m))+b3(f(2)(n)f(2)(m))+R3=k=2N(1)kbk(f(k1)(n)f(k1)(m))+RN
つまり、
k=mnf(k)=mnf(x)dx+12(f(m)+f(n))+k=2N(1)k1bk(f(k1)(n)f(k1)(m))+RN
ここで次のようなBk(x)の母関数F(t,x)=1+k=1Bk(x)tkとおく。すると次の様な計算ができる。
Fx=tFlogF=tx+f(t)F=g(t)etx
また、
{F(t,0)=12t+k=2bktk=g(t)F(t,1)=12t+k=2bktk=g(t)etF(t,1)F(t,0)=t=g(t)(et1)g(t)=tet1F(t,x)=tetxet1
以上から、次の結論を得る。
tetet1=1+12t+k=2bktkt2et+1et1=1+k=2bktk
またt2et+1et1=t2et+1et1より偶数関。ゆえに、b3=b5==0
さらにtについて展開し、tの同じべきで比較すると下記の式を得る。
2(1(2n)!+b2(2n2)!++b2n22!)=1(2n1)!(n=2,3,...)

{b2=112b4=1720

Euler-Maclaurin展開

k=mnf(k)=mnf(x)dx+12(f(m)+f(n))+k=1Nb2k(f(2k1)(n)f(2k1)(m))+R2N

f(x)=xsとすると
k=11ks=1s1+12+k=1Nb2kΓ(s+2k)Γ(s)+R2NR2N=(1)2N1Γ(s+2N)Γ(s)k=101B2N(x)1(x+k)s+2Ndx0(N)
{ζ(s)=1s1+12+k=1b2kΓ(s+2k)Γ(s)(1<s)b2(n1)=1(2n1)!2(1(2n)!+b2(2(n1))!++b2(n2)4!)

望遠鏡和

望遠鏡和

数列{an}nN,{bn}nNan=bnbn+1(nN)を満たすとする。すると下記の事実が成り立つ。
Sn=k=1nan=b1bn+1
さらに、limn|bn|=0ならば下記の式が成り立つ。
limnSn=k=1ak=b1
上記の事実を望遠鏡和という。

n=11an2+bについて、a,bを適切に定めることで望遠鏡和が成り立つようにせよ。
またその値について考察せよ。
n=11an2b2=12bn=1(1anb1an+b)=12b(1ab1a+b  +12ab12a+b   +13ab13a+b    +)
よりna+b=(n+1)abを満たすようにa,bを定めればa=2b。すると
n=114b2n2b2=12b2
すなわち、
n=114n21=12

Kummer変換

数列{an}nN,{bn}nNについてA=n=1an,B=n=1bnγ=limnbnanとする。すると
A=γB+n=1(anγbn)=γB+n=1(1γbnan)an

1(1315)(1719)=12k=11(4k1)(4k+1)=12k=1116k21=12A=π4
また
B=n=114n21=12
このとき、γ=limn116n2114n21=14
Kummer変換を用いると下記の式を得る。
A=1412+n=1(11414n21116n21)116n21=1834n=11(16n21)(4n21)
この式より以下の様な結果を得る。
π=48A=3+6n=11(16n21)(4n21)

Wilf–Zeilberger pair

関数F(n,k),G(n,k)をWZ-pairとすると下記の式を満たす。
F(n+1,k)F(n,k)=G(n,k+1)G(n,k)

k=ab[F(n+1,k)F(n,k)]=G(n,b+1)G(n,a)

k=ab[F(n+1,k)F(n,k)]=k=ab[G(n,k+1)G(n,k)]=G(n,a)+G(n,a+1)  G(n,a+1)+G(n,a+2)            G(n,b)+G(n,b+1)=G(n,b+1)G(n,a)

Gosperのアルゴリズム

Gosperのアルゴリズムについては下記の僕の記事を参照してください。
https://mathlog.info/articles/bHQmk1TsifHrtRUu9vuY

t(n)=nkp(n)とした場合
t(n+1)t(n)=(n+1)knkp(n+1)p(n)
より
{a(n)=(n+1)kb(n)=nkc(n)=p(n)
これからGosper方程式は次の様に書ける。
(n+1)kS(n+1)(n1)kS(n)=p(n)
仮にS(n)=2nであるとすると次の様に書ける。
p(n)=2n(2(n+1)k(n1)k)
ゆえに
u(n)=(n1)kp(n)2nnkp(n)=2n(n(n1))k
これにより次の事実が分かる。
n=0Nnk2n(2(n+1)k(n1)k)=2N+1(N(N+1))k

t(n)=An(2nn)p(n)した場合
t(n+1)t(n)=2A(2n+1)n+1p(n+1)p(n)
より
{a(n)=2A(2n+1)b(n)=n+1c(n)=p(n)
ゆえに、Gosper方程式は次の様に書ける。
2A(2n+1)S(n+1)nS(n)=p(n)
A=12とすると
(2n+1)S(n+1)nS(n)=p(n)
S(n)=1とするとp(n)=n+1
これから、
u(n)=nn+1(2nn)2n(n+1)=n(2nn)2n
ゆえに
n=0N(n+1)(2nn)2n=(N+1)(2N+2N+1)2N+1

t(n)=An(2nn)kp(n)とした場合
t(n+1)t(n)=A2k(2n+1)k(n+1)kp(n+1)p(n)
より
{a(n)=A2k(2n+1)kb(n)=(n+1)kc(n)=p(n)
この事から求めるGosper方程式は次の様になる。
A2k(2n+1)kS(n+1)nkS(n)=p(n)
ここでS(n)=(2n+1)kとすると
p(n)=A22k(n+1)knk
特にA=122kとすると
p(n)=(n+1)knk
u(n)=nkp(n)(2n+1)k(2nn)k22knp(n)=nk(2n+1)k(2nn)k4kn
これから次のような結果を得る
n=0N(2nn)k4kn((n+1)knk)=(N+1)k(2N+3)k(2N+2N+1)k4k(N+1)

超幾何双和

F(n,s1,...,se)s1,...,se,nそれぞれに対して超幾何項であるとする。この時次の様な問題を考える。
k=0γpk(n)S(n+k)=0
ただし、
S(n)=s1seF(n,s1,...,se)

e=1としs1=rと置きなおす。
このとき
k=0γpk(n)f(n+k,r)=g(n,r+1)g(n,r)
であるとすると下記の事実が成り立つ。
k=0γpk(n)S(n+k)=limRg(n,R+1)g(n,0)

参考文献

投稿日:2024103
更新日:2024119
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. 中央ニ項係数って素敵♡
  2. ベータ関数とに項係数を含む和
  3. 逆べき生成変換
  4. q級数とに項係数を含む総和
  5. 超幾何関数と二項係数
  6. 差分と加速
  7. Euler-Maclaurin展開
  8. 望遠鏡和
  9. Wilf–Zeilberger pair
  10. Gosperのアルゴリズム
  11. 超幾何双和
  12. 参考文献