CodeChef AUG14 SIGFIB

　　来自2016集训队作业。
　　CC毒瘤题真是多。。。

---

　　求$\sum_{i+j+k=n}ijkf_if_jf_k\mod m$，其中$f_i$是斐波那契数列。$n\leq 10^{18},m\leq 10^5$。

---

　　答案的生成函数是$(zD\frac{z}{1-z-z^2})^3$，展开后可知答案是$[z^n](z+1/z)^3F(z)^6$，其中$F(z)=\frac{z}{1-z-z^2}$，是$f_i$的生成函数。
　　我们考虑怎么展开$[z^n]F(z)^6$得到其关于$f_n$的表达式。
　　设$\phi,\hat \phi$分别为$z^2-z-1=0$大于0和小于0的两个根。
　　我们知道$F(z)=\frac{1}{\sqrt 5}(\frac{1}{1-\phi z}- \frac{1}{1-\hat \phi z})$，先来试试算个$F(z)^2$。
　　

$$\begin{align*} 　　5F(z)^2&= (\frac{1}{1-\phi z})^2-2F(z)/z+(\frac{1}{1-\hat \phi z})^2 \\ 　　&= \sum_n \binom{n+1}{1}(\phi^n+\hat\phi^n) z^n-2F(z)/z \\ 　　\end{align*}$$
　　注意到
　　 $$\begin{align*} 　　\phi^n+\hat \phi^n&= [z^n](\frac{1}{1-\phi z}+\frac1{1-\hat \phi z})\\ 　　&=[z^n]\frac{2-(\phi+ \hat \phi)z}{(1-\phi z)(1-\hat \phi z)} \\ 　　&=[z^n](2F(z)/z-F(z))\\ 　　&=2f_{n+1}-f_n 　　\end{align*}$$
　　因此 $[z^n]F(z)^2=(2nf_{n+1}-(n+1)f_n)/5$。我们可以用类似的方法算出 $\sum_k kf_kf_{n-k}$， 然后我们算 $F(z)^6$的时候可以先将其二项式展开然后合并成两部分的和，会得到一系列的各种形式的 $F(z)^3$，接着再用 $F(z)^2$和 $\sum_k kf_kf_{n-k}$凑出来。。。
　　于是我们最后得到的大概会是 $[z^n]F(z)^6=\sum_{k=0}^6 C_{n,k} f_{n+k}$的形式，其中 $C_{n,k}$里可能有关于 $\binom{n+O(k)}{O(k)}$的组合数之类。因为下指标很小所以组合数可以暴力算，剩下的 $f_n$直接矩乘算。
　　要注意会有一个 $5^6$的分母，我们把 m <script type="math/tex" id="MathJax-Element-81">m</script>乘上它之后再做，最后再除掉他就好。
　　 （以上都是未经实践的口胡）
　　 其实一开始大力MMA就可以得到系数了
　　用有理生成函数的一般展开定理也可以用来解这个系数，但是得到展开里的多项式也要一点技巧（虽然可以暴力），我打算以后用这个技巧 出题干点奇怪的事所以就不写在这里了。