跳转至

Entringer Number

恩特林格数

恩特林格数(Entringer number,OEIS A008281E(n,k) 是满足下述条件的 0nn+1 个数的置换数目:

  • 首元素是 k
  • 首元素的下一个元素比首元素小,再下一个元素比前一个元素大,再下一个元素比前一个元素小……后面相邻元素的大小关系均满足这样的规则。

恩特林格数的初值有:

E(0,0)=1E(n,0)=0

有递推关系:

E(n,k)=E(n,k1)+E(n1,nk)

Seidel–Entringer–Arnold 三角

恩特林格数的一个适当排列的数字三角,称为 Seidel–Entringer–Arnold 三角(Seidel–Entringer–Arnold triangle,OEIS A008280)。该三角是按照「牛耕」顺序(ox-plowing order)排列的恩特林格数 E(n,k)

E(0,0)E(1,0)E(1,1)E(2,2)E(2,1)E(2,0)E(3,0)E(3,1)E(3,2)E(3,3)E(4,4)E(4,3)E(4,2)E(4,1)E(4,0)

即:

101110012255420

按照这种方式排列的恩特林格数的优势是,与它的递推关系 E(n,k)=E(n,k1)+E(n1,nk) 一致,可以方便记忆和理解。

恩特林格数有一个指数型生成函数:

m=0n=0E(m+n,12(m+n+(1)m+n(nm)))xmm!xnn!=cosx+sinxcos(x+y)

这个生成函数的系数分布事实上是上面的 Seidel–Entringer–Arnold 三角的简单拉伸变形:

E(0,0)E(1,1)E(2,0)E(3,3)E(4,0)E(1,0)E(2,1)E(3,2)E(4,1)E(2,2)E(3,1)E(4,2)E(3,0)E(4,3)E(4,4)

即:

110200122114055

zigzag 置换

一个 zigzag 置换(zigzag permutation)是一个 1n 的排列 c1ci,使得任意一个元素 ci 的大小都不介于 ci1ci+1 之间。

对于 zigzag 置换的个数 ZnOEIS A001250),从 n=0 开始有:

1,1,2,4,10,32,122,544,

例如,前几个 n 的交替置换有:

n=1:{1}n=2:{1,2},{2,1}n=3:{1,3,2},{2,1,3},{2,3,1},{3,1,2}n=4:{1,3,2,4},{1,4,2,3},{2,1,4,3},{2,3,1,4},{2,4,1,3},{3,1,4,2},{3,2,4,1},{3,4,1,2},{4,1,3,2},{4,2,3,1}

交替置换与 zigzag 数

(注意和「错位排列」进行概念上的区分。)

对于大于 1n,每个 zigzag 置换翻转过来仍旧为 zigzag 置换,可以两两配对,所以必然为偶数。

这里再给出一种配对的方法:将 zigzag 置换分为交替置换(alternating permutation)和反交替置换(reverse alternating permutation)。

交替置换的首元素大于第二个元素,大小关系为:

c1>c2<c3>

反交替置换的首元素小于第二个元素,大小关系为:

c1<c2>c3<

如果将 1n 位置互换,2n1 位置互换,以此类推,即可将交替置换与反交替置换两个集合互换。因此,交替置换与反交替置换的个数相等,恰好为 zigzag 置换的一半。

对于大于 1n,记:

An=Zn2

定义初值:

A0=A1=1

这里的 An 称为 zigzag 数(Euler zigzag number,OEIS A000111),从 n=0 开始有:

1,1,1,2,5,16,61,272,

接下来试着求解 An

1n 之中,选取 k 个数构成子集,有 (nk) 种选法。

在这个 k 元子集中,选反交替置换 u,有 Ak 种选法;用全集减掉这个 k 元子集,剩余的 nk 元子集中,选反交替置换 v,有 Ank 种选法。

考虑 n+1 元排列 w,将 u 倒置作为开头,接上 n+1,再接上 v。那么,w 一定是 zigzag 置换,并且任意一个 n+1 元 zigzag 置换,都可以在 n+1 处截断得到对应的反交替置换 uv,并且不同的 n+1 元 zigzag 置换对应的 uv 不同。

因此有递推关系:

2An+1=k=0n(nk)AkAnk2(n+1)An+1(n+1)!=k=0nAkk!Ank(nk)!

n0 时并不满足这个递推式,初值 A0A1 都是 1

可见,这是一个指数型生成函数的卷积。假设 An 的指数型生成函数为 y,就有微分方程:

2dydx=y2+1

等式右面加 1 是为了处理 n0 时的特殊情况。该方程的通解为:

y=tan(12x+C)

代入第 0 项为 1 之后,可以得到特解:

y=tanx+secx

正切函数是奇函数,正割函数是偶函数,两者之和构成 zigzag 数的生成函数。

恩特林格数与 zigzag 数的关系

根据恩特林格数的定义,恩特林格数 E(n,k) 是首元素为 k0n 的交替置换个数。因此恩特林格数与 zigzag 数事实上有关系:

An=E(n,n)

An 称为「zigzag 数」也有原因:记 En 是欧拉数(Euler number),Bn 是伯努利数。

n 为偶数时,偶数项下标的 zigzag 数也称「正割数」Sn 或者「zig 数」。有关系:

An=(1)n/2En

前几项为(OEIS A000364):

1,1,5,61,1385,

n 为奇数时,奇数项下标的 zigzag 数也称「正切数」Tn 或者「zag 数」。有关系:

An=(1)(n1)/22n+1(2n+11)Bn+1n+1

前几项为(OEIS A000182):

1,2,16,272,7936,

于是对于在 x=0 处的泰勒展开,可以给出正割数和正切数:

secx=A0+A2x22!+A4x44!+tanx=A1x+A3x33!+A5x55!+

或者写到一起:

secx+tanx=A0+A1x+A2x22!+A3x33!+A4x44!+A5x55!+

构成 zigzag 数的生成函数。

参考资料与链接

  1. Alternating permutation - Wikipedia