形式陈述
令 A ( n , k ) 采用恰有 $k$ 个下降 理路 Eulerian 数与排列下降 Eulerian number · Eulerian polynomial · 排列下降数 以插入最大元素证明下降数递推,推导Eulerian多项式并用n=4分位置核验十一项。 的零起点约定。在 Q [ x ] (或任意特征零域上的多项式环)中,对 n ≥ 1 ,Worpitzky 恒等式是
x n = ∑ k = 0 n − 1 A ( n , k ) ( x + n − 1 − k n ) = ∑ k = 0 n − 1 A ( n , k ) ( x + k n ) . 此处二项式系数 理路 二项式系数 Binomial coefficient n 元集合的 k 元子集数,记作 C(n,k)。 扩为多项式 ( y n ) = y ( y − 1 ) ⋯ ( y − n + 1 ) / n ! ;第二式由 A ( n , k ) = A ( n , n − 1 − k ) 重排而来。不要把两种方向中的 k 含义随意互换。
它与阶乘基的 x n = ∑ k S ( n , k ) x k ― 同为换基,但基不同:Worpitzky 使用不同平移的 n 次二项式多项式,而 Stirling 使用从零次到 n 次的下降阶乘。两张数表因此不应混用。
直觉
先取正整数 x = m 。左边数所有函数 f : [ n ] → [ m ] 。将标签 i 按 ( f ( i ) , i ) 的字典序排列:先按函数值从小到大;值相同时按标签从小到大。这就是固定了打平规则的稳定排序,得到唯一排列 π 。
沿这份排列,值 a i = f ( π i ) 满足
当 1 ≤ a 1 ≤ ⋯ ≤ a n ≤ m , a i < a i + 1 当 i ∈ Des ( π ) . 若标签出现下降,却让两个函数值相等,就违反了相等时标签递增的规则。反过来,任何满足这些条件的值序列都唯一恢复一个稳定排序为 π 的函数。
设 π 有 k 个下降,令 d i 为位置 i 之前的下降数,取 b i = a i − d i 。则严格增加处的强制一步被消去,得到
1 ≤ b 1 ≤ ⋯ ≤ b n ≤ m − k . 隔板法 理路 隔板法 Stars and bars 把相同对象分入有标号盒子的整数解计数方法。 给这样的弱递增序列数为 ( m − k + n − 1 n ) ;当 m ≤ k 时没有解,且此处上参数在 0 , … , n − 1 ,组合数也为零。按 k 汇总全部稳定排序就得第一式。再由双方多项式在无限多个正整数上一致,得到任意 x 的恒等式。
例子与边界
n = 3 , m = 2 时,
2 3 = ( 4 3 ) + 4 ( 3 3 ) + ( 2 3 ) = 4 + 4 + 0. 没有下降的 123 对应四份弱递增二值序列 111 , 112 , 122 , 222 。一个下降的四份排列各有唯一适合的二值序列。例如 π = 213 要求 a 1 < a 2 ≤ a 3 ,只能是 ( 1 , 2 , 2 ) ,恢复函数值 ( f ( 1 ) , f ( 2 ) , f ( 3 ) ) = ( 2 , 1 , 2 ) 。两次下降的 321 需要三个不同值,二值集合无法提供。
n = 4 , m = 2 时,16 = ( 5 4 ) + 11 ( 4 4 ) = 5 + 11 。这从函数的另一种分类独立核验 A ( 4 , 1 ) = 11 。
必须固定相等值的打平规则。如果同一函数值相等时任意排列标签,一份函数会落入多份 π ,导致重复计数。例如常值函数本应只落到 123 ⋯ n ,不能对它贡献 n ! 份。这里的函数值也有自然顺序;只给一个无序颜色集,须先选定颜色顺序再使用这份分解。
推论与应用
令 A n ( t ) = ∑ k A ( n , k ) t k 。对 n ≥ 1 将上式取整数 m ≥ 0 ,再作普通生成函数 理路 普通生成函数 Ordinary generating function 把序列编码为形式幂级数 Σ a_n x^n。 。由移位二项式级数,
∑ m ≥ 0 m n t m = t A n ( t ) ( 1 − t ) n + 1 . 例如 n = 3 为 t ( 1 + 4 t + t 2 ) / ( 1 − t ) 4 。乘回 ( 1 − t ) n + 1 并取 t k + 1 系数,得到
A ( n , k ) = ∑ j = 0 k + 1 ( − 1 ) j ( n + 1 j ) ( k + 1 − j ) n ( n ≥ 1 , 0 ≤ k < n ) . n = 4 , k = 1 时为 2 4 − 5 ⋅ 1 4 + 10 ⋅ 0 4 = 11 。
还可以完整推导Eulerian 的 EGF 理路 指数生成函数 Exponential generating function · EGF 以 a_n x^n/n! 编码带标号组合对象计数序列的形式幂级数。 。令 E ( z , t ) = ∑ n ≥ 0 A n ( t ) z n / n ! 。把上面的幂级数关系乘 z n / n ! 并对 n ≥ 1 求和,得
t 1 − t [ E ( z 1 − t , t ) − 1 ] = ∑ m ≥ 0 t m ( e m z − 1 ) = 1 1 − t e z − 1 1 − t . 整理、再令 z = ( 1 − t ) w ,得到
E ( w , t ) = 1 − t e ( t − 1 ) w − t . 上面的双和先在 Q [ [ t , z ] ] 中逐系数成立:固定 t m z n 时只对应一个 m 和一个指数项。每个 z n 系数再由前述幂和公式识别为 t 的有理函数,因而可转入 Q ( t ) [ [ z ] ] 作整理与缩放;最终 E 的每个 w n 系数仍是多项式。t = 1 处应按系数代入 A n ( 1 ) = n ! ,得到 E ( w , 1 ) = 1 / ( 1 − w ) ;不能把右边未经消去的 0 / 0 当数值答案。若要在复数点求值,还需另行检查解析收敛。
参考资料