#include #include #include #define llint long long #define mod 998244353 using namespace std; struct Matrix{ int h, w; vector mat; Matrix(){ h = w = 0; } Matrix(int h, int w){ this->h = h, this->w = w; mat.resize(h*w); } llint& at(int i, int j){ return mat[w*(i-1)+(j-1)]; } static Matrix ident(int size){ Matrix ret(size, size); for(int i = 1; i <= size; i++) ret.at(i, i) = 1; return ret; } Matrix operator*(Matrix& ope){ Matrix ret(h, ope.w); for(int i = 1; i <= h; i++){ for(int k = 1; k <= w; k++){ if(at(i, k) == 0) continue; for(int j = 1; j <= ope.w; j++){ ret.at(i, j) += at(i, k) * ope.at(k, j) % mod; ret.at(i, j) %= mod; } } } return ret; } Matrix pow(llint n){ if(n == 0) return ident(h); if(n % 2){ return pow(n-1) * (*this); } else{ Matrix tmp = pow(n/2); return tmp * tmp; } } void print(){ for(int i = 1; i <= h; i++){ for(int j = 1; j <= w; j++){ cout << at(i, j) << " "; } cout << endl; } } }; llint n, x; int main(void) { cin >> n >> x; Matrix mat(x*x*x, x*x*x); llint ans = 0; for(int c = 0; c < x; c++){ //cout << c << endl; for(int i = 1; i <= x*x*x; i++){ for(int j = 1; j <= x*x*x; j++){ mat.at(i, j) = 0; } } for(int i = 0; i < x; i++){ for(int j = 0; j < x; j++){ for(int k = 0; k < x; k++){ int s = i*x*x + j*x + k+1; int ni = i, nj = j, nk = k, ns; ni = (i+1)%x, ns = ni*x*x + nj*x + nk + 1; mat.at(ns, s) = 1; ni = i; nj = (j+1)%x, ns = ni*x*x + nj*x + nk + 1; mat.at(ns, s) = 1; nj = j; int la = i, rc = (c-j+x)%x; nk = k + la*rc%x, nk %= x; ns = ni*x*x + nj*x + nk + 1; mat.at(ns, s) = 1; } } } Matrix pmat = mat.pow(n); for(int i = 0; i < x; i++){ for(int j = 0; j < x; j++){ for(int k = 0; k < x; k++){ if(j != c || k != 0) continue; int s = i*x*x + j*x + k + 1; ans += pmat.at(1, s); ans %= mod; } } } } cout << ans << endl; return 0; }