题目链接: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;
}