牛客挑战赛-----斐波那契数列卷积（矩阵快速幂）

题目链接：https://ac.nowcoder.com/acm/problem/52937

题目大意：太简单了就不解释了。

题解：

先暴力打表：

找出规律来（这一步超难啊啊啊）：

d（n）也就是a（n）。

然后就推出矩阵：

我们要求的答案就是：

直接套矩阵快速幂的板子就可以了。

ac代码：

#include<bits/stdc++.h>
#define Max 10
using namespace std;
typedef long long ll;
const int mod=998244353;
ll n;
ll f[100];
const int N=4;
struct Matrix{///矩阵结构体 
    ll matrix[N][N];
};
void init(Matrix &res)///初始化为A矩阵 
{
    memset(res.matrix,0,sizeof(res.matrix));
    res.matrix[0][0]=2;res.matrix[0][1]=1;
	res.matrix[1][0]=1;res.matrix[1][2]=1;
	res.matrix[2][0]=-2;res.matrix[2][3]=1;
	res.matrix[3][0]=-1;
}
Matrix multiplicative(Matrix a,Matrix b)///矩阵乘法
{
    Matrix res;
    memset(res.matrix,0,sizeof(res.matrix));
    for(int i = 0 ; i < N ; i++){
        for(int j = 0 ; j < N ; j++){
            for(int k = 0 ; k < N ; k++){
                res.matrix[i][j] += a.matrix[i][k]*b.matrix[k][j];
                res.matrix[i][j]=(res.matrix[i][j]+mod)%mod;
            }
        }
    }
    return res;
}
Matrix Pow(Matrix mx,ll m)///矩阵快速幂 
{
    Matrix ans=mx;
    while(m)
    {
        if(m&1)
            ans=multiplicative(ans,mx);
        mx=multiplicative(mx,mx);
        m>>=1;
    }
    return ans;
}

int main(){
	Matrix A,B;
	init(A);
    f[0]=0;
    f[1]=1;
    f[2]=1;
    f[3]=2;
    /*打表找规律：an=2*a(n-1)+a(n-2)-2*a(n-3)-a(n-4) 
    for(int i=2;i<100;i++){
    	f[i]=f[i-1]+f[i-2];
	}
	for(n=1;n<=50;n++){
		ll sum=0; 
    	for(int k=0;k<=n;k++){
		//cout<<"F"<<n-k<<" x F"<<k<<": "<<f[n-k]<<" * "<<f[k]<<endl;
		sum+=(f[n-k]*f[k]);
    	}
	    cout<<"a"<<n<<": "<<sum%mod<<"    f"<<n<<": "<<f[n]<<endl;
	}
	*/
	scanf("%lld",&n);
	if(n==1){
		printf("0\n");
	}
	else if(n==2){
		printf("1\n");
	}
	else if(n==3){
		printf("2\n");
	}
	else{
		B=Pow(A,n-3);
	//	ll ans=0;
	//	ans=ans+(f[0]*map.matrix[0][0]%mod+f[1]*map.matrix[0][1]%mod+f[2]*map.matrix[0][2]%mod+f[3]*map.matrix[0][3]%mod)%mod; 
        printf("%lld\n",B.matrix[0][0]);	
	}
	return 0;
}