题解 LibreOJ #10222. 「一本通 6.5 例 4」佳佳的 Fibonacci

这是一道恶心的数学题，令人谔谔！

---

大家都看过题面了吧？我现在来说说我的思路。

首先我在这篇博客中提出了 $S_{n-2}+1=F_n$ 即 $S_n=F_{n+2}-1$。

然后看题目给出的柿子。把他变成一个三角形：

$F_1$
$F_2{F}_2$
$F_3{F}_3{F}_3$
$F_4{F}_4{F}_4{F}_4$
$F_5{F}_5{F}_5{F}_5{F}_5$
$\text{…………}$

答案就是这个三角形内所有数之和。

我们可以把这个三角形竖着看，发现了啥？这个式子就是：

$$\sum _{i=1}^n(S_n-S_{i-1})$$

把 $S_n$ 提取出来：

$$n{S}_n-\sum _{i=1}^n{S}_{i-1}$$

这个时候，$S_n=F_{n+2}-1$ 派上用场了。把这个柿子往里面套啊！

$$n(F_{n+2}-1)-\sum _{i=1}^n(F_{i+1}-1)$$

发现了啥？$\sum _{i=1}^n(F_{i+1}-1)$ 其实也可以转化为 $S_{n+1}-F_1-n$。

把 $S_n=F_{n+2}-1$ 再往 $S_{n+1}$ 里面套！

答案就变成了：

$n(F_{n+2}-1)-(F_{n+3}-1-F_1-n)$
$=n{F}_{n+2}-n-F_{n+3}+1+F_1+n$
$=n{F}_{n+2}-F_{n+3}+2$

一个矩阵加速就搞定了。

恶心至极啊！

$\text{code:}$

#include<cstdio>
using namespace std;
long long n,m;
struct mat{
    long long a[5][5];
};
inline void clear(mat &z)
{
	for(int i=1;i<=2;i++)
	    for(int j=1;j<=2;j++) z.a[i][j]=0;
}
inline mat operator* (mat A,mat B)
{
	mat c;clear(c);
	for(int i=1;i<=2;i++)
	    for(int j=1;j<=2;j++)
	        for(int k=1;k<=2;k++)
	        c.a[i][j]=(c.a[i][j]+A.a[i][k]*B.a[k][j])%m;
	return c;
}
inline mat operator^ (mat A,long long b)
{
	mat c;clear(c);
	c.a[1][1]=c.a[1][2]=1;
	mat base=A;
	while(b)
	{
		if(b&1) c=c*base;
		base=base*base;
		b/=2;
	}
	return c;
}
int main()
{
	scanf("%lld%lld",&n,&m);
	if(n==1) {printf("%lld\n",1%m);return 0;}
	mat first;clear(first);
	first.a[1][1]=first.a[1][2]=first.a[2][1]=1;
	long long tmp1=n*((first^(n)).a[1][1])%m;
	long long tmp2=((first^(n+1)).a[1][1])%m;
	printf("%lld\n",(tmp1-tmp2+2+m)%m);
	return 0;
}