洛谷P4451 [国家集训队]整数的lqp拆分（生成函数）

最新推荐文章于 2022-10-09 19:08:00 发布

转载最新推荐文章于 2022-10-09 19:08:00 发布 · 103 阅读

CC 4.0 BY-SA版权

原文链接：http://www.cnblogs.com/bztMinamoto/p/10557621.html

博客围绕利用生成函数求解数列通项公式展开。设斐波那契数列生成函数为\(F(x)\)，答案生成函数为\(G(x)\)，通过递推关系得到\(G\)与\(F\)的关系，再对\(F(x)\)进行推导，最终代入求出\(G(x)\)，经因式分解、裂项后得出\(g_n\)的通项公式。

题面

传送门

题解

~~我对生成函数一无所知~~

我们设\(F(x)\)为斐波那契数列的生成函数，\(G(x)\)为答案的生成函数，那么容易得到递推关系

\[g_n=\sum_{i=0}^{n-1}f_ig_{n-i}+f_n\]

其中\(g_0=0,g_1=1\)

那么写成生成函数的形式就是

\[G=FG+F\]

\[G={F\over 1-F}\]

我们考虑一下\(F\)，因为

\[F(x)=\sum_{i=1}^\infty f_ix^i\]

\[xF(x)=\sum_{i=2}^\infty f_{i-1}x^i\]

上面的柿子减去下面的柿子

\[(1-x)F(x)=x+\sum_{i=2}^\infty f_{i-2}x^i\]

即

\[(1-x)F(x)=x+x^2F(x)\]

解得

\[F(x)={x\over 1-x-x^2}\]

代入可解出

\[G(x)={x\over 1-2x-x^2}\]

我们把\(1-2x-x^2\)因式分解一下

\[G(x)={x\over \left(1-(1+\sqrt{2})x\right)\left(1-(1-\sqrt{2})x\right)}\]

然后裂项

\[G(x)={1\over 2\sqrt{2}}{1\over \left(1-(1+\sqrt{2})x\right)}-{1\over 2\sqrt{2}}{1\over \left(1-(1-\sqrt{2})x\right)}\]

那么现在就变成两个等比数列求和的形式了，可以直接求出\(g_n\)的通项公式

\[g_n={(1+\sqrt{2})^n-(1-\sqrt{2})^n\over 2\sqrt{2}}\]

听说大佬们用生成函数只要五行就能写完题解……然而我并看不懂它们在写什么……

//minamoto
#include<bits/stdc++.h>
#define R register
#define fp(i,a,b) for(R int i=(a),I=(b)+1;i<I;++i)
#define fd(i,a,b) for(R int i=(a),I=(b)-1;i>I;--i)
#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)
using namespace std;
const int P=1e9+7,s=59713600;
inline int add(R int x,R int y){return x+y>=P?x+y-P:x+y;}
inline int dec(R int x,R int y){return x-y<0?x-y+P:x-y;}
inline int mul(R int x,R int y){return 1ll*x*y-1ll*x*y/P*P;}
int ksm(R int x,R int y){
    R int res=1;
    for(;y;y>>=1,x=mul(x,x))if(y&1)res=mul(res,x);
    return res;
}
int main(){
//  freopen("testdata.in","r",stdin);
    int n;scanf("%d",&n);
    printf("%d\n",mul(dec(ksm(s+1,n),ksm(P+1-s,n)),ksm(mul(s,2),P-2)));
    return 0;
}

转载于:https://www.cnblogs.com/bztMinamoto/p/10557621.html