矩阵乘法加速线性递推的典型
大概套路就是先构造一个矩阵\( F \)使得另一初始矩阵\( A \)乘以\( F^{x} \)能够得出第n项
跑的飞快
虽然我也不知道那个矩阵要怎么构造
或许就像我使用了瞪眼法和枚举法
#include <cstdio>
#include <algorithm>
#include <cstring>
#define int long long
using namespace std;
const int MOD = ;
int n;
struct Matrix{
static const int MAXN = ;
int alpha[MAXN][MAXN];
int n,m;
void init(void){
for(int i=;i<MAXN;i++)
for(int j=;j<MAXN;j++)
alpha[i][j]=;
n=m=;
}
void init_f2(void){
n=;m=;
alpha[][]=;
alpha[][]=;
alpha[][]=;
alpha[][]=;
}
void init_f(void){
n=;m=;
alpha[][]=;
alpha[][]=;
}
void init_pow(int x){
for(int i=;i<=x;i++)
alpha[i][i]=;
m=n=x;
}
Matrix operator * (Matrix b){
Matrix c;
c.init();
for(int i=;i<=n;i++)
for(int j=;j<=b.m;j++)
for(int k=;k<=m;k++)
c.alpha[i][j]=(c.alpha[i][j]%MOD+alpha[i][k]*b.alpha[k][j]%MOD)%MOD;
c.n=n;
c.m=b.m;
return c;
}
};
Matrix pow(Matrix a,int p){
Matrix ans;
ans.init();
ans.init_pow(a.n);
while(p){
if(p&)
ans=ans*a;
a=a*a;
p>>=;
}
return ans;
}
signed main(){
scanf("%lld",&n);
Matrix f,f2,ans;
f.init();
f.init_f();
f2.init();
f2.init_f2();
ans=pow(f2,n-);
// for(int i=1;i<=ans.n;i++){
// for(int j=1;j<=ans.m;j++)
// printf("%d ",ans.alpha[i][j]);
// printf("\n");
// }
f=f*ans;
printf("%lld",f.alpha[][]%MOD);
return ;
}