RT,一直WA50分

开longlong、改快速幂写法等方法都试过了，但还是不对。。。

#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,d;
void mem(){
    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;
        }
    }
    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();
    d=a;
    b=a^(m-2);
    int ans=-b.mat[1][n+n];
    b=b*d;
    ans+=b.mat[1][n];
    printf("%d\n",ans);
    return 0;
}