五边形数 笔记

博客观赏效果更佳

近日,一国外小哥学习了这个定理,竟然能预处理出整数划分的方案数!快跟小编来看看吧

这个小哥学习的定理,就是小编(我)接下来要讲的五边形数定理

好的,让我们一起来看看这个定理吧

五边形数是啥

百度百科

用图来讲,就是若干个点,排成若干个五边形,需要多少个点。

百度百科上有一个很清楚的图:

它的通项公式为 \(a_n=\dfrac{n(3n-1)}{2}\),前面几项是 \(0, 1, 5, 12, 22, 35, 51, 70, 92, 117, 145, 176, 210...\)。你可以在 OEIS 上找到它,是序列 A000326

我们要研究的是 广义五边形数,是这样一个数列:\(p=a_0,a_{1},a_{-1},a_{2},a_{-2}...\)

其中,如果 \(p\) 的下标是负数,照样代入到上面的通项公式里算即可。容易证明,当 \(n<0\) 时,\(p_n\) 也是正的。前面几项是 \(0, 1, 2, 5, 7, 12, 15, 22, 26, 35, 40, 51, 57, 70, 77, 92, 100, 117...\),同样可以在 OEIS 上找到它是 A001318

它能干啥

欧拉函数(五边形数)

(为了和欧拉数论函数区别开来,这里给它起名叫欧拉函数-五边形数(*/ω\*)

是这样一个函数:\(\varphi(x)=(1-x)(1-x^2)(1-x^3)(1-x^4)...(1-x^\infin)=\prod\limits_{i=1}^{\infin} (1-x^i)\)

这东西有一个很神奇的性质。如果你闲的没事,可以暴力拆开它,发现它等于

\(x^0-x^1-x^2+x^5+x^7-x^{12}-x^{15}+x^{22}+x^{26}...=\sum\limits_{i=1}^{\infin} (-1)^{(i+1)/2} x^{p_i}\)

系数:一个正,两个负,两个正,两个负,两个正,两个负...

指数:广义五边形数

与整数划分问题的联系

然后我们可以来解决整数划分问题。设 \(P(n)\) 表示将 \(n\) 分成若干个整数的无序方案(即,\(a+b\)\(b+a\) 只算一次),同一个数可以用很多次。

比如,\(P(4)=5\),因为 \(4\)

\(=1+1+1+1\)

\(=1+1+2\)

\(=1+3\)

\(=2+2\)

\(=4\)

特殊地,\(P(0)=1\)

我们求出 \(P(n)\) 的生成函数 \(f(x)=\sum\limits_{n=0}^{\infin} P(n)x^n\)

给它换个定义:我们要选出若干个(可以是 \(0\) 个) \(1\),若干个 \(2\),若干个 \(3\)...使得选出来的所有数和为 \(n\)

这样就好考虑了。看这个式子:

\((1+x+x^2+x^3...)(1+x^2+(x^2)^2+(x^2)^3...)(1+x^3+(x^3)^2+(x^3)^3...)...\)

第一个因式表示我们能随便选择用多少个 \(x^1\)

第二个因式表示我们能随便选择用多少个 \(x^2\)

...

乘起来之后,考虑 \(x^n\) 的系数,你会发现它就是 \(P(n)\)!是不是特别巧妙φ(>ω<*)

于是 \(f(x)=(1+x+x^2+x^3...)(1+x^2+(x^2)^2+(x^2)^3...)(1+x^3+(x^3)^2+(x^3)^3...)...=\prod\limits_{n=1}^{\infin} \sum\limits_{i=0}^{\infin} (x^n)^i\)

然后我们用等比数列求和公式化一下,得到:

\(f(x)=\prod\limits_{n=1}^{\infin} \dfrac{1}{1-x^n}\)

我们发现它就是欧拉函数-五边形数的倒数(/≧▽≦/)!

\(\varphi(x)f(x)=1\)

也就是

\((1-x-x^2+x^5+x^7-x^{12}-x^{15}...)(\sum\limits_{n=0}^{\infin} P(n)x^n)=1\)

观察其中 \(x^n\) 项的系数(分别考虑左边和右边的括号出的是几次项),容易发现这个系数应该是 \(P(n)-P(n-1)-P(n-2)+P(n-5)+P(n-7)...\)

\(n=0\) 时,这个式子显然为 \(1\)。那么 \(n>0\) 时, \(x^n\) 乘以它的系数的和就得是 \(0\) 了。

然后我们发现等式右边是 \(1\)。并且原式为恒等式,所以 \(x\) 取任何值都成立。为了让 \(x\) 取任何值时等式都成立,对于 \(n>0\) 时,\(P(n)-P(n-1)-P(n-2)\) 这个可怜的等式只好取 \(0\) 了。

然后就得到了:

\(P(n)-P(n-1)-P(n-2)+P(n-5)+P(n-7)-P(n-12)-P(n-15)...=0\)

于是 \(P(n)=P(n-1)+P(n-2)-P(n-5)-P(n-7)+P(n-12)+P(n-15)...\)

然后 \(\le n\) 的五边形数是 \(\sqrt{n}\) 级别的(注意到五边形数的通项公式是二次的)

于是我们可以递推出 \(P(n)\),是 \(O(n\sqrt{n})\) 的,是不是很厉害呢(✪ω✪)

好了,以上就是这个定理的全部内容了,喜欢记得收藏起来哟

猜你喜欢

转载自www.cnblogs.com/LightningUZ/p/12820073.html
今日推荐