50分求助

#include<cstdio>
#include<cstring>
#define N 110
#define mo 30011
using namespace std;
int n,m;
struct ak{int mat[N][N];ak(){memset(mat,0,sizeof(mat));}}a,b;
void mem(ak &a){
	a.mat[1][2]=a.mat[1][1]=1;
	for(int i=2;i<n;i++)a.mat[i][i+1]=a.mat[i][i]=a.mat[i][i-1]=1;
	a.mat[n][n-1]=a.mat[n][n]=1;
	for(int i=1;i<=n;i++)a.mat[i][i+n]=a.mat[i+n][i]=1;
	return;
}
ak operator * (ak a,ak b){
	ak c;
	for(int i=1;i<=n*2;i++){
		for(int j=1;j<=n*2;j++){
			for(int k=1;k<=n*2;k++)c.mat[i][j]=(c.mat[i][j]+a.mat[i][k]*b.mat[k][j]%mo)%mo;
		}
	}
	return c;
}
ak operator ^ (ak a,int b){
	ak ans;
	for(int i=1;i<=n*2;i++)ans.mat[i][i]=1;
	while(b){
		if(b&1)ans=ans*a;
		a=a*a;b>>=1;
	}
	return ans;
}
int main(){
	scanf("%d%d",&n,&m);
	mem(a);
	b=a^(m-2);
	int ans=-b.mat[1][n*2];
	b=b*a;
	ans+=b.mat[1][n];
	printf("%d\n",ans);
	return 0;
}