inline void prep(int len){
	B[0]=ifac[0]=ifac[1]=fac[0]=fac[1]=1;
	fp(i,2,len) ifac[i]=mul(mod-mod/i,ifac[mod%i]);
	fp(i,2,len) ifac[i]=mul(ifac[i-1],ifac[i]);
	fp(i,2,len) fac[i]=mul(fac[i-1],i);
	get_inv(ifac+1,B,len);
	fp(i,2,len) B[i]=mul(B[i],fac[i]);
}

经过一天的学习，我们发现伯努利数是个非常有用 （个屁） 的数列

定义

但是...伯努利数是什么呢？我们先给伯努利数一个定义：

令 B(i) 表示 伯努利数第 i 项，那么有：

当然 B0=1 ，于是上面的式子就可以一直推下去了...

暴力计算

我们按照上面的式子可以算出：

但是这样的效率嘛 【滑稽

当然是 n2 的啦~ 于是我们试图寻找更高效的办法...

高效计算

我们发现，根据定义有这么个等式：

这里原本 n 是 n+1 的，为了方便我们把 n 带进去...

然后这个式子有什么用呢？当然有用咯~ 这玩意儿可是能化成卷积的形式丫！

于是拆个组合数玩玩，看这里：

我们发现 0!=1

注意上面的 n 仍然不等于 1 ...

然后我们是不是发现这玩意儿左边其实是...卷积呢！1i! 不就是每项系数为 1 的多项式指数函数？

我们把伯努利数除去阶乘当做是一个多项式的系数，那么原来的式子就可以搞生成函数了！

但我们首先还要把 n=1 的情况加进去：

你问我为什么这么写？暴力算一下 n=1 的情况不就发现两边差了多少嘛！

然后转个多项式我们就可以开始生成函数了 QVQ

这里之所以加了 x 是因为上面的 x1 系数多了个 1

然后我们继续推式子：

也就是说，我们只要阶乘算出来整体左移一位然后求个逆就好了

但是这里要注意的是，我们这样求出的是指数型生成函数的伯努利数的多项式，真正的伯努利数列要在我们求出来这个多项式系数的基础上乘上对应的阶乘（也就是在指数型生成函数的多项式中除去的那个阶乘）

code

代码如下...丑的一比...而且还省去了求逆的具体函数...

inline void prep(int len){
	B[0]=ifac[0]=ifac[1]=fac[0]=fac[1]=1;
	fp(i,2,len) ifac[i]=mul(mod-mod/i,ifac[mod%i]);
	fp(i,2,len) ifac[i]=mul(ifac[i-1],ifac[i]);
	fp(i,2,len) fac[i]=mul(fac[i-1],i);
	get_inv(ifac+1,B,len);
	fp(i,2,len) B[i]=mul(B[i],fac[i]);
}

补充

伯努利数其实分两类： B1 这一项了

那么我们重新写一遍伯努利数列就是：

那么另一种不大常见的 B+ 有什么用么？或者说什么特殊的性质？

当然是有的咯，对于下面要说的 S(n,k) 是这么定义的：

那么令 B+ 的第 i 项，则有：

好像并没有什么软用...但是有时候要用到的话可能会让你的推出来的式子更加简洁？【雾

模板

然后这里是板子题：

板子题：51nod1228 序列求和

n2 暴力足以优雅水过~

板子题2：51nod1258 序列求和 V4

要用任意MTT哟~

题解点

应用

伯努利数有个非常良好 （个屁） 的性质：

令：

则：

即：

PS：关于上面等式的证明可以看

然后一系列推倒：

这可不就是卷积的形式嘛！

现在我们就可以在一些 (dl)数学题里面大展拳脚了呢~

比如这道：

转：https://www.cnblogs.com/Judge/p/10722777.html