斐波那契数列卷积（矩阵快速幂+费马小定理）

这个题的名字好高大上哦(:
我打了一下表，然后OEIS到了这个数列；可以用斐波那契数列求的；
自己推还是不太现实，我感觉所以套结论了，嘻嘻；
这个通项为：

其中F(n)是斐波那契数列，所以很明显n为1e18所以上矩阵快速幂；求出分子，然后最后用费马小定理；
什么是费马小定理呢？

就是这个东东，然后就是怎么求F（n）呢？
很简单，可以构造一个矩阵：

然后手写快速幂；
最后取Mod；然后输出答案；
AC代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
ll a[2][2],one[2][2];
ll Mod=998244353;
void Multiply(ll x[][2],ll y[][2])//数组传的地址，所以会改变实参数组
{//乘法
	ll t[2][2];memset(t,0,sizeof(t));
	for(int i=0;i<2;i++){
		  for(int j=0;j<2;j++){
		  	   for(int k=0;k<2;k++){
		  	   	    t[i][j]=(t[i][j]+x[i][k]*y[k][j])%Mod;
				 }
		  }
	}
	for(int i=0;i<2;i++){
		  for(int j=0;j<2;j++){
		  	   x[i][j]=t[i][j];
		  }
	}
}
void Qp(ll x[2][2],ll n){//数组传的地址，所以会改变实参数组
	  while(n){
	  	  if(n&1)Multiply(one,x);
	  	  Multiply(x,x);
	  	  n>>=1;
	  }
}
void init(){//初始化矩阵
	one[0][0]=one[1][1]=1;
	one[0][1]=one[1][0]=0;
     a[0][0]=a[0][1]=1;
     a[1][0]=1;a[1][1]=0;
}
ll QP(ll x,ll n){//指数快速幂
	  ll res=1;
	  while(n){
	  	   if(n&1){
	  	   	  res=(res*x)%Mod;
			 }
		   x=(x*x)%Mod;
			 n>>=1;
	  }
	  return res;
}
int main(){
  ll n;
  scanf("%lld",&n);
  if(n==1)return cout<<0<<endl,0;//n==1时，需要特判一下
   init();
   Qp(a,n-1);
   ll f1=one[0][0];
   //cout<<f1<<endl;
   init();
   Qp(a,n-2);
   ll f2=one[0][0];
   //cout<<f2<<endl;
   ll fenzi=(n-1)%Mod*f1%Mod+2*n%Mod*f2%Mod;
   cout<<fenzi%Mod*QP(5,Mod-2)%Mod<<endl; //费马小定理
	return 0;
}