• 斯特林数学习笔记


    正文开始前,感谢以下几篇博客:

    https://www.cnblogs.com/y2823774827y/p/10700231.html(不认识的)

    https://www.cnblogs.com/p-b-p-b/p/10943453.html(PB的)

    https://www.cnblogs.com/AThousandMoons/p/10949178.html(CDW的)


    第一类斯特林数的定义:$egin{bmatrix}n\mend{bmatrix}$ 表示把 $n$ 个不同的物品组成恰好 $m$ 个没有区别的圆排列的方案数。

    那么有递推公式 $egin{bmatrix}n\mend{bmatrix}=(n-1)egin{bmatrix}n-1\mend{bmatrix}+egin{bmatrix}n-1\m-1end{bmatrix}$。前面是把新的物品塞在一个原来的物品后边,后面是新开一个圆排列。

    初始 $egin{bmatrix}1\1end{bmatrix}=1,egin{bmatrix}n\mend{bmatrix}=0(n<m)$。


    这毒瘤玩意居然没有单个的通项公式

    有个式子:$n^{overline{m}}=sumlimits^{m}_{i=0}egin{bmatrix}m\iend{bmatrix}n^i$。

    组合意义?不会

    理性证明?看PB的去(算了还是搬过来吧)

    $$n^{overline{m+1}}$$

    $$(n+m)n^{overline{m}}$$

    $$(n+m)sumlimits^{m}_{i=0}egin{bmatrix}m\iend{bmatrix}n^i$$

    $$sumlimits^{m+1}_{i=0}(megin{bmatrix}m\iend{bmatrix}+egin{bmatrix}m\i-1end{bmatrix})n^i$$

    $$sumlimits^{m+1}_{i=0}egin{bmatrix}m+1\iend{bmatrix}n^i$$


    求一行的:$x^{overline{n}}=sumlimits^{n}_{i=0}egin{bmatrix}n\iend{bmatrix}x^i$。

    所以只要求出 $x^{overline{n}}$ 这个多项式即可。

    设 $F_n(x)=x^{overline{n}}$。

    当 $n$ 是奇数时 $F_n(x)=(x+n-1)F_{n-1}(x)$,不管。

    当 $n$ 是偶数时 $F_n(x)=F_{frac{n}{2}}(x)F_{frac{n}{2}}(x+frac{n}{2})$。

    (为方便将式子变为 $F_{2n}(x)=F_n(x)F_n(x+n)$)

    假设已经求出 $F_n(x)=sumlimits_{i=0}^na_ix^i$。

    $$F_n(x+n)$$

    $$sumlimits_{i=0}^na_i(x+n)^i$$

    $$sumlimits_{i=0}^na_isumlimits_{j=0}^idbinom{i}{j}x^jn^{i-j}$$

    $$sumlimits_{i=0}^na_ii!sumlimits_{j=0}^idfrac{x^j}{j!}dfrac{n^{i-j}}{(i-j)!}$$

    $$sumlimits_{j=0}^ndfrac{x^j}{j!}sumlimits_{i=j}^n(a_ii!)dfrac{n^{i-j}}{(i-j)!}$$

    令 $A_i=a_ii!,B_i=dfrac{n^{n-i}}{(n-i)!}$,那么 $x^j$ 的系数是 $dfrac{A imes B}{j!}[x^{j+n}]$。

    时间复杂度 $O(nlog n)$。

    求一列的:

    考虑生成函数 $F(x)=sumlimits_{i=1}^{+infty}(i-1)!dfrac{x^i}{i!}=sumlimits_{i=1}^{+infty}dfrac{x^i}{i}$(注意下界是 $1$)。表示把 $n$ 个球形成圆排列(非空)的指数型生成函数(斯特林数中集合没有区别,所以要指数型)。那么 $dfrac{F^m(x)}{m!}$ 就是把 $n$ 个球形成 $m$ 个圆排列且每个圆排列非空的方案数的指数型生成函数。那么 $dfrac{egin{bmatrix}i\mend{bmatrix}}{i!}=dfrac{F^m(x)}{m!}[x^i]$。$O(nlog n)$。


    第二类斯特林数的定义:$egin{Bmatrix}n\mend{Bmatrix}$ 表示把 $n$ 个不同的物品放进 $m$ 个没有区别的集合,且每个集合都非空的方案数。

    那么有递推公式 $egin{Bmatrix}n\mend{Bmatrix}=megin{Bmatrix}n-1\mend{Bmatrix}+egin{Bmatrix}n-1\m-1end{Bmatrix}$。前面是把新的物品随便找个集合塞下,后面是新开一个集合。

    初始 $egin{Bmatrix}1\1end{Bmatrix}=1,egin{Bmatrix}n\mend{Bmatrix}=0(n<m)$。


    有个式子:$n^m=sumlimits^{m}_{i=0}egin{Bmatrix}m\iend{Bmatrix}dbinom{n}{i}i!=sumlimits^m_{i=0}egin{Bmatrix}m\iend{Bmatrix}n^{underline{i}}$。

    组合意义理解一下,左边是 $m$ 个不同物品扔进 $n$ 个不同集合,且可以有空集合的方案数。右边枚举有多少个非空集合,然后把 $m$ 个物品放进这些集合里。注意第二类斯特林数中集合没有区别,所以要乘回一个 $i!$ 使他们有区别。

    二项式反演后可得通项公式 $egin{Bmatrix}n\mend{Bmatrix}=dfrac{1}{m!}sumlimits^m_{k=0}(-1)^kdbinom{m}{k}(m-k)^n$。


    求一行的:

    发现 $egin{Bmatrix}n\iend{Bmatrix}=sumlimits^i_{k=0}dfrac{(-1)^k}{k!}dfrac{(i-k)^n}{(i-k)!}$,是个卷积的形式,$O(nlog n)$ 求解。

    求一列的:

    考虑生成函数 $F(x)=sumlimits_{i=1}^{+infty}dfrac{x^i}{i!}=e^x-1$(注意下界是 $1$)。表示把 $n$ 个球丢入一个集合(非空)的指数型生成函数(斯特林数中集合没有区别,所以要指数型)。类似上面的做法即可。$O(nlog n)$。

    还有一种方法。考虑 $F_m(x)=sumlimits_{i=0}^{+infty}egin{Bmatrix}i\mend{Bmatrix}x^i$。那么显然有 $F_0(x)=1$。

    $$F_m(x)=sumlimits_{i=0}^{+infty}egin{Bmatrix}i\mend{Bmatrix}x^i$$

    $$F_m(x)=sumlimits_{i=0}^{+infty}(megin{Bmatrix}i-1\mend{Bmatrix}+egin{Bmatrix}i-1\m-1end{Bmatrix})x^i$$

    $$F_m(x)=msumlimits_{i=0}^{+infty}egin{Bmatrix}i-1\mend{Bmatrix}x^i+sumlimits_{i=0}^{+infty}egin{Bmatrix}i-1\m-1end{Bmatrix}x^i$$

    $$F_m(x)=msumlimits_{i=0}^{+infty}egin{Bmatrix}i\mend{Bmatrix}x^{i+1}+sumlimits_{i=0}^{+infty}egin{Bmatrix}i\m-1end{Bmatrix}x^{i+1}$$

    $$F_m(x)=mxF_m(x)+xF_{m-1}(x)$$

    $$F_m(x)=dfrac{xF_{m-1}(x)}{1-mx}$$

    $$F_m(x)=dfrac{x^m}{prodlimits_{i=1}^m(1-ix)}$$

    分母实际上是 $prodlimits_{i=1}^m(x-i)$ 的系数倒置。可以用上面的方法。也可以 $O(nlog n)$,常数更小。


    首先由两个很显然的式子:

    $n^{overline{m}}=(-1)^{n+1}(-n)^{underline{m}}$

    $n^{underline{m}}=(-1)^{n+1}(-n)^{overline{m}}$

    总结一下上面提到的性质:(顺便引入一些新的)

    $n!=sumlimits^n_{i=0}egin{bmatrix}n\iend{bmatrix}$(其实就是下面的式子令 $n=1$,$m=n$,也可以用抽象代数的意义证)

    $n^{overline{m}}=sumlimits^m_{i=0}egin{bmatrix}m\iend{bmatrix}n^i$

    $n^{underline{m}}=sumlimits^m_{i=0}(-1)^{m-i}egin{bmatrix}m\iend{bmatrix}n^i$(也可以用数学归纳法证)

    $n^m=sumlimits^m_{i=0}egin{Bmatrix}m\iend{Bmatrix}n^{underline{i}}=sumlimits^m_{i=0}egin{Bmatrix}m\iend{Bmatrix}dbinom{n}{i}i!$

    $n^m=sumlimits^m_{i=0}(-1)^{m-i}egin{Bmatrix}m\iend{Bmatrix}n^{overline{i}}$(下面会讲到怎么证)

    $sumlimits^n_{i=0}i^k=frac{(n+1)^{underline{k+1}}}{k+1}-sumlimits_{i=0}^{k-1}(-1)^{k-i}egin{bmatrix}k\iend{bmatrix}sumlimits_{j=0}^nj^i$

    $sumlimits^n_{i=0}i^k=sumlimits^k_{i=0}egin{Bmatrix}k\iend{Bmatrix}dfrac{(n+1)^{underline{i+1}}}{i+1}$

    关于最后一个式子的证明:(倒数第二个自己来咕了

    $$sumlimits^n_{i=0}sumlimits^k_{j=0}egin{Bmatrix}k\jend{Bmatrix}dbinom{i}{j}j!$$

    $$sumlimits^k_{j=0}j!egin{Bmatrix}k\jend{Bmatrix}sumlimits^n_{i=0}dbinom{i}{j}$$

    $$sumlimits^k_{j=0}j!egin{Bmatrix}k\jend{Bmatrix}dbinom{n+1}{j+1}$$

    $$sumlimits^k_{j=0}j!egin{Bmatrix}k\jend{Bmatrix}dfrac{(n+1)^{underline{j+1}}}{j+1}$$

    中间 $sumlimits^n_{i=0}dbinom{i}{j}=dbinom{n+1}{j+1}$ 就是传说中的杨辉恒等式。右边是从 $n+1$ 个位置中选 $j+1$ 个位置,左边是枚举最右边的位置($i+1$),然后从前 $i$ 个位置中选剩下的 $j$ 个。

    有些题目中会用到普通幂和上升/下降幂转换的方法来做。根据PB的总结,当普通幂转上升/下降幂时用第二类,当上升/下降幂转普通幂时用第一类。


    斯特林反演

    长得和二项式反演很像。

    $$f_n=sumlimits^n_{k=0}(-1)^kegin{bmatrix}n\kend{bmatrix}g_kiff g_n=sumlimits^n_{k=0}(-1)^kegin{Bmatrix}n\kend{Bmatrix}f_k$$

    $$f_n=sumlimits^n_{k=0}egin{bmatrix}n\kend{bmatrix}g_kiff g_n=sumlimits^n_{k=0}(-1)^{n-k}egin{Bmatrix}n\kend{Bmatrix}f_k$$

    用这个式子也可以证明上面的一些式子。

    证明的话需要两个反转公式:

    $$sum_{i=m}^n(-1)^{n-i}egin{bmatrix}n\iend{bmatrix}egin{Bmatrix}i\mend{Bmatrix}=[n=m]$$
    $$sum_{i=m}^n(-1)^{n-i}egin{Bmatrix}n\iend{Bmatrix}egin{bmatrix}i\mend{bmatrix}=[n=m]$$


    两个奇怪的柿子:

    $$egin{Bmatrix}n+1\m+1end{Bmatrix}=sum_{k}inom{n}{k}egin{Bmatrix}k\mend{Bmatrix}$$

    证明:后面那个相当于从 $n$ 个物品中选出一些,并分成 $m$ 个无序非空集合。

    前面那个是把 $n+1$ 个物品分成 $m+1$ 个无序非空集合。

    钦定这 $n+1$ 个物品中有一个关键物品,分成的 $m+1$ 个集合中恰好一个含有关键物品,删掉这个集合,显然和后半部分等价。

    $$egin{bmatrix}n+1\m+1end{bmatrix}=sum_{k}egin{bmatrix}n\kend{bmatrix}inom{k}{m}$$

    证明:后面那个相当于把 $n$ 个物品分成若干个圆排列,再选出 $m$ 个圆排列。

    前面那个是把 $n+1$ 个物品分成 $m+1$ 个圆排列。

    钦定这 $n+1$ 个物品中有一个关键物品,分成的 $m+1$ 个圆排列中恰好一个含有关键物品,将包含关键物品的那个环从关键物品那里断开。

    断开之后会变成一个排列,注意到排列是可以和大小之和一定的圆排列集合一一对应的。

    将这些圆排列和原来没被断掉的 $m$ 个圆排列合在一起看,显然和后半部分等价。


    暂时先咕到这里。

  • 相关阅读:
    JS转义 escape()、encodeURI()、encodeURIComponent()区别详解
    PHP解决搜索时在URL地址栏输入中文字符搜索结果出现乱码
    CMSPRESS-PHP无限级分类2
    CMSPRESS-PHP无限级分类
    HTML5-CSS3-JavaScript(3)
    HTML5-CSS3-JavaScript(2)
    HTML5-CSS3-JavaScript(1)
    CSS3-Hover 效果 展示
    JAVA Socket地址绑定
    JAVA Socket无参构造方法的使用
  • 原文地址:https://www.cnblogs.com/1000Suns/p/11011492.html
Copyright © 2020-2023  润新知