luogu SP31428 FIBONOMIAL - Fibonacci Polynomial

背景：

偶然发现$\text{luogu rank1}$巨佬推荐这道题，就进来做做。
大水题。

题目传送门：

https://www.luogu.org/problemnew/show/SP31428

题意：

给定$n,x$，求$D(n,x)={\sum }_{i=1}^n{f}_i{x}^i$。其中$f$为$\text{Fibonacci}$序列，$f_1=1,f_2=1,f_n=f_{n-1}+f_{n-2}$。

思路：

这种数据范围显然要用矩阵乘法优化或考虑通项式。
我们在生成函数——OGF证明了$\text{Fibonacci}$序列的通项式为：$\frac{1}{\sqrt{5}}\left((\frac{1+\sqrt{5}}{2}{)}^n-(\frac{1-\sqrt{5}}{2}{)}^n\right)$，也就是说：
$$D(n,x)=\sum _{i=1}^n\frac{1}{\sqrt{5}}\left((\frac{1+\sqrt{5}}{2}{)}^i-(\frac{1-\sqrt{5}}{2}{)}^i\right)x^i$$

$$D(n,x)=\frac{1}{\sqrt{5}}\sum _{i=1}^n\left((\frac{1+\sqrt{5}}{2}x{)}^i-(\frac{1-\sqrt{5}}{2}x{)}^i\right)$$

$$D(n,x)=\frac{\sqrt{5}}{5}\left(\sum _{i=1}^n(\frac{x+\sqrt{5}x}{2}{)}^i-\sum _{i=1}^n(\frac{x-\sqrt{5}x}{2}{)}^i\right)$$

注意到两个$\sum$其实就是等比数列求和，剩下的等比数列求和即可。

如何解决$\sqrt{5}$呢？
考虑到$\sqrt{5}$在$\mathrm{m}\mathrm{o}\mathrm{d}(10^9+7)$下没有二次剩余，因此我们可以扩域。
简单来说就是令$z=x+yi$，其中$i=\sqrt{5}$的形式。
然后类似复数的计算，有：
$z1\pm z2=(x1+y1i)\pm (x2+y2i)=(x1\pm x2)+(y1\pm y2)i$；
$z1\cdot z2=(x1+y1i)(x2+y2i)=(x1x2)+(x1y2+x2y1)i+y1y2i^2=(x1x2+5y1y2)+(x1y2+x2y1)i$；
$\frac{1}{z}=\frac{1}{x+yi}=\frac{x-yi}{(x+yi)(x-yi)}=\frac{x}{x^2-5y^2}-\frac{y}{x^2-5y^2}i$

代码：

#include<cstdio>
#include<cstring>
#include<algorithm>
#define LL long long
#define mod 1000000007
using namespace std;
struct comp
{
	LL x,y;
	comp(LL xx=0,LL yy=0):x(xx),y(yy) {}
	friend comp operator+(const comp &a,const comp &b) {return comp((a.x+b.x)%mod,(a.y+b.y)%mod);}
	friend comp operator-(const comp &a,const comp &b) {return comp((a.x-b.x+mod)%mod,(a.y-b.y+mod)%mod);}
	friend comp operator*(const comp &a,const comp &b) {return comp((a.x*b.x%mod+5ll*a.y*b.y%mod)%mod,(a.x*b.y%mod+b.x*a.y%mod)%mod);}
};
LL dg1(LL x,LL k)
{
	if(!k) return 1;
	LL op=dg1(x,k>>1);
	return (k&1)?op*op%mod*x%mod:op*op%mod;
}
LL inv1(LL x)
{
	return dg1(x,mod-2);
}
comp dg2(comp x,LL k)
{
	if(!k) return comp(1,0);
	comp op=dg2(x,k>>1);
	return (k&1)?op*op*x:op*op;
}
comp inv2(comp now)
{
	LL x=now.x,y=now.y;
	LL t=dg1((x*x%mod-5ll*y*y%mod+mod)%mod,mod-2);
    return comp(x*t%mod,(-y*t%mod+mod)%mod);
}
comp work(LL n,comp now)
{
	return (dg2(now,n+1)-now)*inv2(now-comp(1,0));
}
int main()
{
	int T;
	LL n,x;
	LL inv2=inv1(2),inv5=inv1(5);
	scanf("%d",&T);
	while(T--)
	{
		scanf("%lld %lld",&n,&x);
		comp a=comp(inv2*x%mod,inv2*x%mod),b=comp(inv2*x%mod,(-inv2+mod)%mod*x%mod);
		printf("%lld\n",(comp(0,inv5)*(work(n,a)-work(n,b))).x);
	}
}