#include using namespace std; #define rep(i,n) for(int i = 0;i<((int)(n));i++) #define reg(i,a,b) for(int i = ((int)(a));i<=((int)(b));i++) #define irep(i,n) for(int i = ((int)(n)-1);i>=0;i--) #define ireg(i,a,b) for(int i = ((int)(b));i>=((int)(a));i--) typedef long long ll; /* */ const int MAT=125; struct Mat { ll v[MAT][MAT]; }; Mat mulmat(Mat& a,Mat& b,int n=MAT,ll mod=998244353) { Mat r; rep(x,n) rep(y,n) r.v[x][y]=0; rep(x,n) rep(z,n) rep(y,n) { r.v[x][y] += a.v[x][z]*b.v[z][y]; r.v[x][y] %= mod; } rep(x,n) rep(y,n) r.v[x][y]%=mod; return r; } Mat powmat(ll p,Mat a,int n=MAT) { Mat r; rep(x,n) rep(y,n) r.v[x][y]=0; rep(i,n) r.v[i][i]=1; while(p) { if(p%2) r=mulmat(r,a,n); a=mulmat(a,a,n); p>>=1; } return r; } ll n,k,ans=0; ll mod=998244353; Mat mt; void init(){ cin>>n>>k; } int main(void){ init(); rep(i,k*k*k)rep(j,k*k*k)mt.v[i][j]=0; rep(a,k)rep(b,k)rep(c,k){ ll from = a*k*k+b*k+c; mt.v[((a+1)%k)*k*k+b*k+c][from]++; mt.v[a*k*k+((b+a)%k)*k+c][from]++; mt.v[a*k*k+b*k+((c+b)%k)][from]++; } mt=powmat(n,mt,k*k*k); rep(a,k)rep(b,k){ ans=(ans+mt.v[a*k*k+b*k][0])%mod; } cout<